Next Article in Journal
Cutting Edge Research for Exploration of Biomolecules for Gemcitabine-Based Chemo-Resistant Advanced Bile Duct Cancer: From Basic Study to Clinical Trial
Previous Article in Journal
Long-Term Consequences of Placental Vascular Pathology on the Maternal and Offspring Cardiovascular Systems
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Discovering the Ultimate Limits of Protein Secondary Structure Prediction

1
Institute of Bioinformatics and Systems Biology, National Yang Ming Chiao Tung University, Hsinchu 300, Taiwan
2
Department of Biological Science and Technology, National Yang Ming Chiao Tung University, Hsinchu 300, Taiwan
3
The Center for Bioinformatics Research, National Yang Ming Chiao Tung University, Hsinchu 300, Taiwan
*
Author to whom correspondence should be addressed.
Biomolecules 2021, 11(11), 1627; https://doi.org/10.3390/biom11111627
Submission received: 21 September 2021 / Revised: 25 October 2021 / Accepted: 28 October 2021 / Published: 3 November 2021

Abstract

:
Secondary structure prediction (SSP) of proteins is an important structural biology technique with many applications. There have been ~300 algorithms published in the past seven decades with fierce competition in accuracy. In the first 60 years, the accuracy of three-state SSP rose from ~56% to 81%; after that, it has long stayed at 81–86%. In the 1990s, the theoretical limit of three-state SSP accuracy had been estimated to be 88%. Thus, SSP is now generally considered not challenging or too challenging to improve. However, we found that the limit of three-state SSP might be underestimated. Besides, there is still much room for improving segment-based and eight-state SSPs, but the limits of these emerging topics have not been determined. This work performs large-scale sequence and structural analyses to estimate SSP accuracy limits and assess state-of-the-art SSP methods. The limit of three-state SSP is re-estimated to be ~92%, 4–5% higher than previously expected, indicating that SSP is still challenging. The estimated limit of eight-state SSP is 84–87%. Several proposals for improving future SSP algorithms are made based on our results. We hope that these findings will help move forward the development of SSP and all its applications.

1. Introduction

Protein secondary structure prediction (SSP) means to predict the per-residue backbone conformation of a protein based on the amino acid sequence. It is an essential structural biology technique with a variety of applications. New SSP algorithms have been published almost every year for seven decades, and the competition for accuracy has always been fierce [1]. Recently, it is generally believed that the accuracy of current SSP methods has approached the theoretical upper limit of 88%, and it may be difficult to make breakthroughs; thus, the intensity of research and development in SSP seems weakened (see Figure 1 of [1]). Nevertheless, we found that the limit of SSP estimated by previous studies might not be precise enough and seemed underestimated. We believe that precisely determining the limit of SSP will help re-energize this field, set new directions for SSP developments, and ultimately benefit all applications relying on SSP. In this study, we performed exhaustive pairwise sequence and structure alignments on protein structure databases and estimated the theoretical limits of three- and eight-state SSPs. Besides, the experimental results revealed valuable information for future SSP developments.
Theoretically, it is possible to predict the three-dimensional structure and function of a protein through its sequence. However, such prediction is still challenging and time-consuming. A simplified strategy predicts the secondary structure based on sequence and uses the sequence and predicted secondary structural information to predict its 3D structure. Recently, genomes have been rapidly sequenced, and a large number of protein sequences are determined so that the number of proteins with known sequences exceeds the number of proteins with known structures by more than a thousand times. SSP is particularly important in this postgenomics era because we may use it to predict the structure and function of proteins rapidly and achieve many applications. For example, identification of disease-causing genetic mutations [2,3,4], homology detection and multiple alignments of proteins [5,6,7], evolutionary analyses [8,9,10], folding studies [11,12], and the prediction of enzyme target sites [13,14], functional sites [15,16], binding sites [17,18,19], and bioengineering sites [20,21,22,23]. Determinations of epitopes [24,25,26], protein-protein interaction surfaces [27,28], disordered regions [29,30,31,32], and protein subcellular localization [33,34] utilize SSP as well. Predicted secondary structures can be used in de novo protein design [35,36] and drug design [37,38,39]. Perhaps these wide applications are the main reason why, even when SSP accuracy has been considered to approach the limit, there have still been around five methods proposed every year since 2010 [1].
Before carrying out SSP, protein local backbone conformations were first classified into secondary structure elements (SSEs). There are two general sets of SSE: the three-state SSEs that describe protein conformations as helixes, strands, and coils, and the eight-state SSEs defined by the DSSP (Define Secondary Structure of Proteins) algorithm [40]. Hence, the development and assessment of SSP methods also fall into two categories, i.e., three-state predictions with Q3 accuracy or eight-state predictions with Q8 accuracy.
The SSP research field has been initiated since Pauling and Corey proposed the helix and sheet conformations of polypeptide chains roughly 70 years ago (1951) [41,42]. In the 1970s, amino acid propensities and residue or segmental physiochemical properties of polypeptides were used to perform SSP based on statistical approaches and achieved ~60% Q3 accuracy [43,44]. Since the late 1980s, as machine learning techniques were applied, the scale of the SSP feature set had been greatly expanded by adopting more residue properties and wider “window” (segment) sizes, and the Q3 was pushed to ~70% [45]. In 1999, a neural network machine learning SSP method, Psipred, first utilized the position-specific scoring matrix (PSSM) generated by PSI-BLAST [46] to be the feature set and accomplished a Q3 of 76.5% [47]. Since the big jump made by Psipred, the basic procedure of modern SSP methods has been established, which utilizes the PSSM as the main feature set and machine learning as the core algorithm. The main competition for accuracy then shifted to improvements in machine learning. The Q3 of methods developed in the early 2000s was 76–78% in general [48]. Meanwhile, the Q8 of SSpro8 was 62.6% [49]. When Dor and Zhou developed the SPINE using neural networks with a large amount of training data, the Q3 reached 79.5% in 2007 [50]. Soon after that, another neural network-based method, Jpred 3 [51], broke the 80% record of Q3 in 2008. Neural network-derived algorithms hence became the core of most SSP methods. In recent years, the architecture of neural networks used in SSP has become increasingly complex, including the BRNN [52], LSTM-BRNN [53,54], CNF [55], CRNN [56], and many deep learning techniques [57,58,59]. However, the progress of SSP accuracy seems to have encountered a plateau. Compared to the earliest sixty years, the improvement in accuracy during the past decade has been limited. The Q3 of current SSP methods, tested with independent datasets (<25% sequence identity with the training data), is between 81% and 86%, while the accuracy of Q8 is not yet above 77% [54,60].
Interestingly, even when SSP accuracy was still rapidly increasing, there had been studies estimating the upper limit of Q3 to be ~88% [61,62]. The ultimate reason for the existence of an upper limit in SSP is that proteins are not static/rigid bodies but dynamic objects with some parts being more flexible than others [48]. A protein may adopt different conformations under different conditions; besides, conformational changes commonly occur when interacting with ligands, substrates, other proteins, or cellular components. Thus, the conformation of a protein determined with different pH values, temperatures, metal ions, or cofactors would not be identical. Moreover, structures of the same protein solved by different techniques may have noticeable differences due to experimental equipment or conditions. For instance, X-ray diffraction requires crystallization of the protein, a relatively rigid and stable state, while nuclear magnetic resonance (NMR) spectroscopy determines protein structure in the solution phase, a state so flexible that an NMR experiment typically produces an ensemble of models. Inevitably, resolution, errors, and software inconsistencies are responsible for the limit of SSP, too [1,48,62]. Another type of limitation lies in how SSEs are assigned. Automated secondary structure assignment programs define SSEs by hydrogen-bond pattern, geometric features, or expert knowledge, and inconsistencies between programs are often observed [63]. Some technical characteristics in SSP methods may also cause limitations. The general methodology applied by state-of-the-art SSP methods works based on PSSM, implying their core assumption is that homologous proteins have similar structures. Indeed, it has been known that proteins sharing similar sequences have similar 3D structures [64]. However, being similar does not mean being identical. Therefore, in the current homology-based SSP methodology, the fundamental reason for the limit of accuracy is the structural difference or inconsistency between homologs.
Previous studies on the limit of SSP made estimations by measuring the secondary structural (in)consistency between structural homologs. In 1993, Levin et al. analyzed seven structural families (82 proteins sharing >17% sequence identities) from a protein structure database [65] and concluded that those structural homologs had only 88.8% secondary structural consistency [61]. In 1994, Rost et al. analyzed 140 pairs of structures sharing >30% sequence identities and reported that homologous structures differed by 12% in their three-state secondary structures, and thus the upper limit of Q3 should be 88% [62]. They also proposed a new accuracy measure termed SOV (segment overlap), which calculated the secondary structural consistency based on segments instead of residues. The upper limit of the first proposed SOV (v’94) was estimated to be 90 [62], which was then revised downward when an updated SOV (v’99) was proposed (see Table I in [66]). Yang et al. used PSI-BLAST to retrieve homologs from the Protein Data Bank (PDB) for a set of 1199 representative structures and classified those homologs according to sequence identities. They concluded that the limit of Q3 is >88% for structural homologs sharing ≥30% sequence identities [1].
After reviewing these pioneering works, we noticed that the limit of SSP can be more precisely estimated and that several estimates for evaluating current SSP methods are still in need. First of all, due to the amount of data available at that time, the earliest studies only analyzed a small number of proteins. Second, the secondary structural differences or consistencies of structural homologs were measured by sequence alignments, or by structural alignments with algorithms that are relatively primitive compared to today. Large-scale structural alignments were time-consuming, yet if sufficient computing power were available, it would be more appropriate to measure structural homologs’ structural differences/consistencies by structural alignments. Moreover, previous studies had only analyzed Q3 and the first SOV3 (v’94). The limits of Q8 and the new SOV3/8 (v’99) have not been estimated. Despite their high performance in residue-based three-state predictions, current SSP methods are still weak in segment-based and eight-state predictions [55]. The competition of SSP accuracy is about to switch to Q8 and SOV. These estimates shall be very informative for future SSP developments.
To precisely estimate the limit of SSP, large-scale analyses with multiple repeats were conducted in this study. Over ten times the number of PDB structures (~15,000) were used to repeat the experiments of previous studies. Also, exhaustive pairwise alignments were applied to 3192 structural families (~170,000 proteins) of the SCOP (Structural Classification of Proteins) database [67]. Four sequence and three structural alignment algorithms were recruited to measure the secondary structural (in)consistency between homologs. Figure A1 in Appendix A provides a schema of this study and the presented figures. Our results revealed that the theoretical upper limits of Q3 and v’99 SOV3 might be both revised upwards to 90–92%. As for eight-state measures, the upper limit of Q8 is 84–86%, and v’99 SOV8 is 85–87%. However, these theoretical limits of SSP should be challenging to achieve. Considering the current homology-based methodology, we proposed that the practical limits of SSP methods that work with PSI-BLAST PSSMs and PSSM reference datasets of <90% sequence identities would be Q3: 87.0 ± 1.2%, SOV3: 86.6 ± 1.5%, Q8: 78.8 ± 1.6%, and SOV8: 79.6 ± 1.6%. The results are divided into theoretical and practical limits because, although SSP limits can be accurately estimated by analyzing structural differences between structural homologs using structural alignments, in practice, current SSP methods can only rely on the homology detected by sequence alignments (structure of the query protein is unknown in an actual SSP). After further investigation, we found that the restriction caused by this “sequence-only” scenario may be conquerable, and specific suggestions have been made in the Discussion.
This study aimed to precisely estimate SSP limits and provide useful information for improving future SSP methods. Hopefully, our discoveries will help push SSP out of the accuracy plateau and eventually move forward all research and application fields depending on accurate protein secondary structure predictions.

2. Materials and Methods

2.1. Experimental Datasets

2.1.1. Datasets Prepared from the Protein Data Bank

In the review by Yang et al. [1], a nonredundant (nr) query dataset of 1199 proteins and a PSSM reference dataset of all known protein structures were obtained from PDB 2005. In our work, the applied PDB data were released in 2015 Dec. The nr query dataset was prepared by downloading the 100% sequence identity representative PDB entity list of 2015 established by NCBI [68,69] (73,073 polypeptides; termed nrPDB100-2015 in this report) and then reducing the sequence identity from 90%, 80% … to 30% by three algorithms to ensure low redundancy of the final query proteins: CD-HIT [70], USEARCH [71], and MMseqs2 [72]. The produced nr datasets were named nrPDB90-2015, nrPDB80-2015, and so on. In the first experiment of Results, the 30% identity nonredundant dataset nrPDB30-2015 (15,059 protein entities; see Dataset S1) was used as the query dataset, and the whole 2015 Dec. PDB (PDB-2015; 273,920 entities; Dataset S2) was applied as the PSSM reference set. Other nr datasets were utilized to test how the homology of developmental datasets affects the accuracy of SSP methods.

