# A Comprehensive Review of Computation-Based Metal-Binding Prediction Approaches at the Residue Level

Authors: Nan Ye; Feng Zhou; Xingchen Liang; Haiting Chai; Jianwei Fan; Bo Li; Jian Zhang
Journal: BioMed Research International
 (2022)

---

## Abstract

Clear evidence has shown that metal ions strongly connect and delicately tune the dynamic homeostasis in living bodies. They have been proved to be associated with protein structure, stability, regulation, and function. Even small changes in the concentration of metal ions can shift their effects from natural beneficial functions to harmful. This leads to degenerative diseases, malignant tumors, and cancers. Accurate characterizations and predictions of metalloproteins at the residue level promise informative clues to the investigation of intrinsic mechanisms of protein-metal ion interactions. Compared to biophysical or biochemical wet-lab technologies, computational methods provide open web interfaces of high-resolution databases and high-throughput predictors for efficient investigation of metal-binding residues. This review surveys and details 18 public databases of metal-protein binding. We collect a comprehensive set of 44 computation-based methods and classify them into four categories, namely, learning-, docking-, template-, and meta-based methods. We analyze the benchmark datasets, assessment criteria, feature construction, and algorithms. We also compare several methods on two benchmark testing datasets and include a discussion about currently publicly available predictive tools. Finally, we summarize the challenges and underlying limitations of the current studies and propose several prospective directions concerning the future development of the related databases and methods.


---

## Body

## 1. Introduction

Metal ions are certain atom compounds that usually form cations that have (a) positive electric charge(s). Metal ions play pivotal roles in protein structure, function, regulation, and stability [1, 2]. Common metal ions include zinc (Zn2+), calcium (Ca2+), magnesium (Mg2+), manganese (Mn2+), iron (Fe3+ or Fe2+), copper (Cu2+), cobalt (Co2+), sodium (Na+), potassium (K+), and nickel (Ni2+) ions. Recent estimates have shown that approximately 30%-40% of proteins require one or several metal cofactors to together express biological function [3]. The proportion varies in different types of organisms or tissues. For instance, K+ is mostly found inside the cell, while Na+ is abundant outside of the cell [4]. Mn2+ is found accumulated in leafy green plants [5]. In the human body, Ca2+ accounts for approximately 1.5% of total body weight. The bulk of Ca2+ is aggregated in bones and teeth [6].Metal ion binding proteins, i.e., metalloproteins, play critical roles in a biological and chemical process in cellular reactions [7]. Inside the cell, the dynamic homeostasis of the metal ions is strongly connected and delicately tuned [8]. Reinhard et al. claimed that K+ and Na+ are involved in processing cell signaling, intercellular communication, and maintaining tissue electrolyte balance [9]. A small change in the concentration of metal ions may shift the effects of metal ions from natural beneficial to harmful [10]. A recent study pointed out that metalloproteins are associated with degenerative diseases, including Parkinson’s disease and Alzheimer’s disease [11]. For instance, α-synuclein (Cu2+-protein complex) constitutes the main component in Lewy bodies in Parkinson’s disease [11]. Mn2+ and Fe3+ are responsible for inducing tangle pathology in Alzheimer’s disease [8]. The aging of the brain or the development of diseases is associated with the deregulation of the management of metal ions [10]. Particularly, recent evidence indicates that if different types of metalloproteins interact in a certain salt solution, the potential galvanic erosion may dissolute the compound surface and result in inducing tumor formation [12, 13].Elucidated protein-metal ion interactions rely in part on the advancement of various accurate characterizations and predictions of metalloproteins at the residue level. The traditional methods that are used to identify metal-binding conformation or binding residues include biophysics- or biochemistry-related wet-lab experiments, such as mass spectrometry [14], X-ray crystallography [15], and surface Plasmon resonance [16]. Since these technologies need expensive instruments, complex procedures, and elaborate labors, they shall benefit from the recent development of computation-aided methods.We found 12 reviews that focused on the topic of exploring metal-binding residues or proteins in the past decade [7, 10, 11, 17–24]. Mallick et al. shed light on in silico methods including nine predictive tools and discussed the intrinsic mechanisms of metal-protein binding [24]. Thirumoorthy et al. investigated metallothionein isoforms and their role in pathophysiology [17]. They also provided the analysis of how metallothionein impact complex disease scenarios. In [18], the authors focused on structural variability and corresponding mechanisms of polymorphic amyloid oligomers complexed with metal ions. Bal et al. discussed ability constants, dissociation rates, and coordination chemistry of metal-binding residues in albumin [19]. Roohani et al. reviewed the literature related to zinc biochemical and physiological functions, metabolism, and zinc bioavailability in the human body [20]. The authors in [21] summarized the web tools that were proposed to identify metal-binding residues. Liu et al. systematically analyzed the structural features of Zn2+-binding sites and proposed an online predictor [22]. Akcapinar and Sezerman collected and surveyed computational toolboxes designed for the recognition of metal-binding sites or metalloproteins [7]. Quintanar and Kim summarized the research in degenerative diseases related in metal ions [11]. Witkowska and Rowińska-Żyrek overviewed the analytical and biophysical methods utilized for studies on metal-protein interactions [23]. Krzywoszyńska detailed the involvement of metal ions in signaling processes within the cell and its influence in health and disease [10]. Rauer et al. scrutinized computational approaches that are associated with the prediction of protein functional sites and also discussed metal-binding related works [25].Broadly speaking, these reviews discuss some aspects of the predictive methods. Some of them provide sufficient coverage of databases and predictive models and discuss the challenges and limitations of considered approaches. These reviews bring informative clues for the following researchers in this field. From the pertinence of the research, the prediction of metal-binding can be divided into general and specific approaches. The former recognizes metal-binding residues without considering their types, while, the latter is aimed at identifying one or several specific metal-protein interactions. According to the basic design and scheme, we classify these methods into four categories, namely, learning-, docking-, template-, and meta-based methods.This review covers a comprehensive set of 44 computation-based methods, and 25 of them were published in the past three years. Specifically, we survey 32 learning-based, 4 docking-based, 6 template-based, and 2 meta-based methods. Depending on whether the structure of a target protein is known or available, we further divide learn-binding methods into the structure- and sequence-based ones. We discuss their benchmark datasets, features, algorithms, and measurements, respectively. We also detail the docking-, template-, and meta-based methods and point out their advantages and limitations.


## 2. Public Databases for Metal Binding

