pLoc_Deep-mGneg: Predict Subcellular Localization of Gram Negative Bacterial Proteins by Deep Learning

The recent worldwide spreading of pneumonia-causing virus, such as Coronavirus, COVID-19, and H1N1, has been endangering the life of human beings all around the world. In order to really understand the biological process within a cell level and provide useful clues to develop antiviral drugs, information of Gram negative bacterial protein subcellular localization is vitally important. In view of this, a CNN based protein subcellular localization predictor called “pLoc_Deep-mGnet” was developed. The predictor is particularly useful in dealing with the multi-sites systems in which some proteins may simultaneously occur in two or more different organelles that are the current focus of pharmaceutical industry. The global absolute true rate achieved by the new predictor is over 98% and its local accuracy is around 94% 100%. Both are transcending other existing state-of-the-art predictors significantly. To maximize the convenience for most experimental scientists, a user-friendly web-server for the new predictor has been established at http://www.jci-bioinfo.cn/pLoc_Deep-mGneg/, which will become a very useful tool for fighting pandemic coronavirus and save the mankind of this planet.


Introduction
Knowledge of the subcellular localization of proteins is crucially important for fulfilling the following two important goals: 1) revealing the intricate pathways Advances in Bioscience and Biotechnology that regulate biological processes at the cellular level [1] [2]; 2) selecting the right targets [3] for developing new drugs.
With the avalanche of protein sequences in the post-genomic age, we are challenged to develop computational tools for effectively identifying their subcellular localization purely based on the sequence information.
In 2018, a very powerful predictor, called "pLoc_bal-mGneg" [4], was developed for predicting the subcellular localization of Gram negative bacterial proteins based on their sequences information alone. It has the following remarkable advantages. 1) Most existing protein subcellular location prediction methods were developed based on the single-label system in which it was assumed that each constituent protein had one, and only one, subcellular location (see, e.g., [5] [6] [7] and a long list of references cited in a review papers [8]). With more experimental data uncovered, however, the localization of proteins in a cell is actually a multi-label system, where some proteins may simultaneously occur in two or more different location sites. This kind of multiplex proteins often bears some exceptional functions worthy of our special notice [2]. And the pLoc_bal-mGneg predictor [4] can cover this kind of important information missed by most other methods since it was established based on the multi-label benchmark dataset and theory. 2) Although there are a few methods (see, e.g., [9] [10]) that can be used to deal with multi-label subcellular localization for proteins, the prediction quality achieved by pLoc_bal-mGneg [4] is overwhelmingly higher, particularly in the absolute true rate. 3) Although the pLoc_bal-mGneg predictor [4] has the aforementioned merits, it has not been trained at a deeper level yet [11] [12] [13] [14].
The present study was initiated in an attempt to address this problem. As done in pLoc_bal-mGneg [4] as well as many other recent publications in developing new prediction methods (see, e.g., [15] [16]), the guidelines of the 5-step rule [17] are followed. They are about the detailed procedures for 1) benchmark dataset, 2) sample formulation, 3) operation engine or algorithm, 4) cross-validation, and 5) web-server. But here our attentions are focused on the procedures that significantly differ from those in developing the predictor pLoc_bal-mGneg [4].