2.1.2. Homolog Data Obtained from the SCOP Database

Homologous protein structural data of SCOP were obtained from the SCOPe (v.2.07), which hierarchically classified protein structural domains into classes, folds, superfamilies, and families [67,73]. The four major classes were used, a: all-alpha, b: all-beta, c: alpha/beta (alpha-beta mixed), and d: alpha+beta (alpha-beta segregated). After (1) removing proteins with homology classifications only by automated matching procedures but not yet by SCOPe manual verifications, and (2) discarding entities unable to be correctly processed by DSSP, 169,977 domain entities were retrieved from 3192 families, which belonged to 927 folds. We termed this dataset SCOP-2.07 (see Dataset S3). To examine the performance of current SSP methods on proteins of different structural classes or sizes, we prepared a query dataset from the SCOP-2.07 following these steps: (1) reduce the sequence identity of SCOP-2.07 to 30% using the homology reduction methods mentioned above. (2) Classify the 30% identity representative sequences into four subsets according to SCOP classes. (3) From each class subset, randomly select 60 proteins whose sizes are small (<150 residues), medium (150–299 residues), or large (≥300 residues). (4) Collect the selected proteins into a dataset, which hence comprise the same amount of proteins from the four major SCOP classes and the three size groups. The dataset prepared by this procedure contained 720 proteins (i.e., 4 × 60 × 3) and was called the SCOP720 dataset.

2.1.3. The Source of PSSM Reference Sequences

For verifying whether our results properly estimated SSP limits, several state-of-the-art SSP methods were used to perform the tests described in Discussion. In those tests, the source of PSSM reference sequences was the 90% identity nr UniRef90 dataset established by Universal Protein Resource [74] in 2015 (UniRef90-2015; 38.2 million proteins). When necessary to ensure that the homology between the reference and query datasets was lower than some sequence identity, the interdataset homology reduction software CD-HIT-2D [70] was applied. Note that a PSSM reference dataset is typically termed a “target dataset/database” in fields like sequence similarity searches. However, in many protein structural prediction fields, a “target sequence” means the query sequence for prediction. To avoid confusion, we termed a target dataset a reference dataset in this report.

2.2. Assignment of Protein Secondary Structure

The secondary structure of PDB and SCOP structure entities was computed by DSSP [40] (release date: 1 April 2000), which assigned for each residue an SSE code. The DSSP classified protein backbone conformations into α-helix (H), 310-helix (G), π-helix (I), extended strand participating in β ladder (E), isolated β-bridge (B), hydrogen-bonded bend (S), turn (T), and loop/irregular (C). In most SSP works, this eight-state classification was reduced into a three-state version, where the codes H, G, and I were represented by H (helices), E and B by E (extended strands), and the others by C (coils/loops). The three-state SSE codes used in this study were converted from the DSSP eight-state codes following those rules. It is noteworthy that the program of DSSP was modified in 2015 (v2.2.1) such that much more π-helices would be assigned [75]. This study and all the applied SSP methods were established based on the DSSP released before 2015. Because the occurrence of π-helices is very low in protein structures, the influence of the DSSP version on the estimated SSP limits and assessed accuracies of SSP methods is negligible. However, for some specific research topics where the correct assignment of π-helices is critical, attention must be paid to the choice of DSSP version.

2.3. Experimental Procedures

2.3.1. Estimating SSP Limits by Structural Homologs Identified from the PDB Datasets

Following the procedure of [1], homologous protein structures were identified using the PSI-BLAST sequence alignment algorithm to search a PDB reference dataset against a set of query sequences for homolog pairs whose E-value was lower than a cutoff of 10–3. The reference dataset used in this study was PDB-2015, and the query set was nrPDB30-2015. Since nrPDB30-2015 was more than 10 times larger than the 1199-protein query set applied in [1], we were able to perform multiple repeats in this experiment. The nrPDB30-2015 was randomly divided into 10 subsets, each with 1500 proteins (remnants discarded). These query sets were used to perform PSI-BLAST searches from PDB-2015 for query-reference homolog pairs, which were then classified according to the sequence identity computed by PSI-BLAST into a series of identity levels, 0, 5, 10, …, and 100%. For each level, in addition to PSI-BLAST, homolog pairs were aligned by the structure alignment algorithm FAST [76]. Based on these sequence and structure alignments, the percentages of the secondary structural consistency between all pairs were averaged. Standard deviations of the average percent consistencies of the 10 repeats were also computed. The secondary structural consistency between homologs was calculated in two ways: the three-state and eight-state consistencies, equivalent to the Q3 and Q8 accuracies widely used in SSP research (see Equation (1) in Section 2.8.1).

2.3.2. Estimating SSP Limits by Structural Homologs Determined by the SCOP Database

Proteins classified in the same SCOP fold possess similar composition of secondary structures, structural architecture, and topology; proteins in the same family are closely related homologs with the same functions and clear evidence for their common evolutionary origin [67,77]. There were 927 folds comprising 3192 families in the SCOP-2.07 dataset, and we randomly divided them into 10 subsets for multiple repeats of the experiment. We performed all-against-all pairwise sequence and structure alignments for each family by several algorithms. According to the alignments, the percent secondary structural consistency between each pair was computed. The secondary structural consistencies of homolog pairs belonging to the same fold were weighted averaged according to protein sizes to represent the fold. The secondary structural consistencies of all folds belonging to the same SCOP-2.07 subset were finally averaged to represent the subset. After applying the same procedure to all subsets, the standard deviation of these 10 repeats was calculated. The largest families in the SCOP-2.07 were composed of 3000 to 6000 protein entities. All-against-all alignments were very time-consuming. For instance, ~38.3 million alignments were made for the C1 set domain family (SCOPe ID: b.1.1.1; 6190 entities). For efficiency, the distributed computation system we developed for the iSARST protein structural similarity search system was utilized [78].

2.4. Applied Sequence Alignment Methods

The sequence alignment methods utilized in this study included local alignment algorithms PSI-BLAST (v.2.3.0) [46], the classic BLAST (v.2.2.13) [79], and Water [80], as well as the global sequence alignment algorithm Stretcher [81]. In the experiment on PDB, the PSI-BLAST was applied using its default parameters except for the E-value cutoff (set as 10−3, stricter than the default value 10). To ensure all applied algorithms were executed on an equivalent basis, their gap opening and extension penalties were uniformly set to be 11 and 2, respectively. Water was an implementation of the Smith-Waterman algorithm [82], and Stretcher was developed based on the Needleman–Wunsch algorithm [83]; both programs were obtained from the EMBOSS server (v.EMBOSS:6.6.0.0) [80,81].

2.5. Applied Structure Alignment Methods

Structure alignment algorithms including FAST [76], TM-align [84], and SARST [85] were applied. The programs of these algorithms accepted single-chain and single-model protein structure entity files in the PDB format. To apply them to the structure entry files obtained from the PDB database, when an entry file contained multiple chains, the chain of interest was extracted and saved in an individual chain entity file. When a chain entity had multiple models, only the first was used. As for domain entities from the SCOP database, the preprocessing was more complicated because a SCOP domain might comprise polypeptide fragments from multiple chains. Before a SCOP domain entity file containing multiple fragments could be accepted by these alignment algorithms, (1) the chain identifiers of the fragments were unified as “A,” (2) all residues of the entity were renumbered according to the sequence they appeared in the structure file, and (3) all atoms of the entity were also renumbered.

2.6. Applied Secondary Structure Prediction Methods

Several state-of-the-art SSP methods were applied to help verify the estimated upper limits of SSP, inclusive of three-state algorithms Psipred v3.3 [47], SpineX v2.0 [86], Scorpion v1.0 [87] and Spider2 v2.0 [88], and eight-state algorithms RaptorX v1.0 [89], SSpro8 v5.2 [90] and DeepCNF v1.02 [55]. The utilized versions of these methods were all trained and released before 2016. Their original pipelines relied on different editions and settings of PSI-BLAST to generate PSSM. To make all methods applied on an equal basis, we modified their pipeline scripts to unify the edition and settings of the PSI-BLAST engine (v2.3.0) [46]. The unified PSI-BLAST settings were based on the common parameter settings used by most of these SSP methods and the defaults recommended by PSI-BLAST (see Table S1 for the original and modified settings). To ensure that we correctly applied these methods and the modified settings of PSI-BLAST would not hinder their performance, we preliminarily conducted SSP using these state-of-the-art methods with UniRef90-2015 as the PSSM reference dataset. The query sets were the TS115 and CASP12 prepared by Yang et al. [1], which were consisted of novel PDB structures released later than 1 January 2016, meaning that the machine learning core of the applied SSP methods had not learned those proteins. In particular, the CASP12 was obtained from the competition of the 12th biannual meeting of Critical Assessment of Structure Prediction techniques [91]. Results listed in Table 1 showed that the accuracy values computed using our pipelines were all close to the values obtained from previous reports of these SSP methods assessed with equivalent datasets [1,55].

2.7. Implemented Secondary Structure Prediction Models

In implementing the general SSP methodology (see Section 4.3 of Discussion), the PSI-BLAST-generated PSSM utilized by current SSP methods and the amino acid information used in many recent SSP works [55,92,93,94,95] were applied as the predictive feature set. For each residue, the scale of the feature set was expanded by referring to the features of its nearby residues based on the “window” concept proposed in the Psipred algorithm [47]. The window size was set at 5. The SSP models were trained with the integrated machine learning system we developed [21,22]. The integrated algorithms included artificial neural networks, decision trees, support vector machines, and bootstrap sampling. Settings of these algorithms were according to the defaults of this system [22].

2.8. Computation of Accuracy Measures

2.8.1. The Q Accuracy Measure and the Percent Secondary Structural Consistency between Homologous Protein Pairs

For measuring the secondary structural consistency between a pair of homologous proteins, the secondary structures of the two proteins were first represented by their per residue SSE sequences generated by DSSP. Next, according to the alignment between the proteins, the percent consistency (C3 or C8) in SSE sequences was calculated for all aligned equivalent residues using the following equation:
C s = n c n e × 100   % ,     s = 3 ,   for   three - class   SSE   codes 8 ,   for   eight - class   SSE   codes
where ne stands for the number of equivalent residue pairs based on the alignment, and nc represents the number of equivalent residue pairs having the same SSE codes. This definition of percent SSE consistency is equivalent to the SSP accuracy measure Q, defined as the number of correctly predicted residues divided by the total number of predicted residues. In this report, the C3/C8 measures are simply written as Q3/Q8.

2.8.2. The SOV Measure

In addition to the Q measure, we also estimated the limit of SSP accuracy using the SOV v’99 edition [66]. SOV was developed for evaluating the consistency between the predicted and actual SSE sequences of a protein based on secondary structure segments instead of residues. It has been generally regarded as a more critical accuracy measure than Q for its capability of capturing the overall consistency of two SSE sequences and reducing noises from individual residues [62,66]. Most studies calculated SOV based on the three-state SSEs, i.e., the SOV3. In this work, the eight-state SOV (SOV8) was also calculated.

2.8.3. The Weighted Average of Accuracy Measures

Protein pairs may have quite different sizes. If the classic arithmetic mean were applied to compute the average Q or SOV (i.e., sum/n, where sum is the summation of measure values from all pairs and n denotes the number of pairs), the influence of small proteins might be overestimated while that of large proteins be underestimated. In this study, the average Q and SOV were calculated as the weighted average, with the number of aligned residues for each protein pair used as the weight, that is:
Q ¯ = i = 1 N n e , i × Q i i = 1 N n e , i
SOV ¯ = i = 1 N n e , i × S O V i i = 1 N n e , i
where Q ¯ and SOV ¯ represent the weighted average of Q and SOV, N denotes the number of protein pairs, ne,i means the number of aligned equivalent residue pairs in the ith protein pair, and Qi and SOVi represent the measure values of the ith protein pair.

3. Results

3.1. Experiments Based on PDB Datasets