The development in biochemistry and biophysics leads to a fast increasing number of protein-metal ion binding complexes. Figure1 draws the top 10 metal-binding annotations in PDB. Our survey reveals that Zn2+, Ca2+, and Mg2+ occupy the top three prevalent metal ions. The Zn2+ is currently the best-explored and described metal ion [26]. Zn2+ participates in many biological processes, such as metabolism, immune system, neurotransmission, hormone secretion, and signaling [27]. According to a rough statistic, approximately 10% of eukaryotic proteins bind Zn2+ [28]. Ca2+ is mainly aggregated in bones and teeth vertebrates [29]. It helps form solid support structures through biomineralization [6]. Mg2+ is usually associated with solvent water molecules, which endow it with a good capability of binding affinity with proteins and movement. The solvation state of Mg2+ usually serves as the enzyme in which Mg2+ acts as a coenzyme [6].Figure 1
Fraction of top 10 metal-binding interactions that stored in PDB (date: December 20, 2021).Besides RCSB PDB (https://www.rcsb.org/) [30], recent years have witnessed several specific databases that collect, categorize, and store these metal-protein interactions. Table 1 summarizes the publication year, considered metal ions, size of the database, web link, citations, and availability for the recently released database. We use citations as a one direct and good way to quantify the impact of these resources within the community [31]. The citation counts were collected from Google Scholar (https://scholar.google.com/) on December 20, 2021.Table 1
Summary of recently released database of metal ion binding interactions.
NameYearConsidered metal ionsNumber of sitesWeb linkRef.CitationAvailabilityInterMetalDB2021All metal ion binding6,423https://intermetaldb.biotech.uni.wroc.pl/[26]N/AYesMeLAD2020All metal ion bindingN/Ahttps://melad.ddtmlab.org/[33]9YesZincBindDB2019Zn24,992https://github.com/samirelanduk/ZincBindDB[49]23YesMetalPDB (v2)2018All metal ion bindingN/Ahttp://metalweb.cerm.unifi.it[34]90NoBioLiP2013All metal ion binding146,969https://zhanggroup.org/BioLiP/[36]446YesZiFDB (v2)2013ZnN/Ahttp://bindr.gdcb.iastate.edu/ZiFDB[37]25NoMetalPDB (v1)2013All metal ion bindingN/Ahttp://metalweb.cerm.unifi.it[35]108NoBioMe2012All metal ion binding20,307http://metals.zesoi.fer.hr[39]30NoMetLigDB2011Zn, Mn, Fe, Ni, mg, cu, co, Mo732http://silver.sejong.ac.kr/MetLigDB[40]13YesMIPS2010All metal ion bindingN/Ahttp://dicsoft2.physics.iisc.ernet.in/mips/[41]28YesMEDB2010All metal ion bindingN/Ahttp://www.uohyd.ernet.in/anambs/[42]14NoZiFDB (v1)2009ZnN/Ahttp://bindr.gdcb.iastate.edu/ZiFDB[38]87NoMetalMine2009All metal ion binding412http://metalmine.naist.jp[43]3NoMetal-MACiE2009All metal ion bindingN/Ahttps://www.ebi.ac.uk/thornton-srv/databases/Metal_MACiE/home.html[44]60YesZifBASE2009ZnN/Ahttps://web.iitd.ac.in/~sundar/zifbase/[45]35YesMESPEUS2008Na, mg, K, ca, Mn, Fe, co, Ni, cu, Zn34,896http://eduliss.bch.ed.ac.uk/MESPEUS/[46]102NoMSDsite2005All metal ion bindingN/Ahttp://www.ebi.ac.uk/msd-srv/msdsite[47]122YesMDB2002All metal ion bindingN/Ahttp://metallo.scripps.edu/[48]276No1We estimate the availability on December 1st, 10th, and 20th of 2021, respectively.Specifically, InterMetalDB collects and presents metal ion binding proteins from RCSB PDB. It uses MMseq2 [32] to cluster the structure chains with the 50% sequence identity. Then, it groups similar binding sites and selects the best-resolution structure as a representative. MeLAD is a metalloenzyme-ligand association database, which contains structural data, metal-binding pharmacophores, and ligand chemical similarity of metalloenzyme-ligand interactions [33]. MetalPDB details the local environment, three-dimensional (3D) structure, secondary structure, and solvent accessibility of the metal ion binding sites [34, 35]. BioLiP is a semimanually curated database, which includes protein-peptide, protein-nucleic acid, and protein-ligand annotations [36]. BioLiP stores and periodically updates all types of metal ion binding information from PDB. ZiFDB is a database that collects information about individual zinc fingers, engineered zinc-finger arrays, and related target sequences [37, 38]. BioMe provides a web interface for biologists to capture coordination numbers, distances, geometry, and percentage of monodentate and bidentate bound aspartic acid and glutamic acid carboxyl groups [39]. MetLigDB is specially designed to select chelating groups or chemical moieties that might be presented in the inhibitor of a metalloprotein [40]. MIPS stores the geometric information, macromolecular function, different chemical behavior of metals, and metalloproteins [41]. MEDB presents quantitative information on metal-binding sites in protein structures and can be used for the identification of trends or patterns in the metal-binding sites [42]. MetalMine automatically collects and categorizes different types of metal-binding sites that derived from the structures of protein-metal-ion complexes [43]. Metal-MACiE gathers all available metalloenzymes and includes structural and functional information of metal ions in the context of the catalytic mechanisms of these metalloenzymes [44]. ZifBASE deposits engineered and natural zinc finger proteins and provides sequences and structural features and associated potential target sites of these proteins [45]. MESPEUS [46] focuses on the geometry of metal sites in proteins at resolution≤2.5 Å. It provides an open web interface for further identifying and displaying the metal sites. MSDsite deposits computation-based metal-binding geometries and residues [47]. MDB offers quantitative information about metalloproteins [48]. It provides functions to analyze the binding attributes such as metal-ligand bond distances and side-chain torsion angles in metal sites.We show that twelve source databases are designed for all metal ion binding data. Two databases, namely, MetLigDB and MESPEUS, consider several types of metal ions. There are four specific zinc-binding-related databases. Our survey also reveals that BioLiP is the most favored database, given the fact that its citations are average about 56 (446/8≈56) per year. Moreover, we notice that only half of the databases are available. Thus, we recommend that future databases shall be chronically maintained, periodically updated, and easy expanded.


## 3. Method Development of Metal-Binding Prediction

Figure2 illustrates the flowchart of computation-based methods for the prediction of metal-binding residues. Generally, based on the basic design and scheme, these methods can be categorized into four groups. The learning-based methods regard the identification of metal-binding residues as a typical classification problem and attempt to use machine learning or deep learning algorithms to construct prediction models. The docking-based approaches are aimed at finding proper binding conformation as well as the appropriate target binding residues by scanning protein surface. The scoring functions are introduced to assess the selected pockets and quantify the strength of binding affinity. The template-based methods are designed to select the optimal template structures for a given unknown protein. Then, they map and transfer the binding annotations from similar spatial conformation to the target protein. By contrast, the meta-based methods focus on combing the predictions from other methods in order to build more accurate predictors.Figure 2
The flowchart of computation-based methods for prediction of metal-binding residues.


### 3.1. Benchmark Datasets

The sequences and structures of protein-metal ion complexes are available in public databases for the end-users to customize the benchmark datasets. As shown in Table2, the considered methods use various numbers of sequences/chains, ranging from several dozens to thousands. Besides that, protein complexes with high resolution indicate relatively more comprehensive and accurate annotations of protein-metal ion interactions. According to our survey, 12 out of 23 sequence-based and 6 out of 9 structure-based methods filter the candidate complexes using high resolution with ≤3 Å. Some methods [50–57] remove the sequences/chains whose lengths are less than 50 residues (or 45 residues [58]) since they might be potential segments or peptides. To build an unbiased dataset, it is necessary to remove homologous or redundant proteins. The cutoff threshold which researchers choose varies from minimal 25% to maximal 90%. Generally, a higher identity means a higher chance in local alignments [59]. The literatures in [60, 61] point out that if a pair of proteins have a sequence identity lower than 30%, they have little chance to share the same biological processes. Three tools, namely, BLASTclust [62], PISCES [63], and CD-HIT [64], are mainly used to cluster homologous proteins.Table 2
Summary of learning-based methods.
TypeMethod1Ref.YearMetal ion binding2Dataset3ResolutionSequence similarity (tool)4Prediction model5Cross-validationIndependent testMeasurements6Sequence-basedLiu et al.[56]2020Zn, Cu, Fe, Co, Mn, Ca, Mg, Na, K5,340≤3 Å30% (CD-HIT)RF5-fold√SN, SP, ACC, MCCMIonSite[90]2019Zn, Ca, Mg, Mn, Fe, Cu, Fe, Co, Na, K, Cd, Ni7,676N/A30% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCC, AUCMPLs-Pred[91]2019General metal ions1,492N/A30% (CD-HIT)RF10-fold√SN, SP, ACC, MCCSXGBsite[92]2019Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)GBM5-fold√SN, SP, ACC, MCC, AUCWang et al.[50]2019Zn, Cu, Fe, Mn, Ca, Mg, Na, K5,146≤3 Å30% (N/A)SVM, SMO5-fold√SN, SP, ACC, MCCznMachine[51]2019Zn2,043≤3 Å30% (BLASTclust)SVM, NN5-fold√SN, SP, ACC, MCC, PRE, AUCSSWPNN[93]2019Zn213≤2.5 Å70% (N/A)SVM, NN5-fold√SN, SP, PRE, F1, MCC, ACCZinCaps[94]2019Zn738≤3 ÅN/A (N/A)CN5-fold√SN, SP, ACC, MCC, AUCHaberal and Oğul[65]2018General metal ions2,727N/AN/A (N/A)CNN5-fold√SN, ACC, PRE, F1ZincBinder[87]2018Zn738≤2.5 Å30% (PISCES)SVM5-fold√SN, SP, ACC, MCC, AUCEC-RUS[95]2017Ca, Mg, Mn, Fe, Zn4,421N/A40% (PISCES)WSRC5-fold√SN, SP, ACC, MCC, AUCCao et al.[52]2017Zn, Cu, Fe, Co, Mn, Ca, Mg, K, Na5,340≤3 Å30% (CD-HIT)SVM5-fold√SN, SP, ACC, MCCKumar[96]2017Cu, Ca, Co, Fe, Mg, Mn, Ni, Zn3,922N/A50% (CD-HIT)RF10-fold√SN, SP, ACC, MCCDeepMBS[97]2017General metal ions2,727≤3 ÅN/A (N/A)CNN5-fold√SN, PRE, F1Qiao et al.[98]2017Ca2,239N/A30% (CD-HIT)SVM5-fold√SN, ACC, PRE, MCC, AUCIonCom[99]2016Zn, Cu, Fe, Ca, Mg, Mn, Na, K1,374N/A30% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCCJiang et al.[77]2016Ca1,885≤3 Å25% (N/A)SVM5-fold√SN, SP, ACC, MCCTargetCom[53]2016Cu, Fe, Zn1,373≤3 Å40% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCCOSML[100]2015Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)SVM5-fold√SN, SP, ACC, MCCTargetS[101]2013Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)SVM, AdaBoost5-fold√SN, SP, ACC, MCC, AUCETMB-RBF[102]2013General metal ions55N/A20% (BLASTclust)RBFN10-fold√SN, SP, ACC, MCCZincExplorer[103]2013Zn392≤3 ÅN/A (N/A)SVM5-fold√SN, SP, PRE, MCC, AUPRCHorst et al.[58]2010Ca635≤2.1 Å35% (N/A)LR10-fold√MCC,AUC,AUPRCStructure-basedNguyen et al.[104]2021Mn, Fe, Co, Ni, Cu, Zn9,955≤2.5 Å90% (N/A)RF5-fold×ACCTMP-MIBS[54]2021General metal ions427N/A40% (CD-HIT)RF10-fold√SN, SP, ACC, MCC, AUCZincbindpredict[105]2021ZnN/A≤ 2 Å40% (CD-HIT)RF5-fold√SN, PRE, F1, MCCWang et al.[81]2021Zn, Cu, Fe, Ca, Mg, Mn, Na, K, Co5,340≤3 Å30% (N/A)MLP,SVM5-fold√SN, SP, ACC, MCCDELIA[80]2020Ca, Mn, Mg3,966N/A30% (CD-HIT)CNN5-fold√SN, PRE, MCC, AUCHu et al.[57]2020Zn, Cu, Fe, Co, Mn, Ca, Mg, Na, K5,340≤3 Å30% (CD-HIT)GBM5-fold√SN, SP, FPR, ACC, MCCMetalExplorer[79]2017Ca, Co, Cu, Fe, Ni, Mg, Mn, Zn3,192≤2.5 Å30% (CD-HIT)RF5-fold√SN, FPR, PRE, AUC, AURPCFINDSITE-metal[55]2011Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn860N/A35% (PISCES)SVM2-fold√ACC, SPC, PPVZinc identifier[78]2011Zn1,103≤2.5 ÅN/A (N/A)RF5-fold√SN, PRE, SP, FPR, AUC, AUPRC1The name of each method is provided in either the publication or the last name of its first author. 2General metal ions mean that the related predictor does not differentiate the types of metal ion binding. Otherwise, we list the specific types of metal-binding in detail. 3The number represents the size of the benchmark dataset. 4The value reveals the protein similarity threshold in the benchmark dataset. The content in the blanket indicates the tool that is used for clustering proteins. 5SMO: sequential minimal optimization; SVM: support vector machine; WSRC: weighted sparse representation based classifier; NN: neural network; CN: capsule network; CNN: convolutional neural networks; RF: random forest; GBM: gradient boosting machine; RBFN: radial basis function networks; LR: logistic regression; MLP: multilayer perceptron. 6SN: sensitivity/recall; SP: specificity; ACC: accuracy; MCC: Matthew’s correlation coefficient; PRE: precision; F1: F1-score; AUC: area under the ROC curve; AUPRC: area under the precision recall curve; FPR: false positive rate (FPR=1-SP).


### 3.2. The Validation and Evaluation Metrics

#### 3.2.1. Cross-Validation and Independent Test

To construct a predictor with high accuracy and decent generalization ability, it is necessary to avoid potential overfitting. In practice, cross-validation and independent test are two popular ways (Table2) to evaluate the proposed models [31]. Specifically, k-fold cross-validation is usually adopted on the training dataset when building the prediction model and optimizing the related parameters [61]. First, the training dataset is equally divided into k parts. The division can be done at residue level or protein level. Next, k-1 subsets are used to train the model, and the last one subset is used for testing. The procedure repeats k times until every subset is been predicted. The performance of the model is usually evaluated by averaging the results of the k repeats.


#### 3.2.2. Performance Measures

According to Table2, the measures that used to evaluate the performance of the predictors can be divided into binary value-based and propensity score-based ones. The former needs preset thresholds to compute the number of putative binding residues and nonbinding residues. These measures include sensitivity (SN)/recall/true positive rate (TPR), specificity (SP), false positive rate (FPR, FPR=1-SP) precision (PRE), accuracy (ACC), F1-score (F1), and Matthew’s correlation coefficient (MCC). They are defined as follows:
(1)SN=TPR=TPTP+FN,SP=1−FPR=TNTN+FP,PRE=TPTP+FP,F1=2TP2TP+FP+FN,ACC=TP+TNTP+FN+TN+FP,MCC=TP×TN−FN×FPTP+FN×TP+FP×TN+FP×TN+FN,where TP (true positive) indicates the number of correctly recognized metal-binding residues, FP (false positive) means the number of non-metal-binding residues that are incorrectly predicted as metal-binding residues, TN (true negative) stands for the number of correctly predicted non-metal-binding residues, and FN (false negative) is the number of metal-binding residues that are incorrectly predicted as non-metal-binding residues.The prediction of metal-binding residues is a typical imbalanced classification problem. That is, the number of metal-binding residues is much less than that of the non-metal-binding ones. Therefore, F1-score and MCC are regarded as key criteria since they are featured by assessing the prediction performance for both metal-binding and non-metal ion binding residues.The propensity score-based measures include receiver operating characteristic curve (ROC curve) and precision-recall curve (PR curve). The ROC curve draws the TPR (true positive rate) against the FPR (false positive rate) at various thresholds. The AUC computes the area under the ROC curve and can be used to quantify the ROC curve. The PR curve plots PRE values on they-axis and recalls values on the x-axis, and the AUPRC estimates the area under the PR curve.


