Next Article in Journal
Semi-Supervised Learning Using Hierarchical Mixture Models: Gene Essentiality Case Study
Next Article in Special Issue
Optimization of Power Generation Grids: A Case of Study in Eastern Mexico
Previous Article in Journal
Operational Risk Reverse Stress Testing: Optimal Solutions
Previous Article in Special Issue
Modeling and Optimizing the Multi-Objective Portfolio Optimization Problem with Trapezoidal Fuzzy Parameters
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing

by
Juan P. Sánchez-Hernández
1,†,
Juan Frausto-Solís
2,*,†,
Juan J. González-Barbosa
2,
Diego A. Soto-Monterrubio
2,
Fanny G. Maldonado-Nava
2 and
Guadalupe Castilla-Valdez
2
1
Dirección de Informática, Electrónica y Telecomunicaciones, Universidad Politécnica del Estado de Morelos, Boulevard Cuauhnáhuac 566, Jiutepec 62574, Mexico
2
Graduate Program Division, Tecnológico Nacional de México/Instituto Tecnológico de Ciudad Madero, Cd. Madero 89440, Mexico
*
Author to whom correspondence should be addressed.
These authors contributed equally to the development of this paper.
Math. Comput. Appl. 2021, 26(2), 39; https://doi.org/10.3390/mca26020039
Submission received: 23 February 2021 / Revised: 23 April 2021 / Accepted: 27 April 2021 / Published: 29 April 2021
(This article belongs to the Special Issue Numerical and Evolutionary Optimization 2020)

Abstract

:
The Protein Folding Problem (PFP) is a big challenge that has remained unsolved for more than fifty years. This problem consists of obtaining the tertiary structure or Native Structure (NS) of a protein knowing its amino acid sequence. The computational methodologies applied to this problem are classified into two groups, known as Template-Based Modeling (TBM) and ab initio models. In the latter methodology, only information from the primary structure of the target protein is used. In the literature, Hybrid Simulated Annealing (HSA) algorithms are among the best ab initio algorithms for PFP; Golden Ratio Simulated Annealing (GRSA) is a PFP family of these algorithms designed for peptides. Moreover, for the algorithms designed with TBM, they use information from a target protein’s primary structure and information from similar or analog proteins. This paper presents GRSA-SSP methodology that implements a secondary structure prediction to build an initial model and refine it with HSA algorithms. Additionally, we compare the performance of the GRSAX-SSP algorithms versus its corresponding GRSAX. Finally, our best algorithm GRSAX-SSP is compared with PEP-FOLD3, I-TASSER, QUARK, and Rosetta, showing that it competes in small peptides except when predicting the largest peptides.

1. Introduction

Proteins or polypeptides are macromolecules built from amino acids (aa) and are mainly responsible for living beings’ functionality. Proteins are essentials elements because every protein has a specific function related to its unique three-dimensional structure named Native Structure (NS). All the proteins consist of a polymer chain of aa; the junctions with a small number of them are named peptides. The peptides have significant importance in the science community because of their multiple applications, for instance, in pharmaceutical research [1,2,3,4], drug design [5,6,7], diagnosis [8,9,10], and therapy [11,12]. To obtain the NS of proteins from an amino acid sequence could bring benefits to human beings.
The PFP has been identified as an important problem since Kendrew and Perutz’s research teams obtained the myoglobin and hemoglobin molecules’ tertiary structure, respectively [13,14]. These studies established the relation between function and structure. PFP consists of obtaining the three-dimensional structure of a protein with the lowest Gibbs free energy, thermodynamically stable three-dimensional conformation [15].
The PFP is considered an NP-hard problem [16]. Thus, presumably, none of the known exact algorithms can solve it in polynomial time. In other words, the execution time grows exponentially when using them. In contrast, any protein passes from the aa sequence to its NS three-dimensional structure very rapidly in nature. The latter issue is known as the Levinthal Paradox [17].
Several algorithms have been applied to solve the PFP successfully, and one of the most effective algorithms has been the Simulated Annealing algorithm (SA). The SA is commonly hybridized with other methods; the combination algorithms are called Hybrid Simulated Annealing algorithms (HSA). These algorithms successfully applied to peptides are the following:
(a)
The classical Monte Carlo Method, or SA, was applied to the PFP [18,19]. Additionally, an analytical tuning method to SA was proposed [20].
(b)
Golden Ratio Simulated Annealing (GRSA) family: Original GRSA proposing a cooling strategy [21], Evolutionary Golden Ratio SA (EGRSA) using genetic operators [22], and GRSA2, which is hybridization with the GRSA and Chemical Reaction Optimization algorithm (CRO) [23].
(c)
Metropolis and multiobjective optimization methods were applied in the previous CASP competitions. The approaches that traditionally have obtained the best results were Rossetta [24], QUARK [25], and I-TASSER [26]. However, deep learning applied by the Alphafold algorithm [27] achieved the best score in the CASP13 and CASP14.
(d)
PEP-FOLD3 algorithm, which uses secondary structure information and a Monte Carlo method, and is very successful for small peptides (5 to 50 aa) [28].
The HSA algorithms previously mentioned obtained excellent results for small proteins or peptides. However, when the number of aa increases, the variables (torsional angle of aa) are also increased, the computational time for exploring the solution space is considerable. As a result, the PFP area needs new approaches to obtaining better solutions for large peptides or proteins.
This paper proposes the methodology GRSA-SSP that combines GRSA algorithms with the Secondary Structure Prediction (SSP). For a given chain of aa representing a peptide or a protein, the GRSA-SSP performs two processes:
(a)
To obtain the first protein prediction from the secondary structure of the amino-acids sequence.
(b)
To refine the previous protein prediction by using GRSA family algorithms.
These two processes are performed in several steps described in this paper. The algorithms used in the second phase of GRSA-SSP can be one of the GRSA family algorithms. This paper named these hybrid algorithms GRSAX-SSP, where X is used to distinguish the GRSA algorithm. We evaluate our methodology using RMSD and TM-score metrics [29]. Additionally, experimentation is performed with a set of forty-five instances of peptides and a set of six mini proteins, which are compared with the most popular algorithms in the literature, such as PEP-FOLD3 [28], I-TASSER [30,31], Rosetta [24,32], and QUARK [25,33].
The paper’s organization is as follows: first, we present the introduction to PFP and HSA algorithms. Then, in the Background section, we review the Protein Folding Problem definition and some relevant research in the literature, and we explain the GRSA family of algorithms. In the next section, we describe the GRSA-SSP methodology. In the Results section, we present experimentation comparing the GRSA algorithms with those of the literature; also, we analyze the presented methodology’s performance. Finally, the conclusions of this research are presented.

2. Background

The PFP is a significant multidisciplinary problem that has been investigated for over half a century [34]. Different scientific areas have been studied, for example, computer science, bioinformatics, and molecular biology, concerning this problem, and three questions in particular need to be answered [34].
  • Which is the physical code in which an amino-acids sequence dictates an NS?
  • Why in nature do proteins fold very quickly while in silicon they fold relatively slower?
  • Is there an algorithm that predicts the protein structure from the amino-acids sequence?
This paper is related to the last question. We propose different strategies to obtain the NS tertiary structure using GRSA family algorithms and secondary structure prediction. As we mentioned before, finding new algorithms for PFP is significant not only because of its potential applications but also because it is an NP-hard problem [16], and the number of combinations that determine which algorithms must be explored in a very large solution space.

2.1. Definition of Ab-Initio and Force Fields