The procedure Yang et al. estimated the upper limit of SSP [1] was repeated to verify whether we correctly applied all algorithms and check whether the conclusion of [1] was still applicable as the amount of protein structure data increased rapidly. A set of query structures was used to retrieve homologs from the entire PDB by PSI-BLAST sequence alignment searches [46]. The secondary structural consistency between homologs was computed for a given sequence identity level (individual) or identity threshold (cumulative). Since the query set we utilized was over 10 times larger than previous works’, 10 repeats of the experiment were thus allowed by random sampling.
As Figure 1 shows, both the curve for individual identity levels and the cumulative curve of our results were similar to those reported in [1], with minor differences at identities ~70%. The main reason for the differences was that Yang et al. discarded structures solved not by X-ray or with low resolution, but we did not. By applying the same restrictions, our curves got closer to theirs. However, to make conclusions applicable to a broader scope, we decided to keep using all available structures without restrictions on the method of structural determination or resolution. The data of this figure are listed in Table 2. When the identities between homologous structures were ≥30%, their secondary structural consistency was ≥87.8%, agreed well with [1,62]. The consistency was also analyzed by FAST structural alignment [76]. The results indicated that the limit of SSP could be higher than that estimated by sequence alignments, especially for homologs with low identities. Judging from the structure-alignment data, the limit of SSP was >90% as the identity of homologs was ≥30% (see identity level 30–35 for FAST in Figure 1 and Table 2).
Previous studies estimated SSP limits using three-state SSEs. We also used eight-state SSEs and estimated the Q8 accuracy, a more challenging assessment. In general, the Q8s achieved by state-of-the-art SSP methods are ≤77% [54,60]. Figure 1 and Table 3 demonstrate that when the sequence identities between homologs were >30%, the upper limit of Q8 was >79% as estimated by PSI-BLAST sequence alignments; when estimated by FAST structure alignments, the limit of Q8 was >83%. These data revealed that there is still much space left for improvement in eight-state SSP.
The three-state secondary structural consistency between a pair of proteins is analogous to the Q3 accuracy of SSP, if one protein were considered as the query and the other were the answer. The full three-/eight-state results of this experiment can be found in Table S2.
To analyze the secondary structural consistency of homologs based on secondary structure segments, we computed the SOV for these PDB pairs. Computed by sequence alignments, the SOV3 was a little lower than Q3. Contrarily, the SOV3 computed by structural alignments was slightly higher than Q3 (see Table S2). When it was first proposed, the upper limit of SOV3 (v’94) was estimated to be higher than that of Q3 [62]. However, the upper limit of the updated SOV3 (v’99) has not been determined until this work. Considering structural alignment is better than sequence alignment at detecting residue-residue correspondence between proteins [61,96], it was supposed that the FAST alignments were more reliable than PSI-BLAST alignments in estimating SSP limits, inclusive of SOV.

3.2. Limit of SSP Accuracy Estimated by Protein Sequence Alignments between SCOP-Determined Structural Homologs

In the previous experiment, the quality of homologs was controlled by the E-value cutoff of PSI-BLAST as it searched PDB ([1] and Materials and Methods). Now we modified the experimental design by using the structure families classified by SCOP [67] to define homologs. For each family, all-against-all pairwise sequence alignments were performed by PSI-BLAST [46], classic BLAST [79], Water, and Stretcher [80,81]. Based on these alignments, the secondary structural consistencies of homologs were computed. Figure 2a,b shows that no matter using homologs automatically identified by PSI-BLAST or semi-manually verified by structural biologists like SCOP, as long as the identities were ≥30%, the secondary structural consistencies between homologs were similar. As for identities <30%, the consistency of SCOP family homologs was obviously higher than that of PSI-BLAST identified homologs, especially when measured by PSI-BLAST, classic BLAST, or Water (see Table S3 for raw data).
Following previous works [1,62], we calculated the secondary structural consistency for homologs with identities higher than certain thresholds (e.g., Figure 1). However, preparing protein datasets with identities higher than a threshold is uncommon in structural biology applications. The typical way is making nonredundant datasets with identities lower than certain cutoff(s), just like the nrPDB30-2015, where any two proteins shared identities <30%. We supposed that the information about secondary structural consistency of homologs sharing identities below specific cutoffs would help estimate the practical limit of SSP methods. For example, results of such cutoffs shown in Figure 2c,d and Table S3 demonstrate that the highest average three-state secondary structural consistencies were 87.3% and 80.9% (registered by Water) for homologs with identities <90% and <30%, respectively. These consistencies suggested that when the sequence identity redundancy of the datasets utilized to develop SSP algorithms was <90% (or <30%), the operational upper limit of Q3 accuracy would be around 87% (or 81%). Speculated similarly, the upper limit of Q8 for SSP algorithms would be ~79% or ~70% if the sequence identity redundancy of developmental datasets were <90% or <30%, respectively (refer to the dotted lines in Figure 2c,d). See Discussion for more information.

3.3. Limit of SSP Accuracy Estimated by Protein Structure Alignments between SCOP-Determined Structural Homologs

Since the fundamental reason for the upper limit of SSP is the structural difference among structural homologs, we supposed that it should be more reasonable to estimate SSP limits using structural methods than using sequence-based methods. The above experiment on SCOP families was repeated here using structural alignment methods, FAST [76], TM-align [84], and SARST [85]. As Figure 3 shows, the secondary structural consistencies computed by structural alignments were higher than those by sequence alignments, especially when the identity was <30%. Speculating based on the results of these structural alignments, if the sequence identity redundancy of the developmental datasets for SSP was <90% or <30%, the maximally achievable average Q3 might be 90% or 86% (by FAST), and the maximally achievable average Q8 might be 83% or 78%, respectively. These upper limits estimated by structural alignments were 3–8% higher than those by sequence alignments, compared with Figure 2 (see also Table S3 for the results of SOV3 and SOV8).

3.4. Limit of Accuracy for Various Protein Structure Classes and Sizes

We divided the homologs from SCOP into four structural classes according to SCOP classifications [67] and observed the secondary structural consistency between homologs in each class. Besides, the homologs were divided into three size groups by sequence lengths. For simplicity, only the consistencies at 90% and 30% sequence identity cutoffs are displayed in Figure 4 (see Tables S4 and S5 for complete results). No matter whether categorized by structural class or size, SOVs were generally higher than Q consistencies. Among the four classes, all-beta proteins had the lowest consistency. Among the size groups, large proteins (>300 residues) exhibited the highest consistency. In these plots, both the secondary structural consistencies computed by structure alignment (representative algorithm: FAST) and sequence alignment (representative algorithm: PSI-BLAST) were provided, demonstrating the gap between them, which we supposed was the difference between the theoretical limit of accuracy in SSP and the practical limit of accuracy for SSP algorithms.

3.5. The Lower Bound of Secondary Structure Prediction Accuracy

The experiments described above focused on the upper limit. In this part, we determined the lower bound of SSP. All domain entities of SCOP-2.07 were transformed into secondary structural strings composed of three- or eight-state SSE codes. A million protein pairs were randomly selected from this SSE string dataset, and the secondary structural consistencies (Qs and SOVs) between the random pairs were computed. Distributions of these measures are shown in Figure 5, and the detailed results and the background occurrence frequency of each three-/eight-state SSE are listed in Table S6. The averages of Q3, Q8, SOV3, and SOV8 were 34.9%, 21.9%, 32.3%, and 22.0%, respectively.

4. Discussion

4.1. Difference between Using Homologs Identified by Programs and Those Determined by SCOP

Previous studies mainly analyzed the limit of SSP using homologous protein pairs identified by sequence alignment searches [1,62]. Rost et al. took some representative structures prepared by [97] as the core protein dataset and retrieved homologs sharing >30% identities to those proteins from PDB by sequence alignment [62]. Yang et al. used a set of 1199 representative structures as the core dataset and retrieved the homologs meeting a strict sequence alignment quality cutoff, i.e., E-value <10–3, by PSI-BLAST from PDB [1]. Both studies concluded that the limit of SSP would be >88% for homologs sharing >30% sequence identities. In the present study, homologs identified by PSI-BLAST from PDB based on a much larger amount of representative core proteins were applied and reached the same conclusion (Figure 1). However, there were some shortcomings in this experimental design. First, the logic behind this experiment was to measure the secondary structural differences of homologous “structures.” Determining whether proteins were homologous structures by sequence alignments seemed contradictory and might not be as precise as structural alignments. Second, the PSI-BLAST is a local sequence alignment method, meaning that even when two proteins were not related in evolution, it might still report some aligned fragments between them. PSI-BLAST provided a statistical significance index, the E-value, to help judge the quality of alignments; however, a strict E-value setting would only reduce the number and length of aligned fragments rather than ensure the structural equivalence between the fragments.
Using well-classified protein structural families can avoid these problems. This idea was also implied by Rost et al. [62] and had been implemented by Levin et al. in [61], where seven protein structural families prepared by [65] were used to quantify SSP accuracy. The present study utilized the structural families prepared by SCOP, which used automated screening and manual curation to classify proteins based on their structural, sequence, functional, and evolutionary similarities/relationships [67,77]. Entities in the same SCOP family are closely related global structural homologs with sufficient supporting data for their common origin [67,77]. Estimating SSP limits using SCOP family entities avoids the ambiguity in identifying structural homologs by sequence alignment search. The limitation of local alignments can be ignored because family entities are global homologs. As shown in Figure 2, for identities <30%, the SSP accuracies estimated using SCOP data were much higher than those estimated using PSI-BLAST identified homologs. We supposed that manually-verified structural homologs from a reliable source fit well the purpose of estimating the limit of SSP according to the structural differences of structural homologs and thus used the SCOP database to perform most experiments.

4.2. The Different Meanings between the Upper Limits of SSP Accuracy Estimated by Sequence and Structural Alignments

SSP is limited by the “structural differences” of homologous structures. Previous studies noted that sequence alignment had inherent inaccuracies compared to structure alignment when it was applied to calculating the secondary structural (in)consistency between homologous structures [62], but they still utilized sequence alignments. We speculated that this contradiction was because of some technical limitations. In the 1990s, the speed and precision of protein structural alignment algorithms were not as high as algorithms developed after the mid-2000s. Despite the high precision of algorithms nowadays, performing structural alignments is still much more time-consuming than sequence alignments. We did apply sequence alignments to compute the secondary structural consistency of homolog structures, and the results were very close to previous reports’ (Figure 1 and Figure 2; Table 2). Nevertheless, it seemed more appropriate that the structural differences of structures should be computed by structural alignments. For the preciseness of estimation, we utilized structural alignments. The heavy computing tasks were accomplished by an efficient distributed computation system [78]. Since structural alignment could determine the residue-residue equivalence in 3D space much better than sequence alignment, we expected that the SSP upper limit estimated by structure alignment would be higher than sequence alignment. Our results agreed well with the expectation (Table 2, Table 3, Tables S2 and S3). Take the 30–35% sequence identity level for example, the highest Q3 estimated by structural alignments was 89–90%, while estimated by sequence alignments was only 87–88% (Table 2 and Table S3). We inferred that this difference was not an indicator that structural alignment is superior to sequence alignment in estimating SSP limits, but it implied something interesting:
  • Accuracies estimated by structure alignment should be the ultimate theoretical upper limit of SSP because the way they got computed strictly followed the logic that structural homologs’ structural differences were measured by structural alignments.
  • Accuracies estimated by sequence alignment is the practical upper limit of SSP since their computation folloTable Swed the procedure of the current SSP methodology.
In practice, only sequence alignment, but not structural alignment, is applicable in predicting the secondary structure for a query sequence. In other words, sequence alignments could not estimate the ultimate limit of “SSP” but could estimate the operational limit of “SSP methods.”

4.3. Information Revealed by the Estimated Limits of SSP under Different Sequence Identity Cutoffs