### 3.3. Learning-Based Methods

Learning-based methods treat the recognition of metal-binding residues as a typical pattern recognition problem. Specifically, the metal-binding residues and nonbinding ones are encoded by using mathematical descriptors, i.e. features. According to the information that used to compute the features, the learning-based methods can be further categorized into sequence-based and structure-based methods. The former only needs simple protein sequences to extract features when encoding the binding residues. These features include sequence directly derived, evolutionary profile-based, and putative structure-based features, while the latter uses both sequence and native structure data to mathematically describe a binding residue. We make a comprehensive literature search and collect 23 sequence- and 9 structure-based methods that were published after the year 2010.


#### 3.3.1. Feature Construction

(1) Sequence Directly Derived Features. We define sequence directly derived features as the ones that are computed from protein primary sequences without using any other information. In Figure 3, 14 out of 32 considered methods consider amino acid composition [50], which quantifies the relative difference in abundance of a given amino acid type [65, 66]. Amino acid pairs, or dipeptides, are based on the observation that amino acid pairs show different propensities in protein structure and function. For instance, pairs of lysine are found present in close spatial vicinity [67]. Moreover, the concept of k-spaced amino acid pairs is introduced in [68]. It calculates the amino acid pairs with k spaces between two residues. Our survey also shows that the majority of studies use physicochemical properties to describe the local environment of the metal-binding residues. The basic physicochemical environment of a metal-protein binding interface is reflected by the specific roles the metal plays in biostructural chemistry and protein function. These properties are crucial since they underpin many of the functional roles of metal ions. These properties include aliphatic [69], sulphur [70], aromatic [71], hydrophobic [72], charge [73], polar [74], positive [73], acidic [75], and hydroxylic [76]. The position-related features mainly consider the influence of the specifically located residues, such as autocross covariance [77] and sequence length [78, 79].Figure 3
Summary of the feature construction and selection for learning-based methods. The light green cells indicate sequence directly derived features. The light blue cells stand for profile-based features. The light red cells mean putative structure-based features. The light grey cells are native structure-based features.1In the amino acid pair column, the cells without annotations indicate original amino acid pairs; the cells annotated using “k-spaced” means k-spaced amino acid pairs. 2PWM: position weight matrix; PSSM: position specific scoring matrix; PMS: position matrix scoring; PAM: point accepted mutation; EMS: evolutionary matrix scoring. 3TS: topology structure; DM: distance matrix; GTN: graph theoretic network; RRCG: residue-residue contact graphs. 4FFS: forward feature selection; EB: experience-based; BA: Boruta algorithm; mRMR: minimum-redundancy maximum-relevancy; MDGI: mean decrease Gini index [50–51, 53–58, 65, 77–81, 87, 90–105].(2) Evolutionary Profile-Based Features. Recent studies [54, 56, 57, 80, 81] pointed out that functional or structural important residues tend to show higher evolutionary conservation. The conserved residues are usually involved in enzyme activity, ligand binding, or protein structural stability [82]. The conserved residues and regions can be identified by multiple sequence alignment [83]. These multiple sequence alignments, also named conservation profiles, include aligning families of homologous sequences and having knowledge of their evolutionary relationships [84]. For an unknown protein, although its accurate function is not available, it is expected that we can use its homologous proteins to speculate the function since they share the similar evolutionary profile [85]. Many studies use position-specific scoring matrix (PSSM), which is computed from PSI-BLAST [62], to quantify the evolutionary conservation. PSSM scores the substitution probability of each residue in the protein being substituted by other types of amino acids. Liu et al. [56] and Hu et al. [57] set different weights according to the positions of considered residues within the window and construct position weight matrix (PWM). Wang et al. proposed a customized position matrix scoring (PMS) algorithm, which uses known sequence patterns to describe the composition of amino acids at different positions [50]. Haberal and Oğul introduced a point accepted mutation (PAM) scoring matrix, which measures the rate at which point mutations that substitute one residue for another during evolution [65]. Jiang et al. adopted evolutionary matrix scoring (EMS) algorithm to extract the position conservation of amino acid residues from segments with low dimension feature parameters [77].(3) Putative Structure-Based Features. For an unknown protein, although the accurate function is not available, it is expected that we can use its homologous proteins or template structures to speculate the structure. The secondary structure mainly involves α-helix, β-sheet, and coil, which are fundamental elements of protein tertiary structure [86]. Natively disordered or unstructured regions are proved to be associated with molecular assembly, protein translation, modification, and molecular recognition [78, 79, 87]. Previous studies [79, 87] indicate that disordered regions are strongly correlated with local solvent accessibility areas. Figure 3 reveals that 16 methods introduce secondary structure features and 3 approaches use disorder features, respectively. The secondary structure can be obtained from the primary sequence by using PSIPRED [88]. Putative intrinsic disorder data can be computed by using DISOPRED [89].(4) Structure-Based Features. The structure-based features include descriptors that are computed from protein 3D structure. These features include solvent exposure, B-factor, spatial cluster properties, and native secondary structure. Compared with the abovementioned putative structure-based features, the native structure-based features are more accurate since they are directly computed by using residue coordinate data. Besides that, a residue contact network is also considered by some literature. In [79], two residues are defined as being in contact if the distance of their Cα atoms is less than a predefined cutoff distance of 6.5 Å. These features include clustering coefficient, degree, density, distance, topology structure, and graph theoretic network [55, 79, 80].


#### 3.3.2. Sliding Window Optimization and Feature Selection

As shown in Figure3, many methods adopt a sliding window scheme when they construct different types of features. It is because residues in proteins are influenced by adjacent residues. Besides that, binding residues tend to cluster together. If a central residue is a native-binding residue, its adjacent residues usually have a relatively higher chance to bind the same ligands. Usually, the residues with a long distance away have a lower impact on the central residues when compared with the residues with short distance. Figure 3 summarizes that 19 out of 23 methods use the sliding window scheme. The size of the shortest window is 3 [105], while the size for the longest one is 25 [94]. Some studies [50, 52, 53, 56, 57, 81, 90, 105] use more than one type of window because they consider different types of metal-binding residues. A long window means the introduction of more features.However, a bigger number of features do not absolutely mean a better prediction performance [106, 107]. The existence of potential “bad” features may interfere with the classifiers and cause unpredictable consequences [108]. The so-called “bad” features include irrelevant and redundant ones. To avoid their terrible influences, it is necessary to perform feature selection before training the model [109]. Figure 3 reveals that 6 out of 32 methods adopt feature selection before training the model. These feature selection approaches include forward feature selection [79, 102], experience-based [104], Boruta algorithm [57, 81], minimum-redundancy maximum-relevancy [79], and mean decrease Gini index [78].


#### 3.3.3. Prediction Algorithms

Learning-based methods use machine-learning or deep-learning-based algorithms to train the model and perform predictions [110]. As shown in Table 2, a variety of algorithms are introduced for solving the problem of correctly recognizing metal-binding residues. Support vector machine (SVM) is a popular machine learning algorithm in bioinformatical research. It is aimed at finding a hyperplane or decision boundary that can segregate a high-dimensional space [111]. Particularly, it uses kernel functions to reduce computation time to avoid strapping into dimension disaster [112]. Sequential minimal optimization (SMO) is an algorithm that is specially used for training support vector machines [113]. The procedure of training large data by SVM usually leads to a complex quadratic programming optimization problem [114]. SMO breaks large programming optimization problems into small ones, which endows SVM a good generalization on large data [113]. The idea of a neural network (NN) comes from the work system of neurons in the biological brain [115]. It learns the correlations between inputs and outputs, making generalizations and build models [116]. The NN algorithm assigns and adjusts different weights for neurons and edges as learning proceeds. The radial basis function network (RBFN) is a variant of the original NN [117]. It adopts radial basis functions as activation functions, which can be used for accelerating learning speed due to their universal approximation [118]. The multilayer perceptron (MLP) algorithm is an improved back propagation NN [119]. It mainly includes three procedures: forward propagation, error evaluation, and error backpropagation [120]. The MLP is featured by its strong generalization and fault tolerance [121]. Therefore, it is proved to be an efficient classification algorithm. The logistic regression (LR) adopts a logistic function to model the probability of an unknown sample being a certain class [122].Our survey also reveals that the ensemble algorithms are favored by eight studies. The random forest (RF) aggregates the predictions of all the decision trees and performs decisions by most trees [123]. RF can be used for classification, regression, and optimization problems [124]. Adaptive boosting (AdaBoost) is aimed at combining weak learners with strong ones [125, 126]. The key point of AdaBoost is to ensure the diversities of individual learners, which makes it a good generalization ability [90, 127]. The gradient boosting machine (GBM) is another popular ensemble algorithm. During the iterative process, GBM dynamically increases the weight of wrong recognitions and reduces that of the correct ones [128–131]. It should be noted that GBM focuses on the sample residual of the previous iteration instead of the sample itself [132].Besides machine-learning algorithms, recent studies also use deep-learning methods in this research field. The convolutional neural network (CNN) is one of the most prevalent algorithms that is widely used in bioinformatics [133]. The CNN consists of three main layers, which are the convolutional layer, pooling layer, and fully connected layer [134–136]. Although the CNN is proved to be powerful in dealing with a variety of problems, it performs badly when facing samples with different sizes and orientations [137, 138]. To overcome this shortcoming, the capsule network (CN) is proposed to estimate features of objects by incorporating dynamic routing algorithms [139, 140]. Our review finds two studies use CNN [65, 97] and one uses CN [94].


### 3.4. Docking-Based Methods