The ab initio modeling can be defined as an optimization problem where the Gibbs free energy is the objective function f(n), and this has to be minimized. Thus, this problem is defined as follows: let there be a sequence of amino acids: n = a1, a2, …, an; every amino acid has associated with it a set of angles σ1, σ2, …, σm where m represents a particular dihedral angle; then, minimizing the energy function f(σ|1, σ2, …, σm) provides the best tertiary structure or NS. The energy functions (force fields) are used for determining the energy of a protein structure [35], and some examples of these are AMBER [36], CHARMM [37], ECEPP/2, and ECEPP/3 [38]. The potential energy of ECEPP/2 is given by Equation (1), which is calculated in vacuo for only intramolecular energies, and this is the energy function to be minimized [39].
E t o t a l = j > i ( A i j r i j 12 B i j r i j 6 ) + 332 j > i q i q j ε r i j + j > i ( C i j r i j 12 D i j r i j 10 ) + n U n ( 1 ± c o s ( k n φ n ) )
where: rij is the distance in Å (angstroms) between the atoms i and j; Aij, Bij, Cij, and Dij are the parameters of the empirical potentials; qi and qj are the partial charges in the atoms i and j, respectively; ε is the dielectric constant; Un is the energetic torsion barrier of rotation about the bond n; kn is the multiplicity of the torsion angle φn.
In this paper, we use the potential energy of ECEPP/2 as an objective function because we explore the conformational space, and when the energy of the protein structure is minimized, then the protein structure is accepted.

2.2. Computational Approaches for PFP

The CASP organization has classified PFP models into two main groups:
Group 1: Template-based modeling (TBM). In this group, we find algorithms that use biological information obtained from the secondary structure of the target protein, homology, and fragments of other proteins. These algorithms have achieved good results for predicting protein structures in the CASP [32,40,41]. TBM involves several strategies; some of the most common are homology [42,43], threading [44], and fragment assembly [30,45].
Group 2: Ab initio. This prediction approach classically refers to the determination of the NS using only the aa sequence information. Unfortunately, ab initio algorithms have achieved good PFP results but only for small proteins with less than 120 residues [46]. The Ab initio modeling is the most challenging approach because it uses the amino acids’ sequence as unique information. Finding an optimal solution with ab initio is very difficult for big proteins because the solution space is enormous.
These two groups can be applied to small proteins or peptides (between 5 to 50 aa) [28,47]. There are successful studies applied to protein prediction using SA [48,49,50] or Monte Carlo algorithms with Metropolis-Hasting [26,27]. The Monte Carlo algorithms are also applied to the inverse protein folding problem, which objective function is to find a sequence given a structure [51,52]. This paper focuses on the classical PFP that consists of finding the functional structure given a sequence aa.
The Rosetta is a protein structure prediction or de novo approach that performs models for the tertiary structure using the primary and secondary structure predictions. The algorithm generates a local sequence to produce local structures (fragments) that form a target protein template. Additionally, the fragments are then assembled by randomly using a Monte Carlo simulated annealing algorithm. Finally, the fitness of individual conformation interactions is evaluated based on a scoring function derived from known protein structures. However, only peptides longer than 27 aa can be provided as input [32].
Another PFP approach is I-TASSER (Iterative Threading ASSEmbly Refinement). It has four principal parts: generating a template using a multi-threading method, fragments’ assembly method, refinement process, final model selection, and annotation tools. The I-TASSER applies an alignment of the target sequence and divides it into aligned using LOMETS [53,54] and nonaligned regions using the Monte Carlo algorithm. In the last step, annotation of functions is performed based on the structural models obtained using the BioLIP [55] database of ligand-protein interactions. Finally, the I-TASSER predicts protein structures from 10 to 1500 amino acids [31].
PEP-FOLD3 has a framework to predict the tertiary structure of peptides using de novo structure modeling. The process of predicting structure consists of three stages. Firstly, for a peptide amino acid sequence, a support vector machine is applied to predict the structural alphabet of fragments. Secondly, several models are generated using series of states and refined by a Monte Carlo algorithm. Finally, the five best conformations are selected [28].
Another approach is QUARK [33], in which an ab initio strategy is used to predict protein structures in ranges of 20 to 200 aa. Additionally, an assembly process of fragments with small structures is carefully selected and applied in the target sequence using a Monte Carlo algorithm.
SAINT2 is a fragment-based de novo structure prediction approach that has been successfully compared with the CASP12 approaches [56], which consists of a sequence-to-structure pipeline divided into four principal sections: (a) the secondary structure prediction where PSI-PRED [57] is applied, (b) the torsion angles prediction using SPINE-X [58], (c) a fragment library with the Flib package, and (d) the residue-residue contact prediction applying metaPSICOV [59]. Finally, the highest-scoring model is selected. In our methodology, sections (a) and (b) are applied, and they are shown in Figure 1.

The GRSA Family Algorithms

The SA algorithm is inspired by the physical annealing process of metals [60,61]. The algorithm has been applied with success in many NP-hard problems [20], including the PFP. SA employs the Metropolis algorithm to efficiently explore the solution space and obtain a good solution to optimization problems. We show the pseudocode of SA in Algorithm 1. Ti and Tf parameters define the initial and final temperatures, respectively; the α parameter represents the cooling factor. In the Metropolis cycle, new solutions are generated by a perturbation function. Finally, to accept or reject a new solution, an acceptance criterion based on Boltzmann distribution is applied (lines 11–14). The SA algorithm is executed until the final temperature, Tf, is reached. The SA algorithm source code is available at https://github.com/DrJuanFraustoSolis/SimulatedAnnealing.git (accessed on 28 April 2021).
Algorithm 1. SA algorithm Procedure.
Mca 26 00039 i001
However, when the solution space is very large, the algorithm’s exploration takes a long time to obtain optimal solutions. Thus, new algorithms are necessary. The GRSA algorithm was proposed, which has been successfully applied in different NP problems [62,63], including the PFP [18]. The main characteristics of GRSA are the cooling scheme that decreases according to Tfp temperature cuts calculated by the golden number (ɸ) and then a stop criterion that reduces the cost of exploration (Algorithm 2). GRSA has a similar structure to the SA algorithm (lines 4 to 16). The difference with SA is that the GRSA calculates Tfp temperature cuts (five cuts are recommended), and in each cut, an α parameter in the range [0.7, 1] is associated (the common higher value is 0.95); the intermediate α values in this range are determined with an increment δ which represent the α increment since the lowest until the highest α value (in this case, δ = 0.05 ). These alpha values are associated with each temperature cut (line 17). The algorithm reduces the temperature cooling speed; thus, the execution time, corresponding to lines 18 to 23, decreases. Finally, to reduce wasting time in low temperatures, where the quality of the result is not improved, a stop criterion was implemented using the least-squares method (lines 24 to 29). This stop criterion detects the stochastic equilibrium for some i Metropolis cycles. We measure the slope (m is a global variable) of the linear regression of the energy of these cycles. In this regression, we used the coordinates ( E i ,   i ); where i is in the range [2, κmax]. In our case, we used κmax = 5. The equilibrium is found when m is close to zero, calculated by (2).
m = κ i = 2 κ i E i ( i = 2 κ i ) ( i = 2 κ E i ) κ i = 2 κ i 2 ( i = 2 κ i ) 2
The Equation (2) can be written as follows (3):
m = 12 i = 2 κ i E i 6 ( κ 1 ) ( i = 1 κ E i ) κ 3 κ
where: κ is the number of metropolis cycles for measuring the slope, i is the iteration of every metropolis cycle, and Ei the energy in each iteration.
The evaluation of m in Equation (2) does not imply a significative execution time; the summations on Equation (3) are only cumulative operations in Algorithm 3. This algorithm determines the equilibrium with this Equation (3). The GRSA algorithm source code is available at https://github.com/DrJuanFraustoSolis/GRSA.git (accessed on 28 April 2021).
Algorithm 2. GRSA algorithm Procedure.
Mca 26 00039 i002
Algorithm 3. Equilibrium Function.
Mca 26 00039 i003
The EGRSA (Algorithm 4) is an algorithm integrated by the hybridization of GRSA with evolutionary techniques. This algorithm has an evolutionary perturbation (EGRSApert) in the GRSA phase (line 7), where a genetic algorithm is used. The EGRSA algorithm starts with a set of individuals generated for determining the initial solution designed as Si. Then in the Metropolis Cycle, the Si is perturbated by EGRSApert to generate new solutions. Next, the best individual generated Sj solution is selected of the population (lines 9 and 10). EGRSA is similar to GRSA, and both applied a stop criterion (see Algorithm 2.1) by the least-squares method [64,65] (lines 24–29). Algorithm 5 presents EGRSApert function, where one individual is a set of dihedral angles [ɸ1, Ψ1, Χ1, ω1, ɸ2, Ψ2, Χ2, ω2, …, ɸn, Ψn, Χn, ωn] and a population is a set of individuals. Then crossover and mutation operators are applied to generate new solutions by the perturbation function. Finally, when the number of generations is reached, the best individual of the population is selected. The EGRSA algorithm source code is available at https://github.com/DrJuanFraustoSolis/EGRSA.git (accessed on 28 April 2021).
Algorithm 4. EGRSA algorithm Procedure.
Mca 26 00039 i004
Algorithm 5. EGRSApert Function.
Mca 26 00039 i005
The GRSA2 algorithm [23] is a hybridization of GRSA with the CRO algorithm [66]. GRSA2 (Algorithm 6) is an enhancement of GRSA. It has the same structure as the previous algorithms revised in this paper. Specifically, GRSA2 has two principal differences in the perturbation phase, applying decomposition and soft collision (line 8) and the acceptance criterion (lines 10 to 14). In Algorithm 7, we show the perturbation process implemented in the GRSA2pert function. In GRSA2, two soft collisions are used (unimolecular, Intermolecular). This algorithm has been applied only in the PFP with a set of 19 peptides and compared with I-TASSER and PEP-FOLD3 approaches obtaining outstanding results in the case of peptides [23]. The GRSA2 algorithm source code is available at https://github.com/DrJuanFraustoSolis/GRSA2.git (accessed on 28 April 2021).
Algorithm 6. GRSA2 algorithm Procedure.
Mca 26 00039 i006
Algorithm 7. GRSA2pert Function.
Mca 26 00039 i007