Previous studies on SSP limits made estimates for homolog structures with sequence identities above given thresholds (e.g., >30%). To our knowledge, the present study is the first to provide the estimates for homologs with identities below given cutoffs (e.g., <90% or <30%), supposedly helpful for knowing the limit of SSP under situations where only homologs sharing identities lower than some cutoffs can be used to make SSPs. Based on this supposition and Figure 2c, it was expected that the SSP limit estimated by PSI-BLAST at the 90% identity cutoff would define the upper bound of PSI-BLAST-driven SSP methods in which the sequence identities between (and within) the query and reference sets for training and testing were <90%. Similarly, if the sequence identities between/within the datasets were <30%, the estimates by PSI-BLAST at the 30% identity cutoff might define the upper bound of SSP accuracy.
We utilized three query datasets and seven state-of-the-art SSP algorithms to perform actual SSPs to verify the supposition. The datasets included the representative datasets prepared by Yang et al., the TS115 and CASP12, comprising structures deposited in PDB 2016 and highly nonredundant from those deposited before 2016 [1], and the SCOP720 prepared in this work. TS115 and CASP12 were particularly suitable because all the applied SSP algorithms were developed before 2016, meaning that their prediction models were not familiar with these 2016 novel proteins. The SCOP720 consisted of an equivalent number of small, medium, and large proteins from the four major structural classes (see Materials and Methods) and should achieve nonbiased evaluations. The source of PSSM reference proteins was the 90% identity nr dataset UniRef90-2015. Before the experiment, homology reductions were performed on UniRef90-2015 by CD-HIT-2D [70] to ensure that the reference dataset shared <90% identities with any protein in the query sets. The SSP methods were then applied using the same version of PSI-BLAST as the search engine. Figure 6 demonstrated that, in this 90% identity cutoff test, there was not yet any SSP method reaching the upper limits estimated by PSI-BLAST sequence alignments.
The homology between/within the datasets for “testing” was known in this test, but not the homology between/within those used to “train” the prediction models. Because of the unavailability of the source code of those SSP methods, we were unable to retrain their models. To further examine whether the secondary structural consistencies of homologs at various identity cutoffs can be used to estimate the limit of SSP methods trained/operated with protein datasets with different levels of homology reduction, we performed the following experiment. First, nr sets of sequence identity cutoffs 90%, 80%, …, and 30% were reduced from the PDB-2015 as described in Materials and Methods. Second, at each identity cutoff, from the corresponding nr dataset, 250 proteins were randomly selected as the training query set, another 250 as the testing query set, and 10,000 as the PSSM reference dataset. The reference dataset was small because the lowest identity dataset (i.e., nrPDB30-2015) contained only 15,059 proteins. Third, at each identity cutoff, a machine learning prediction model was trained with a PSSM-based feature set generated by performing PSI-BLAST searches against the reference dataset for all query sequences of the training set. Finally, by predicting the secondary structure of proteins in the testing set, accuracies were measured. As expected, results shown in Figure 7 revealed that as the homology of developmental datasets lowered, the accuracy of the established SSP model decreased, with a trend similar to the decreasing SSP accuracy limits estimated by PSI-BLAST. These results supported our idea that the accuracy estimated at different identity cutoffs may serve as a good indicator of SSP’s practical limit in research or applications where SSP is operated with homology-reduced datasets.

4.4. On the SSP Accuracies for Proteins of Different Structural Classes or Sizes and Residues with Different Physical Properties

To see how current SSP methods perform for different types of proteins and examine in detail our proposal that estimates made by structural and sequence alignments respectively stand for the theoretical and practical limits of SSP, we observed the accuracies of state-of-the-art SSP methods on different protein structural classes or sizes. The SCOP720 query set and UniRef90-2015 reference set were used. Figure 8 and Figure S1 show that the accuracies of all tested methods were lower than the limits estimated by PSI-BLAST sequence alignment. All methods performed relatively low for all-beta proteins. Previous studies observed that among helices, strands (β-ladders and β-bridges), and coils, the accuracy of SSP methods is usually lowest in strands [1,88], perhaps due to difficulty in catching long-range interactions of residues forming β-sheet strands [1,48]. By estimating the upper limits of SSP for every eight-state SSE and examining for each SSE the prediction accuracy of state-of-the-art SSP methods, we clearly saw that, compared with helices, the performance of current methods is low for β-ladders and extremely low for β-bridges (see Data S1 for figures and Table S7 for detailed data). Since all-beta proteins mainly comprise β-strands, this class’s low accuracy seemed reasonable. Moreover, the Q3/8 of these methods were high for all-alpha proteins, but their SOVs were significantly lower than Qs. The difference between Q and SOV for all-alpha proteins was the largest among these classes. Considering SOV was designed to assess SSP by secondary structural segments [62,66], there seemed to be some weakness of current SSP methods in predicting the pattern of segments for proteins composed mainly of helices. The monotonic SSE composition of such proteins might be where the difficulty lies. For instance, between an actual SSE segment HHHCHHHCHHHC and a predicted HHHHHHHCHHHC, with the difference of just one position, the Q3 is 91.7% while the SOV3 is only 70.2%. Unlike structural classes, the SSP methods’ performances for different sizes were even, although the accuracy slightly decreased as the size increased. These results suggested that overcoming the difficulties brought about by long-range interacting β-sheets and monotonic helices codes shall be the key to improving future SSP algorithms.
In addition to protein-level classifications, we have also classified the aligned residues between SCOP family homologs according to some physical properties, for instance, solvent-accessible residues (≥25% RSA, relative solvent accessibility) vs. inaccessible ones (otherwise) based on the RSA calculated by the NACCESS [98], and residues having high vs. low B-factors by the normalization and threshold applied in [99]. As shown in Data S1, we observed that solvent-inaccessible residues (buried region) have higher SSP limits than accessible ones, and residues with low B-factors (relatively inflexible) have higher SSP limits than those with high B-factors. Interestingly, judging by the estimated upper limits, the average accuracy of modern SSP algorithms is relatively lower for solvent inaccessible residues or residues with low B-factors than for solvent-accessible or flexible residues, revealing the direction toward improvements for future SSP algorithms.

4.5. Light Shed from the Discovery of This Study

Although structure alignments drew a higher upper bound on SSP accuracy than sequence alignment, the difference between the limits estimated by them was not significant until the identity of homologs became lower than 30–35% (Figure 3), implying that the alignment algorithm’s capability of determining the residue-residue equivalence between homologs with low sequence identity would be critical to improving SSP. This possibility might be supported by the fact that, for low-identity homologs, the secondary structural consistency measured by Water [80] was higher than that by PSI-BLAST or the classic BLAST (Figure 2 and Table S3). Water is a full implementation of the Smith-Waterman dynamic programming [80,82], whereas (PSI-)BLAST are heuristic implementations much faster than the Smith-Waterman algorithm without guarantee of finding the optimal alignment [79]. To further test this possibility of SSP improvement, we repeated the sequence-alignment-based experiment of Figure 2 using (PSI-)BLAST, which provided a word size parameter to control the nucleation of alignments. A short word size could find more exact local matches between sequences, which might increase the amount of correctly detected residue-residue equivalences between structural homologs. The (PSI-)BLAST would most mimic the full Smith-Waterman dynamic programming with the word size 1 (default = 3). By decreasing the word size, we found the secondary structural consistencies reported by (PSI-)BLAST all raised (Figure 9), especially when the identity between homologs was low. These results supported that the accuracy of future SSP methods could be enhanced by improving the sequence alignment engine in determining the residue-residue equivalence between low-identity homologs. Drawbacks of this idea may include (1) a small word size increases the time cost, and (2) the shortest word size 1 is accepted only by the classic BLAST running in the pairwise mode but not PSI-BLAST nor BLAST in the batch searching mode, without which the alignment search against a large dataset like UniRef90 will be extremely time-consuming. Reprogramming accurate pairwise alignment methods like Water to support batch searching or modifying PSI-BLAST to accept the shortest word size may be part of the solution.
The difference between the secondary structural consistency measured by sequence and structure alignments implied that the nature of the methodology might restrict the operational limit of SSP. By “nature of the methodology,” we mean that the core procedure of current SSP methods is to generate the PSSM-based feature set by sequence alignment rather than structural alignment because the structure of the query is unknown. The ways out of this restriction may include, (1) updating the fundamental SSP feature set, for the PSSM-based set has lasted for 20 years without significant modifications (see [94] for a feature set proposed recently), (2) enhancing the balance of SSEs in the training data, for we observed that the SSEs with low occurrence were barely or never predicted by current SSP method, such as β-bridges and π-helices (Data S1 and Table S7; see also [100] for similar observations), (3) redesigning the algorithm to be staged, for instance, by predicting the structural class of the query in the first stage and making the final prediction with a tailored model established for the class, (4) refining per-residue SSE predictions according to predicted physical properties of residues like solvent accessibility and backbone flexibility (see Data S1 for supporting data), and (5) introducing sequence-structure relationship information generated by statistical techniques (e.g., the sequence structural entropy [101], and the per-residue flexibility encoded by local structural alphabets [100,102]) into the methodology to overcome the sequence-only scenario.

5. Conclusions

This study applied PDB and SCOP data, sequence and structural alignments, and random pairings to estimate the upper and lower limits of SSP accuracy, as summarized in Figure 10. Estimated by sequence alignments between homologous structures, the upper limit of Q3 is 87.0 ± 1.2 (%), agreed with previous studies. The upper limit observed in this situation is supposed to be the practical limit of the current SSP methodology, which works based on sequence homology encoded in the PSSM. Estimated by structural alignments between homologous structures, the upper limit of Q3 is 91.4 ± 0.8 (%). The upper limit obtained through this procedure is supposed the ultimate theoretical limit of SSP. Besides, the practical and theoretical limit of Q8 was estimated to be 78.8 ± 1.6 (%) and 85.0 ± 0.9 (%), respectively. The Q3 upper limit is ~4% higher than previously expected, indicating that developing SSP algorithms is still challenging. The Q8 accuracy of state-of-the-art SSP methods is ≤77%; hence, there is still 7–9% space to develop eight-state SSP algorithms. The gap between the practical and theoretical limits may be resolvable, and several proposals for breaking through the practical limits have been made in this report. We hope that these accuracy estimates can help set clear goals for future algorithm developments, and our discoveries will stimulate the SSP field and benefit all research and applications where protein secondary structure prediction plays a key role.

Supplementary Materials

The following are available online at https://www.mdpi.com/article/10.3390/biom11111627/s1. Figure S1: The eight-state prediction accuracy of state-of-the-art SSP methods for proteins of different structural classes or sizes, Table S1: The version and settings of PSI-BLAST for the tested SSP methods, Table S2: Secondary structural consistency measured by sequence and structural alignments for homologous structures identified from the PDB, Table S3: Secondary structural consistency measured by sequence and structural alignments for SCOP-verified homologous protein structures sharing different sequence identities, Table S4: Secondary structural consistency for homologous protein structures of different structural classes, Table S5: Secondary structural consistency for homologous protein structures of different sizes, Table S6: Results of random and monotonic SSPs, Table S7: Consistency of each three- and eight-state secondary structure element between homologous protein structures, Dataset S1: The nrPDB30-2015 dataset, Dataset S2: The PDB-2015 dataset, Dataset S3: The SCOP classification and entities of the SCOP-2.07 dataset, Data S1: Additional supporting data.

Author Contributions

C.-T.H. performed experiments and data analyses and helped draft the manuscript. Y.-W.H. performed experiments, participated in data analyses, and helped prepare software packages. T.-R.C. prepared software packages and established the experimental datasets. C.-H.L. performed experiments and helped proofread the manuscript. W.-C.L. conceived of the study, prepared experimental datasets, performed experiments, analyzed results, and was the major contributor in writing the manuscript. All authors have read and agreed to the published version of the manuscript.

Funding

This work was funded by the Ministry of Science and Technology (MOST), Taiwan, with grant number NSC 101-2311-B-009-006-MY2 to WCL. The funders had no role in the design of the study, the collection, analysis, and interpretation of data, or in writing the manuscript.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data generated or analyzed during this study are included in this published article and its supplementary information files.

Acknowledgments

We would like to thank Sheng-Hung Juan, a student of WCL, for helping with the literature survey and dataset preparations. The progress of this work was greatly accelerated due to the computing power offered by Jinn-Moon Yang and Jenn-Kang Hwang at National Chiao Tung University and Ping-Chiang Lyu at National Tsing Hua University, Taiwan.

Conflicts of Interest

The authors declare no conflict of interest.

Appendix A

Figure A1. A general schema of the methodology. Sequence and structural alignments between protein structural homologs were performed to estimate the upper limits of SSP accuracy, while random pairing between protein structures was conducted to draw the lower bound of SSP. Note that only aligned residues were considered in estimating SSP limits. The effects of disordered regions, which were likely presented in the unaligned segments, might thus be ignored.
Figure A1. A general schema of the methodology. Sequence and structural alignments between protein structural homologs were performed to estimate the upper limits of SSP accuracy, while random pairing between protein structures was conducted to draw the lower bound of SSP. Note that only aligned residues were considered in estimating SSP limits. The effects of disordered regions, which were likely presented in the unaligned segments, might thus be ignored.
Biomolecules 11 01627 g0a1