Benchmark Dataset
The benchmark dataset used in this study is exactly the same as that in pLoc_bal-mGneg [4]; i.e., where 1  only contains the Gram-negative bacteria protein samples from the "Cell inner membrane" organelle (cf.

Proteins Sample Formulation
Now let us consider the 2 nd step of the 5-step rule [17]; i.e., how to formulate the biological sequence samples with an effective mathematical expression that can truly reflect their essential correlation with the target concerned. Given a protein sequence P, its most straightforward expression is where L denotes the protein's length or the number of its constituent amino acid residues, 1 R is the 1 st residue, 2 R the 2 nd residue, 3 R the 3 rd residue, and so forth. Since all the existing machine-learning algorithms} can only handle vectors as elaborated in [3], one has to convert a protein sample from its sequential expression (Equation (2)) to a vector. But a vector defined in a discrete model might completely miss all the sequence-order or pattern information. To deal with this problem, the Pseudo Amino Acid Composition [18] or PseAAC [19] has been proposed. Ever since then, the concept of "Pseudo Amino Acid Composition" has been widely used in nearly all the areas of computational proteomics with the aim to grasp various different sequence patterns that are essential to the targets investigated (see, e.g., [20] [21] as well as a long list of references cited in [22]). Meanwhile, because it has been widely and increasingly used, four powerful open access soft-wares, called "PseAAC" [23], "PseAAC-Builder" [24], "propy" [25], and "PseAAC-General" [26], were established: the former three are for generating various modes of special PseAAC [27]; while the 4 th one for those of general PseAAC [17], including not only all the special modes of feature vectors for proteins but also the higher level feature vectors such as "Functional Domain" mode, "Gene Ontology" mode, and "Sequential Evolution" or "PSSM" mode. Encouraged by the successes of using PseAAC to deal with protein/peptide sequences, its idea and approach were extended to PseKNC (Pseudo K-tuple Nucleotide Composition) to generate various feature vectors for DNA/RNA sequences [28] that have proved very successful as well (see, e.g., [29] [30]).
According to the concept of general PseAAC [17], any protein sequence can be formulated as a PseAAC vector given by [ ] where T is a transpose operator, while the integer Ω is a parameter and its value as well as the components will depend on how to extract the desired information from the amino acid sequence of P, as elaborated in [4]. Thus, by following exactly the same procedures as described in the Sec-Advances in Bioscience and Biotechnology tion 2.2 of [4], each of the protein samples in the benchmark dataset can be uniquely defined as an 8-D numerical vector as given in Supporting Information S2, which can also be directly downloaded at http://www.jci-bioinfo.cn/pLoc_bal-mGneg/Supp2.pdf.

Installing Deep-Learning for Three Deeper Levels
In this study, we use multilayer perceptron neural network model, which consists of 3 fully connected layers and was used to predict subcellular localization of multi-label human proteins, as illustrated in Figure 1. "accuracy". The batch size is set as 28, and the epochs is 100. The predicted results were decided by the output of the threshold θ. If the output is greater than 0.5, the outcome was true; otherwise, false. For more information about this, see [11], where the details have been clearly elaborated and hence there is no need to repeat here.
The new predictor developed via the above procedures is called "pLoc_ Deep-mGneg", where "pLoc_Deep" stands for "predict subcellular localization by deep learning", and "mGneg" for "multi-label Gram negative proteins".

Results and Discussion
According to the 5-step rules [17], one of the important procedures in developing Figure 1. An illustration to show a dense neural network with 3 fully connected layers. Adapted from [11] with permission.

A Set of Five Metrics for Multi-Label Systems
Different from the metrics used to measure the prediction quality of single-label systems, the metrics for the multi-label systems are much more complicated [31]. To make them more intuitive and easier to understand for most experimental scientists, here we use the following intuitive Chou's five metrics [32] where q N is the total number of query proteins or tested proteins, M is the total number of different labels for the investigated system (for the current study it is cell 8 L = ), means the operator acting on the set therein to count the number of its elements,  means the symbol for the "union" in the set theory,  denotes the symbol for the "intersection", k  denotes the subset that contains all the labels observed by experiments for the k-th tested sample, * k  represents the subset that contains all the labels predicted for the k-th sample, and ( ) * * 1, if all the labels in are identical to those in Δ , 0, otherwise In Equation (4), the first four metrics with an upper arrow ↑ are called positive metrics, meaning that the larger the rate is the better the prediction quality will be; the 5 th metrics with a down arrow ↓ is called negative metrics, implying just the opposite meaning.
From Equation (4) we can see the following: 1) the "Aiming" defined by the 1 st sub-equation is for checking the rate or percentage of the correctly predicted la-bels over the practically predicted labels; 2) the "Coverage" defined in the 2 nd sub-equation is for checking the rate of the correctly predicted labels over the actual labels in the system concerned; 3) the "Accuracy" in the 3 rd sub-equation is for checking the average ratio of correctly predicted labels over the total labels including correctly and incorrectly predicted labels as well as those real labels but are missed in the prediction; 4) the "Absolute true" in the 4 th sub-equation is for checking the ratio of the perfectly or completely correct prediction events over the total prediction events; 5) the "Absolute false" in the 5 th sub-equation is for checking the ratio of the completely wrong prediction over the total prediction events.