3. GRSA-SSP Methodology

In this section, we present the GRSA-SSP methodology (Figure 1). This methodology has two main processes:
(a)
The prediction of the torsion angles (initial solution) from the secondary structure; that corresponds to stages 1 to 4 in Figure 1.
(b)
The refinement of the solution obtained from the secondary structure. This is performed with GRSA algorithms showed in stage four (Figure 1).
The GRSA-SSP methodology has an input (amino acid sequence), an output (tertiary structure prediction), and four stages: (1) secondary structure prediction, (2) torsion angles prediction, (3) template construction, and (4) refinement by GRSAX algorithms. Next, we explain each of these stages:
Input (Amino acid sequence). The amino acid sequences are taken as input.
(1)
Secondary structure prediction. This secondary structure, which corresponds to the amino acid sequence and is predicted using PSI-PRED [57]. This algorithm generates a sequence profile with PSI-BLAST [67] and performs the prediction of the stage, such as the helix (H), strand (E), and coil (C). PSI-PRED calculates the probability of each possible state and defines the most likely structure.
(2)
Torsion angles prediction. The secondary structure’s prediction is essential for this stage, where SPINE-X is used to obtain the torsion angles (ɸ, Ψ, and ω) of each amino acid. This process is realized through the Position-Specific Score Matrix and Physical Parameters [58]. SPINE-X applies artificial neural networks to obtain the best predictions of the target’s proteins.
(3)
Model construction. In this stage, the torsion angles or variables are used to construct a template as initial solution Si = [ɸ1, Ψ1, Χ1, ω1, ɸ2, Ψ2, Χ2, ω2, …, ɸn, Ψn, Χn, ωn] that is represented by amino acids subscript 1 to n and the same form by the following amino acids up to n; n is dependent on the size of an amino acid sequence of the target protein. The torsion angles represent the base column of the peptide on which the refinement will be performed.
(4)
Refinement by GRSAX. When the previous stages construct the peptide template, we can apply a GRSAX algorithm such as GRSA (renamed GRSA1), EGRSA (renamed GRSAE), and GRSA2, as well as the classical SA (GRSA0). The GRSAX algorithms are tested individually for comparison, which obtains a better tertiary structure of the target peptide. Moreover, once the energy and three-dimensional structure is obtained, the structure is evaluated with the RMSD and TM-score [29] metrics.
Output. The GRSAX-SSP algorithm obtains the tertiary structure prediction.

4. Results