The investigation on the protein-metal complex helps biologists to understand the mechanism of protein-metal interactions. Protein-ligand docking approaches are always based on molecular structure and are used to explore biomolecular interactions and mechanisms [141]. It can be adopted to predict binding conformation as well as the appropriate target binding residues [142, 143].As shown in Figure2, the docking-based methods mainly include three steps: searching algorithm, scoring function, and docking assessment [141]. The searching algorithm focuses on creating an optimum number of configurations that properly include the determined binding modes [144]. To reduce computation time, it is necessary to make a balance between the computational expense and the searching space. The scoring function includes a series of mathematical functions that quantify the strength of binding affinity [145]. The energy-based scoring functions are always introduced to score the potential interactions between the protein and the corresponding ligands [141]. The frequently used functions include empirical-based, knowledge-based, and consensus-based ones. Finally, the putative docking can be evaluated by using docking accuracy and the correlation between putative and native docking scores [145]. Figure 4 illustrates the structure of a calmodulin (PDB: 4HEX) that is secreted by Escherichia coli in Mus musculus [146]. Calmodulin is one of the most prevalent EF-hand calcium sensor proteins in eukaryotic cells [147]. It is a highly conserved and soluble protein, which activates enzymes and regulates many cellular functions. 4HEX has three Ca2+-binding and two Zn2+-binding sites. Ca2+-binding causes a change in calmodulin conformation opening both globular domains and exposing hydrophobic surfaces that form binding sites for the target enzymes. Figure 4 shows that these three Ca2+ are in the pockets. The binding pockets are half-closed and buried, which substantially limits the capability of Ca2+ to escape. Two Zn2+-binding sites are surrounded by a shell of hydrophilic groups that are embedded into a larger shell of hydrophobic groups. The amino acid side chains providing ligands to Zn2+ in these structures often form hydrogen bonds with other residues [147].Figure 4
Ribbon and surface model of X-ray structure of Ca2+- and Zn2+-bound calmodulin (PDB: 4HEX) in Mus musculus. Red sphere represents bound zinc ion; green one indicates calcium ion; the spatial adjacent residues participating its coordination are shown by the stick model.
(a)(b)(c)In [148], He et al. proposed a docking-based predictor named mFASD. It first explored the local biochemical environment of potential functional atoms and then measured the distances between the atoms and bound metal. mFASD also claimed that it can differentiate different types of metal-binding sites. Zhou et al. improved the FEATURE-based calcium model and used the grid scan algorithm to recognize binding sites [149]. GaudiMM [150] adopted a multiobjective genetic algorithm to search metal-binding sites in biological scaffolds. BioMetAll focused on the conformation of the potential metal-binding site, associated with the geometric organization of the protein backbone [151]. It was also proved to have good performance on the applications including the modulation and mutation of the metal-binding residues. Table 3 summarizes the key notes of the abovementioned 4 docking-based methods.Table 3
Summary of docking-based, template-based, and meta-based methods.
TypeMethodYearNotesDocking-basedmFASD [148]2015Capture the characteristics of metal-binding sites and discriminate most types of these sitesZhou et al. [149]2015Use a FEATURE-based calcium model and convert high scoring regions into specific site predictionsGaudiMM [150]2019Find poses that satisfy metal-derived geometrical rules and use post optimizationsBioMetAll [151]2020Predict metal-binding sites with particular motifs, determine transient sites in structures, and predict potential mutations to generate convenient sitesTemplate-basedDeng et al. [152]2006Use a graph theory algorithm to find oxygen clusters of the protein (high potential for calcium binding)Goyal et al. [153]2008Describe generation of 3D-structural motifs for metal-binding sites from the known metalloproteinsLevy et al. [154]2009Analyze whether structural models based on remote homology are effective in predicting 3D metal binding sitesFunFOLD [155]2011Use an automated method for ligand clustering and identification of binding residuesFunFOLDQA [156]2012Use a fully automated agglomerative clustering approach for both ligand identification and residue selectionFunFOLD2 [157]2013Propose a method that include protein-ligand binding prediction and quality assessment protocolMeta-basedLi et al. [158]2017Integrate the results of ZincExplorer [103], zincFinder [159], and zincPred [160]IBayes_Zinc [161]2019Adopt Bayesian method and combine the predictions from ZincExplorer [103], zincFinder [159], and zincPred [160]


### 3.5. Template-Based Methods

It is well known that protein structure determines function, and similar interface conformation indicates similar bound regions [162]. The template-based methods are based on the abovementioned hypothesis. Therefore, the most important thing for template-based methods is to find and validate proper structural templates. The fold recognition algorithms, which quantify the best matches from candidate templates, are commonly used to select the optimal template structures [163]. Next, the selected templates are used to map onto the target protein given the alignments with the template structures [164].As shown in Table3, Goyal and Mande analyzed the metal-binding sites by using structure templates and designing 3D motifs for several types of metal-binding interactions [153]. In [154], the authors analyzed whether structural models based on remote homology are effective in recognizing structural metal-binding residues based on simple protein primary sequences. Deng et al. applied a graph theory algorithm to identify, predict, and analyze calcium-binding residues [152]. However, it should be noted that this strategy produces good prediction performance when a decent complex is available as a template. If the template structure information is not available, this strategy might have poor predictions [164]. The FunFOLD was an automatic method that uses protein structure superposition of distantly related templates to a modelled protein for the clustering of ligands and prediction of metal binding residues [155]. The FunFOLDQA [156] approach determined the reliability of our FunFOLD [155] by assigning the quality assessment scores. FunFOLD2 was a web server that integrated cutting edge function and putative 3D structures to identify metal-binding residues [157].


### 3.6. Meta-Based Methods

The meta-based methods use a meta-learning strategy from fewer samples than traditional machine learning models. Since meta-based methods can only use limited data, they must ensure that the data is featured with high accuracy. As a result, a meta-based approach always directly combines the predictions of other methods. It uses weights or voting strategy on the available propensity scores or binary values. Thus, the meta-based method promises a robust accurate prediction on the metal-binding residues. In [158], Li et al. collected the predictions from ZincExplorer [103], ZincFinder [159], and ZincPred [160] (Table 3). Then, they built a linear regression model and optimized corresponding parameters on the training dataset. They claimed that the meta-model, which was named meta-zincPrediction, improves the AUPRC by about 2%~8%. IBayes_Zinc [161] was another meta-based predictor for the identification of zinc-binding residues (Table 3). It firstly computed the predictions of zinc-binding probabilities from ZincExplorer [103], ZincFinder [159], and ZincPred [160]. Next, IBayes_Zinc processed the missing attribute values and adopts Bayesian theory [165] to construct a meta-based model. The performance on the independent dataset proved that the MCC value of IBayes_Zinc was about 5~13% higher than the considered three predictors.


### 3.7. Prediction Results

This review surveys 44 computation-based methods. It is necessary to make a consensus comparison for these methods. However, since there is no standard benchmark dataset and some methods are not currently available, we use two datasets that are used by some methods to perform the evaluations. The first dataset is compiled by Yu et al. in [101], which includes five types of metal ions binding annotations. The second dataset is obtained from [52], consisting of ten types of metal ions binding annotations.Figure5 illustrates the predictive performance on two benchmark test datasets. Details are provided in Table S1 and Table S2 in Supplementary Materials, respectively; the corresponding results are sourced from [56, 57, 81, 92]. We notice that the predictors show relatively big differences in recognizing various types of metal-binding residues. On Yu et al.’s dataset [101], TargetS shows the best results in predicting Ca2+-, Zn2+-, and Mn2+-binding residues; EC-RUS [95] performs best in recognizing Mg2+-binding residues; OSML [100] achieves the highest MCC on Fe3+-binding predictions. Besides that, Figure 5(a) indicates that all five methods show a decent performance on recognizing Fe3+-binding residues (MCC values close or higher than 0.4), compared with MCC close or less than 0.2 on Ca2+ binding residues. Figure 5(b) draws the bars of AUC values for SXGBsite [92], EC-RUS [95], and TargetS [101], respectively. These three predictors all achieve high AUC scores (close or higher than 0.9) on Zn2+- and Fe3+-binding residues. Figure 5(c) summarizes the results of ten metal ions binding residues on Cao et al.’s dataset [52]. Among these predictors, Liu et al. [56] performs the best on Zn2+, Fe3+, and Cu2+, compared to [81], Wang et al. shows best on Zn2+ and Cu2+, and Hu et al. [57] achieves the highest on Fe2+. Interestingly, the binding residues associated with relatively inactive metal (Zn2+, Fe3+, and Cu2+) ions show relatively better results compared to that of the active metal ions (Na+ and K+). Particularly, four methods all give better results on Fe3+-binding residues than that on Fe2+-binding residues, which keep consistent with our observations as mentioned above.Figure 5
Comparative assessment of several predictors on two benchmark dataset. (a) and (c) indicate the MCC bar charts for considered methods on different metal ion binding residues on Yu et al.’s and Cao et al.’s testing datasets, respectively. (b) draws the AUC values of three predictors on corresponding metal ion binding residues.
(a)(b)(c)


### 3.8. Publicly Available Tools

The publicly available standalone software or web server that implements the proposed approach provides convenience for biologists and researchers [79, 105, 122]. These tools help the community to repeat the results and build a platform for easy understanding and improvement. Table 4 summarizes the public availability of implementations for the considered methods. These 28 predictors are implemented as standalone software or web servers. Among these predictive tools, 16 (or 57%) of them are currently publicly available. Standalone software requires the biologists to build the same running environment. By contrast, the web server provides the most convenient since the users only need to submit their queries via the browser, and the server helps to do the computations. Three methods, namely, znMachine [51], ZincBinder [103], and FunFOLD [155], provide both web server and standalone software. TMP_MIBS [54] is designed to predict general metal-binding residues and deployed using the Python language. DELIA [80] requires PDB-formatted 3D coordinates input and produces both binary prediction and putative probability of a residue being potential specific metal-binding. BioMetAll uses a docking-based strategy to scan specific motifs, putative mutations, and binding residues. Another available docking-based method is mFASD, which distinguishes different types of metal-binding sites according to the interaction distances. MPLs-Pred [91], SXGBsite [92], MIonSite [90], OSML [100], EC-RUS [95], and TargetS [101] are all sequence-based predictive tools, which accepts FASTA-formatted input and produced the results of putative metal-binding residues. ZinCaps [94], SSWPNN [93], and ZincBinder [103] are specially designed for the identification of zinc-binding residues. FunFOLD [155], FunFOLDQA [156], and FunFOLD2 [157] are a series of template-based methods.Table 4
A breakdown of predictive tools of metal-binding residues.
MethodYearPlatform1Web linkAvailability2TMP-MIBS [54]2021SShttps://github.com/QuJing785464/TMP_MIBSYesWang et al. [50]2021WShttp://39.104.77.103:8081/lsb/HomePage/HomePage.htmlNoZincbindpredict [105]2021WShttps://zincbind.bioinf.org.uk/predict/NoDELIA [80]2020WShttp://www.csbio.sjtu.edu.cn/bioinf/delia/YesBioMetAll [151]2020SShttps://github.com/insilichem/biometallYesMPLs-Pred [91]2019WShttp://icdtools.nenu.edu.cn/YesSXGBsite [92]2019SShttps://github.com/Lightness7/SXGBsiteYesMIonSite [90]2019SShttps://github.com/LiangQiaoGu/MIonSite.gitYesznMachine [51]2019WS&SShttp://bioinformatics.fzu.edu.cn/znMachine.htmlNoZinCaps [94]2019SShttps://github.com/clemEssien/ActiveSitePredictionYesEC-RUS [95]2017SShttps://github.com/6gbluewind/protein_ligand_binding_siteYesMetalExplorer [79]2017WShttp://metalexplorer.erc.monash.edu.au/NoCao et al. [52]2017WShttp://60.31.198.140:8081/metal/HomePage/HomePage.htmlNoZincBinder [103]2017WS&SShttp://proteininformatics.org/mkumar/znbinder/YesSSWPNN [93]2017SShttp://net.jitsec.cn:88/UploadedImages/SSWPNN.rarYesJiang et al. [77]2016WShttp://202.207.29.245/NoTargetCom [53]2016SShttp://dase.ecnu.edu.cn/qwdong/TargetCom/TargetCom_standalone.tar.gzNoOSML [100]2015WShttp://www.csbio.sjtu.edu.cn/OSML/YesmFASD [148]2015SShttp://staff.ustc.edu.cn/liangzhi/mfasd/YesFunFOLD2 [157]2013WShttp://www.reading.ac.uk/bioinf/FunFOLD/FunFOLD_form_2_0.htmlYesZincExplorer [103]2013WShttp://protein.cau.edu.cn/ZincExplorerNoTargetS [101]2013WShttp://www.csbio.sjtu.edu.cn/TargetS/YesFunFOLDQA [156]2012SShttp://www.reading.ac.uk/bioinf/downloads/YesZincidentifier [78]2012WShttp://protein.cau.edu.cn/zincidentifier/NoFINDSITE-metal [55]2011WShttp://cssb.biology.gatech.edu/findsite-metal/NoFunFOLD [155]2011WS&SShttp://www.reading.ac.uk/bioinf/FunFOLD/YesGoyal et al. [153]2008WShttp://sunserver.cdfd.org.in:8080/protease/PAR_3D/index.htmlNoDeng et al. [152]2006SShttp://chemistry.gsu.edu/faculty/Yang/GG.htmNo1WS: web server; SS: standalone software. 2The availability was estimated on Dec 1st, 10th, and 20th of 2021, respectively.