References

  1. Yang, Y.; Gao, J.; Wang, J.; Heffernan, R.; Hanson, J.; Paliwal, K.; Zhou, Y. Sixty-five years of the long march in protein secondary structure prediction: The final stretch? Brief. Bioinform. 2018, 19, 482–494. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  2. Li, B.; Krishnan, V.G.; Mort, M.E.; Xin, F.; Kamati, K.K.; Cooper, D.N.; Mooney, S.D.; Radivojac, P. Automated inference of molecular mechanisms of disease from amino acid substitutions. Bioinformatics 2009, 25, 2744–2750. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  3. Folkman, L.; Yang, Y.; Li, Z.; Stantic, B.; Sattar, A.; Mort, M.; Cooper, D.N.; Liu, Y.; Zhou, Y. DDIG-in: Detecting disease-causing genetic variations due to frameshifting indels and nonsense mutations employing sequence and structural properties at nucleotide and protein levels. Bioinformatics 2015, 31, 1599–1606. [Google Scholar] [CrossRef] [PubMed]
  4. Zhao, H.; Yang, Y.; Lin, H.; Zhang, X.; Mort, M.; Cooper, D.N.; Liu, Y.; Zhou, Y. DDIG-in: Discriminating between disease-associated and neutral non-frameshifting micro-indels. Genome Biol. 2013, 14, R23. [Google Scholar] [CrossRef] [Green Version]
  5. Do, C.B.; Mahabhashyam, M.S.; Brudno, M.; Batzoglou, S. ProbCons: Probabilistic consistency-based multiple sequence alignment. Genome Res. 2005, 15, 330–340. [Google Scholar] [CrossRef] [Green Version]
  6. Pei, J.M.; Kim, B.H.; Grishin, N.V. PROMALS3D: A tool for multiple protein sequence and structure alignments. Nucleic Acids Res. 2008, 36, 2295–2300. [Google Scholar] [CrossRef]
  7. Soding, J.; Biegert, A.; Lupas, A.N. The HHpred interactive server for protein homology detection and structure prediction. Nucleic Acids Res. 2005, 33, W244–W248. [Google Scholar] [CrossRef] [Green Version]
  8. Cuthbertson, L.; Mainprize, I.L.; Naismith, J.H.; Whitfield, C. Pivotal roles of the outer membrane polysaccharide export and polysaccharide copolymerase protein families in export of extracellular polysaccharides in gram-negative bacteria. Microbiol. Mol. Biol. Rev. 2009, 73, 155–177. [Google Scholar] [CrossRef] [Green Version]
  9. Ambrosi, C.; Gassmann, O.; Pranskevich, J.N.; Boassa, D.; Smock, A.; Wang, J.; Dahl, G.; Steinem, C.; Sosinsky, G.E. Pannexin1 and Pannexin2 channels show quaternary similarities to connexons and different oligomerization numbers from each other. J. Biol. Chem. 2010, 285, 24420–24431. [Google Scholar] [CrossRef] [Green Version]
  10. Makarova, K.S.; Aravind, L.; Wolf, Y.I.; Koonin, E.V. Unification of Cas protein families and a simple scenario for the origin and evolution of CRISPR-Cas systems. Biol. Direct 2011, 6, 38. [Google Scholar] [CrossRef] [Green Version]
  11. Kifer, I.; Nussinov, R.; Wolfson, H.J. Constructing templates for protein structure prediction by simulation of protein folding pathways. Proteins 2008, 73, 380–394. [Google Scholar] [CrossRef] [Green Version]
  12. Nalini, V.; Bax, B.; Driessen, H.; Moss, D.S.; Lindley, P.F.; Slingsby, C. Close packing of an oligomeric eye lens beta-crystallin induces loss of symmetry and ordering of sequence extensions. J. Mol. Biol. 1994, 236, 1250–1258. [Google Scholar] [CrossRef]
  13. Song, J.N.; Tan, H.; Perry, A.J.; Akutsu, T.; Webb, G.I.; Whisstock, J.C.; Pike, R.N. PROSPER: An Integrated Feature-Based Tool for Predicting Protease Substrate Cleavage Sites. PLoS ONE 2012, 7, e50300. [Google Scholar] [CrossRef] [Green Version]
  14. Song, J.N.; Tan, H.; Shen, H.B.; Mahmood, K.; Boyd, S.E.; Webb, G.I.; Akutsu, T.; Whisstock, J.C. Cascleave: Towards more accurate prediction of caspase substrate cleavage sites. Bioinformatics 2010, 26, 752–760. [Google Scholar] [CrossRef] [Green Version]
  15. Iwakura, M.; Nakamura, T.; Yamane, C.; Maki, K. Systematic circular permutation of an entire protein reveals essential folding elements. Nat. Struct. Biol. 2000, 7, 580–585. [Google Scholar] [CrossRef]
  16. Wright, G.; Basak, A.K.; Wieligmann, K.; Mayr, E.M.; Slingsby, C. Circular permutation of betaB2-crystallin changes the hierarchy of domain assembly. Protein Sci 1998, 7, 1280–1285. [Google Scholar] [CrossRef]
  17. Fiser, A. Template-based protein structure modeling. Methods Mol. Biol. 2010, 673, 73–94. [Google Scholar] [CrossRef] [Green Version]
  18. Madhusudhan, M.S.; Marti-Renom, M.A.; Sanchez, R.; Sali, A. Variable gap penalty for protein sequence-structure alignment. Protein Eng. Des. Sel. 2006, 19, 129–133. [Google Scholar] [CrossRef] [Green Version]
  19. Vakser, I.A. Protein-protein docking: From interaction to interactome. Biophys. J. 2014, 107, 1785–1793. [Google Scholar] [CrossRef] [Green Version]
  20. Lee, Y.Z.; Lo, W.C.; Sue, S.C. Computational Prediction of New Intein Split Sites. Methods Mol. Biol. 2017, 1495, 259–268. [Google Scholar] [CrossRef]
  21. Lo, W.C.; Wang, L.F.; Liu, Y.Y.; Dai, T.; Hwang, J.K.; Lyu, P.C. CPred: A web server for predicting viable circular permutations in proteins. Nucleic Acids Res. 2012, 40, W232–W237. [Google Scholar] [CrossRef] [Green Version]
  22. Lo, W.C.; Dai, T.; Liu, Y.Y.; Wang, L.F.; Hwang, J.K.; Lyu, P.C. Deciphering the preference and predicting the viability of circular permutations in proteins. PLoS ONE 2012, 7, e31791. [Google Scholar] [CrossRef] [Green Version]
  23. Lee, Y.T.; Su, T.H.; Lo, W.C.; Lyu, P.C.; Sue, S.C. Circular permutation prediction reveals a viable backbone disconnection for split proteins: An approach in identifying a new functional split intein. PLoS ONE 2012, 7, e43820. [Google Scholar] [CrossRef] [Green Version]
  24. Pellequer, J.L.; Westhof, E.; Vanregenmortel, M.H.V. Correlation between the Location of Antigenic Sites and the Prediction of Turns in Proteins. Immunol. Lett. 1993, 36, 83–100. [Google Scholar] [CrossRef]
  25. Gao, J.; Faraggi, E.; Zhou, Y.; Ruan, J.; Kurgan, L. BEST: Improved prediction of B-cell epitopes from antigen sequences. PLoS ONE 2012, 7, e40104. [Google Scholar] [CrossRef] [Green Version]
  26. Li, Y.; Liu, X.; Zhu, Y.; Zhou, X.; Cao, C.; Hu, X.; Ma, H.; Wen, H.; Ma, X.; Ding, J.B. Bioinformatic prediction of epitopes in the Emy162 antigen of Echinococcus multilocularis. Exp. Ther. Med. 2013, 6, 335–340. [Google Scholar] [CrossRef] [Green Version]
  27. Zhou, H.X.; Shan, Y. Prediction of protein interaction sites from sequence profile and residue neighbor list. Proteins 2001, 44, 336–343. [Google Scholar] [CrossRef]
  28. Mukherjee, S.; Zhang, Y. Protein-Protein Complex Structure Predictions by Multimeric Threading and Template Recombination. Structure 2011, 19, 955–966. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  29. Ward, J.J.; Sodhi, J.S.; McGuffin, L.J.; Buxton, B.F.; Jones, D.T. Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J. Mol. Biol. 2004, 337, 635–645. [Google Scholar] [CrossRef]
  30. Deng, X.; Eickholt, J.; Cheng, J.L. PreDisorder: Ab initio sequence-based prediction of protein disordered regions. BMC Bioinform. 2009, 10, 436. [Google Scholar] [CrossRef] [Green Version]
  31. Xue, B.; Dunbrack, R.L.; Williams, R.W.; Dunker, A.K.; Uversky, V.N. PONDR-FIT: A meta-predictor of intrinsically disordered amino acids. Biochim. Biophys. Acta 2010, 1804, 996–1010. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Zhang, T.; Faraggi, E.; Xue, B.; Dunker, A.K.; Uversky, V.N.; Zhou, Y.Q. SPINE-D: Accurate Prediction of Short and Long Disordered Regions by a Single Neural-Network Based Method. J. Biomol. Struct. Dyn. 2012, 29, 799–813. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  33. Tardif, M.; Atteia, A.; Specht, M.; Cogne, G.; Rolland, N.; Brugiere, S.; Hippler, M.; Ferro, M.; Bruley, C.; Peltier, G.; et al. PredAlgo: A New Subcellular Localization Prediction Tool Dedicated to Green Algae. Mol. Biol. Evo.L 2012, 29, 3625–3639. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  34. Yu, C.S.; Chen, Y.C.; Lu, C.H.; Hwang, J.K. Prediction of protein subcellular localization. Proteins 2006, 64, 643–651. [Google Scholar] [CrossRef]
  35. Li, Z.X.; Yang, Y.D.; Zhan, J.; Dai, L.; Zhou, Y.Q. Energy Functions in De Novo Protein Design: Current Challenges and Future Prospects. Annu Rev. Biophys. 2013, 42, 315–335. [Google Scholar] [CrossRef] [Green Version]
  36. Xiong, P.; Wang, M.; Zhou, X.; Zhang, T.; Zhang, J.; Chen, Q.; Liu, H. Protein design with a comprehensive statistical energy function and boosted by experimental selection for foldability. Nat. Commun. 2014, 5, 5330. [Google Scholar] [CrossRef]
  37. Gebhard, L.G.; Risso, V.A.; Santos, J.; Ferreyra, R.G.; Noguera, M.E.; Ermacora, M.R. Mapping the distribution of conformational information throughout a protein sequence. J. Mol. Biol. 2006, 358, 280–288. [Google Scholar] [CrossRef]
  38. Michalsky, E.; Goede, A.; Preissner, R. Loops In Proteins (LIP)—A comprehensive loop database for homology modelling. Protein Eng. 2003, 16, 979–985. [Google Scholar] [CrossRef]
  39. Hu, X.; Wang, H.; Ke, H.; Kuhlman, B. High-resolution design of a protein loop. Proc. Natl. Acad. Sci. USA 2007, 104, 17668–17673. [Google Scholar] [CrossRef] [Green Version]
  40. Kabsch, W.; Sander, C. Dictionary of protein secondary structure: Pattern recognition of hydrogen-bonded and geometrical features. Biopolymers 1983, 22, 2577–2637. [Google Scholar] [CrossRef]
  41. Pauling, L.; Corey, R.B. Configurations of Polypeptide Chains with Favored Orientations Around Single Bonds: Two New Pleated Sheets. Proc. Natl. Acad. Sci. USA 1951, 37, 729–740. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  42. Pauling, L.; Corey, R.B.; Branson, H.R. The structure of proteins; two hydrogen-bonded helical configurations of the polypeptide chain. Proc. Natl. Acad. Sci. USA 1951, 37, 205–211. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  43. Chou, P.Y.; Fasman, G.D. Prediction of protein conformation. Biochemistry 1974, 13, 222–245. [Google Scholar] [CrossRef] [PubMed]
  44. Garnier, J.; Osguthorpe, D.J.; Robson, B. Analysis of the accuracy and implications of simple methods for predicting the secondary structure of globular proteins. J. Mol. Biol. 1978, 120, 97–120. [Google Scholar] [CrossRef]
  45. Rost, B.; Sander, C. Improved Prediction of Protein Secondary Structure by Use of Sequence Profiles and Neural Networks. Proc. Natl. Acad. Sci. USA 1993, 90, 7558–7562. [Google Scholar] [CrossRef] [Green Version]
  46. Altschul, S.F.; Madden, T.L.; Schaffer, A.A.; Zhang, J.H.; Zhang, Z.; Miller, W.; Lipman, D.J. Gapped BLAST and PSI-BLAST: A new generation of protein database search programs. Nucleic Acids Res. 1997, 25, 3389–3402. [Google Scholar] [CrossRef] [Green Version]
  47. Jones, D.T. Protein secondary structure prediction based on position-specific scoring matrices. J. Mol. Biol. 1999, 292, 195–202. [Google Scholar] [CrossRef] [Green Version]
  48. Rost, B. Review: Protein secondary structure prediction continues to rise. J. Struct. Biol. 2001, 134, 204–218. [Google Scholar] [CrossRef]
  49. Pollastri, G.; Przybylski, D.; Rost, B.; Baldi, P. Improving the prediction of protein secondary structure in three and eight classes using recurrent neural networks and profiles. Proteins 2002, 47, 228–235. [Google Scholar] [CrossRef]
  50. Dor, O.; Zhou, Y. Achieving 80% ten-fold cross-validated accuracy for secondary structure prediction by large-scale training. Proteins 2007, 66, 838–845. [Google Scholar] [CrossRef]
  51. Cole, C.; Barber, J.D.; Barton, G.J. The Jpred 3 secondary structure prediction server. Nucleic Acids Res. 2008, 36, W197–W201. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  52. Mirabello, C.; Pollastri, G. Porter, PaleAle 4.0: High-accuracy prediction of protein secondary structure and relative solvent accessibility. Bioinformatics 2013, 29, 2056–2058. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  53. Heffernan, R.; Yang, Y.; Paliwal, K.; Zhou, Y. Capturing non-local interactions by long short-term memory bidirectional recurrent neural networks for improving prediction of protein secondary structure, backbone angles, contact numbers and solvent accessibility. Bioinformatics 2017, 33, 2842–2849. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  54. Hanson, J.; Paliwal, K.; Litfin, T.; Yang, Y.; Zhou, Y. Improving prediction of protein secondary structure, backbone angles, solvent accessibility and contact numbers by using predicted contact maps and an ensemble of recurrent and residual convolutional neural networks. Bioinformatics 2019, 35, 2403–2410. [Google Scholar] [CrossRef]
  55. Wang, S.; Peng, J.; Ma, J.Z.; Xu, J.B. Protein Secondary Structure Prediction Using Deep Convolutional Neural Fields. Sci. Rep.-UK 2016, 6, 18962. [Google Scholar] [CrossRef] [Green Version]
  56. Torrisi, M.; Kaleel, M.; Pollastri, G. Deeper Profiles and Cascaded Recurrent and Convolutional Neural Networks for state-of-the-art Protein Secondary Structure Prediction. Sci. Rep. 2019, 9, 12374. [Google Scholar] [CrossRef] [Green Version]
  57. Heffernan, R.; Paliwal, K.; Lyons, J.; Singh, J.; Yang, Y.; Zhou, Y. Single-sequence-based prediction of protein secondary structures and solvent accessibility by deep whole-sequence learning. J. Comput. Chem. 2018, 39, 2210–2216. [Google Scholar] [CrossRef]
  58. Fang, C.; Shang, Y.; Xu, D. MUFOLD-SS: New deep inception-inside-inception networks for protein secondary structure prediction. Proteins 2018, 86, 592–598. [Google Scholar] [CrossRef]
  59. Klausen, M.S.; Jespersen, M.C.; Nielsen, H.; Jensen, K.K.; Jurtz, V.I.; Sonderby, C.K.; Sommer, M.O.A.; Winther, O.; Nielsen, M.; Petersen, B.; et al. NetSurfP-2.0: Improved prediction of protein structural features by integrated deep learning. Proteins 2019, 87, 520–527. [Google Scholar] [CrossRef] [Green Version]
  60. Zhou, J.; Wang, H.; Zhao, Z.; Xu, R.; Lu, Q. CNNH_PSS: Protein 8-class secondary structure prediction by convolutional neural network with highway. BMC Bioinform. 2018, 19, 60. [Google Scholar] [CrossRef] [Green Version]
  61. Levin, J.M.; Pascarella, S.; Argos, P.; Garnier, J. Quantification of secondary structure prediction improvement using multiple alignments. Protein Eng. 1993, 6, 849–854. [Google Scholar] [CrossRef]
  62. Rost, B.; Sander, C.; Schneider, R. Redefining the goals of protein secondary structure prediction. J. Mol. Biol. 1994, 235, 13–26. [Google Scholar] [CrossRef]
  63. Zhang, W.; Dunker, A.K.; Zhou, Y. Assessing secondary structure assignment of protein structures by using pairwise sequence-alignment benchmarks. Proteins 2008, 71, 61–67. [Google Scholar] [CrossRef]
  64. Kuziemko, A.; Honig, B.; Petrey, D. Using structure to explore the sequence alignment space of remote homologs. PLoS Comput. Biol. 2011, 7, e1002175. [Google Scholar] [CrossRef]
  65. Pascarella, S.; Argos, P. A data bank merging related protein structures and sequences. Protein Eng. 1992, 5, 121–137. [Google Scholar] [CrossRef]
  66. Zemla, A.; Venclovas, C.; Fidelis, K.; Rost, B. A modified definition of Sov, a segment-based measure for protein secondary structure prediction assessment. Proteins 1999, 34, 220–223. [Google Scholar] [CrossRef]
  67. Murzin, A.G.; Brenner, S.E.; Hubbard, T.; Chothia, C. SCOP: A structural classification of proteins database for the investigation of sequences and structures. J. Mol. Biol. 1995, 247, 536–540. [Google Scholar] [CrossRef]
  68. Madej, T.; Lanczycki, C.J.; Zhang, D.; Thiessen, P.A.; Geer, R.C.; Marchler-Bauer, A.; Bryant, S.H. MMDB and VAST+: Tracking structural similarities between macromolecular complexes. Nucleic Acids Res. 2014, 42, D297–D303. [Google Scholar] [CrossRef] [Green Version]
  69. NCBI nr-PDB: Non-Redundant PDB Data Set for VAST. Available online: https://www.ncbi.nlm.nih.gov/Structure/VAST/vast.shtml (accessed on 21 September 2021).
  70. Fu, L.; Niu, B.; Zhu, Z.; Wu, S.; Li, W. CD-HIT: Accelerated for clustering the next-generation sequencing data. Bioinformatics 2012, 28, 3150–3152. [Google Scholar] [CrossRef]
  71. Edgar, R.C. Search and clustering orders of magnitude faster than BLAST. Bioinformatics 2010, 26, 2460–2461. [Google Scholar] [CrossRef] [Green Version]
  72. Steinegger, M.; Soding, J. MMseqs2 enables sensitive protein sequence searching for the analysis of massive data sets. Nat. Biotechnol. 2017, 35, 1026–1028. [Google Scholar] [CrossRef] [Green Version]
  73. Fox, N.K.; Brenner, S.E.; Chandonia, J.M. SCOPe: Structural Classification of Proteins--extended, integrating SCOP and ASTRAL data and classification of new structures. Nucleic Acids Res. 2014, 42, D304–D309. [Google Scholar] [CrossRef]
  74. UniProt, C. Activities at the Universal Protein Resource (UniProt). Nucleic Acids Res. 2014, 42, D191–D198. [Google Scholar] [CrossRef] [Green Version]
  75. Touw, W.G.; Baakman, C.; Black, J.; te Beek, T.A.; Krieger, E.; Joosten, R.P.; Vriend, G. A series of PDB-related databanks for everyday needs. Nucleic Acids Res. 2015, 43, D364–D368. [Google Scholar] [CrossRef]
  76. Zhu, J.; Weng, Z. FAST: A novel protein structure alignment algorithm. Proteins 2005, 58, 618–627. [Google Scholar] [CrossRef] [Green Version]
  77. Lo Conte, L.; Ailey, B.; Hubbard, T.J.; Brenner, S.E.; Murzin, A.G.; Chothia, C. SCOP: A structural classification of proteins database. Nucleic Acids Res. 2000, 28, 257–259. [Google Scholar] [CrossRef] [Green Version]
  78. Lo, W.C.; Lee, C.Y.; Lee, C.C.; Lyu, P.C. iSARST: An integrated SARST web server for rapid protein structural similarity searches. Nucleic Acids Res. 2009, 37, W545–W551. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  79. Altschul, S.F.; Gish, W.; Miller, W.; Myers, E.W.; Lipman, D.J. Basic local alignment search tool. J. Mol. Biol. 1990, 215, 403–410. [Google Scholar] [CrossRef]
  80. Rice, P.; Longden, I.; Bleasby, A. EMBOSS: The European Molecular Biology Open Software Suite. Trends Genet. TIG 2000, 16, 276–277. [Google Scholar] [CrossRef]
  81. Myers, E.W.; Miller, W. Optimal alignments in linear space. Comput. Appl. Biosci. 1988, 4, 11–17. [Google Scholar] [CrossRef] [PubMed]
  82. Smith, T.F.; Waterman, M.S. Identification of common molecular subsequences. J. Mol. Biol. 1981, 147, 195–197. [Google Scholar] [CrossRef]
  83. Needleman, S.B.; Wunsch, C.D. A general method applicable to the search for similarities in the amino acid sequence of two proteins. J. Mol. Biol. 1970, 48, 443–453. [Google Scholar] [CrossRef]
  84. Zhang, Y.; Skolnick, J. TM-align: A protein structure alignment algorithm based on the TM-score. Nucleic Acids Res. 2005, 33, 2302–2309. [Google Scholar] [CrossRef]
  85. Lo, W.C.; Huang, P.J.; Chang, C.H.; Lyu, P.C. Protein structural similarity search by Ramachandran codes. BMC Bioinform. 2007, 8, 307. [Google Scholar] [CrossRef] [Green Version]
  86. Faraggi, E.; Zhang, T.; Yang, Y.; Kurgan, L.; Zhou, Y. SPINE X: Improving protein secondary structure prediction by multistep learning coupled with prediction of solvent accessible surface area and backbone torsion angles. J. Comput. Chem. 2012, 33, 259–267. [Google Scholar] [CrossRef] [Green Version]
  87. Yaseen, A.; Li, Y. Context-based features enhance protein secondary structure prediction accuracy. J. Chem. Inf. Model. 2014, 54, 992–1002. [Google Scholar] [CrossRef]
  88. Heffernan, R.; Paliwal, K.; Lyons, J.; Dehzangi, A.; Sharma, A.; Wang, J.; Sattar, A.; Yang, Y.; Zhou, Y. Improving prediction of secondary structure, local backbone angles, and solvent accessible surface area of proteins by iterative deep learning. Sci. Rep. 2015, 5, 11476. [Google Scholar] [CrossRef] [Green Version]
  89. Wang, Z.; Zhao, F.; Peng, J.; Xu, J. Protein 8-class secondary structure prediction using conditional neural fields. Proteomics 2011, 11, 3786–3792. [Google Scholar] [CrossRef] [Green Version]
  90. Magnan, C.N.; Baldi, P. SSpro/ACCpro 5: Almost perfect prediction of protein secondary structure and relative solvent accessibility using profiles, machine learning and structural similarity. Bioinformatics 2014, 30, 2592–2597. [Google Scholar] [CrossRef] [Green Version]
  91. Moult, J.; Fidelis, K.; Kryshtafovych, A.; Schwede, T.; Tramontano, A. Critical assessment of methods of protein structure prediction (CASP)-Round XII. Proteins 2018, 86 (Suppl. 1), 7–15. [Google Scholar] [CrossRef]
  92. Zhou, J.; Troyanskaya, O.G. Deep supervised and convolutional generative stochastic network for protein secondary structure prediction. In Proceedings of the 31st International Conference on International Conference on Machine Learning, Beijing, China, 21–26 June 2014; Volume 32, pp. 745–753. [Google Scholar]
  93. Zhang, B.; Li, J.; Lu, Q. Prediction of 8-state protein secondary structures by a novel deep learning architecture. BMC Bioinform. 2018, 19, 1–13. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  94. Chen, T.R.; Juan, S.H.; Huang, Y.W.; Lin, Y.C.; Lo, W.C. A secondary structure-based position-specific scoring matrix applied to the improvement in protein secondary structure prediction. PLoS ONE 2021, 16, e0255076. [Google Scholar] [CrossRef]
  95. Chen, T.R.; Lo, C.H.; Juan, S.H.; Lo, W.C. The influence of dataset homology and a rigorous evaluation strategy on protein secondary structure prediction. PLoS ONE 2021, 16, e0254555. [Google Scholar] [CrossRef]
  96. Wilson, C.A.; Kreychman, J.; Gerstein, M. Assessing annotation transfer for genomics: Quantifying the relations between protein sequence, structure and function through traditional and probabilistic scores. J. Mol. Biol. 2000, 297, 233–249. [Google Scholar] [CrossRef] [Green Version]
  97. Hobohm, U.; Scharf, M.; Schneider, R.; Sander, C. Selection of representative protein data sets. Protein Sci. 1992, 1, 409–417. [Google Scholar] [CrossRef] [Green Version]
  98. Hubbard, S.J.; Thornton, J.M. NACCESS V2.1.1. Computer Program, Department of Biochemistry and Molecular Biology, University College London. 1993. Available online: http://www.bioinf.manchester.ac.uk/naccess/ (accessed on 21 September 2021).
  99. Vander Meersche, Y.; Cretin, G.; de Brevern, A.G.; Gelly, J.C.; Galochkina, T. MEDUSA: Prediction of Protein Flexibility from Sequence. J. Mol. Biol. 2021, 433, 166882. [Google Scholar] [CrossRef]
  100. de Brevern, A.G. Impact of protein dynamics on secondary structure prediction. Biochimie 2020, 179, 14–22. [Google Scholar] [CrossRef]
  101. Shih, C.H.; Chang, C.M.; Lin, Y.S.; Lo, W.C.; Hwang, J.K. Evolutionary information hidden in a single protein structure. Proteins 2012, 80, 1647–1657. [Google Scholar] [CrossRef]
  102. Joseph, A.P.; Agarwal, G.; Mahajan, S.; Gelly, J.C.; Swapna, L.S.; Offmann, B.; Cadet, F.; Bornot, A.; Tyagi, M.; Valadie, H.; et al. A short survey on protein blocks. Biophys. Rev. 2010, 2, 137–147. [Google Scholar] [CrossRef]