We performed the next GRSAX-SSP algorithms with the proposed methodology: (a) GRSA0-SSP using classical SA [19], (b) GRSA1-SSP using original GRSA [21], (c) GRSAE-SSP using EGRSA [22], and (d) GRSA2-SSP using GRSA2 [23]. For all of them, we used the methodology presented in Figure 1. The peptides in this experimentation have 9 to 49 amino acids. The number of variables (torsion angles) for each peptide in this data set is in the range [49, 304]. We chose this set because these instances (peptides) were used before in the literature. This set was also useful for comparing the GRSA2-SSP algorithm with the top-performing approaches of the CASP, which can be used for small peptides. We compared the last algorithm with I-TASSER, PEP-FOLD3, QUARK, and Rosetta, which are among the best algorithms in the CASP competition. We noted a difference between the GRSAX-SSP algorithms and the one that only applies ab initio by naming it GRSAX. Table 1 presents the set of 45 instances sorted by the number of variables taken from [23,28,68,69] and a PDB code represents each peptide.
In the experimentation, the GRSAX-SSP algorithms were executed 30 times to validate the results. The energy function ECEPP/2 is determined with SMMP framework [38]; it is the objective function of our optimization algorithms. An analytical tuning [20] was performed to obtain the initial and final temperature for each instance. In GRSA0-SSP the α value is 0.95, and the temperature range has zero golden sections. For GRSA1-SSP, GRSAE-SSP, and GRSA2-SSP algorithms, the same cooling scheme was used, using the α parameter with values from 0.75 to 0.95 with five golden ratio sections, which was determined by experimentation [21,22,23]. The GRSAX-SSP algorithms were executed in one of the terminals of the Ehecatl cluster in TecNM/IT Ciudad Madero, and it has the following characteristics: Intel® Xeon® processor at 2.30 GHz, Memory: 64 GB (4 × 16 GB) ddr4-2133, Linux CentOS operating system, and Fortran language.
We used the minimum energy quality values, the RMSD, and TM-score to evaluate the results, which are two metrics of the structural quality used for PFP algorithms. The RMSD is a structural measure between the native structure and the one predicted by the GRSAX-SSP and classical SA named here as GRSA0:
(a)
If the RMSD has a value close to zero, the quality of the structure is considered excellent. On the contrary, the quality is worse.
(b)
The TM-score is also used to measure the similarity between two structures. When the TM-score is greater than 0.5, it indicates that there is a good similarity between the two structures, and the tested one has the same fold. Otherwise, as the TM-score is lower than 0.5, the target peptide has a different fold [29].
The TM-score metrics can be calculated using the TM-align [70] (an algorithm to obtain the best structural alignment between two proteins) or in a classical formulation [29]. In this paper, we use the classical formulation of TM-score.
GRSAX-SSP algorithms use a model determined by the secondary structure, and then it is refined for obtaining a better prediction. The results are compared with the GRSAX based on ab initio that only uses the amino acid sequence as information. Figure 2, Figure 3, Figure 4 and Figure 5 show average results related to energy (kcal/mol), RMSD, and TM-score for each peptide. The numbers in the x-axis, represent the instances or peptides of Table 1, and each instance is a set of torsional angles X = [ɸ1, Ψ1, Χ1, ω1, ɸ2, Ψ2, Χ2, ω2, …, ɸn, Ψn, Χn, ωn] associated to each amino acid. We averaged the results of 30 executions for comparison.
Figure 2 shows that GRSA0-SSP has better behavior than GRSA0 or classical SA. Note that in all the peptides, GRSA0-SSP obtained the lowest energy. In other cases, the RMSD is more stable with small instances (1–16), and in the next instances, the behavior is equal. Additionally, when we compared with TM-score, the behavior, in general, is similar. In conclusion, by implementing this methodology in GRSA0-SSP with these instances, we obtained slightly improved results.
Figure 3 presents the comparison of the GRSA1-SSP versus GRSA1 with the same metrics; we observed the behavior with the 45 instances evaluated. In terms of energy, RMSD, and TM-score, the performance of GRSA1-SSP is equivalent to GRSA1.
Figure 4 shows the behavior of GRSAE-SSP, and we compared it with the original GRSAE algorithm. In this figure, we can appreciate that the results are equivalent in all cases when energy, RMSD, and TM-score are used for comparison.
In Figure 5, we present the comparison of GRSA2 versus GRSA2-SSP. Note that the results obtained in every instance are very remarkable, and the superiority of GRSA2-SSP uses the metrics of energy, RMSD, and TM-Score. In this case, we applied the methodology GRSA-SSP to improve the behavior of the classical GRSA2 algorithm.
Finally, in Figure 6, we present the comparison of the GRSAX-SSP family algorithms. We observe that GRSA2-SSP has the best values in several instances against the other algorithms, being higher than the others. Therefore, the best behavior of the algorithms with secondary structure prediction is GRSA2-SSP.
Furthermore, Figure 7 presents the computational time of the GRSAX-SSP family algorithms. The GRSA2-SSP has the best behavior in time with low values in most of the instances compared to the other algorithms.
Table 2 presents the results obtained by GRSA2-SSP. For each instance, we show the best TM-score and their RMSD. Additionally, we calculated the average of the RMSD and TM-score for the five best predictions. Complementing the results, we determined the standard deviation (std) of the RMSD and TM-score for the five best predictions and included the best type of secondary structure: A (mainly alpha), B (mainly beta), and N (mainly none). This classification as A, B, and N is based on the secondary structure predominating in each peptide [27,68,69,71,72]. We sort Table 2 by the number of amino acids for comparing the best results obtained by GRSA2-SSP with the best algorithms of the literature. This comparison is presented in Figures 9–11.
Figure 8 shows the GRSA2-SSP algorithm performance with instances classified by secondary structure. We show that the GRSA2-SSP algorithm has the best behavior in alpha structure instances evaluated with TM-score in Figure 8a and RMSD metrics in Figure 8b. The values in Figure 8 are the best obtained using TM-score and their RMSD. In Figure 8c,d, we present the TM-score average for the five best predictions and their RMSD average.
In Figure 9, Figure 10 and Figure 11, we present the behavior of the GRSA2-SSP algorithm, and we compare it with the results obtained from the approaches PEP-FOLD3, I-TASSER, QUARK, and Rosetta. We divided the dataset of Table 1 into three groups of 15 instances; groups 1, 2, and 3 have instances 1–15, 16–30, and 31–45. We compared these groups using the metrics RMSD, TM-score, GDT-TS [73], and TM-score (classical), and we present the best TM-score, the average of the five best predictions of the TM-score, and their RMSD. Additionally, we present the GDT-TS average and TM-score average.
In Figure 9, we introduced the comparison of the first group, and we observed that GRSA2-SSP behaves similarly to I-TASSER and PEP-FOLD3, but in this group of small peptides, PEP-FOLD3 is slightly better than our algorithm and I-TASSER when GDT-TS is compared (Figure 9e). Furthermore, we observed that our algorithm is competitive in this group. In this comparison, Rossetta and QUARK were not added because the minimal number of amino acids predicted are 27 and 20, respectively.
Figure 10 compares the second group of 16 to 30 amino acids with the best and the five best obtained using the TM-score metric and their RMSD, and the GDT-TS average. In this comparison, we added the second group of instances’ results of QUARK; Rosetta was omitted because it is unable to predict most of the instances of this group.
In Figure 10a we observe very similar behavior among GRSA2-SSP, PEP-FOLD3, I-TASSER, and Rosetta. Note in this figure, GRSA2-SSP and PEP-FOLD3 obtain the best prediction. In Figure 10c, when the best five predictions are compared, I-TASSER obtains the best results, followed by PEPFOLD3 and GRSA2-SSP. Additionally, when the RMSD average is compared (Figure 10d), I-TASSER is the best, followed by PEP-FOLD3 and GRSA2-SSP. Finally, in Figure 10e, when GDT-TS is compared, GRSA2-SSP has a similar performance to PEP-FOLD3, I-TASSER, and QUARK. According to this figure, GRSA2-SSP and I-TASSER obtained a similar average.
Figure 11 compares the third group of 31 to 49 amino acids with the five best results obtained using the TM-score metric and their RMSD y GDT-TS. This comparison added the Rosetta approach because it can process the number of aa in this group. As we observe, the best algorithm is I-TASSER, followed by Rosetta, QUARK, PEP-FOLD3, and finally GRSA2-SSP.
The 45 instances evaluated in the below experimentation show the application of the secondary structure results and refine them with the GRSAX algorithms, enhancing the performance in energy, RMSD, and TM-score. Specifically, when GRSA2-SSP is compared with PEP-FOLD3, I-TASSER, QUARK, and Rosetta, we observed that our algorithm performs well in small instances (Group 1 and 2). Nevertheless, in the largest instances, our algorithm is not the best, but it is competitive.
We carried out a second experimentation with six mini-proteins (5wll, 5lo2, 5up5, 5uoi, 2ki0, and 2kik) presented in Table 3. The mini-proteins come from the de novo protein design field [74,75,76,77,78]. This data set was proposed to observe the behavior of our best algorithm in these kinds of instances.
We applied the same evaluation of all the algorithms, as in the first experimentation, using RMSD, TM-score, and GDT-TS metrics. Table 4 shows the results of all the algorithms in this data set. Evaluating them with TM-score and GDT-TS, we observe that the best algorithms were Rosetta, I-TASSER, and GRSA2-SSP, where the number of times the best results were achieved 3, 2, and 1, respectively. Additionally, evaluating with the RMSD, the best algorithms were again Rosseta, I-TASSER, and GRSA2-SSP, but this time they obtained the best results in two instances, which were (5uoi, 2kik), (2ki0, 5up5), and (5wll, 5lo2), respectively. As a result, we can say that Rosetta is the best algorithm, followed by I-TASSER, and GRSA2-SSP.

5. Conclusions

In this paper, we present the methodology GRSA-SSP for Protein Folding Problem applied to peptides. The objective of this problem is to predict the functional tridimensional protein structure. The algorithms developed with this methodology are GRSA0-SSP, GRSA1-SSP, GRSAE-SSP, and GRSA2-SSP. The main relevance of the algorithm GRSA2-SSP, developed with this methodology, is that it produces very good results in the case of peptides; specifically, it is similar or better than the algorithms Rosetta, PEP-FOLD3, QUARK, and I-TASSER for the small and medium peptides, according to the experimentation presented. The last algorithms have traditionally been among the best of the CASP competition; besides, they use modern machine learning techniques like artificial neural networks.
We compared the algorithms developed with the original algorithms GRSA0, GRSA1, GRSAE, and GRSA2; we used a data set of 45 instances for this comparison. We showed that the hybrid algorithms produced with the GRSA-SSP methodology outperform the original ones. For this comparison, we used the metrics Energy, RMSD, TM-score, and execution time. We observed that the best of all these algorithms is GRSA2-SSP formulated with the proposed methodology.
We made a second evaluation comparing the GRSA2-SSP algorithm with the best state-of-the-art algorithms (we used the same data set of 45 instances). We selected for this comparison PEP-FOLD3, I-TASSER, QUARK, and Rosetta. We used a data set of forty-five instances divided into three groups, from small to large peptides. The experimentation shows that for groups 1 and 2, GRSA2-SSP performs as well as these algorithms. We observe that for the first group PEP-FOLD3 was the best, followed by GRSA2-SSP, while in the second group, the best algorithm was I-TASSER followed by GRSA2-SSP and PEP-FOLD3. Finally, in the third group, the best algorithm was Rosseta, followed by I-TASSER. Additionally, we present an analysis of GRSA2-SSP results for each type of secondary structure, obtaining a better behavior with alpha structures.
Furthermore, we assessed GRSA2-SSP with a second data set of six instances named mini proteins. The GRSA2-SSP results were compared with PEP-FOLD3, I-TASSER, QUARK, and Rosetta. The best algorithms in this data set were Rosetta, I-TASSER, and GRSA2-SSP because the number of times the best TM-score and GDT-TS were 3, 2, and 1, respectively. However, each of the three achieved two times the first place when RMSD was evaluated. As a result, the best of these algorithms for this data set is Rosetta, followed by I-TASSER and GRSA2-SSP.
We conclude that GRSAX-SSP algorithms enhance the original GRSA algorithms. The best of them is GRSA2-SSP which achieves very good results, surpassing the best state-of-art for peptides up to thirty amino acids. Finally, we note that the main advantage of our methodology is that it is simpler than the most powerful approaches of the literature.