## 3.1. Benchmark Datasets

The sequences and structures of protein-metal ion complexes are available in public databases for the end-users to customize the benchmark datasets. As shown in Table2, the considered methods use various numbers of sequences/chains, ranging from several dozens to thousands. Besides that, protein complexes with high resolution indicate relatively more comprehensive and accurate annotations of protein-metal ion interactions. According to our survey, 12 out of 23 sequence-based and 6 out of 9 structure-based methods filter the candidate complexes using high resolution with ≤3 Å. Some methods [50–57] remove the sequences/chains whose lengths are less than 50 residues (or 45 residues [58]) since they might be potential segments or peptides. To build an unbiased dataset, it is necessary to remove homologous or redundant proteins. The cutoff threshold which researchers choose varies from minimal 25% to maximal 90%. Generally, a higher identity means a higher chance in local alignments [59]. The literatures in [60, 61] point out that if a pair of proteins have a sequence identity lower than 30%, they have little chance to share the same biological processes. Three tools, namely, BLASTclust [62], PISCES [63], and CD-HIT [64], are mainly used to cluster homologous proteins.Table 2
Summary of learning-based methods.
TypeMethod1Ref.YearMetal ion binding2Dataset3ResolutionSequence similarity (tool)4Prediction model5Cross-validationIndependent testMeasurements6Sequence-basedLiu et al.[56]2020Zn, Cu, Fe, Co, Mn, Ca, Mg, Na, K5,340≤3 Å30% (CD-HIT)RF5-fold√SN, SP, ACC, MCCMIonSite[90]2019Zn, Ca, Mg, Mn, Fe, Cu, Fe, Co, Na, K, Cd, Ni7,676N/A30% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCC, AUCMPLs-Pred[91]2019General metal ions1,492N/A30% (CD-HIT)RF10-fold√SN, SP, ACC, MCCSXGBsite[92]2019Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)GBM5-fold√SN, SP, ACC, MCC, AUCWang et al.[50]2019Zn, Cu, Fe, Mn, Ca, Mg, Na, K5,146≤3 Å30% (N/A)SVM, SMO5-fold√SN, SP, ACC, MCCznMachine[51]2019Zn2,043≤3 Å30% (BLASTclust)SVM, NN5-fold√SN, SP, ACC, MCC, PRE, AUCSSWPNN[93]2019Zn213≤2.5 Å70% (N/A)SVM, NN5-fold√SN, SP, PRE, F1, MCC, ACCZinCaps[94]2019Zn738≤3 ÅN/A (N/A)CN5-fold√SN, SP, ACC, MCC, AUCHaberal and Oğul[65]2018General metal ions2,727N/AN/A (N/A)CNN5-fold√SN, ACC, PRE, F1ZincBinder[87]2018Zn738≤2.5 Å30% (PISCES)SVM5-fold√SN, SP, ACC, MCC, AUCEC-RUS[95]2017Ca, Mg, Mn, Fe, Zn4,421N/A40% (PISCES)WSRC5-fold√SN, SP, ACC, MCC, AUCCao et al.[52]2017Zn, Cu, Fe, Co, Mn, Ca, Mg, K, Na5,340≤3 Å30% (CD-HIT)SVM5-fold√SN, SP, ACC, MCCKumar[96]2017Cu, Ca, Co, Fe, Mg, Mn, Ni, Zn3,922N/A50% (CD-HIT)RF10-fold√SN, SP, ACC, MCCDeepMBS[97]2017General metal ions2,727≤3 ÅN/A (N/A)CNN5-fold√SN, PRE, F1Qiao et al.[98]2017Ca2,239N/A30% (CD-HIT)SVM5-fold√SN, ACC, PRE, MCC, AUCIonCom[99]2016Zn, Cu, Fe, Ca, Mg, Mn, Na, K1,374N/A30% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCCJiang et al.[77]2016Ca1,885≤3 Å25% (N/A)SVM5-fold√SN, SP, ACC, MCCTargetCom[53]2016Cu, Fe, Zn1,373≤3 Å40% (CD-HIT)SVM, AdaBoost5-fold√SN, SP, ACC, MCCOSML[100]2015Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)SVM5-fold√SN, SP, ACC, MCCTargetS[101]2013Ca, Zn, Mg, Mn, Fe4,421N/A40% (PISCES)SVM, AdaBoost5-fold√SN, SP, ACC, MCC, AUCETMB-RBF[102]2013General metal ions55N/A20% (BLASTclust)RBFN10-fold√SN, SP, ACC, MCCZincExplorer[103]2013Zn392≤3 ÅN/A (N/A)SVM5-fold√SN, SP, PRE, MCC, AUPRCHorst et al.[58]2010Ca635≤2.1 Å35% (N/A)LR10-fold√MCC,AUC,AUPRCStructure-basedNguyen et al.[104]2021Mn, Fe, Co, Ni, Cu, Zn9,955≤2.5 Å90% (N/A)RF5-fold×ACCTMP-MIBS[54]2021General metal ions427N/A40% (CD-HIT)RF10-fold√SN, SP, ACC, MCC, AUCZincbindpredict[105]2021ZnN/A≤ 2 Å40% (CD-HIT)RF5-fold√SN, PRE, F1, MCCWang et al.[81]2021Zn, Cu, Fe, Ca, Mg, Mn, Na, K, Co5,340≤3 Å30% (N/A)MLP,SVM5-fold√SN, SP, ACC, MCCDELIA[80]2020Ca, Mn, Mg3,966N/A30% (CD-HIT)CNN5-fold√SN, PRE, MCC, AUCHu et al.[57]2020Zn, Cu, Fe, Co, Mn, Ca, Mg, Na, K5,340≤3 Å30% (CD-HIT)GBM5-fold√SN, SP, FPR, ACC, MCCMetalExplorer[79]2017Ca, Co, Cu, Fe, Ni, Mg, Mn, Zn3,192≤2.5 Å30% (CD-HIT)RF5-fold√SN, FPR, PRE, AUC, AURPCFINDSITE-metal[55]2011Ca, Co, Cu, Fe, Mg, Mn, Ni, Zn860N/A35% (PISCES)SVM2-fold√ACC, SPC, PPVZinc identifier[78]2011Zn1,103≤2.5 ÅN/A (N/A)RF5-fold√SN, PRE, SP, FPR, AUC, AUPRC1The name of each method is provided in either the publication or the last name of its first author. 2General metal ions mean that the related predictor does not differentiate the types of metal ion binding. Otherwise, we list the specific types of metal-binding in detail. 3The number represents the size of the benchmark dataset. 4The value reveals the protein similarity threshold in the benchmark dataset. The content in the blanket indicates the tool that is used for clustering proteins. 5SMO: sequential minimal optimization; SVM: support vector machine; WSRC: weighted sparse representation based classifier; NN: neural network; CN: capsule network; CNN: convolutional neural networks; RF: random forest; GBM: gradient boosting machine; RBFN: radial basis function networks; LR: logistic regression; MLP: multilayer perceptron. 6SN: sensitivity/recall; SP: specificity; ACC: accuracy; MCC: Matthew’s correlation coefficient; PRE: precision; F1: F1-score; AUC: area under the ROC curve; AUPRC: area under the precision recall curve; FPR: false positive rate (FPR=1-SP).


## 3.2. The Validation and Evaluation Metrics

### 3.2.1. Cross-Validation and Independent Test

To construct a predictor with high accuracy and decent generalization ability, it is necessary to avoid potential overfitting. In practice, cross-validation and independent test are two popular ways (Table2) to evaluate the proposed models [31]. Specifically, k-fold cross-validation is usually adopted on the training dataset when building the prediction model and optimizing the related parameters [61]. First, the training dataset is equally divided into k parts. The division can be done at residue level or protein level. Next, k-1 subsets are used to train the model, and the last one subset is used for testing. The procedure repeats k times until every subset is been predicted. The performance of the model is usually evaluated by averaging the results of the k repeats.


### 3.2.2. Performance Measures

According to Table2, the measures that used to evaluate the performance of the predictors can be divided into binary value-based and propensity score-based ones. The former needs preset thresholds to compute the number of putative binding residues and nonbinding residues. These measures include sensitivity (SN)/recall/true positive rate (TPR), specificity (SP), false positive rate (FPR, FPR=1-SP) precision (PRE), accuracy (ACC), F1-score (F1), and Matthew’s correlation coefficient (MCC). They are defined as follows:
(1)SN=TPR=TPTP+FN,SP=1−FPR=TNTN+FP,PRE=TPTP+FP,F1=2TP2TP+FP+FN,ACC=TP+TNTP+FN+TN+FP,MCC=TP×TN−FN×FPTP+FN×TP+FP×TN+FP×TN+FN,where TP (true positive) indicates the number of correctly recognized metal-binding residues, FP (false positive) means the number of non-metal-binding residues that are incorrectly predicted as metal-binding residues, TN (true negative) stands for the number of correctly predicted non-metal-binding residues, and FN (false negative) is the number of metal-binding residues that are incorrectly predicted as non-metal-binding residues.The prediction of metal-binding residues is a typical imbalanced classification problem. That is, the number of metal-binding residues is much less than that of the non-metal-binding ones. Therefore, F1-score and MCC are regarded as key criteria since they are featured by assessing the prediction performance for both metal-binding and non-metal ion binding residues.The propensity score-based measures include receiver operating characteristic curve (ROC curve) and precision-recall curve (PR curve). The ROC curve draws the TPR (true positive rate) against the FPR (false positive rate) at various thresholds. The AUC computes the area under the ROC curve and can be used to quantify the ROC curve. The PR curve plots PRE values on they-axis and recalls values on the x-axis, and the AUPRC estimates the area under the PR curve.


## 3.2.1. Cross-Validation and Independent Test

To construct a predictor with high accuracy and decent generalization ability, it is necessary to avoid potential overfitting. In practice, cross-validation and independent test are two popular ways (Table2) to evaluate the proposed models [31]. Specifically, k-fold cross-validation is usually adopted on the training dataset when building the prediction model and optimizing the related parameters [61]. First, the training dataset is equally divided into k parts. The division can be done at residue level or protein level. Next, k-1 subsets are used to train the model, and the last one subset is used for testing. The procedure repeats k times until every subset is been predicted. The performance of the model is usually evaluated by averaging the results of the k repeats.


## 3.2.2. Performance Measures