Comparison with the State-of-the-Art Predictor
Listed in Table 1 are the rates achieved by the current pLoc_Deep-mGneg predictor via the cross validations on the same experiment-confirmed dataset as used in [4]. For facilitating comparison, listed there are also the corresponding results obtained by the pLoc_bal-mGneg [4], the existing most powerful predictor for identifying the subcellular localization of Gram negative proteins with both single and multiple location sites. As shown in Table 1, the newly proposed predictor pLoc_Deep-mGneg is remarkably superior to the existing state-of-the-art predictor pLoc_bal-mGneg in all the five metrics. Particularly, it can be seen from the table that the absolute true rate achieved by the new predictor is over 98%, which is far beyond the reach of any other existing methods (see, e.g., [35] [36] [37]). This is because it is extremely difficult to enhance the absolute true rate of a prediction method for a multi-label system as clearly elucidated in [4].
Actually, to avoid embarrassment, many investigators even chose not to mention the metrics of absolute true rate in dealing with multi-label systems (see, e.g., [38] [39]).
Moreover, to in-depth examine the prediction quality of the new predictor for the proteins in each of the subcellular locations concerned (cf. Table 2), we used a set of four intuitive metrics that were derived in [40] based on the Chou's symbols introduced for studying protein signal peptides [41] and that have ever since been widely concurred or justified (see, e.g., [42] [43]). For the current study, the set of metrics can be formulated as:  [4], where the reported metrics rates were obtained by the jackknife test on the benchmark dataset of Supporting Information S1 that contains experiment-confirmed proteins only. c The proposed predictor; to assure that the test was performed on exactly the same experimental data as reported in [4] for pLoc_bal-mGneg. Advances in Bioscience and Biotechnology where Sn, Sp, Acc, and MCC represent the sensitivity, specificity, accuracy, and Mathew's correlation coefficient, respectively, and i denotes the i-th subcellular location (or subset) in the benchmark dataset. Listed in Table 2 are the results achieved by pLoc_Deep-mGneg for the Gram negative proteins in each of 8 subcellular locations. As we can see from the table, nearly all the success rates achieved by the new predictor for the Gram negative proteins in each of the 8 subcellular locations are within the range of 90% -100%, which is once again far beyond the reach of any of its counterparts.
Meanwhile, as a byproduct, the present paper has also stimulated a series of somewhat provocative but quite interesting papers (see, e.g., [44]- [49]).

Web Server and User Guide
As pointed out in [50], user-friendly and publicly accessible web-servers represent the future direction for developing practically more useful predictors.
Actually, user-friendly web-servers will significantly enhance the impacts of theoretical work because they can attract the broad experimental scientists [22].
In view of this, the web-server of the current pLoc_Deep-mGneg predictor has also been established at http://www.jci-bioinfo.cn/pLoc_Deep-mGneg/, by which users can easily get their desired data without the need to go thru the mathematical details.

Conclusion
It is anticipated that the pLoc_Deep-mGneg predictor holds very high potential to become a useful high throughput tool in identifying the subcellular localization of Gram negative bacterial proteins, particularly for finding multi-target drugs that is currently a very hot trend in drug development. Most important is that the predictor will become a very useful tool for fighting against the coronavirus to save mankind on this planet.