Author Contributions

J.F.-S. and J.P.S.-H. contributed equally to the development of this paper. Conceptualization, J.P.S.-H., D.A.S.-M. and J.F.-S.; methodology J.F.-S., D.A.S.-M., J.P.S.-H., and J.J.G.-B.; Software J.P.S.-H., D.A.S.-M. and F.G.M.-N.; validation, J.P.S.-H. and J.F.-S.; formal analysis, D.A.S.-M., F.G.M.-N., J.J.G.-B., and G.C.-V.; writing—original draft J.F.-S., J.P.S.-H., and D.A.S.-M.; writing—review and editing, J.F.-S., D.A.S.-M. and J.P.S.-H. All authors have read and agreed to the published version of the manuscript.

Funding

This research received no external funding.

Acknowledgments

The authors would like to acknowledge with appreciation and gratitude CONACYT and TecNM/Instituto Tecnológico de Ciudad Madero. Also, we acknowledge Laboratorio Nacional de Tecnologías de la Información (LaNTI) for the access to the cluster.

Conflicts of Interest

The authors declare that they have no competing interests.

References

  1. Uhlig, T.; Kyprianou, T.; Martinelli, F.G.; Oppici, C.A.; Heiligers, D.; Hills, D.; Calvo, X.R.; Verhaert, P. The emergence of peptides in the pharmaceutical business: From exploration to exploitation. EuPA Open Proteom. 2014, 4, 58–69. [Google Scholar] [CrossRef] [Green Version]
  2. Agyei, D.; Danquah, M.K. Industrial-scale manufacturing of pharmaceutical-grade bioactive peptides. Biotechnol. Adv. 2011, 29, 272–277. [Google Scholar] [CrossRef]
  3. Patel, L.N.; Zaro, J.L.; Shen, W.-C. Cell Penetrating Peptides: Intracellular Pathways and Pharmaceutical Perspectives. Pharm. Res. 2007, 24, 1977–1992. [Google Scholar] [CrossRef]
  4. Danquah, M.; Agyei, D. Pharmaceutical applications of bioactive peptides. OA Biotechnol. 2012, 1. [Google Scholar] [CrossRef] [Green Version]
  5. Fosgerau, K.; Hoffmann, T. Peptide therapeutics: Current status and future directions. Drug Discov. Today 2015, 20, 122–128. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  6. Vetter, I.; Davis, J.L.; Rash, L.D.; Anangi, R.; Mobli, M.; Alewood, P.F.; Lewis, R.J.; King, G.F. Venomics: A new paradigm for natural products-based drug discovery. Amino Acids 2010, 40, 15–28. [Google Scholar] [CrossRef] [PubMed]
  7. Craik, D.J.; Fairlie, D.P.; Liras, S.; Price, D. The Future of Peptide-based Drugs. Chem. Biol. Drug Des. 2013, 81, 136–147. [Google Scholar] [CrossRef] [PubMed]
  8. Stalmach, A.; Johnsson, H.; McInnes, I.B.; Husi, H.; Klein, J.; Dakna, M. Identification of urinary peptide biomarkers associated with rheumatoid arthritis. PLoS ONE 2014, 9, e104625. [Google Scholar]
  9. Gautam, A.; Kapoor, P.; Chaudhary, K.; Kumar, R.; Raghava, G.P. Tumor homingpeptides as molecular probes for cancer therapeutics, diagnostics and theranostics. Curr. Med. Chem. 2014, 21, 2367–2391. [Google Scholar] [CrossRef] [PubMed]
  10. Li, Z.J.; Cho, C.H. Peptides as targeting probes against tumor vasculature for diagnosis and drug delivery. J. Transl. Med. 2012, 10 (Suppl. S1). [Google Scholar] [CrossRef] [Green Version]
  11. Lau, J.L.; Dunn, M.K. Therapeutic peptides: Historical perspectives, current development trends, and future directions. Bioorgan. Med. Chem. 2018, 26, 2700–2707. [Google Scholar] [CrossRef] [PubMed]
  12. Vlieghe, P.; Lisowski, V.; Martinez, J.; Khrestchatisky, M. Synthetic therapeutic peptides: Science and market. Drug Discov. Today 2010, 15, 40–56. [Google Scholar] [CrossRef] [PubMed]
  13. Kendrew, J.C.; Bodo, G.; Dintzis, H.M.; Parrish, R.G.; Wyckoff, H.; Phillips, D.C. A three-dimensional model of the myoglobin molecule obtained by X-ray analysis. Nature 1958, 181, 662–666. [Google Scholar] [CrossRef] [PubMed]
  14. Perutz, M.F.; Rossmann, M.G.; Cullis, A.F.; Muirhead, H.I.L.A.R.Y.; Will, G.; North, A.C.T. Structure of hemoglobin. Brookhaven. Symp. Biol. 1960, 13, 165–183. [Google Scholar] [PubMed]
  15. Anfinsen, C.B. Principles that Govern the Folding of Protein Chains. Science 1973, 181, 223–230. [Google Scholar] [CrossRef] [Green Version]
  16. Hart, W.E.; Istrail, S. Robust Proofs of NP-Hardness for Protein Folding: General Lattices and Energy Potentials. J. Comput. Biol. 1997, 4, 1–22. [Google Scholar] [CrossRef] [PubMed]
  17. Levinthal, C. Are There Pathways for Protein Folding. J. Chim. Phys. 1968, 65, 44–45. [Google Scholar] [CrossRef]
  18. Li, Z.; Scheraga, H.A. Monte Carlo-minimization approach to the multiple-minima problem in protein folding. Proc. Natl. Acad. Sci. USA 1987, 84, 6611–6615. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  19. Morales, L.B.; Garduño-Juárez, R.; Romero, D. Applications of Simulated Annealing to the Multiple-Minima Problem in Small Peptides. J. Biomol. Struct. Dyn. 1991, 8, 721–735. [Google Scholar] [CrossRef]
  20. Frausto, J.; Román, E.F.; Romero, D.; Soberon, X.; Liñán, E. Analytically Tuned Simulated Annealing Applied to the Protein Folding Problem. In Proceedings of the 7th International Conference on Computational Science, Beijing, China, 27–30 May 2007; Springer: Berlin/Heidelberg, Germany, 2007; Volume 4488, pp. 370–377. [Google Scholar]
  21. Frausto, J.; Sánchez, J.P.; Sánchez, M.; García, E.L. Golden Ratio Simulated Annealing for Protein Folding Problem. Int. J. Comput. Methods 2015, 12, 1550037. [Google Scholar] [CrossRef]
  22. Maldonado, F.; Frausto, J.; Sánchez, J.; González, J.; Liñán, E.; Castilla, G. Evolutionary GRSA for Protein Structure Prediction. Int. J. Comb. Optim. Probl. Inform. 2016, 7, 75–86. [Google Scholar]
  23. Frausto, J.; Sánchez, J.P.; Maldonado, F.; González, J.J. GRSA Enhanced for Protein Folding Problem in the Case of Peptides. Axioms 2019, 8, 136. [Google Scholar] [CrossRef] [Green Version]
  24. Hiranuma, N.; Park, H.; Baek, M.; Anishchenko, I.; Dauparas, J.; Baker, D. Improved protein structure refinement guided by deep learning based accuracy estimation. Nat. Commun. 2021, 12, 1–11. [Google Scholar] [CrossRef]
  25. Xu, D.; Zhang, Y. Toward optimal fragment generations for ab initio protein structure assembly. Proteins 2012, 81, 229–239. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Wang, D.; Geng, L.; Zhao, Y.-J.; Yang, Y.; Huang, Y.; Zhang, Y.; Shen, H.-B. Artificial intelligence-based multi-objective optimization protocol for protein structure refinement. Bioinformatics 2020, 36, 437–448. [Google Scholar] [CrossRef] [PubMed]
  27. Senior, A.W.; Evans, R.; Jumper, J.; Kirkpatrick, J.; Sifre, L.; Green, T.; Qin, C.; Žídek, A.; Nelson, A.W.R.; Bridgland, A.; et al. Protein structure prediction using multiple deep neural networks in the 13th Critical Assessment of Protein Structure Prediction (CASP13). Proteins Struct. Funct. Bioinform. 2019, 87, 1141–1148. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  28. Lamiable, A.; Thévenet, P.; Rey, J.; Vavrusa, M.; Derreumaux, P.; Tufféry, P. PEP-FOLD3: Faster de Novo Structure Prediction for Linear Peptides in Solution and in Complex. Nucleic Acids Res. 2016, 44, W449–W454. [Google Scholar] [CrossRef] [Green Version]
  29. Zhang, Y.; Skolnick, J. Scoring function for automated assessment of protein structure template quality. Proteins 2004, 57, 702–710. [Google Scholar] [CrossRef]
  30. Yang, J.; Yan, R.; Roy, A.; Xu, D.; Poisson, J.; Zhang, Y. The I-TASSER Suite: Protein structure and function prediction. Nat. Methods 2015, 12, 7–8. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  31. Yang, J.; Zhang, Y. I-TASSER server: New development for protein structure and function predictions. Nucleic Acids Res. 2015, 43, W174–W181. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  32. Rohl, C.A.; Strauss, C.E.; Misura, K.M.; Baker, D. Protein Structure Prediction Using Rosetta. Oncogene Tech. 2004, 383, 66–93. [Google Scholar] [CrossRef]
  33. Xu, D.; Zhang, Y. Ab initio protein structure assembly using continuous structure fragments and optimized knowledge-based force field. Proteins 2012, 80, 1715–1735. [Google Scholar] [CrossRef] [Green Version]
  34. Dill, K.A.; Maccallum, J.L. The Protein-Folding Problem, 50 Years On. Science 2012, 338, 1042–1046. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  35. Dill, K.A. Dominant forces in protein folding. Biochemistry 1990, 29, 7133–7155. [Google Scholar] [CrossRef] [PubMed]
  36. Ponder, J.W.; Case, D.A. Force Fields for Protein Simulations. Accessory Fold. Proteins 2003, 66, 27–85. [Google Scholar] [CrossRef]
  37. Brooks, B.R.; Bruccoleri, R.E.; Olafson, B.D.; States, D.J.; Swaminathan, S.; Karplus, M. CHARMM: A program for macromolecular energy, minimization, and dynamics calculations. J. Comput. Chem. 1983, 4, 187–217. [Google Scholar] [CrossRef]
  38. Momany, F.A.; McGuire, R.F.; Burgess, A.W.; Scheraga, H.A. Energy parameters in polypeptides. VII. Geometric parameters, partial atomic charges, nonbonded interactions, hydrogen bond interactions, and intrinsic torsional potentials for the naturally occurring amino acids. J. Phys. Chem. 1975, 79, 2361–2381. [Google Scholar] [CrossRef]
  39. Eisenmenger, F.; Hansmann, U.H.; Hayryan, S.; Hu, C.-K. [SMMP] A modern package for simulation of proteins. Comput. Phys. Commun. 2001, 138, 192–212. [Google Scholar] [CrossRef]
  40. Jiang, P.; Xu, J. RaptorX: Exploiting structure information for protein alignment by statistical inference. Proteins 2011, 79, 161–171. [Google Scholar]
  41. Zhou, H.; Pandit, S.B.; Skolnick, J. Performance of the Pro-sp3-TASSER server in CASP8. Proteins 2009, 77, 123–127. [Google Scholar] [CrossRef] [Green Version]
  42. Konstantin, A.; Lorenza, B.; Jürgen, K.; Torsten, S. The SWISS-MODEL workspace: A web-based environment for protein structure homology modelling. Bioinformatics 2006, 22, 195–201. [Google Scholar]
  43. Schmitt, S.; Kuhn, D.; Klebe, G. A New Method to Detect Related Function among Proteins Independent of Sequence and Fold Homology. J. Mol. Biol. 2002, 323, 387–406. [Google Scholar] [CrossRef]
  44. Lemer, C.M.-R.; Rooman, M.J.; Wodak, S.J. Protein structure prediction by threading methods: Evaluation of current techniques. Proteins 1995, 23, 337–355. [Google Scholar] [CrossRef]
  45. Dorn, M.; e Silva, M.B.; Buriol, L.S.; Lamb, L.C. Three-dimensional protein structure prediction: Methods and computational strategies. Comput. Biol. Chem. 2014, 53, 251–276. [Google Scholar] [CrossRef]
  46. Zhang, Y. Interplay of I-TASSER and QUARK for template-based and ab initio protein structure prediction in CASP10. Proteins 2013, 82, 175–187. [Google Scholar] [CrossRef] [Green Version]
  47. Bhardwaj, G.; Mulligan, V.K.; Bahl, C.D.; Gilmore, J.M.; Harvey, P.J.; Cheneval, O.; Buchko, G.W.; Pulavarti, S.V.S.R.K.; Kaas, Q.; Eletsky, A.; et al. Accurate de novo design of hyperstable constrained peptides. Nature 2016, 538, 329–335. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  48. Harada, R.; Nakamura, T.; Shigeta, Y. A Fast Convergent Simulated Annealing Algorithm for Protein-Folding: Simulated Annealing Outlier FLOODing (SA-OFLOOD) Method. Bull. Chem. Soc. Jpn. 2016, 89, 1361–1367. [Google Scholar] [CrossRef]
  49. Zhang, L.; Ma, H.; Qian, W.; Li, H. Protein structure optimization using improved simulated annealing algorithm on a three-dimensional AB off-lattice model. Comput. Biol. Chem. 2020, 85, 107237. [Google Scholar] [CrossRef] [PubMed]
  50. Zhang, L.; Ma, H.; Qian, W.; Li, H. Sequence-based protein structure optimization using enhanced simulated annealing al-gorithm on a coarse-grained model. J. Mol. Model. 2020, 26, 1–13. [Google Scholar] [CrossRef]
  51. Mitra, P.; Shultis, D.; Brender, J.R.; Czajka, J.; Marsh, D.; Gray, F.; Cierpicki, T.; Zhang, Y. An Evolution-Based Approach to De Novo Protein Design and Case Study on Mycobacterium tuberculosis. PLoS Comput. Biol. 2013, 9, e1003298. [Google Scholar] [CrossRef] [Green Version]
  52. Banerjee, A.; Pal, K.; Mitra, P. An evolutionary profile guided greedy parallel replica-exchange Monte Carlo search algorithm for rapid convergence in protein design. IEEE/ACM Trans. Comput. Biol. Bioinform. 2019, 18, 489–499. [Google Scholar] [CrossRef] [PubMed]
  53. Wu, S.; Zhang, Y. LOMETS: A local meta-threading-server for protein structure prediction. Nucleic Acids Res. 2007, 35, 3375–3382. [Google Scholar] [CrossRef] [Green Version]
  54. Zheng, W.; Zhang, C.; Wuyun, Q.; Pearce, R.; Li, Y.; Zhang, Y. LOMETS2: Improved meta-threading server for fold-recognition and structure-based function annotation for distant-homology proteins. Nucleic Acids Res. 2019, 47, W429–W436. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  55. Yang, J.; Roy, A.; Zhang, Y. BioLiP: A semi-manually curated database for biologically relevant ligand–protein interactions. Nucleic Acids Res. 2012, 41, D1096–D1103. [Google Scholar] [CrossRef] [Green Version]
  56. De Oliveira, S.; Law, E.C.; Shi, J.; Deane, C.M. Sequential search leads to faster, more efficient fragment-based de novo protein structure prediction. Bioinformatics 2017, 34, 1132–1140. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  57. Jones, D.T. Protein secondary structure prediction based on position-specific scoring matrices. J. Mol. Biol. 1999, 292, 195–202. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  58. Faraggi, E.; Kloczkowski, A. Accurate Prediction of One-Dimensional Protein Structure Features Using SPINE-X. In Prediction of Protein Secondary Structure; Methods in Molecular Biology; Humana Press: New York, NY, USA, 2017; Volume 1484, pp. 45–53. [Google Scholar]
  59. Jones, D.T.; Singh, T.; Kosciolek, T.; Tetchner, S. MetaPSICOV: Combining coevolution methods for accurate prediction of contacts and long range hydrogen bonding in proteins. Bioinformatics 2015, 31, 999–1006. [Google Scholar] [CrossRef] [PubMed]
  60. Kirkpatrick, S.; Gelatt, C.D.; Vecchi, M.P. Optimization by Simulated Annealing. Science 1983, 220, 671–680. [Google Scholar] [CrossRef]
  61. Černý, V. Thermodynamical approach to the traveling salesman problem: An efficient simulation algorithm. J. Optim. Theory Appl. 1985, 45, 41–51. [Google Scholar] [CrossRef]
  62. Frausto, J.; Martinez, F. Golden Ratio Annealing for Satisfiability Problems Using Dynamically Cooling Schemes. In Foundations of Intelligent Systems; Springer: Berlin/Heidelberg, Germany, 2008; Volume 4994, pp. 215–224. [Google Scholar]
  63. Frausto, J.; Martinez, F. Golden annealing method for job shop scheduling problem. In MACMESE’08: Proceedings of the 10th WSEAS International Conference on Mathematical and Computational Methods in Science and Engineering; World Scientific and Engineering Academy and Society (WSEAS): Stevens Point, WI, USA, 2008; pp. 374–379. [Google Scholar]
  64. Frausto, J.; Liñán, E.; Sánchez, J.P.; González, J.J.; González, C.; Castilla, G. Multiphase Simulated Annealing Based on Boltzmann and Bose–Einstein Distribution Applied to Protein Folding Problem. Adv. Bioinform. 2016, 2016, 7357123. [Google Scholar]
  65. Martinez, F.; Frausto, J. A simulated annealing algorithm for the satisfiability problem using dynamic Markov chains with linear regression equilibrium. Simulated Annealing. InTechOpen 2012, 21, 281–285. [Google Scholar]
  66. Lam, A.Y.S.; Li, V.O.K. Chemical Reaction Optimization: A tutorial. Memetic Comput. 2012, 4, 3–17. [Google Scholar] [CrossRef] [Green Version]
  67. Altschul, S.F.; Madden, T.L.; Schäffer, A.A.; Zhang, J.; 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]
  68. Maupetit, J.; Derreumaux, P.; Tuffery, P. PEP-FOLD: An online resource for de novo peptide structure prediction. Nucleic Acids Res. 2009, 37 (Suppl. S2), W498–W503. [Google Scholar] [CrossRef] [Green Version]
  69. Shen, Y.; Maupetit, J.; Derreumaux, P.; Tufféry, P. Improved PEP-FOLD approach for peptide and miniprotein structure pre-diction. J. Chem. Theory Comput. 2014, 10, 4745–4758. [Google Scholar] [CrossRef] [PubMed]
  70. 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] [PubMed]
  71. Munte, C.E.; Vilela, L.; Kalbitzer, H.R.; Garratt, R.C. Solution structure of human proinsulin C-peptide. FEBS J. 2005, 272, 4284–4293. [Google Scholar] [CrossRef]
  72. Luitz, M.P.; Bomblies, R.; Zacharias, M. Comparative Molecular Dynamics Analysis of RNase-S Complex Formation. Biophys. J. 2017, 113, 1466–1474. [Google Scholar] [CrossRef] [Green Version]
  73. Zemla, A.; Moult, J.; Fidelis, K. Processing and evaluation of predictions in CASP4. Proteins 2001, 45, 13–21. [Google Scholar] [CrossRef]
  74. Lombardi, A.; Pirro, F.; Maglio, O.; Chino, M.; DeGrado, W.F. De Novo Design of Four-Helix Bundle Metalloproteins: One Scaffold, Diverse Reactivities. Accounts Chem. Res. 2019, 52, 1148–1159. [Google Scholar] [CrossRef]
  75. Liang, H.; Chen, H.; Fan, K.; Wei, P.; Guo, X.; Jin, C.; Zeng, C.; Tang, C.; Lai, L. De novo design of a beta alpha beta motif. Angew. Chem. Int. Ed. Engl. 2009, 48, 3301–3303. [Google Scholar] [CrossRef] [PubMed]
  76. Baker, E.G.; Bartlett, G.J.; Goff, K.L.P.; Woolfson, D.N. Miniprotein Design: Past, Present, and Prospects. Accounts Chem. Res. 2017, 50, 2085–2092. [Google Scholar] [CrossRef] [PubMed]
  77. Rocklin, G.J.; Chidyausiku, T.M.; Goreshnik, I.; Ford, A.; Houliston, S.; Lemak, A.; Carter, L.; Ravichandran, R.; Mulligan, V.K.; Chevalier, A.; et al. Global analysis of protein folding using massively parallel design, synthesis, and testing. Science 2017, 357, 168–175. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  78. Zhang, S.-Q.; Chino, M.; Liu, L.; Tang, Y.; Hu, X.; DeGrado, W.F.; Lombardi, A. De Novo Design of Tetranuclear Transition Metal Clusters Stabilized by Hydrogen-Bonded Networks in Helical Bundles. J. Am. Chem. Soc. 2018, 140, 1294–1304. [Google Scholar] [CrossRef] [PubMed] [Green Version]