According to Table2, the measures that used to evaluate the performance of the predictors can be divided into binary value-based and propensity score-based ones. The former needs preset thresholds to compute the number of putative binding residues and nonbinding residues. These measures include sensitivity (SN)/recall/true positive rate (TPR), specificity (SP), false positive rate (FPR, FPR=1-SP) precision (PRE), accuracy (ACC), F1-score (F1), and Matthew’s correlation coefficient (MCC). They are defined as follows:
(1)SN=TPR=TPTP+FN,SP=1−FPR=TNTN+FP,PRE=TPTP+FP,F1=2TP2TP+FP+FN,ACC=TP+TNTP+FN+TN+FP,MCC=TP×TN−FN×FPTP+FN×TP+FP×TN+FP×TN+FN,where TP (true positive) indicates the number of correctly recognized metal-binding residues, FP (false positive) means the number of non-metal-binding residues that are incorrectly predicted as metal-binding residues, TN (true negative) stands for the number of correctly predicted non-metal-binding residues, and FN (false negative) is the number of metal-binding residues that are incorrectly predicted as non-metal-binding residues.The prediction of metal-binding residues is a typical imbalanced classification problem. That is, the number of metal-binding residues is much less than that of the non-metal-binding ones. Therefore, F1-score and MCC are regarded as key criteria since they are featured by assessing the prediction performance for both metal-binding and non-metal ion binding residues.The propensity score-based measures include receiver operating characteristic curve (ROC curve) and precision-recall curve (PR curve). The ROC curve draws the TPR (true positive rate) against the FPR (false positive rate) at various thresholds. The AUC computes the area under the ROC curve and can be used to quantify the ROC curve. The PR curve plots PRE values on they-axis and recalls values on the x-axis, and the AUPRC estimates the area under the PR curve.


## 3.3. Learning-Based Methods

Learning-based methods treat the recognition of metal-binding residues as a typical pattern recognition problem. Specifically, the metal-binding residues and nonbinding ones are encoded by using mathematical descriptors, i.e. features. According to the information that used to compute the features, the learning-based methods can be further categorized into sequence-based and structure-based methods. The former only needs simple protein sequences to extract features when encoding the binding residues. These features include sequence directly derived, evolutionary profile-based, and putative structure-based features, while the latter uses both sequence and native structure data to mathematically describe a binding residue. We make a comprehensive literature search and collect 23 sequence- and 9 structure-based methods that were published after the year 2010.


### 3.3.1. Feature Construction

(1) Sequence Directly Derived Features. We define sequence directly derived features as the ones that are computed from protein primary sequences without using any other information. In Figure 3, 14 out of 32 considered methods consider amino acid composition [50], which quantifies the relative difference in abundance of a given amino acid type [65, 66]. Amino acid pairs, or dipeptides, are based on the observation that amino acid pairs show different propensities in protein structure and function. For instance, pairs of lysine are found present in close spatial vicinity [67]. Moreover, the concept of k-spaced amino acid pairs is introduced in [68]. It calculates the amino acid pairs with k spaces between two residues. Our survey also shows that the majority of studies use physicochemical properties to describe the local environment of the metal-binding residues. The basic physicochemical environment of a metal-protein binding interface is reflected by the specific roles the metal plays in biostructural chemistry and protein function. These properties are crucial since they underpin many of the functional roles of metal ions. These properties include aliphatic [69], sulphur [70], aromatic [71], hydrophobic [72], charge [73], polar [74], positive [73], acidic [75], and hydroxylic [76]. The position-related features mainly consider the influence of the specifically located residues, such as autocross covariance [77] and sequence length [78, 79].Figure 3
Summary of the feature construction and selection for learning-based methods. The light green cells indicate sequence directly derived features. The light blue cells stand for profile-based features. The light red cells mean putative structure-based features. The light grey cells are native structure-based features.1In the amino acid pair column, the cells without annotations indicate original amino acid pairs; the cells annotated using “k-spaced” means k-spaced amino acid pairs. 2PWM: position weight matrix; PSSM: position specific scoring matrix; PMS: position matrix scoring; PAM: point accepted mutation; EMS: evolutionary matrix scoring. 3TS: topology structure; DM: distance matrix; GTN: graph theoretic network; RRCG: residue-residue contact graphs. 4FFS: forward feature selection; EB: experience-based; BA: Boruta algorithm; mRMR: minimum-redundancy maximum-relevancy; MDGI: mean decrease Gini index [50–51, 53–58, 65, 77–81, 87, 90–105].(2) Evolutionary Profile-Based Features. Recent studies [54, 56, 57, 80, 81] pointed out that functional or structural important residues tend to show higher evolutionary conservation. The conserved residues are usually involved in enzyme activity, ligand binding, or protein structural stability [82]. The conserved residues and regions can be identified by multiple sequence alignment [83]. These multiple sequence alignments, also named conservation profiles, include aligning families of homologous sequences and having knowledge of their evolutionary relationships [84]. For an unknown protein, although its accurate function is not available, it is expected that we can use its homologous proteins to speculate the function since they share the similar evolutionary profile [85]. Many studies use position-specific scoring matrix (PSSM), which is computed from PSI-BLAST [62], to quantify the evolutionary conservation. PSSM scores the substitution probability of each residue in the protein being substituted by other types of amino acids. Liu et al. [56] and Hu et al. [57] set different weights according to the positions of considered residues within the window and construct position weight matrix (PWM). Wang et al. proposed a customized position matrix scoring (PMS) algorithm, which uses known sequence patterns to describe the composition of amino acids at different positions [50]. Haberal and Oğul introduced a point accepted mutation (PAM) scoring matrix, which measures the rate at which point mutations that substitute one residue for another during evolution [65]. Jiang et al. adopted evolutionary matrix scoring (EMS) algorithm to extract the position conservation of amino acid residues from segments with low dimension feature parameters [77].(3) Putative Structure-Based Features. For an unknown protein, although the accurate function is not available, it is expected that we can use its homologous proteins or template structures to speculate the structure. The secondary structure mainly involves α-helix, β-sheet, and coil, which are fundamental elements of protein tertiary structure [86]. Natively disordered or unstructured regions are proved to be associated with molecular assembly, protein translation, modification, and molecular recognition [78, 79, 87]. Previous studies [79, 87] indicate that disordered regions are strongly correlated with local solvent accessibility areas. Figure 3 reveals that 16 methods introduce secondary structure features and 3 approaches use disorder features, respectively. The secondary structure can be obtained from the primary sequence by using PSIPRED [88]. Putative intrinsic disorder data can be computed by using DISOPRED [89].(4) Structure-Based Features. The structure-based features include descriptors that are computed from protein 3D structure. These features include solvent exposure, B-factor, spatial cluster properties, and native secondary structure. Compared with the abovementioned putative structure-based features, the native structure-based features are more accurate since they are directly computed by using residue coordinate data. Besides that, a residue contact network is also considered by some literature. In [79], two residues are defined as being in contact if the distance of their Cα atoms is less than a predefined cutoff distance of 6.5 Å. These features include clustering coefficient, degree, density, distance, topology structure, and graph theoretic network [55, 79, 80].


### 3.3.2. Sliding Window Optimization and Feature Selection

As shown in Figure3, many methods adopt a sliding window scheme when they construct different types of features. It is because residues in proteins are influenced by adjacent residues. Besides that, binding residues tend to cluster together. If a central residue is a native-binding residue, its adjacent residues usually have a relatively higher chance to bind the same ligands. Usually, the residues with a long distance away have a lower impact on the central residues when compared with the residues with short distance. Figure 3 summarizes that 19 out of 23 methods use the sliding window scheme. The size of the shortest window is 3 [105], while the size for the longest one is 25 [94]. Some studies [50, 52, 53, 56, 57, 81, 90, 105] use more than one type of window because they consider different types of metal-binding residues. A long window means the introduction of more features.However, a bigger number of features do not absolutely mean a better prediction performance [106, 107]. The existence of potential “bad” features may interfere with the classifiers and cause unpredictable consequences [108]. The so-called “bad” features include irrelevant and redundant ones. To avoid their terrible influences, it is necessary to perform feature selection before training the model [109]. Figure 3 reveals that 6 out of 32 methods adopt feature selection before training the model. These feature selection approaches include forward feature selection [79, 102], experience-based [104], Boruta algorithm [57, 81], minimum-redundancy maximum-relevancy [79], and mean decrease Gini index [78].


### 3.3.3. Prediction Algorithms

Learning-based methods use machine-learning or deep-learning-based algorithms to train the model and perform predictions [110]. As shown in Table 2, a variety of algorithms are introduced for solving the problem of correctly recognizing metal-binding residues. Support vector machine (SVM) is a popular machine learning algorithm in bioinformatical research. It is aimed at finding a hyperplane or decision boundary that can segregate a high-dimensional space [111]. Particularly, it uses kernel functions to reduce computation time to avoid strapping into dimension disaster [112]. Sequential minimal optimization (SMO) is an algorithm that is specially used for training support vector machines [113]. The procedure of training large data by SVM usually leads to a complex quadratic programming optimization problem [114]. SMO breaks large programming optimization problems into small ones, which endows SVM a good generalization on large data [113]. The idea of a neural network (NN) comes from the work system of neurons in the biological brain [115]. It learns the correlations between inputs and outputs, making generalizations and build models [116]. The NN algorithm assigns and adjusts different weights for neurons and edges as learning proceeds. The radial basis function network (RBFN) is a variant of the original NN [117]. It adopts radial basis functions as activation functions, which can be used for accelerating learning speed due to their universal approximation [118]. The multilayer perceptron (MLP) algorithm is an improved back propagation NN [119]. It mainly includes three procedures: forward propagation, error evaluation, and error backpropagation [120]. The MLP is featured by its strong generalization and fault tolerance [121]. Therefore, it is proved to be an efficient classification algorithm. The logistic regression (LR) adopts a logistic function to model the probability of an unknown sample being a certain class [122].Our survey also reveals that the ensemble algorithms are favored by eight studies. The random forest (RF) aggregates the predictions of all the decision trees and performs decisions by most trees [123]. RF can be used for classification, regression, and optimization problems [124]. Adaptive boosting (AdaBoost) is aimed at combining weak learners with strong ones [125, 126]. The key point of AdaBoost is to ensure the diversities of individual learners, which makes it a good generalization ability [90, 127]. The gradient boosting machine (GBM) is another popular ensemble algorithm. During the iterative process, GBM dynamically increases the weight of wrong recognitions and reduces that of the correct ones [128–131]. It should be noted that GBM focuses on the sample residual of the previous iteration instead of the sample itself [132].Besides machine-learning algorithms, recent studies also use deep-learning methods in this research field. The convolutional neural network (CNN) is one of the most prevalent algorithms that is widely used in bioinformatics [133]. The CNN consists of three main layers, which are the convolutional layer, pooling layer, and fully connected layer [134–136]. Although the CNN is proved to be powerful in dealing with a variety of problems, it performs badly when facing samples with different sizes and orientations [137, 138]. To overcome this shortcoming, the capsule network (CN) is proposed to estimate features of objects by incorporating dynamic routing algorithms [139, 140]. Our review finds two studies use CNN [65, 97] and one uses CN [94].


## 3.3.1. Feature Construction