Figure 1. Secondary structural consistency between homologs identified by PSI-BLAST from PDB. (a) Secondary structural consistency at different sequence identity levels. The horizontal axis indicates the sequence identity level where, for instance, 30 means that the sequence identities of homologous structures fall between ≥30% and <35%. These levels were made according to Yang et al. [1]. The consistency dropped as the identity decreased, and our 10-repeat experimental results agreed well with [1] (the solid purple curve of Q3 versus the red). Yang et al. only analyzed high resolutions (<3 Å) X-ray structures. Following the same filtering criteria, our results (dotted purple curves) were closer to theirs. For easy observation, we enlarged the vertical axis of Q3 with identities ≥50. The secondary structural consistency revealed the structural difference between homologs, which is the ultimate factor limiting SSP accuracy. For instance, at the 30% identity level, the averaged consistency is 87.8%, meaning an 11.2% difference in secondary structure, which indicates that the SSP accuracy of such homologs can be at most ~88%. Previous studies measured the secondary structure consistency by sequence alignments [1,62]. We also performed structural alignments (black lines). (b) Secondary structural consistency at different sequence identity thresholds. The horizontal axis indicates the threshold where, for instance, 30 means that the sequence identities of homologous structures are ≥30%. Previous studies estimated the limit of SSP on the basis of three-state SSEs. We also made estimations based on the eight-state SSEs.
Figure 1. Secondary structural consistency between homologs identified by PSI-BLAST from PDB. (a) Secondary structural consistency at different sequence identity levels. The horizontal axis indicates the sequence identity level where, for instance, 30 means that the sequence identities of homologous structures fall between ≥30% and <35%. These levels were made according to Yang et al. [1]. The consistency dropped as the identity decreased, and our 10-repeat experimental results agreed well with [1] (the solid purple curve of Q3 versus the red). Yang et al. only analyzed high resolutions (<3 Å) X-ray structures. Following the same filtering criteria, our results (dotted purple curves) were closer to theirs. For easy observation, we enlarged the vertical axis of Q3 with identities ≥50. The secondary structural consistency revealed the structural difference between homologs, which is the ultimate factor limiting SSP accuracy. For instance, at the 30% identity level, the averaged consistency is 87.8%, meaning an 11.2% difference in secondary structure, which indicates that the SSP accuracy of such homologs can be at most ~88%. Previous studies measured the secondary structure consistency by sequence alignments [1,62]. We also performed structural alignments (black lines). (b) Secondary structural consistency at different sequence identity thresholds. The horizontal axis indicates the threshold where, for instance, 30 means that the sequence identities of homologous structures are ≥30%. Previous studies estimated the limit of SSP on the basis of three-state SSEs. We also made estimations based on the eight-state SSEs.
Biomolecules 11 01627 g001
Figure 2. Secondary structural consistency measured by sequence alignments between structural homologs verified by SCOP. (a) Three-state secondary structural consistency (Q3) at various sequence identity levels. Using the structural homologs determined by SCOP to repeat the experiment of previous studies, the Q3 data of homologs sharing ≥30% identity were very similar to previous results. The purple curve of PSI-BLAST and PDB is obtained from Figure 1a for comparison. The Q3 differed significantly at low identities, with SCOP-determined homologs being much higher than PSI-BLAST-identified ones. Note that the difference was mainly caused by how homologs were defined (manual curations vs. programmatic searches) rather than by the source of data (SCOP vs. PDB). The Q3 was measured by PSI-BLAST, classic BLAST, Water, and Stretcher. The classic BLAST was omitted in this figure for clarity because its curve lay closely to that of PSI-BLAST (see Table S3 for full data, with the classic BLAST included). (b) Eight-state secondary structural consistency (Q8) at various sequence identity levels. (c) Q3 at various sequence identity cutoffs. The horizontal axis indicates the cutoff where, for instance, 30 and 90 mean that the identities of homologous structures are lower than 30% and 90%, respectively. Since protein datasets are typically prepared as nonredundant sets of identities lower than given cutoffs in practical research, in all figures after this, we only display the consistency computed under given identity cutoffs instead of thresholds. (d) Q8 at various sequence identity cutoffs.
Figure 2. Secondary structural consistency measured by sequence alignments between structural homologs verified by SCOP. (a) Three-state secondary structural consistency (Q3) at various sequence identity levels. Using the structural homologs determined by SCOP to repeat the experiment of previous studies, the Q3 data of homologs sharing ≥30% identity were very similar to previous results. The purple curve of PSI-BLAST and PDB is obtained from Figure 1a for comparison. The Q3 differed significantly at low identities, with SCOP-determined homologs being much higher than PSI-BLAST-identified ones. Note that the difference was mainly caused by how homologs were defined (manual curations vs. programmatic searches) rather than by the source of data (SCOP vs. PDB). The Q3 was measured by PSI-BLAST, classic BLAST, Water, and Stretcher. The classic BLAST was omitted in this figure for clarity because its curve lay closely to that of PSI-BLAST (see Table S3 for full data, with the classic BLAST included). (b) Eight-state secondary structural consistency (Q8) at various sequence identity levels. (c) Q3 at various sequence identity cutoffs. The horizontal axis indicates the cutoff where, for instance, 30 and 90 mean that the identities of homologous structures are lower than 30% and 90%, respectively. Since protein datasets are typically prepared as nonredundant sets of identities lower than given cutoffs in practical research, in all figures after this, we only display the consistency computed under given identity cutoffs instead of thresholds. (d) Q8 at various sequence identity cutoffs.
Biomolecules 11 01627 g002
Figure 3. Secondary structural consistency measured by structural alignments between structural homologs determined by SCOP. (a) Three-state secondary structural consistency (Q3) at various sequence identity levels. (b) Eight-state secondary structural consistency (Q8) at various sequence identity levels. (c) Q3 at different various identity cutoffs. (d) Q8 at various sequence identity cutoffs. To make comparisons easier in these plots, the results of Water, which marked the highest Q3 and Q8 measured by sequence alignments in Figure 2, are shown again for reference. Structure alignment methods applied in this experiment were FAST, TM-align, and SARST. All the Q3/8 values obtained with them are higher than those obtained with sequence alignments, especially at low identities. Compared with Q3, the differences between structural and sequence alignments were much more significant in Q8. Structural alignments had a better capability of detecting the residue-residue equivalences between structural homologs than sequence alignments, even when they had distant evolutionary relationships. Thus, structural alignments should be more suitable than sequence alignments for estimating the limit of SSP.
Figure 3. Secondary structural consistency measured by structural alignments between structural homologs determined by SCOP. (a) Three-state secondary structural consistency (Q3) at various sequence identity levels. (b) Eight-state secondary structural consistency (Q8) at various sequence identity levels. (c) Q3 at different various identity cutoffs. (d) Q8 at various sequence identity cutoffs. To make comparisons easier in these plots, the results of Water, which marked the highest Q3 and Q8 measured by sequence alignments in Figure 2, are shown again for reference. Structure alignment methods applied in this experiment were FAST, TM-align, and SARST. All the Q3/8 values obtained with them are higher than those obtained with sequence alignments, especially at low identities. Compared with Q3, the differences between structural and sequence alignments were much more significant in Q8. Structural alignments had a better capability of detecting the residue-residue equivalences between structural homologs than sequence alignments, even when they had distant evolutionary relationships. Thus, structural alignments should be more suitable than sequence alignments for estimating the limit of SSP.
Biomolecules 11 01627 g003
Figure 4. Secondary structural consistency for structural homologs of different structural classes or sizes. (a) The secondary structural consistency for homologs belonging to different structural classes and sharing <90% identities. (b) The consistency for homologs having different sizes and sharing <90% identities. (c) The consistency for homologs belonging to different structural classes and sharing <30% identities. (d) The consistency for homologs having different sizes and sharing <30% identities. In these plots, the secondary structural consistencies computed by sequence alignments (method: PSI-BLAST) and structural alignments (method: FAST) are drawn as blue and black bars, respectively; because the latter are all higher than the former, only black caps are visible. PSI-BLAST represents sequence alignment methods because most SSP algorithms work depending on it. FAST represents structure alignment methods because it reported the highest Q3/8 between SCOP homologs (Figure 3).
Figure 4. Secondary structural consistency for structural homologs of different structural classes or sizes. (a) The secondary structural consistency for homologs belonging to different structural classes and sharing <90% identities. (b) The consistency for homologs having different sizes and sharing <90% identities. (c) The consistency for homologs belonging to different structural classes and sharing <30% identities. (d) The consistency for homologs having different sizes and sharing <30% identities. In these plots, the secondary structural consistencies computed by sequence alignments (method: PSI-BLAST) and structural alignments (method: FAST) are drawn as blue and black bars, respectively; because the latter are all higher than the former, only black caps are visible. PSI-BLAST represents sequence alignment methods because most SSP algorithms work depending on it. FAST represents structure alignment methods because it reported the highest Q3/8 between SCOP homologs (Figure 3).
Biomolecules 11 01627 g004
Figure 5. Distribution of secondary structural consistency between random protein pairs. One million rounds of random pairing were performed to compute the distribution of secondary structural consistencies between proteins to mark the lower bound of SSP accuracy. Two proteins were randomly selected from the SCOP-2.07 dataset in each round, and the secondary structural consistency between them was measured with the Q and SOV scores. The average SOV3 of these random pairs was lower than the average Q3 (SOV3 = 32% versus Q3 = 35%), different from the results of the first paper of SOV (SOV3 = 37% versus Q3 = 35%) [62], because the algorithm of SOV applied in this study was an updated version (v’99) [66], the value of which is lower than the first defined SOV (v’94; see Table I of [66] for comparisons).
Figure 5. Distribution of secondary structural consistency between random protein pairs. One million rounds of random pairing were performed to compute the distribution of secondary structural consistencies between proteins to mark the lower bound of SSP accuracy. Two proteins were randomly selected from the SCOP-2.07 dataset in each round, and the secondary structural consistency between them was measured with the Q and SOV scores. The average SOV3 of these random pairs was lower than the average Q3 (SOV3 = 32% versus Q3 = 35%), different from the results of the first paper of SOV (SOV3 = 37% versus Q3 = 35%) [62], because the algorithm of SOV applied in this study was an updated version (v’99) [66], the value of which is lower than the first defined SOV (v’94; see Table I of [66] for comparisons).
Biomolecules 11 01627 g005
Figure 6. Accuracies of state-of-the-art SSP methods evaluated with datasets sharing <90% sequence identities. This test utilized the independent query sets prepared by Yang et al. for SSP evaluations [1] and the SCOP720 query set. The source of reference sequences for PSSM generation was UniRef90-2015. The blue and black bars indicate the SSP upper limits estimated at a <90% sequence identity cutoff by sequence and structural alignments, respectively. The sequence alignment algorithm used to make those estimates was PSI-BLAST, the PSSM generator of all the tested SSP methods. These results reveal that, when the homology of the SSP query and reference sequences was <90% identity, the accuracy of current SSP methods has not reached the limit estimated by sequence alignments of homologs sharing <90% sequence identities.
Figure 6. Accuracies of state-of-the-art SSP methods evaluated with datasets sharing <90% sequence identities. This test utilized the independent query sets prepared by Yang et al. for SSP evaluations [1] and the SCOP720 query set. The source of reference sequences for PSSM generation was UniRef90-2015. The blue and black bars indicate the SSP upper limits estimated at a <90% sequence identity cutoff by sequence and structural alignments, respectively. The sequence alignment algorithm used to make those estimates was PSI-BLAST, the PSSM generator of all the tested SSP methods. These results reveal that, when the homology of the SSP query and reference sequences was <90% identity, the accuracy of current SSP methods has not reached the limit estimated by sequence alignments of homologs sharing <90% sequence identities.
Biomolecules 11 01627 g006
Figure 7. The accuracy of the current SSP methodology applied with restricted homology between developmental and operational datasets. For avoiding the side effects of sequence redundancy, which might cause information leakage and overfitting for a predictor, it is typical for computational biology studies to restrict protein datasets with some sequence homology cutoff. We speculated that the SSP limits estimated based on homologs with different homology cutoffs (blue curves, adopted from the result of PSI-BLAST in Figure 2c and Table S3) would draw the accuracy upper bound for SSP methods trained and operated with datasets where the homology of proteins is restricted. The protein materials used here were prepared so that the sequence identities within and between the query and PSSM reference sets were lower than the given cutoffs. As expected, as the homology decreased, the accuracy of SSP models trained and operated with such proteins decreased. This experiment was repeated 10 times by random sampling.
Figure 7. The accuracy of the current SSP methodology applied with restricted homology between developmental and operational datasets. For avoiding the side effects of sequence redundancy, which might cause information leakage and overfitting for a predictor, it is typical for computational biology studies to restrict protein datasets with some sequence homology cutoff. We speculated that the SSP limits estimated based on homologs with different homology cutoffs (blue curves, adopted from the result of PSI-BLAST in Figure 2c and Table S3) would draw the accuracy upper bound for SSP methods trained and operated with datasets where the homology of proteins is restricted. The protein materials used here were prepared so that the sequence identities within and between the query and PSSM reference sets were lower than the given cutoffs. As expected, as the homology decreased, the accuracy of SSP models trained and operated with such proteins decreased. This experiment was repeated 10 times by random sampling.
Biomolecules 11 01627 g007
Figure 8. Accuracy of state-of-the-art SSP methods for proteins of different (a) structural classes or (b) sizes. The upper limit of SSP accuracy estimated by sequence and structural alignments at 90% identity cutoff of homologs are respectively indicated by blue and black bars. No method exceeded the estimated limits in any class or size group, and most methods met the greatest challenge in predicting all-beta proteins. SSP accuracy decreased a little as the protein size increased. Since the same tendencies were observed in three- and eight-state predictions, this figure only displays Q3 and SOV3. See Figure S1 for Q8 and SOV8 results.
Figure 8. Accuracy of state-of-the-art SSP methods for proteins of different (a) structural classes or (b) sizes. The upper limit of SSP accuracy estimated by sequence and structural alignments at 90% identity cutoff of homologs are respectively indicated by blue and black bars. No method exceeded the estimated limits in any class or size group, and most methods met the greatest challenge in predicting all-beta proteins. SSP accuracy decreased a little as the protein size increased. Since the same tendencies were observed in three- and eight-state predictions, this figure only displays Q3 and SOV3. See Figure S1 for Q8 and SOV8 results.
Biomolecules 11 01627 g008
Figure 9. Secondary structural consistency measured by PSI-BLAST and BLAST with different word sizes. (a) Secondary structural consistencies between SCOP family homologs measured by PSI-BLAST at different identity levels. As illustrated by the color codes, word sizes 3 and 2 were applied. The default word size of PSI-BLAST was 3, but changing it to 2 increased the measured secondary structural consistency between homologs, especially at low identities. (b) Secondary structural consistencies between SCOP family homologs measured by the classic BLAST at different identity levels. The BLAST accepted word sizes 3 and 2 in the database searching mode (blastall) and word size 1 in the pairwise alignment mode (bl2seq). A small word size notably increased the measured secondary structural consistency. (c), (d) Secondary structural consistencies between SCOP family homologs measured by (PSI-)BLAST at different identity cutoffs. Results obtained according to identity cutoffs make the effects of word size more observable.
Figure 9. Secondary structural consistency measured by PSI-BLAST and BLAST with different word sizes. (a) Secondary structural consistencies between SCOP family homologs measured by PSI-BLAST at different identity levels. As illustrated by the color codes, word sizes 3 and 2 were applied. The default word size of PSI-BLAST was 3, but changing it to 2 increased the measured secondary structural consistency between homologs, especially at low identities. (b) Secondary structural consistencies between SCOP family homologs measured by the classic BLAST at different identity levels. The BLAST accepted word sizes 3 and 2 in the database searching mode (blastall) and word size 1 in the pairwise alignment mode (bl2seq). A small word size notably increased the measured secondary structural consistency. (c), (d) Secondary structural consistencies between SCOP family homologs measured by (PSI-)BLAST at different identity cutoffs. Results obtained according to identity cutoffs make the effects of word size more observable.
Biomolecules 11 01627 g009
Figure 10. Summary of the estimated limits of SSP accuracy. Both the lower limits (light blue strips) and upper limits of three-/eight-state SSP accuracies are illustrated. The practical upper limits of accuracies (blue strips to the right) were estimated by PSI-BLAST sequence alignments with the word size 2 (see Figure 9 and Table S3), and the theoretical limits (black strips to the right) were estimated by structural alignments (Figure 3 and Table S3).
Figure 10. Summary of the estimated limits of SSP accuracy. Both the lower limits (light blue strips) and upper limits of three-/eight-state SSP accuracies are illustrated. The practical upper limits of accuracies (blue strips to the right) were estimated by PSI-BLAST sequence alignments with the word size 2 (see Figure 9 and Table S3), and the theoretical limits (black strips to the right) were estimated by structural alignments (Figure 3 and Table S3).
Biomolecules 11 01627 g010
Table 1. Performance pretest results of several state-of-the-art SSP methods.
Table 1. Performance pretest results of several state-of-the-art SSP methods.
SSE SetMethodDatasets: TS115 vs. UniRef90-2015Datasets: CASP12 vs. UniRef90-2015
Three-stateMeasureQ3 (reported)Q3Q3 (reported)Q3
Scorpion 10.8170.8150.8050.823
Spider2 10.8190.8190.7980.812
SpineX 10.8010.8000.7690.783
Psipred 10.8020.8010.7800.783
DeepCNF 10.8230.8190.8210.828
RaptorX 20.8120.8070.7910.793
SSpro820.7950.7880.7760.779
Eight-stateMeasureQ8 (reported)Q8Q8 (reported)Q8
DeepCNF 10.7200.7030.7300.728
RaptorX 20.6970.6980.6510.694
SSpro8 10.6800.6710.6900.675
1 Obtained from [1], where the query datasets were TS115 and CASP12, and the PSSM reference dataset was UniRef90-2015. 2 Obtained from [55], where the query datasets were CullPDB and CASP11, and the PSSM reference dataset was UniRef90-2015.
Table 2. Three-state secondary structural consistency measured by sequence/structure alignments for homologs identified from PDB.
Table 2. Three-state secondary structural consistency measured by sequence/structure alignments for homologs identified from PDB.
Identity
(%)
Yang et al. [1]PSI-BLAST
(Avg ± Std)
FAST
(Avg ± Std)
Identity
Threshold (%)
Yang et al. [1]PSI-BLAST
(Avg ± Std)
FAST
(Avg ± Std)
10096.897.1 ± 0.397.1 ± 0.310096.897.1 ± 0.397.1 ± 0.3
90‒10096.295.6 ± 0.495.7 ± 0.4≥9096.896.1 ± 0.396.2 ± 0.4
80‒9095.193.5 ± 1.394.0 ± 1.2≥8096.695.6 ± 0.595.8 ± 0.5
70‒8095.892.5 ± 1.693.2 ± 1.3≥7096.595.4 ± 0.695.6 ± 0.5
60‒7094.792.9 ± 0.893.4 ± 0.7≥6096.395.2 ± 0.695.5 ± 0.6
55‒6094.292.6 ± 0.793.1 ± 0.8≥5596.295.0 ± 0.695.3 ± 0.6
50‒5593.192.3 ± 0.692.8 ± 0.6≥5096.094.8 ± 0.695.1 ± 0.6
45‒5092.691.8 ± 0.992.6 ± 0.8≥4595.794.6 ± 0.694.9 ± 0.6
40‒4591.090.1 ± 1.391.2 ± 1.3≥4095.194.1 ± 0.794.5 ± 0.7
35‒4089.989.6 ± 0.890.8 ± 0.9≥3594.693.6 ± 0.794.1 ± 0.7
30‒3587.787.8 ± 0.890.2 ± 0.4≥3093.492.4 ± 0.793.3 ± 0.6
25‒3084.284.2 ± 0.589.0 ± 0.4≥2590.789.1 ± 0.691.6 ± 0.5
20‒2581.181.4 ± 0.388.2 ± 0.2≥2087.086.0 ± 0.590.2 ± 0.4
15‒2075.977.5 ± 0.487.2 ± 0.2≥1584.083.8 ± 0.489.4 ± 0.3
<10‒1571.472.1 ± 0.985.7 ± 0.3≥1083.780.8 ± 0.588.5 ± 0.3
Table 3. Eight-state secondary structural consistency measured by sequence/structure alignments for homologs identified from PDB.
Table 3. Eight-state secondary structural consistency measured by sequence/structure alignments for homologs identified from PDB.
Identity
(%)
PSI-BLAST
(Avg ± Std)
FAST
(Avg ± Std)
Identity
Threshold (%)
PSI-BLAST
(Avg ± Std)
FAST
(Avg ± Std)
10094.8 ± 0.494.9 ± 0.410094.8 ± 0.494.9 ± 0.4
90‒10092.2 ± 0.792.6 ± 0.7≥9093.1 ± 0.593.4 ± 0.5
80‒9088.5 ± 2.089.6 ± 1.7≥8092.3 ± 0.792.7 ± 0.7
70‒8086.7 ± 2.888.1 ± 2.2≥7092.0 ± 0.992.4 ± 0.8
60‒7087.2 ± 1.388.4 ± 1.1≥6091.6 ± 0.992.1 ± 0.9
55‒6086.2 ± 1.887.3 ± 1.8≥5591.2 ± 1.091.7 ± 0.9
50‒5585.8 ± 0.587.0 ± 0.5≥5090.8 ± 1.091.4 ± 0.9
45‒5085.3 ± 1.286.8 ± 1.1≥4590.4 ± 1.091.1 ± 0.9
40‒4582.6 ± 2.084.6 ± 1.9≥4089.6 ± 1.190.4 ± 1.0
35‒4081.6 ± 1.283.9 ± 1.3≥3588.6 ± 1.189.6 ± 1.1
30‒3578.9 ± 1.182.7 ± 0.5≥3086.7 ± 1.188.2 ± 0.9
25‒3074.0 ± 0.780.9 ± 0.6≥2581.6 ± 0.985.3 ± 0.8
20‒2570.7 ± 0.479.7 ± 0.3≥2077.2 ± 0.783.0 ± 0.6
15‒2065.9 ± 0.478.3 ± 0.3≥1574.2 ± 0.681.8 ± 0.5
<10‒1559.7 ± 1.076.5 ± 0.5≥1070.5 ± 0.780.4 ± 0.5
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Ho, C.-T.; Huang, Y.-W.; Chen, T.-R.; Lo, C.-H.; Lo, W.-C. Discovering the Ultimate Limits of Protein Secondary Structure Prediction. Biomolecules 2021, 11, 1627. https://doi.org/10.3390/biom11111627

AMA Style

Ho C-T, Huang Y-W, Chen T-R, Lo C-H, Lo W-C. Discovering the Ultimate Limits of Protein Secondary Structure Prediction. Biomolecules. 2021; 11(11):1627. https://doi.org/10.3390/biom11111627

Chicago/Turabian Style

Ho, Chia-Tzu, Yu-Wei Huang, Teng-Ruei Chen, Chia-Hua Lo, and Wei-Cheng Lo. 2021. "Discovering the Ultimate Limits of Protein Secondary Structure Prediction" Biomolecules 11, no. 11: 1627. https://doi.org/10.3390/biom11111627

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

Article Metrics

Back to TopTop