Figure 1. Methodology GRSA-SSP for peptide prediction.
Figure 1. Methodology GRSA-SSP for peptide prediction.
Mca 26 00039 g001
Figure 2. Comparison of GRSA0 versus GRSA0-SSP.
Figure 2. Comparison of GRSA0 versus GRSA0-SSP.
Mca 26 00039 g002
Figure 3. Comparison of GRSA1 versus GRSA1-SSP.
Figure 3. Comparison of GRSA1 versus GRSA1-SSP.
Mca 26 00039 g003
Figure 4. Comparison of GRSAE versus GRSAE-SSP.
Figure 4. Comparison of GRSAE versus GRSAE-SSP.
Mca 26 00039 g004
Figure 5. Comparison of GRSA2 versus GRSA2-SSP.
Figure 5. Comparison of GRSA2 versus GRSA2-SSP.
Mca 26 00039 g005
Figure 6. Comparison of GRSAX-SSP algorithms.
Figure 6. Comparison of GRSAX-SSP algorithms.
Mca 26 00039 g006
Figure 7. Comparison of the average time of the 30 execution of GRSAX-SSP algorithms.
Figure 7. Comparison of the average time of the 30 execution of GRSAX-SSP algorithms.
Mca 26 00039 g007
Figure 8. GRSA2-SSP according to the type of secondary structure.
Figure 8. GRSA2-SSP according to the type of secondary structure.
Mca 26 00039 g008
Figure 9. Comparison of GRSA2-SSP, PEP-FOLD3, and I-TASSER by RMSD (up to 15 amino acids). Figure 9 (a) best TM-score and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD average of the five best predictions, (e) GDT-TS average.
Figure 9. Comparison of GRSA2-SSP, PEP-FOLD3, and I-TASSER by RMSD (up to 15 amino acids). Figure 9 (a) best TM-score and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD average of the five best predictions, (e) GDT-TS average.
Mca 26 00039 g009
Figure 10. Comparison of GRSA2-TBM, PEP-FOLD3, and I-TASSER by TM-score (16 to 30 amino acids). Figure 10 (a) best TM-score, and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD of the five best predictions, and (e) GDT-TS average of the five best predictions.
Figure 10. Comparison of GRSA2-TBM, PEP-FOLD3, and I-TASSER by TM-score (16 to 30 amino acids). Figure 10 (a) best TM-score, and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD of the five best predictions, and (e) GDT-TS average of the five best predictions.
Mca 26 00039 g010
Figure 11. Comparison of GRSA2-SSP, PEP-FOLD3, I-TASSER, QUARK, and Rosetta by TM-Score (31 to 49 amino acids). Figure 11 (a) best TM-score, and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD average of the five best predictions, and (e) GDT-TS average of the five best predictions.
Figure 11. Comparison of GRSA2-SSP, PEP-FOLD3, I-TASSER, QUARK, and Rosetta by TM-Score (31 to 49 amino acids). Figure 11 (a) best TM-score, and (b) their RMSD, (c) TM-score average of the five best predictions, (d) RMSD average of the five best predictions, and (e) GDT-TS average of the five best predictions.
Mca 26 00039 g011
Table 1. Set of instances (peptides).
Table 1. Set of instances (peptides).
PDB-CodeaaNumber of Variables
(Torsion Angles)
PDB-CodeaaNumber of Variables
(Torsion Angles)
11uao1047241wz423123
21egs949251rpv17124
31eg41361261pef18124
41l3q1262271du120134
52evq1266281pei22143
61le11269291yyb27160
71in31274301t0c31163
83bu31474311by027193
91gjf1479322bn633200
101rnu1381331wr436206
111lcx1381341yiu37206
121k431484352ysh40213
131a131485361bhi38216
141nkf1686371i6c39218
151le31691381wr741222
161pgbF1693392dmv43229
171dep1594401bwx39242
181niz1697411f4i45276
192bta15100421dv047279
201l2y20100431ify49290
211e0q17109442p8144295
221b0318109451pgy47304
231wbr17120----
Table 2. Results obtained by GRSA2-SSP.
Table 2. Results obtained by GRSA2-SSP.
PDB CodeaaSSRMSD RMSD AveRMSD stdTM1 BestTM1 AveTM1 stdPDB
Code
aaSSRMSD RMSD AveRMSD stdTM1 BestTM1 AveTM1 std
11egs9N1.470.7280.7370.4110.36300.043241pef18A1.50.7060.4680.6860.6610.014
21uao10B0.711.2140.8280.4010.3750.022251l2y20A0.772.2680.9140.2580.2430.008
31l3q12N1.551.4860.7270.2710.2520.025261du120A1.131.620.4630.2660.2660.001
42evq12B2.431.2741.0200.3820.3180.031271pei22A2.021.430.3660.3790.3640.010
51le112B0.381.3561.2080.3160.3010.011281wz423A2.662.660.4240.2720.2650.015
61in312A1.071.0540.3410.3950.3870.007291yyb27A1.471.750.3060.3970.3950.002
71eg413N1.591.6320.3970.3390.3300.006301by027A1.161.440.2170.4130.4080.003
81rnu13A0.260.2880.0330.6280.6160.010311t0c31N2.733.040.3440.2160.20.009
91lcx13N1.081.4120.4220.3340.3230.009322bn633A2.172.330.220.3290.3190.010
103bu314N1.021.1220.470.2940.2630.019331wr436B3.183.090.550.2430.210.018
111gjf14A1.370.8740.4610.5610.5470.040341yiu37B3.013.170.4550.2210.2020.011
121k4314B2.921.4880.9160.3030.2610.027351bhi38N2.762.7360.7940.3060.2960.007
131a1314N1.381.290.1260.3130.3020.007361i6c39B4.293.510.5050.2050.1910.010
141dep15A0.980.7620.3520.6410.6030.023371bwx39A2.982.580.2820.4510.4430.005
152bta15N2.471.7160.4550.2270.1960.018382ysh40B3.213.460.4930.2430.2220.016
161nkf16A3.031.8380.8420.2870.2780.009391wr741B3.713.550.1460.2230.2080.011
171le316B1.021.250.770.2240.2150.007402dmv43B3.273.4020.60.2170.2010.013
181pgbF16B1.542.030.4090.2290.2090.018412p8144A3.523.210.4760.1850.1780.007
191niz16B2.41.770.5720.2350.2140.016421f4i45A3.133.460.2210.310.3020.006
201e0q17B0.791.4940.5360.2260.2210.008431dv047A2.652.940.4370.3030.2830.011
211wbr17N1.681.310.3630.2950.27160.016441pgy47A3.222.620.460.3450.3360.006
221rpv17A0.810.710.0960.4690.4630.005451ify49A2.562.770.40.3110.2970.008
231b0318B3.042.3560.6290.21430.2080.004----------
Note: PDB code (Instance), number of amino acids (aa), SS is the predominant secondary structure type: beta strand (B), alpha-helix (A) and none (N), TM1 = TM-score.
Table 3. Mini-proteins.
Table 3. Mini-proteins.
Instances
PDB Codeaa Number of Variables
(Torsion Angles)
SS
15wll26174A
25lo234192A
32ki036214N
45up540266N
55uoi43282A
62kik48306A
Note: alpha-helix (A) and none (N) for secondary structure.
Table 4. Average metrics results of Mini-proteins.
Table 4. Average metrics results of Mini-proteins.
ApproachesInstances
5wll5lo22ki0
RMSDTM-ScoreGDT-TSRMSDTM-ScoreGDT-TSRMSDTM-ScoreGDT-TS
GRSA2-SSP0.656 *0.642 *0.944 *1.504 *0.5010.6492.1720.3540.504
PEP-FOLD31.0740.5260.8921.9220.5320.7692.4220.4660.697
I-TASSER0.8230.5300.7371.7340.6080.7760.620 *0.899 *0.986 *
QUARK0.8970.5650.7881.8480.5270.7132.2280.4500.688
RosettaN/AN/AN/A1.5520.694 *0.849 *2.1460.4600.710
ApproachesInstances
5up55uoi2kik
RMSDTM-ScoreGDT-TSRMSDTM-ScoreGDT-TSRMSDTM-ScoreGDT-TS
GRSA2-SSP2.2340.2770.4033.1940.1920.3402.7560.3390.508
PEP-FOLD32.5120.3720.5412.5160.4810.6292.2820.3950.597
I-TASSER1.390 *0.782 *0.900 *2.5650.5120.6642.1870.4480.557
QUARK1.8800.6140.7782.0220.6330.7772.0280.4620.627
Rosetta1.7160.6920.8381.642 *0.753 *0.871 *1.968 *0.665 *0.785 *
Note: The asterisk (*) represents the best result in each column.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Sánchez-Hernández, J.P.; Frausto-Solís, J.; González-Barbosa, J.J.; Soto-Monterrubio, D.A.; Maldonado-Nava, F.G.; Castilla-Valdez, G. A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing. Math. Comput. Appl. 2021, 26, 39. https://doi.org/10.3390/mca26020039

AMA Style

Sánchez-Hernández JP, Frausto-Solís J, González-Barbosa JJ, Soto-Monterrubio DA, Maldonado-Nava FG, Castilla-Valdez G. A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing. Mathematical and Computational Applications. 2021; 26(2):39. https://doi.org/10.3390/mca26020039

Chicago/Turabian Style

Sánchez-Hernández, Juan P., Juan Frausto-Solís, Juan J. González-Barbosa, Diego A. Soto-Monterrubio, Fanny G. Maldonado-Nava, and Guadalupe Castilla-Valdez. 2021. "A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing" Mathematical and Computational Applications 26, no. 2: 39. https://doi.org/10.3390/mca26020039

APA Style

Sánchez-Hernández, J. P., Frausto-Solís, J., González-Barbosa, J. J., Soto-Monterrubio, D. A., Maldonado-Nava, F. G., & Castilla-Valdez, G. (2021). A Peptides Prediction Methodology for Tertiary Structure Based on Simulated Annealing. Mathematical and Computational Applications, 26(2), 39. https://doi.org/10.3390/mca26020039

Article Metrics

Back to TopTop