(1) Sequence Directly Derived Features. We define sequence directly derived features as the ones that are computed from protein primary sequences without using any other information. In Figure 3, 14 out of 32 considered methods consider amino acid composition [50], which quantifies the relative difference in abundance of a given amino acid type [65, 66]. Amino acid pairs, or dipeptides, are based on the observation that amino acid pairs show different propensities in protein structure and function. For instance, pairs of lysine are found present in close spatial vicinity [67]. Moreover, the concept of k-spaced amino acid pairs is introduced in [68]. It calculates the amino acid pairs with k spaces between two residues. Our survey also shows that the majority of studies use physicochemical properties to describe the local environment of the metal-binding residues. The basic physicochemical environment of a metal-protein binding interface is reflected by the specific roles the metal plays in biostructural chemistry and protein function. These properties are crucial since they underpin many of the functional roles of metal ions. These properties include aliphatic [69], sulphur [70], aromatic [71], hydrophobic [72], charge [73], polar [74], positive [73], acidic [75], and hydroxylic [76]. The position-related features mainly consider the influence of the specifically located residues, such as autocross covariance [77] and sequence length [78, 79].Figure 3
Summary of the feature construction and selection for learning-based methods. The light green cells indicate sequence directly derived features. The light blue cells stand for profile-based features. The light red cells mean putative structure-based features. The light grey cells are native structure-based features.1In the amino acid pair column, the cells without annotations indicate original amino acid pairs; the cells annotated using “k-spaced” means k-spaced amino acid pairs. 2PWM: position weight matrix; PSSM: position specific scoring matrix; PMS: position matrix scoring; PAM: point accepted mutation; EMS: evolutionary matrix scoring. 3TS: topology structure; DM: distance matrix; GTN: graph theoretic network; RRCG: residue-residue contact graphs. 4FFS: forward feature selection; EB: experience-based; BA: Boruta algorithm; mRMR: minimum-redundancy maximum-relevancy; MDGI: mean decrease Gini index [50–51, 53–58, 65, 77–81, 87, 90–105].(2) Evolutionary Profile-Based Features. Recent studies [54, 56, 57, 80, 81] pointed out that functional or structural important residues tend to show higher evolutionary conservation. The conserved residues are usually involved in enzyme activity, ligand binding, or protein structural stability [82]. The conserved residues and regions can be identified by multiple sequence alignment [83]. These multiple sequence alignments, also named conservation profiles, include aligning families of homologous sequences and having knowledge of their evolutionary relationships [84]. For an unknown protein, although its accurate function is not available, it is expected that we can use its homologous proteins to speculate the function since they share the similar evolutionary profile [85]. Many studies use position-specific scoring matrix (PSSM), which is computed from PSI-BLAST [62], to quantify the evolutionary conservation. PSSM scores the substitution probability of each residue in the protein being substituted by other types of amino acids. Liu et al. [56] and Hu et al. [57] set different weights according to the positions of considered residues within the window and construct position weight matrix (PWM). Wang et al. proposed a customized position matrix scoring (PMS) algorithm, which uses known sequence patterns to describe the composition of amino acids at different positions [50]. Haberal and Oğul introduced a point accepted mutation (PAM) scoring matrix, which measures the rate at which point mutations that substitute one residue for another during evolution [65]. Jiang et al. adopted evolutionary matrix scoring (EMS) algorithm to extract the position conservation of amino acid residues from segments with low dimension feature parameters [77].(3) Putative Structure-Based Features. For an unknown protein, although the accurate function is not available, it is expected that we can use its homologous proteins or template structures to speculate the structure. The secondary structure mainly involves α-helix, β-sheet, and coil, which are fundamental elements of protein tertiary structure [86]. Natively disordered or unstructured regions are proved to be associated with molecular assembly, protein translation, modification, and molecular recognition [78, 79, 87]. Previous studies [79, 87] indicate that disordered regions are strongly correlated with local solvent accessibility areas. Figure 3 reveals that 16 methods introduce secondary structure features and 3 approaches use disorder features, respectively. The secondary structure can be obtained from the primary sequence by using PSIPRED [88]. Putative intrinsic disorder data can be computed by using DISOPRED [89].(4) Structure-Based Features. The structure-based features include descriptors that are computed from protein 3D structure. These features include solvent exposure, B-factor, spatial cluster properties, and native secondary structure. Compared with the abovementioned putative structure-based features, the native structure-based features are more accurate since they are directly computed by using residue coordinate data. Besides that, a residue contact network is also considered by some literature. In [79], two residues are defined as being in contact if the distance of their Cα atoms is less than a predefined cutoff distance of 6.5 Å. These features include clustering coefficient, degree, density, distance, topology structure, and graph theoretic network [55, 79, 80].


## 3.3.2. Sliding Window Optimization and Feature Selection

As shown in Figure3, many methods adopt a sliding window scheme when they construct different types of features. It is because residues in proteins are influenced by adjacent residues. Besides that, binding residues tend to cluster together. If a central residue is a native-binding residue, its adjacent residues usually have a relatively higher chance to bind the same ligands. Usually, the residues with a long distance away have a lower impact on the central residues when compared with the residues with short distance. Figure 3 summarizes that 19 out of 23 methods use the sliding window scheme. The size of the shortest window is 3 [105], while the size for the longest one is 25 [94]. Some studies [50, 52, 53, 56, 57, 81, 90, 105] use more than one type of window because they consider different types of metal-binding residues. A long window means the introduction of more features.However, a bigger number of features do not absolutely mean a better prediction performance [106, 107]. The existence of potential “bad” features may interfere with the classifiers and cause unpredictable consequences [108]. The so-called “bad” features include irrelevant and redundant ones. To avoid their terrible influences, it is necessary to perform feature selection before training the model [109]. Figure 3 reveals that 6 out of 32 methods adopt feature selection before training the model. These feature selection approaches include forward feature selection [79, 102], experience-based [104], Boruta algorithm [57, 81], minimum-redundancy maximum-relevancy [79], and mean decrease Gini index [78].


## 3.3.3. Prediction Algorithms

Learning-based methods use machine-learning or deep-learning-based algorithms to train the model and perform predictions [110]. As shown in Table 2, a variety of algorithms are introduced for solving the problem of correctly recognizing metal-binding residues. Support vector machine (SVM) is a popular machine learning algorithm in bioinformatical research. It is aimed at finding a hyperplane or decision boundary that can segregate a high-dimensional space [111]. Particularly, it uses kernel functions to reduce computation time to avoid strapping into dimension disaster [112]. Sequential minimal optimization (SMO) is an algorithm that is specially used for training support vector machines [113]. The procedure of training large data by SVM usually leads to a complex quadratic programming optimization problem [114]. SMO breaks large programming optimization problems into small ones, which endows SVM a good generalization on large data [113]. The idea of a neural network (NN) comes from the work system of neurons in the biological brain [115]. It learns the correlations between inputs and outputs, making generalizations and build models [116]. The NN algorithm assigns and adjusts different weights for neurons and edges as learning proceeds. The radial basis function network (RBFN) is a variant of the original NN [117]. It adopts radial basis functions as activation functions, which can be used for accelerating learning speed due to their universal approximation [118]. The multilayer perceptron (MLP) algorithm is an improved back propagation NN [119]. It mainly includes three procedures: forward propagation, error evaluation, and error backpropagation [120]. The MLP is featured by its strong generalization and fault tolerance [121]. Therefore, it is proved to be an efficient classification algorithm. The logistic regression (LR) adopts a logistic function to model the probability of an unknown sample being a certain class [122].Our survey also reveals that the ensemble algorithms are favored by eight studies. The random forest (RF) aggregates the predictions of all the decision trees and performs decisions by most trees [123]. RF can be used for classification, regression, and optimization problems [124]. Adaptive boosting (AdaBoost) is aimed at combining weak learners with strong ones [125, 126]. The key point of AdaBoost is to ensure the diversities of individual learners, which makes it a good generalization ability [90, 127]. The gradient boosting machine (GBM) is another popular ensemble algorithm. During the iterative process, GBM dynamically increases the weight of wrong recognitions and reduces that of the correct ones [128–131]. It should be noted that GBM focuses on the sample residual of the previous iteration instead of the sample itself [132].Besides machine-learning algorithms, recent studies also use deep-learning methods in this research field. The convolutional neural network (CNN) is one of the most prevalent algorithms that is widely used in bioinformatics [133]. The CNN consists of three main layers, which are the convolutional layer, pooling layer, and fully connected layer [134–136]. Although the CNN is proved to be powerful in dealing with a variety of problems, it performs badly when facing samples with different sizes and orientations [137, 138]. To overcome this shortcoming, the capsule network (CN) is proposed to estimate features of objects by incorporating dynamic routing algorithms [139, 140]. Our review finds two studies use CNN [65, 97] and one uses CN [94].


## 3.4. Docking-Based Methods

The investigation on the protein-metal complex helps biologists to understand the mechanism of protein-metal interactions. Protein-ligand docking approaches are always based on molecular structure and are used to explore biomolecular interactions and mechanisms [141]. It can be adopted to predict binding conformation as well as the appropriate target binding residues [142, 143].As shown in Figure2, the docking-based methods mainly include three steps: searching algorithm, scoring function, and docking assessment [141]. The searching algorithm focuses on creating an optimum number of configurations that properly include the determined binding modes [144]. To reduce computation time, it is necessary to make a balance between the computational expense and the searching space. The scoring function includes a series of mathematical functions that quantify the strength of binding affinity [145]. The energy-based scoring functions are always introduced to score the potential interactions between the protein and the corresponding ligands [141]. The frequently used functions include empirical-based, knowledge-based, and consensus-based ones. Finally, the putative docking can be evaluated by using docking accuracy and the correlation between putative and native docking scores [145]. Figure 4 illustrates the structure of a calmodulin (PDB: 4HEX) that is secreted by Escherichia coli in Mus musculus [146]. Calmodulin is one of the most prevalent EF-hand calcium sensor proteins in eukaryotic cells [147]. It is a highly conserved and soluble protein, which activates enzymes and regulates many cellular functions. 4HEX has three Ca2+-binding and two Zn2+-binding sites. Ca2+-binding causes a change in calmodulin conformation opening both globular domains and exposing hydrophobic surfaces that form binding sites for the target enzymes. Figure 4 shows that these three Ca2+ are in the pockets. The binding pockets are half-closed and buried, which substantially limits the capability of Ca2+ to escape. Two Zn2+-binding sites are surrounded by a shell of hydrophilic groups that are embedded into a larger shell of hydrophobic groups. The amino acid side chains providing ligands to Zn2+ in these structures often form hydrogen bonds with other residues [147].Figure 4
Ribbon and surface model of X-ray structure of Ca2+- and Zn2+-bound calmodulin (PDB: 4HEX) in Mus musculus. Red sphere represents bound zinc ion; green one indicates calcium ion; the spatial adjacent residues participating its coordination are shown by the stick model.
(a)(b)(c)In [148], He et al. proposed a docking-based predictor named mFASD. It first explored the local biochemical environment of potential functional atoms and then measured the distances between the atoms and bound metal. mFASD also claimed that it can differentiate different types of metal-binding sites. Zhou et al. improved the FEATURE-based calcium model and used the grid scan algorithm to recognize binding sites [149]. GaudiMM [150] adopted a multiobjective genetic algorithm to search metal-binding sites in biological scaffolds. BioMetAll focused on the conformation of the potential metal-binding site, associated with the geometric organization of the protein backbone [151]. It was also proved to have good performance on the applications including the modulation and mutation of the metal-binding residues. Table 3 summarizes the key notes of the abovementioned 4 docking-based methods.Table 3
Summary of docking-based, template-based, and meta-based methods.
TypeMethodYearNotesDocking-basedmFASD [148]2015Capture the characteristics of metal-binding sites and discriminate most types of these sitesZhou et al. [149]2015Use a FEATURE-based calcium model and convert high scoring regions into specific site predictionsGaudiMM [150]2019Find poses that satisfy metal-derived geometrical rules and use post optimizationsBioMetAll [151]2020Predict metal-binding sites with particular motifs, determine transient sites in structures, and predict potential mutations to generate convenient sitesTemplate-basedDeng et al. [152]2006Use a graph theory algorithm to find oxygen clusters of the protein (high potential for calcium binding)Goyal et al. [153]2008Describe generation of 3D-structural motifs for metal-binding sites from the known metalloproteinsLevy et al. [154]2009Analyze whether structural models based on remote homology are effective in predicting 3D metal binding sitesFunFOLD [155]2011Use an automated method for ligand clustering and identification of binding residuesFunFOLDQA [156]2012Use a fully automated agglomerative clustering approach for both ligand identification and residue selectionFunFOLD2 [157]2013Propose a method that include protein-ligand binding prediction and quality assessment protocolMeta-basedLi et al. [158]2017Integrate the results of ZincExplorer [103], zincFinder [159], and zincPred [160]IBayes_Zinc [161]2019Adopt Bayesian method and combine the predictions from ZincExplorer [103], zincFinder [159], and zincPred [160]


## 3.5. Template-Based Methods

It is well known that protein structure determines function, and similar interface conformation indicates similar bound regions [162]. The template-based methods are based on the abovementioned hypothesis. Therefore, the most important thing for template-based methods is to find and validate proper structural templates. The fold recognition algorithms, which quantify the best matches from candidate templates, are commonly used to select the optimal template structures [163]. Next, the selected templates are used to map onto the target protein given the alignments with the template structures [164].As shown in Table3, Goyal and Mande analyzed the metal-binding sites by using structure templates and designing 3D motifs for several types of metal-binding interactions [153]. In [154], the authors analyzed whether structural models based on remote homology are effective in recognizing structural metal-binding residues based on simple protein primary sequences. Deng et al. applied a graph theory algorithm to identify, predict, and analyze calcium-binding residues [152]. However, it should be noted that this strategy produces good prediction performance when a decent complex is available as a template. If the template structure information is not available, this strategy might have poor predictions [164]. The FunFOLD was an automatic method that uses protein structure superposition of distantly related templates to a modelled protein for the clustering of ligands and prediction of metal binding residues [155]. The FunFOLDQA [156] approach determined the reliability of our FunFOLD [155] by assigning the quality assessment scores. FunFOLD2 was a web server that integrated cutting edge function and putative 3D structures to identify metal-binding residues [157].


## 3.6. Meta-Based Methods

The meta-based methods use a meta-learning strategy from fewer samples than traditional machine learning models. Since meta-based methods can only use limited data, they must ensure that the data is featured with high accuracy. As a result, a meta-based approach always directly combines the predictions of other methods. It uses weights or voting strategy on the available propensity scores or binary values. Thus, the meta-based method promises a robust accurate prediction on the metal-binding residues. In [158], Li et al. collected the predictions from ZincExplorer [103], ZincFinder [159], and ZincPred [160] (Table 3). Then, they built a linear regression model and optimized corresponding parameters on the training dataset. They claimed that the meta-model, which was named meta-zincPrediction, improves the AUPRC by about 2%~8%. IBayes_Zinc [161] was another meta-based predictor for the identification of zinc-binding residues (Table 3). It firstly computed the predictions of zinc-binding probabilities from ZincExplorer [103], ZincFinder [159], and ZincPred [160]. Next, IBayes_Zinc processed the missing attribute values and adopts Bayesian theory [165] to construct a meta-based model. The performance on the independent dataset proved that the MCC value of IBayes_Zinc was about 5~13% higher than the considered three predictors.


## 3.7. Prediction Results

This review surveys 44 computation-based methods. It is necessary to make a consensus comparison for these methods. However, since there is no standard benchmark dataset and some methods are not currently available, we use two datasets that are used by some methods to perform the evaluations. The first dataset is compiled by Yu et al. in [101], which includes five types of metal ions binding annotations. The second dataset is obtained from [52], consisting of ten types of metal ions binding annotations.Figure5 illustrates the predictive performance on two benchmark test datasets. Details are provided in Table S1 and Table S2 in Supplementary Materials, respectively; the corresponding results are sourced from [56, 57, 81, 92]. We notice that the predictors show relatively big differences in recognizing various types of metal-binding residues. On Yu et al.’s dataset [101], TargetS shows the best results in predicting Ca2+-, Zn2+-, and Mn2+-binding residues; EC-RUS [95] performs best in recognizing Mg2+-binding residues; OSML [100] achieves the highest MCC on Fe3+-binding predictions. Besides that, Figure 5(a) indicates that all five methods show a decent performance on recognizing Fe3+-binding residues (MCC values close or higher than 0.4), compared with MCC close or less than 0.2 on Ca2+ binding residues. Figure 5(b) draws the bars of AUC values for SXGBsite [92], EC-RUS [95], and TargetS [101], respectively. These three predictors all achieve high AUC scores (close or higher than 0.9) on Zn2+- and Fe3+-binding residues. Figure 5(c) summarizes the results of ten metal ions binding residues on Cao et al.’s dataset [52]. Among these predictors, Liu et al. [56] performs the best on Zn2+, Fe3+, and Cu2+, compared to [81], Wang et al. shows best on Zn2+ and Cu2+, and Hu et al. [57] achieves the highest on Fe2+. Interestingly, the binding residues associated with relatively inactive metal (Zn2+, Fe3+, and Cu2+) ions show relatively better results compared to that of the active metal ions (Na+ and K+). Particularly, four methods all give better results on Fe3+-binding residues than that on Fe2+-binding residues, which keep consistent with our observations as mentioned above.Figure 5
Comparative assessment of several predictors on two benchmark dataset. (a) and (c) indicate the MCC bar charts for considered methods on different metal ion binding residues on Yu et al.’s and Cao et al.’s testing datasets, respectively. (b) draws the AUC values of three predictors on corresponding metal ion binding residues.
(a)(b)(c)


## 3.8. Publicly Available Tools

The publicly available standalone software or web server that implements the proposed approach provides convenience for biologists and researchers [79, 105, 122]. These tools help the community to repeat the results and build a platform for easy understanding and improvement. Table 4 summarizes the public availability of implementations for the considered methods. These 28 predictors are implemented as standalone software or web servers. Among these predictive tools, 16 (or 57%) of them are currently publicly available. Standalone software requires the biologists to build the same running environment. By contrast, the web server provides the most convenient since the users only need to submit their queries via the browser, and the server helps to do the computations. Three methods, namely, znMachine [51], ZincBinder [103], and FunFOLD [155], provide both web server and standalone software. TMP_MIBS [54] is designed to predict general metal-binding residues and deployed using the Python language. DELIA [80] requires PDB-formatted 3D coordinates input and produces both binary prediction and putative probability of a residue being potential specific metal-binding. BioMetAll uses a docking-based strategy to scan specific motifs, putative mutations, and binding residues. Another available docking-based method is mFASD, which distinguishes different types of metal-binding sites according to the interaction distances. MPLs-Pred [91], SXGBsite [92], MIonSite [90], OSML [100], EC-RUS [95], and TargetS [101] are all sequence-based predictive tools, which accepts FASTA-formatted input and produced the results of putative metal-binding residues. ZinCaps [94], SSWPNN [93], and ZincBinder [103] are specially designed for the identification of zinc-binding residues. FunFOLD [155], FunFOLDQA [156], and FunFOLD2 [157] are a series of template-based methods.Table 4
A breakdown of predictive tools of metal-binding residues.
MethodYearPlatform1Web linkAvailability2TMP-MIBS [54]2021SShttps://github.com/QuJing785464/TMP_MIBSYesWang et al. [50]2021WShttp://39.104.77.103:8081/lsb/HomePage/HomePage.htmlNoZincbindpredict [105]2021WShttps://zincbind.bioinf.org.uk/predict/NoDELIA [80]2020WShttp://www.csbio.sjtu.edu.cn/bioinf/delia/YesBioMetAll [151]2020SShttps://github.com/insilichem/biometallYesMPLs-Pred [91]2019WShttp://icdtools.nenu.edu.cn/YesSXGBsite [92]2019SShttps://github.com/Lightness7/SXGBsiteYesMIonSite [90]2019SShttps://github.com/LiangQiaoGu/MIonSite.gitYesznMachine [51]2019WS&SShttp://bioinformatics.fzu.edu.cn/znMachine.htmlNoZinCaps [94]2019SShttps://github.com/clemEssien/ActiveSitePredictionYesEC-RUS [95]2017SShttps://github.com/6gbluewind/protein_ligand_binding_siteYesMetalExplorer [79]2017WShttp://metalexplorer.erc.monash.edu.au/NoCao et al. [52]2017WShttp://60.31.198.140:8081/metal/HomePage/HomePage.htmlNoZincBinder [103]2017WS&SShttp://proteininformatics.org/mkumar/znbinder/YesSSWPNN [93]2017SShttp://net.jitsec.cn:88/UploadedImages/SSWPNN.rarYesJiang et al. [77]2016WShttp://202.207.29.245/NoTargetCom [53]2016SShttp://dase.ecnu.edu.cn/qwdong/TargetCom/TargetCom_standalone.tar.gzNoOSML [100]2015WShttp://www.csbio.sjtu.edu.cn/OSML/YesmFASD [148]2015SShttp://staff.ustc.edu.cn/liangzhi/mfasd/YesFunFOLD2 [157]2013WShttp://www.reading.ac.uk/bioinf/FunFOLD/FunFOLD_form_2_0.htmlYesZincExplorer [103]2013WShttp://protein.cau.edu.cn/ZincExplorerNoTargetS [101]2013WShttp://www.csbio.sjtu.edu.cn/TargetS/YesFunFOLDQA [156]2012SShttp://www.reading.ac.uk/bioinf/downloads/YesZincidentifier [78]2012WShttp://protein.cau.edu.cn/zincidentifier/NoFINDSITE-metal [55]2011WShttp://cssb.biology.gatech.edu/findsite-metal/NoFunFOLD [155]2011WS&SShttp://www.reading.ac.uk/bioinf/FunFOLD/YesGoyal et al. [153]2008WShttp://sunserver.cdfd.org.in:8080/protease/PAR_3D/index.htmlNoDeng et al. [152]2006SShttp://chemistry.gsu.edu/faculty/Yang/GG.htmNo1WS: web server; SS: standalone software. 2The availability was estimated on Dec 1st, 10th, and 20th of 2021, respectively.


## 4. Conclusions and Future Perspectives

This review summarizes the public database of metal ions binding interactions, discusses the architectures of computation-based methods for identifying binding residues, and comparatively evaluates four types of methods. Based on the observations made in this work, we propose a few recommendations for future research in this field:First, the researchers should maintain and update the database regularly. This will significantly improve effectiveness and completeness for these databases and provide convenience for the computation-based methods, which depend on the accurate internal database. We expect a high-quality metal ion binding-related database with an advanced searching engine, high-speed download service, complete annotation information, etc. Particularly, a decent database should be designed to open for easy expanding and improvement. Second, standard benchmark datasets that related to general or ligand-specific metal-binding residues should be periodically compiled and made available. This will ensure consistent evaluation and comparative analysis of the performance of the existing and novel methods. Third, these predictors are expected to use delicate architectures and powerful algorithms. Since the differences between different types are quite small, the novel predictors shall not only correctly identify metal-binding residues but also distinguish different types of metal ions. Fourth, the authors of the metal-binding predictors are suggested to make their approaches publicly available, preferably as both webservers and standalone software. Particularly, high-throughput predictors promise a wide application among the research community since they can be used to perform large-scale computations, such as proteome-level predictions.



---

*Source: 8965712-2022-03-31.xml*