Next Article in Journal
Supporting Structure of Steel Corrugated Plate-Mold Bag Concrete and Its Application in a Circular Shaft
Previous Article in Journal
Multiport Energy Management System Design for a 150 kW Range-Extended Towing Vessel
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

Analysis of Structure–Activity Relationships of Food-Derived DPP IV-Inhibitory Di- and Tripeptides Using Interpretable Descriptors

by
Monika Hrynkiewicz
1,†,
Anna Iwaniak
1,*,†,
Piotr Minkiewicz
1,
Małgorzata Darewicz
1 and
Wojciech Płonka
2
1
Chair of Food Biochemistry, Faculty of Food Science, University of Warmia and Mazury in Olsztyn, Pl. Cieszyński 1, 10-726 Olsztyn-Kortowo, Poland
2
FQS Poland (Fujitsu Group), Parkowa 11, 30-438 Kraków, Poland
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Appl. Sci. 2023, 13(23), 12935; https://doi.org/10.3390/app132312935
Submission received: 24 October 2023 / Revised: 29 November 2023 / Accepted: 30 November 2023 / Published: 4 December 2023
(This article belongs to the Section Food Science and Technology)

Abstract

:
This study aimed to analyze the structural requirements for di- and tripeptides exhibiting a DPP IV-inhibitory effect. The sequences of 46 di- and 33 tripeptides, including their bioactivity (IC50; μM), were implemented from the BIOPEP-UWM database, whereas modeling was performed using SCIGRESS Explorer: Version FJ 3.5.1 software. Models included 336 (dipeptide dataset) and 184 descriptors (tripeptide dataset). The values of the determination coefficient (R2) defining model reliability were 0.782 and 0.829 for di- and tripeptides, respectively. Based on the implemented descriptors, it was concluded that increased numbers of nitrogen atoms, as well as the methyl groups, are required for dipeptides to enhance the DPP IV-inhibitory effect. This was indicated by the presence of amino acids with an aliphatic side chain (e.g., Leu, Val, Ile) and an aromatic ring (Trp). In the case of tripeptides, a correlation was found between their molecular weight (MW) and studied bioactivity. A tripeptide with a molecular weight of up to 500 Da was found suitable for the sequence to act as the DPP IV inhibitor. Although there is still a gap in explaining the relations between the structural nature and the DPP IV-inhibitory activity of peptides, and certain issues related to this topic still remain unknown, the results are in line with those reported by other authors. Additionally, the suitability of the SCIGRESS tool in the QSAR analysis of peptides derived from foods can be confirmed. Interpretable descriptors enabled the achievement of more unequivocal results concerning the main structural factors affecting the DPP IV inhibition of di- and tripeptides.

1. Introduction

It has recently been highlighted that food, as well as its main function, i.e., providing essential nutrients that regulate bodily functions, is also expected to prevent, or even cure, diseases. Its health potential results from the presence of bioactive compounds like, e.g., vitamins, phytochemicals, and peptides [1]. The latter ones have attracted the attention of scientists as non-toxic and multi-active molecules [2]. They were proven to act as, i.a., glucose level reductants, blood hypotensors, cholesterol level reductants, anti-oxidants [3,4], anti-inflammatory, and anti-obesity agents [5].
Reduction of the blood glucose level is related to, i.a., the action of peptides that inhibit dipeptidyl peptidase IV (DPP IV). DPP IV itself (EC 3.4.14.5) is an enzyme involved in the degradation of GLP 1(glucagon-like peptide 1) and GIP (glucose-dependent insulinotropic peptide). These two molecules are released into the gastrointestinal tract and are responsible for glucose homeostasis. The action of DPP IV briefly relies on the cleavage of N-terminal dipeptides from GLP 1 and GIP to produce inactive incretins; hence, it triggers the loss of insulinotropic activity of GLP 1 and GIP. Thus, peptides that inhibit DPP IV action display the anti-diabetic function [6,7].
DPP IV inhibitors were identified in, e.g., milk and dairy products, rice bran, salmon, tuna, and azuki beans [8]. The hydrolysis of proteins followed by the determination of biological activity of a hydrolysate and the identification and/or synthesis of peptides is the classical procedure of studying bioactive peptides derived from foods. Another research strategy includes an in silico analysis of biopeptides. It involves the application of, i.a., bioinformatic and cheminformatic databases as sources of information about peptide sequences, tools for the molecular docking of peptides to predict the possible mechanism of their action, and the prediction of peptide bioactivity based on its structural nature [9]. The latter is called a QSAR, i.e., the quantitative structure–activity relationship, and, so far, it has been applied in the study of peptides showing a one [9,10,11] and dual biological effect [12,13].
A literature search using the Scopus [14] platform (accession date: October 2023) revealed dozens of records (i.e., 51, including reviews and original papers as well as conference abstracts) dedicated to QSAR of DPP IV inhibitors using “QSAR and dipeptidyl peptidase 4 and inhibitors” as a query. Further analysis of these records showed that none of them concerned peptides with such an activity. The articles mostly concerned the QSAR analysis of non-peptide DPP IV inhibitors like, e.g., 4-fluoropyrrolidine-2-carbonitrile and octahydrocyclopenta[b]pyrrole-2-carbonitrile [15], tetralin-sulfonamide derivatives [16], and polyphenols [17]. This literature search via Scopus has shown that there is paucity of literature data related to the QSAR of DPP IV peptide inhibitors, especially those derived from food proteins. The query “QSAR and dipeptidyl peptidase-IV and inhibition or inhibitors and peptides” revealed 15 records (as of October 2023), but barely 2 of them were related to the food-originating DPP IV inhibitors. They included an original article concerning the identification of DPP IV inhibitors from camel milk hydrolysates [18] and a review of the efficiency of the QSAR of food-derived bioactive peptides [19]. One of the latest articles concerned the impact of certain structural features on the DPP IV-inhibitory activity of dipeptides. However, the models were based on several attributes describing protein structures (long-chain molecules) which are not applicable to dipeptides (short-chain compounds) [20]. Other articles were related to DPP IV inhibition described by means of programs useful for reporting the DPP IV inhibition of compounds like, e.g., iDPPIV-SCM and StackDPPIV [21,22]. Hence, even though there is ample literature data on the potential of bioactive peptides, including DPP IV inhibitors, there is a still a need to revise and/or systemize knowledge about DPP IV inhibitors. This knowledge gap has stemmed from: (a) the not-fully elucidated mechanism of DPP IV inhibition induced by peptides [6]; (b) the not-fully explained structure–activity relationship of these peptides [6]; and c) an increase in the number of new biopeptides (including DPP IV inhibitors) deposited in databases [9]. Therefore, taking into account the latter of the above-mentioned issues in particular, this study aimed to predict the DPP IV-inhibitory activity of di- and tripeptides based on their structural nature using interpretable descriptors.

2. Materials and Methods

2.1. Sequences of DPP IV Inhibitors

The sequences of di- and tripeptides with DPP IV-inhibitory activity were acquired from the BIOPEP-UWM database of protein and bioactive peptide sequences [23], available at: https://biochemia.uwm.edu.pl/en/biopep-uwm-2/ (accessed on 12 June 2023). The complete list of DPP IV inhibitors, including their bioactivity measure expressed as the concentration of peptides corresponding to their half-maximal activity (IC50; μM), is shown in Results and Discussion (see the appropriate Tables). To summarize, 46 dipeptides and 33 tripeptides were used for QSAR modeling. Data on peptides available in the BIOPEP-UWM database were acquired in June 2023.

2.2. Analysis of Structure–Activity Relationship of DPP IV Inhibitors

The procedure of QSAR modeling of di- and tripeptide DPP IV inhibitors using BIOPEP-UWM database (see Section 2.1) and SCIGRESS Explorer: Version FJ 3.5.1 is shown in Figure 1.
Briefly, QSAR models calculated using SCIGRESS software were based on MLR. All the structures of di- (46 sequences) and tripeptides (33 sequences) were drawn using the above-mentioned software by means of the following tabs available in the applied program: Action → Protein → Sequence. They were protonated to obtain the zwitterion by exemplifying the occurrence of peptides in a solution at a neutral pH. The geometries of peptides were optimized using the SCIGRESS option Experiment → Run by selecting “chemical sample” and the procedure called “MO-G PM 6 in water”. QSAR analyses were initially computed using a method called “Complete Topological QSAR” (see Table 1), taking into account the function called the Enhanced Replacement Method (ERM) (denoted as the ERM model) [24] implemented in the SCIGRESS. A parameter called Descriptor Count was set to 5 in order to obtain an acceptable ratio of the number of descriptors to the number of experimental points. The initial models were refined using the “Custom QSAR” method in order to manually remove descriptors which are difficult to interpret mechanistically while retaining satisfying statistical performance parameters. Finally, the QSAR model created for dipeptides was built up from 46 sequences and 336 potentially usable descriptors, whereas 33 sequences and 184 descriptors were used for tripeptides modeling. The list of descriptors, as well as the complete reports concerning the models, may be found in the Supplementary Materials. The IC50 values (µM), defining the DDP-IV-inhibitory activity of peptides, were converted to log10. The software applied in our study automatically adjusts the measure of bioactivity to produce optimized results.

3. Results and Discussion

The structure–analysis relationship approach relies on the construction of mathematical models that enable computing and finding (if possible) the correlation between the structure of a compound and its bioactivity [6]. In the case of peptides, the structure is understood as their amino acid sequences, whereas bioactivity refers to the measure of a molecule activity [25]. Usually, it is the concentration of a compound that corresponds to its-half maximal effect (like, e.g., IC50). Moreover, such a measure is usually transformed into logarithms [26].
The calculations were made using SCIGRESS software, as the QSAR model development in SCIGRESS is GUI-driven (i.e., graphical user interface), that is, by a user-interactive process with the possibility of manual removal of hard-to-mechanistically-interpret descriptors. This is an advantage of this software over other tools, like, e.g., RDKit/Scikit-learn Python scripts [27,28], where the feature (descriptor) selection process is mostly automated and based on statistical parameters, with no simple way of taking descriptor mechanistic interpretability.
Analyses concerned 46 di- and 33 tripeptide DPP IV inhibitors. According to Tropsha [29], the number of samples (i.e., cases, objects) in the datasets should not be too small or too big. The low number of samples results in false-positive or false-negative results, which may lead to uncertain conclusions in biological and medical sciences [30]. The upper limit of the samples is often defined by the computer and the method by which the QSAR model was constructed [29]. The analysis using SCIGRESS was based on multilinear regression (MLR), i.e., a traditional approach in QSAR [31]. According to Cargnelutti Filho et al. [32], in MLR, the sample size depends on effect size, whereas the number of independent variables and recommendations for the number of chemicals are associated with the different criteria adopted by each researcher. Based on information from the theoretical model selection, Jenkins and Quintana-Asencio [30] tried to recommend the minimum N (understood as the number of samples) to correctly match a model to a data shape in regression models. In the words by Jenkins and Quintana-Asencio [30], “statistical limbo may be better avoided and reproducibility may be improved if research based on regressions and meta-regressions uses N ≥ 25” [30]. Apart from the number of molecules (samples), the number of descriptors (variables) in a model is important. To conclude, the number of peptides taken for QSAR analysis in the present study is in line with the above statement. However, due to a relatively small dataset, there was no possibility of performing a typical train/test dataset split. Instead, leave-one-out validation was used to measure the performance of the models (see discussion below).
Table 1 shows the statistical results of analyses obtained for the di- and tripeptide DPP IV-inhibitory models. For both models, similar statistics were obtained in terms of the R2 values, which were R2 = 0.782 and R2 = 0.829 for di- and tripeptides, respectively. R2 is a statistical parameter that shows which percent of endogenous variable is explained by exogenous variable(s) [33]. Endogenous and exogenous variables refer to the bioactivity of peptides and descriptors expressing their properties, respectively [14]. On the other hand, the parameter called adjusted R2 (see Table 1) includes degrees of freedom, and its value is usually smaller than that of R2 [16] (0.755 and 0.798, respectively (see Table 1)).
The validity of the model was calculated by default using the leave-one-out (LOO) cross-validation procedure. The LOO procedure is the one of the simplest statistical procedures for model validation and relies on excluding each sample once and constructing the model without this sample to predict the value of a dependent variable (bioactivity) [34]. Moreover, the LOO procedure enables nearly unbiased estimation of the true generalization performance [35]. The statistical parameters of LOO cross-validation are represented by CVR2 (correlation coefficient of leave-one-out cross-validation, also called Q2), RMSEC (root mean square error of calibration), and RMSECV (root mean square error of cross-validation). For di- and tripeptides, the CVR2 values were 0.817 and 0.722, respectively (see Table 1), whereas the respective RMSECV values were 0.377 and 0.387. Usually, CVR2 > 0.5 is found sufficient to prove model accuracy [36]. Kiralj and Fereira [34] reported that relations of the above-mentioned statistics when using the LOO procedure should be as follows: R2 > Q2 and RMSEC < RMSECV. The statistics present in Table 1 meet these requirements. Moreover, according to Kiralj and Fereira [34], the sufficient requirements for regression models include the following conditions: CVR2 > 0.5 and R2 > 0.6. The difference between R2 and CVR2 exceeding 0.2–0.3 indicates the model overfitting. Based on these findings, the models presented in this study could be considered statistically acceptable (see Table 1). The statistical significance of the estimated QSAR models was analyzed based on F-statistics (α = 0.05), the values of which reached 28.773 and 26.228 for di- and tripeptides, respectively. The p-value was lower than α. To summarize, the results of computations made for CVR2, F-statistics, and p-values indicate some significant relationships between descriptors and bioactivity (logIC50) of both DPP IV-inhibitory QSAR models and enable developing equations (see Figure 2a,b) for predicting the relationship between the structure and DPP IV-inhibitory activity of dipeptides/tripeptides. The names of the descriptors found in these equations are provided in the exact words as are given in the program used.
The impact of individual descriptors on the DPP IV activity of di- and tripeptides is shown in Table 2. The respective plots referring to the correlations between the observed and the predicted activity of di- and tripeptide DPP IV inhibitors are shown in Figure 3 and Figure 4, respectively. The predicted and observed values of logIC50 for di- and tripeptides (tripeptides) are provided in Table 3 and Table 4, respectively.
Statistical parameters present in Table 2 can be explained as the weights showing the contribution of each descriptor as a factor influencing the activity of DPP IV inhibitors. Partial-F indicates the importance of a descriptor in the SAR prediction. The higher its value is, the higher the significance of the descriptor. The relative weights of the normalized coefficients are the measures of the impact of the specific property on the peptide structure. The negative values (bold in Table 2) of normalized coefficients suggest that the property improves the bioactivity of peptides, whereas the positive values mean that the property should rather be minimized (and/or eliminated) from the structure of a molecule in order to “improve” its biological function [37]. To recapitulate, the advantageous impact on dipeptide DPP IV-inhibitory activity was correlated with the descriptors defined as “nitrogen count2” and “methyl count2” (having been squared). Disadvantageous impact was correlated with the “ring count 5” (square value). In the case of tripeptides, the descriptors improving their bioactivity were: “methyl count2” (having been squared) and “sqrt(molecular weight)”. Based on the normalized coefficient value, the importance of the first descriptor can be ignored as similar to that of the descriptor called “Csp3 bonded to 3C” (the higher its value, the lesser its impact). Some descriptors summarized in Table 2 are counts of simple structure features, possess only integer values, and may thus be understood as molecular quantum numbers according to a definition proposed by Nguyen et al. [38]. Descriptors proposed here include, i.a., nitrogen count, methyl count, ring count 5 member, and Csp3 bonded to 3 C. Among them, the ring count 5 member, i.e., the number of rings containing five atoms, is one of the molecular quantum numbers proposed by Nguyen and coworkers [38]. The above authors have divided nitrogen atoms into two classes: acyclic and cyclic nitrogens. In this work, it was enough to include the total count of nitrogen atoms in a molecule. The set of descriptors proposed by Nguyen et al. [38] is designed to enable the classification of all classes of low-molecular compounds, whereas the set analyzed in this study is adapted to work with a limited number of short peptides, with limited variability of backbones and a limited number of possible side chains. Peptide fragments taken into account while constructing a set of descriptors, summarized in Table 2, are presented in Figure 5 using four peptides as examples.
Focusing on the descriptors enhancing the DPP IV-inhibitory activity of dipeptides, it can be concluded that peptides rich in methyl groups as well as nitrogen might be greater enhancers of DPP IV inhibition. Data presented in Table 3 show 23 dipeptides out of the 46 in total whose predicted (ERM variant (see Section 2)) and observed values of logIC50 were comparable. These peptides were as follows: GL, GP, KP, ML, MP, NH, RP, SL, TP, VL, VP, VR, WC, WI, WM, WN, WP, WQ, WR, WS, WT, WW, and WY. When looking at their sequences, some of them contained an amino acid with a non-polar side chain (e.g., G, V, L, I), which can be represented as the presence of methyl groups (the higher their number, the better their activity) in a dipeptide chain. Several dipeptides were composed of tryptophan (W; 11 peptides out of 23), asparagine (N), and/or glutamine (Q). Compared to other amino acids, these three residues contain “additional” nitrogen. Thus, the presence of one of the above-mentioned amino acids promoting the DPP IV-inhibitory activity can be explained by the descriptor defined as “nitrogen count2” (the higher the number, the more favorable the impact on dipeptide bioactivity). Three dipeptides (RP, VR, WR) were composed of arginine (R). Although the presence of guanidine (R) increases the number of nitrogens in a dipeptide chain, the statistical parameters (see Table 3), as well as the equation (Figure 2a), indicated that the number of guanidine groups should be reduced (or eliminated) to improve the DPP IV-inhibitory effect of dipeptide. Moreover, RP and WR peptides—numbers 22 and 40 in the plot (see Figure 3)—were located out of the confidence interval line, whereas VR (number 31) was found in the confidence interval line (dashed line in Figure 3). Some of the dipeptides were built up of proline (P), a residue containing a five-atom aliphatic ring. According to the statistics (see Table 3), the lesser the number of such rings in the peptide chain, the stronger the impact on the DPP IV-inhibitory effect.
The analysis of tripeptide DPP IV inhibitors showed that the following 16 tripeptides out of the 33 sequences analyzed had comparable values of the observed and predicted logIC50 (ERM variant): APG, GPA, PPG, PPL, WRF, WRG, WRH, WRI, WRL, WRP, WRQ, WRR, WRS, WRT, WRY, and WWW (see Table 4 and Figure 4). As was mentioned above, due to the minor impact of such descriptors as the number of -CH3 groups squared (“methyl count2”), as well as the number of sp3 carbon atoms bound to the other three carbon atoms (“Csp3 bonded to 3 C”), they were not further considered as factors affecting the DPP IV-inhibitory activity of tripeptides. Finally, the bioactivity of tripeptides was found to be related to their molecular weight (MW). Moreover, based on the normalized coefficients (see Table 3), it can be concluded that the higher molecular weight of a tripeptide is more beneficial to its DPP IV-inhibitory effect. When looking first at the sequences of tripeptides with comparable values of the predicted and observed logIC50, the majority of them contained tryptophans (W), prolines (P), and/or tyrosines (Y), which made them heavier than the other tripeptide sequences. Nevertheless, the next step of the research was to create the non-linear function of three descriptors related to molecular weight (see Figure 6a), with the plot (see Figure 6b) showing the relationship between the predicted logIC50 and the molecular weight of a tripeptide. This function was applied for the first time to analyze the impact of specific descriptors of DPP IV inhibition by peptides. The plot enables the distinguishing of the minimum and the maximum of the function. The minima and maxima of the functions have a long history in terms of their practical applications in life to solve the problems of the so-called “optimizations” [39]. According to the plot, the minimum of the function was reached by peptides IPI and LPL (number 6 and 8, respectively). It can be concluded that the minimal structural requirements for a tripeptide to act as a DPP IV inhibitor is for it to be composed of hydrophobic amino acids (Leu, Ile) and the ones with rings (Pro). According to Guasch et al. [40], the presence of hydrophobic residues may enhance the interactions between the ligand (peptide) and the enzyme (DPP IV). Moreover, the S1 subsite of DPP IV was composed of hydrophobic pockets, which are important in enzyme inhibition by peptides [40]. Liu et al. [6] reported that DPP IV possessed two pockets, namely, S1, consisting of Tyr631, Val656, Trp659, Tyr662, Tyr666, and Val711; and S2, which was charged due to the presence of Arg125, Glu205, Glu206, Phe357, Ser209, and Arg358. According to other scientific data, three pockets can be distinguished in DPP IV: S1 (Tyr547, Ser630, Tyr631, Val656, Trp659, Tyr662, Asn710, Val711, His740), S2 (Glu205, Glu206, and Tyr662), and S3 (Ser 209, Arg358, Phe357) [6]. The tripeptides from the studied dataset that reached the maximum on the plot were WRD, WRI, WRL, and WRN (numbers 15, 20, 22, and 24 in the plot, respectively). The impact of Trp and/or aliphatic amino acids, like, e.g., Ala or Ile, has already been discussed (see above). In the case of hydrophilic residues, like Arg (R) and/or Asn (N), their role in DPP IV inhibition remains unknown [21]; but, as it could be seen above, Arg is a residue forming the charged pocket of DPP IV [6].
The molecular weight of the tripeptides that reached the minimum of the function (see Figure 6b) was ca. 345 Da, whereas the MW of the tripeptides reaching the maximum of the function ranged from ca. 473 to 476 Da. According to Liu et al. [6], the molecular weight of the peptides is one of the factors affecting, e.g., the bioavailability of peptides, mode of action, and/or bioactivity. Liu et al. [6] analyzed the molecular weight of 222 peptides with different chain lengths for their IC50 values (μM). It was found that peptide DPP IV inhibitors with MW lesser than 500 Da had IC50 lower than 1000 μM (logIC50 = 3). These findings are in line with the results obtained in the present study, indicating that predicted logarithmic values of IC50 of tripeptides located in the minimum and maximum of the plot (Figure 6b) ranged from 1.45 (IPI; minimum of the function) to 2.82 (WRI; maximum of the function).
Among the results presented in Figure 6, four peptides with the lowest molecular weight (MW) had the following predicted log IC50 values: 1.45 (IPI and LPL), 1.47 (VPL), and 2.41 (LQP). There are no significant differences in predicted log IC50 (MW) between peptides with a molecular weight lower or higher than 500 Da. The molecular weight of 500 Da is the upper limit for compounds recommended as orally delivered drugs according to the Rule of Five [41]. This rule may also be applied to the bioactive components of food. Recommendation takes into account the fact that molecules with too high a molecular weight and containing too many groups able to form hydrogen bonds have lower bioavailability than the smaller ones. On the other hand, peptides too long to fulfill the Rule of Five are considered to be compounds with relatively high bioavailability [42,43]. Peptides with molecular weight exceeding 500 Da may thus be considered objects of interest as bioactive compounds together with the smaller ones.
Comparison of Figure 4 and Figure 7 shows that the prediction of log IC50 based on a wider variety of descriptors provides better correlation between predicted and observed values than the calculation taking into account predictors associated with molecular weight only. Log IC50 (ERM) reveals a better correlation with experimental results than log IC50 (MW).
When discussing the results, it is worth mentioning that the SCIGRESS program was applied for the first time to study the SAR of DPP IV-inhibitory peptides. Its advantage is the possibility of refining the models by removing the descriptors which are difficult to interpret in order to produce reliable statistics (see Section 2). Nevertheless, the outputs obtained using SCIGRESS and the implemented descriptors are comparable with the results of Nongonierma et al. [44], who confirmed that some characteristic residues occurred in different locations of peptide DPP IV inhibitors. Dipeptide DPP IV inhibitors were often rich in Trp, which was also confirmed in this study. Additionally, amino acids—like Ala, Leu, and/or His—were typical for dipeptides capable of DPP IV inhibition [44]. This regularity was also found in the present study, in which the presence of the above-mentioned amino acids in a dipeptide sequence was defined by descriptors related to the number of methyl groups and the number of nitrogens (the higher this number is, the better the impact on bioactivity).
The analyses of tripeptide SAR were focused on the impact of their molecular weight on DPP IV inhibition. Nongonierma et al. [44] demonstrated more the diversified composition of tripeptides acting as DPP IV inhibitors, but some specific residues were indicated as well, including, Gly, Ala, Leu, Ile, Gln, Met, Phe, and Pro, which were found in different locations of the sequence. The results achieved using SCIGRESS software showed a more specific focus on the presence of Trp in the tripeptide chain to enhance the DPP IV-inhibiting effect.
To summarize, QSAR studies of compounds have been known since 1962 and are still being developed with increasingly sophisticated methods of computation to become a standard approach in drug discovery due to their efficiency in prediction and low cost of analyses [45]. In the area of the QSAR of peptides, particularly those of food origin, studies were undertaken on the sequences showing ACE-inhibitory [46], bitter-tasting [47], and anti-oxidative [11] effects. To date, there has been little research on QSAR of peptide DPP IV inhibitors. This may be due to the fact that although new peptides with such an activity are being identified, it is still not enough to merit creating a dataset for QSAR analysis. The results obtained using SCIGRESS bring more insight into the structural nature of peptides acting as DPP IV inhibitors.

4. Conclusions

The applied descriptors proved useful in revealing some structural requirements for food-derived di- and tripeptides to act as DPP IV inhibitors. Compared to the findings reported by other authors, this study results indicate more unequivocal factors affecting this particular bioactivity of peptides. For example, an increasing number of nitrogen atoms, as well as methyl groups, enhanced the DPP IV inhibition of dipeptides, whereas molecular weight turned out to be crucial in the case of tripeptides. Based on the latter, the requirement for tripeptide DPP IV inhibitor was a molecular weight up to 500 Da.

Supplementary Materials

The following supporting information can be downloaded at https://www.mdpi.com/article/10.3390/app132312935/s1. The file contains the list of descriptors as well as the complete reports concerning the QSAR models of di- and tripeptide DPP IV inhibitors.

Author Contributions

Conceptualization, A.I. and M.H.; methodology, W.P.; validation, W.P.; formal analysis, A.I. and M.H.; investigation, A.I. and M.H.; resources, P.M. and M.D.; data curation, P.M.; writing—original draft preparation, A.I., P.M. and M.H.; writing—review and editing, A.I. and M.H.; visualization, P.M. and M.H.; supervision, A.I.; funding acquisition, M.D. All authors have read and agreed to the published version of the manuscript.

Funding

This project was financially supported by Minister of Education and Science in the range of the program entitled “Regional Initiative of Excellence” for the years 2019–2023, Project No. 010/RID/2018/19 (amount of funding: 12.000.000 PLN).

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

All data are presented in the main text and in Supplementary file.

Conflicts of Interest

The authors declare no conflict of interest. Wojciech Płonka works in FQS Poland (Fujitsu Group).The company was not involved in the study design, collection, analysis, interpretation of data, the writing of this article or the decision to submit it for publication.

Abbreviations

Csp3 bonded to 3 CNumber of carbon atoms bound to three carbon atoms with sp3 hybridization (number of tertiary carbon atoms)
DPP IVDipeptidyl peptidase IV (EC 3.4.14.5)
ERMEnhanced Replacement Method
IC50Concentration corresponding to half-maximal inhibition of DPP IV
log IC50Logarithm of IC50
MWMolecular weight
Predicted log IC50 (ERM)Logarithm of IC50 predicted using model based on enhanced replacement method
Predicted log IC50 (MW)Logarithm of IC50 predicted using model based on molecular weight
QSARQuantitative structure–activity relationship
Ring count 5 memberNumber of rings containing five atoms in peptide molecule
UWMUniversity of Warmia and Mazury

References

  1. Tadesse, S.A.; Emire, S.A. Production and processing of antioxidant bioactive peptides: A driving force for the functional food market. Heliyon 2020, 6, e04765. [Google Scholar] [CrossRef] [PubMed]
  2. He, R.; Wang, Y.; Yang, Y.; Wang, Z.; Ju, X.; Yuan, J. Rapeseed protein-derived ACE inhibitory peptides LY, RALP and GHS show antioxidant and anti-inflammatory effects on spontaneously hypertensive rats. J. Funct. Foods 2019, 55, 211–219. [Google Scholar] [CrossRef]
  3. Barati, M.; Javanmardi, F.; Jazayeri, S.M.H.M.; Jabbari, M.; Rahmani, J.; Barati, F.; Nickho, H.; Davoodi, S.H.; Roshanravan, N.; Khaneghah, A.M. Techniques, perspectives, and challenges of bioactive peptide generation: A comprehensive systematic review. Compr. Rev. Food Sci. Food Saf. 2020, 19, 1488–1520. [Google Scholar] [CrossRef] [PubMed]
  4. Dhar, H.; Verma, S.; Dogra, S.; Katoch, S.; Vij, R.; Singh, G.; Sharma, M. Functional attributes of bioactive peptides of bovine milk origin and application of in silico approaches for peptide prediction and functional annotations. Crit. Rev. Food Sci. Nutr. 2023. [Google Scholar] [CrossRef]
  5. Jakubczyk, A.; Szymanowska, U.; Karaś, M.; Złotek, U.; Kowalczyk, D. Potential anti-inflammatory and lipase inhibitory peptides generated by in vitro gastrointestinal hydrolysis of heat treated millet grains. CyTA J. Food 2019, 17, 324–333. [Google Scholar] [CrossRef]
  6. Liu, R.; Cheng, J.; Wu, H. Discovery of food-derived dipeptidyl peptidase IV inhibitory peptides: A review. Int. J. Mol. Sci. 2019, 20, 463. [Google Scholar] [CrossRef]
  7. Zhang, M.; Zhu, L.; Wu, G.; Liu, T.; Qi, X.; Zhang, H. Food-derived dipeptidyl peptidase IV inhibitory peptides: Production, identification, structure-activity relationship, and their potential role in glycemic regulation. Crit. Rev. Food Sci. Nutr. 2023. [Google Scholar] [CrossRef]
  8. Power, O.; Nongonierma, A.; Jakeman, P.; FitzGerald, R. Food protein hydrolysates as a source of dipeptidyl peptidase IV inhibitory peptides for the management of type 2 diabetes. Proc. Nutr. Soc. 2014, 73, 34–46. [Google Scholar] [CrossRef]
  9. Du, Z.; Comer, J.; Li, Y. Bioinformatics approaches to discovering food-derived bioactive peptides: Reviews and perspectives. Trends Anal. Chem. 2023, 162, 117051. [Google Scholar] [CrossRef]
  10. Wu, J.; Aluko, R.E.; Nakai, S. Structural requirements of angiotensin I-converting enzyme inhibitory peptides: Quantitative structure-activity relationship study of di- and tripeptides. J. Agric. Food Chem. 2006, 54, 732–738. [Google Scholar] [CrossRef]
  11. Udenigwe, C.C.; Aluko, R.E. Chemometric analysis of the amino acid requirements of antioxidant food protein hydrolysates. Int. J. Mol. Sci. 2011, 12, 3148–3161. [Google Scholar] [CrossRef] [PubMed]
  12. Pripp, A.H.; Ardö, S. Modelling relationship between angiotensin-(I)-converting enzyme inhibition and the bitter taste of peptides. Food Chem. 2007, 102, 880–888. [Google Scholar] [CrossRef]
  13. Hrynkiewicz, M.; Iwaniak, A.; Bucholska, J.; Minkiewicz, P.; Darewicz, M. Structure-activity prediction of ACE inhibitory/bitter dipeptides-a chemometric approach based on stepwise regression. Molecules 2019, 24, 950. [Google Scholar] [CrossRef] [PubMed]
  14. Available online: https://scopus.com (accessed on 20 October 2023).
  15. Arba, M.; Ruslin, R.; Akib, N.; Yamin, Y.; Ombe, S.; Jessi, J.; Muzakkar, M.Z.; Tjahjono, D.R. Molecular modeling of 4-fluoropyrrolidine-2-carbonitrile and octahydrocyclopenta[b]pyrrole-2-carbonitrile as a Dipeptidyl Peptidase IV (DPP4) Inhibitor. J. App. Pharm. Sci. 2018, 8, 1–7. [Google Scholar] [CrossRef]
  16. Abd El-Karim, S.S.; Anwar, M.M.; Syam, Y.M.; Nael, M.A.; Ali, H.F.; Motaleb, M.A. Rational design and synthesis of new tetralin-sulfonamide derivatives as potent anti-diabetics and DPP-4 inhibitors: 2D &3D QSAR, in vivo radiolabeling and bio distribution studies. Bioorg. Chem. 2018, 81, 481–493. [Google Scholar] [CrossRef] [PubMed]
  17. Riyaphan, J.; Pham, D.C.; Leong, M.K.; Weng, C.F. In Silico Approaches to Identify Polyphenol Compounds as α-Glucosidase and α-Amylase Inhibitors against Type-II Diabetes. Biomolecules 2021, 11, 1877. [Google Scholar] [CrossRef] [PubMed]
  18. Nongonierma, A.B.; Paolella, S.; Mudgil, P.; Maqsood, S.; FitzGerald, R.J. Identification of novel dipeptidyl peptidase IV (DPP-IV) inhibitory peptides in camel milk protein hydrolysates. Food Chem. 2018, 244, 340–348. [Google Scholar] [CrossRef]
  19. Nongonierma, A.; FitzGerald, R. Learnings from quantitative structure activity relationship (QSAR) studies with respect to food protein-derived bioactive peptides: A review. RSC Adv. 2016, 79, 6. [Google Scholar] [CrossRef]
  20. Kęska, P.; Stadnik, J. Structure-activity relationships study on biological activity of peptides as dipeptidyl peptidase IV inhibitors by chemometric modeling. Chem. Biol. Drug Des. 2020, 95, 291–301. [Google Scholar] [CrossRef]
  21. Charoenkwan, P.; Kanthawong, K.; Nantasenamat, C.; Hasan, M.; Shoombuatong, W. iDPPIV-SCM: A sequence-based predictor for identifying and analyzing dipeptidyl peptidase IV (DPP-IV) inhibitory peptides using a scoring card method. J. Proteome Res. 2020, 19, 4125–4136. [Google Scholar] [CrossRef]
  22. Charoenkwan, P.; Nantasenamat, C.; Hasan, M.; Moni, M.A.; Lio’, P.; Manavalan, B.; Shoombuatong, W. StackDPPIV: A novel computational approach for accurate prediction of dipeptidyl peptidase IV (DPP-IV) inhibitory peptides. Methods 2022, 204, 189–198. [Google Scholar] [CrossRef] [PubMed]
  23. Minkiewicz, P.; Iwaniak, A.; Darewicz, M. BIOPEP-UWM database of bioactive peptides: Current opportunities. Int. J. Mol. Sci. 2019, 20, 5978. [Google Scholar] [CrossRef] [PubMed]
  24. Mercader, A.G.; Duchowicz, P.R.; Fernández, F.M.; Castro, E.A. Replacement method and enhanced replacement method versus the genetic algorithm approach for the selection of molecular descriptors in QSPR/QSAR theories. J. Chem. Inf. Model. 2010, 50, 1542–1548. [Google Scholar] [CrossRef] [PubMed]
  25. Rivero-Pino, F.; Millan-Linares, M.C.; Montserrat-de-la-Paz, S. Strengths and limitations of in silico tools to assess physicochemical properties, bioactivity, and bioavailability of food-derived peptides. Trends Food Sci. Technol. 2023, 138, 433–440. [Google Scholar] [CrossRef]
  26. Franke, R.; Gruska, A. General Introduction to QSAR. In Quantitative Structure-Activity Relationship (QSAR) Models of Mutagens and Carcinogens; Benigni, R., Ed.; CRC Press: Boca Raton, FL, USA; London, UK; New York, NY, USA; Washington, DC, USA, 2003; pp. 1–40. [Google Scholar]
  27. Sivaraman, G.; Jackson, N.E.; Sanchez, B.; Vázquez-Mayagoitia, Á.; Aspuru-Guzik, A.; Vishwanath, V.; de Pablo, J. A machine learning workflow for molecular analysis: Application to melting points. Mach. LearnSci. Technol. 2020, 1, 025015. [Google Scholar] [CrossRef]
  28. Scikit-Learn Integration Package for Apache Spark. Available online: https://github.com/databricks/spark-sklearn (accessed on 20 October 2023).
  29. Tropsha, A. Best Practices for QSAR Model Development, Validation, and Exploitation. Mol. Inform. 2010, 29, 476–488. [Google Scholar] [CrossRef] [PubMed]
  30. Jenkins, D.G.; Quintana-Ascencio, P.F. A solution to minimum sample size for regressions. PLoS ONE 2020, 15, e0229345. [Google Scholar] [CrossRef]
  31. Piir, G.; Kahn, I.; García-Sosa, A.T.; Sild, S.; Ahte, P.; Maran, U. Best Practices for QSAR Model Reporting: Physical and Chemical Properties, Ecotoxicity, Environmental Fate, Human Health, and Toxicokinetics Endpoints. Environ. Health Perspect. 2018, 126, 126001. [Google Scholar] [CrossRef]
  32. Cargnelutti Filho, A.; Toebe, M. Reference sample size for multiple regression in corn. Pesqui. Agropec. Bras. 2020, 55, e01400. [Google Scholar] [CrossRef]
  33. Yadav, D.K.; Khan, F.; Negi, A.S. Pharmacophore modeling, molecular docking, QSAR, and in silico ADMET studies of gallic acid derivatives for immunomodulatory activity. J. Mol. Model. 2012, 18, 2513–2525. [Google Scholar] [CrossRef]
  34. Kiralj, R.; Ferreira, M. Basic Validation Procedures for Regression Models in QSAR and QSPR Studies: Theory and Application. J. Braz. Chem. Soc. 2008, 20, 770–787. [Google Scholar] [CrossRef]
  35. Cawley, G.; Talbot, N. On Over-fitting in Model Selection and Subsequent Selection Bias in Performance Evaluation. J. Mach. Learn. Res. 2010, 11, 2079–2107. [Google Scholar]
  36. Hosseini, R.; Newlands, N.K.; Dean, C.B.; Takemura, A. Statistical modeling of soil moisture, integrating satellite remote-sensing (SAR) and ground-based data. Remote Sens. 2015, 7, 2752–2780. [Google Scholar] [CrossRef]
  37. Paneth, A.; Płonka, W.; Paneth, P. Assessment of nonnucleoside inhibitors binding to HIV-1 reverse transcriptase using HYDE scoring. Pharmaceuticals 2019, 12, 64. [Google Scholar] [CrossRef] [PubMed]
  38. Nguyen, K.T.; Blum, L.C.; van Deursen, R.; Reymond, J.-L. Classification of organic molecules by Molecular Quantum Numbers. ChemMedChem 2009, 4, 1803–1805. [Google Scholar] [CrossRef]
  39. Andrescu, T.; Mushkarov, O.; Stoyanov, L. Geometric Problems in Maxima and Minima; Birkhäuser: Boston, MA, USA; Basel, Switzerland; Berlin, Germany, 2006; pp. vii–ixi. Available online: https://www.isinj.com/mt-usamo/Geometric%20Problems%20on%20Maxima%20and%20Minima%20-%20Titu%20Andreescu,%20Oleg%20Mushkarov,%20Luchezar%20Stoyanov.pdf (accessed on 18 October 2023).
  40. Guasch, L.; Ojeda, M.J.; González-Abuín, N.; Sala, E.; Cereto-Massagué, A.; Mulero, M.; Valls, C.; Pinent, M. Identification of novel human dipeptidyl peptidase-IV inhibitors of natural origin (part I): Virtual screening and activity assays. PLoS ONE 2012, 7, e44971. [Google Scholar] [CrossRef]
  41. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Delivery Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  42. Santos, G.B.; Ganesan, A.; Emery, F.S. Oral administration of peptide-based drugs: Beyond Lipinski’s Rule. ChemMedChem 2016, 11, 2245–2251. [Google Scholar] [CrossRef]
  43. Capecchi, A.; Awale, M.; Probst, D.; Reymond, J.-L. PubChem and ChEMBL beyond Lipinski. Mol. Inf. 2019, 38, 1900016. [Google Scholar] [CrossRef]
  44. Nongonierma, A.B.; FitzGerald, R.J. Features of dipeptidyl peptidase IV (DPP-IV) inhibitory peptides from dietary proteins. J. Food Biochem. 2019, 43, e12451. [Google Scholar] [CrossRef]
  45. Soares, T.A.; Nunes-Alves, A.; Mazollari, A.; Ruggiu, F.; Wei, G.-W.; Merz, K. The (re)-evolution of Quantitative Structure–Activity Relationship (QSAR) studies propelled by the surge of machine learning methods. J. Chem. Inf. Model. 2022, 62, 5317–5320. [Google Scholar] [CrossRef] [PubMed]
  46. Tong, J.-B.; Wang, T.-H.; Feng, Y.; Jiang, G.-Y. QSAR studies of angiotensin converting enzyme inhibitors using CoMFA, CoMSIA and molecular docking: Scientific paper. J. Serb. Chem. Soc. 2021, 86, 469–482. [Google Scholar] [CrossRef]
  47. Iwaniak, A.; Hrynkiewicz, M.; Bucholska, J.; Minkiewicz, P.; Darewicz, M. Understanding the nature of bitter-taste di- and tripeptides derived from food proteins based on chemometric analysis. J. Food Biochem. 2019, 43, e12500. [Google Scholar] [CrossRef] [PubMed]
Figure 1. Procedure of QSAR methodology for di- and tripeptide DPP IV inhibitors.
Figure 1. Procedure of QSAR methodology for di- and tripeptide DPP IV inhibitors.
Applsci 13 12935 g001
Figure 2. Equations describing structure–activity relationships of dipeptide (a) and tripeptide (b) DPP IV inhibitors.
Figure 2. Equations describing structure–activity relationships of dipeptide (a) and tripeptide (b) DPP IV inhibitors.
Applsci 13 12935 g002
Figure 3. Plot of the observed and predicted logIC50 (ERM model) values for dipeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 3. CI—confidence interval.
Figure 3. Plot of the observed and predicted logIC50 (ERM model) values for dipeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 3. CI—confidence interval.
Applsci 13 12935 g003
Figure 4. Plot of the observed and predicted logIC50 (ERM model) values for tripeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 4. CI—confidence interval.
Figure 4. Plot of the observed and predicted logIC50 (ERM model) values for tripeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 4. CI—confidence interval.
Applsci 13 12935 g004
Figure 5. Backbones and structural features of peptides, summarized in Table 2: (1) peptide backbone; (2) hydrogen donors (hydrogen bond donors); (3) guanidine group; (4) nitrogen atoms; (5) methyl groups; (6) rings containing 5 atoms; (7) Csp3 bonded to 3 C atoms (tertiary carbon atoms).
Figure 5. Backbones and structural features of peptides, summarized in Table 2: (1) peptide backbone; (2) hydrogen donors (hydrogen bond donors); (3) guanidine group; (4) nitrogen atoms; (5) methyl groups; (6) rings containing 5 atoms; (7) Csp3 bonded to 3 C atoms (tertiary carbon atoms).
Applsci 13 12935 g005
Figure 6. Relationships between molecular weight and predicted tripeptide DPP IV-inhibitory activity (model based on MW). (a,b) equation and plot, being non-linear functions of the DPP IV-inhibiting activity and adequate plot, respectively.
Figure 6. Relationships between molecular weight and predicted tripeptide DPP IV-inhibitory activity (model based on MW). (a,b) equation and plot, being non-linear functions of the DPP IV-inhibiting activity and adequate plot, respectively.
Applsci 13 12935 g006
Figure 7. Plot of the observed and predicted logIC50 (MW) values for tripeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 4. CI—confidence interval.
Figure 7. Plot of the observed and predicted logIC50 (MW) values for tripeptide DPP IV-inhibitory activity. The numbers in the plot refer to the peptide sequence provided in Table 4. CI—confidence interval.
Applsci 13 12935 g007
Table 1. Statistical parameters of QSAR models obtained for di- and tripeptides with DPP IV-inhibitory activity (understood as logIC50).
Table 1. Statistical parameters of QSAR models obtained for di- and tripeptides with DPP IV-inhibitory activity (understood as logIC50).
ParameterDPP IV Inhibitors
DipeptidesTripeptides
R20.7820.829
Degrees-of-freedom adjusted R20.7550.798
CVR20.7370.710
Average error0.2550.209
RMSEC0.3430.296
RMSECV0.3770.386
F-statistics28.77326.228
p-value0.00000.0000
α0.050.05
R2—determination coefficient; CVR2—correlation coefficient of leave-one-out cross-validation; RMSEC—root mean square error of calibration; RMSECV—root mean square error of cross-validation.
Table 2. The relative weight of each descriptor contributing to the structure–activity prediction of DPP IV-inhibitory di- and tripeptides (the statistical parameters affecting the improvement of bioactivity are in bold).
Table 2. The relative weight of each descriptor contributing to the structure–activity prediction of DPP IV-inhibitory di- and tripeptides (the statistical parameters affecting the improvement of bioactivity are in bold).
DPP IV inhibitorsDipeptidesDescriptorNormalized
Coefficient
Partial-F
hydrogen donor partial surface area/total accessible surface area/MW0.755951.3355
guanidine count/MW0.508614.5636
nitrogen count2−1.000053.3368
methyl count2−0.34149.8481
ring count 5 member20.607028.3351
TripeptidesCsp3 bonded to 3 C0.00533.4655
methyl count2−0.009111.0810
ln(molecular weight)1.000020.7152
1.0/molecular weight0.365724.0262
sqrt(molecular weight)−0.641319.3041
Table 3. Observed and predicted logIC50 values obtained for dipeptide DPP IV inhibitors.
Table 3. Observed and predicted logIC50 values obtained for dipeptide DPP IV inhibitors.
No.1SequenceID 2IC50Log IC50Predicted
log IC50 3
1AA86379400.003.973.77
2AL8559882.132.952.72
3AP31777950.003.903.71
4EK85583216.753.513.00
5FL8555399.582.602.33
6FP8506363.002.562.99
7GL85612615.033.423.37
8GP31699690.003.994.12
9HL8557143.192.161.71
10HP85202820.003.453.14
11IP8501410.002.612.96
12KA31746270.003.803.53
13KP85192540.003.413.54
14LP31802370.003.382.98
15LW8688993.403.002.13
16ML883291.001.962.11
17MM883393.001.972.33
18MP3171870.002.943.02
19MW86901691.403.232.25
20NH884469.001.841.85
21PP31705860.003.774.01
22RP85182240.003.353.44
23SL85602517.083.403.16
24SP85055980.003.783.97
25TH890249.001.692.21
26TP85032370.003.383.52
27TW891384.001.922.55
28VA3172168.242.232.72
29VL892274.001.871.87
30VP3181880.002.943.07
31VR8594826.102.922.83
32WC8684420.002.622.71
33WI8679138.702.142.16
34WK867640.601.612.29
35WL867743.601.642.13
36WM8682243.102.392.30
37WN8680148.502.172.25
38WP85044530.003.663.49
39WQ8678120.302.082.15
40WR867537.801.581.57
41WS8687643.502.812.76
42WT8685482.102.682.60
43WV855665.691.822.15
44WW8686554.802.742.69
45WY8683281.002.452.40
46YP85213170.003.503.16
1 Numbers used in Figure 3. 2 ID in the BIOPEP-UWM database of bioactive peptides (see Section 2). 3 Denoted as Log IC50 (ERM).
Table 4. Observed and predicted logIC50 values obtained for tripeptide DPP IV inhibitors.
Table 4. Observed and predicted logIC50 values obtained for tripeptide DPP IV inhibitors.
No 1SequenceID 2IC50Log IC50Predicted
log IC50 3
1APG850040,0004.604.50
2GPA852240,0004.604.50
3GPM9117417.902.622.43
4GPV9116794.802.903.26
5IPA830449.001.692.03
6IPI31677.400.871.45
7IPM923369.501.841.81
8LPL8616241.402.381.45
9LPQ933982.001.912.41
10LQP86891181.103.072.41
11PPG86532252.683.353.32
12PPL8652390.142.592.42
13VPL834715.801.201.47
14WRA86086902.842.58
15WRD85963762.582.79
16WRE85953502.542.79
17WRF85994132.622.76
18WRG86004732.682.64
19WRH86066702.832.78
20WRI86107302.862.82
21WRK85984062.612.79
22WRL86129032.962.82
23WRM86076732.832.67
24WRN85974032.612.79
25WRP86117802.892.77
26WRQ86097202.862.79
27WRR86045702.762.72
28WRS86014832.682.75
29WRT86035262.722.66
30WRW86024872.692.57
31WRY86056402.812.70
32WWW86812162.332.32
33YPY8617243.702.392.73
1 Numbers used in Figure 4 and Figure 6. 2 ID in the BIOPEP-UWM database of bioactive peptides (see Section 2). 3 Denoted as Log IC50 (ERM).
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.

Share and Cite

MDPI and ACS Style

Hrynkiewicz, M.; Iwaniak, A.; Minkiewicz, P.; Darewicz, M.; Płonka, W. Analysis of Structure–Activity Relationships of Food-Derived DPP IV-Inhibitory Di- and Tripeptides Using Interpretable Descriptors. Appl. Sci. 2023, 13, 12935. https://doi.org/10.3390/app132312935

AMA Style

Hrynkiewicz M, Iwaniak A, Minkiewicz P, Darewicz M, Płonka W. Analysis of Structure–Activity Relationships of Food-Derived DPP IV-Inhibitory Di- and Tripeptides Using Interpretable Descriptors. Applied Sciences. 2023; 13(23):12935. https://doi.org/10.3390/app132312935

Chicago/Turabian Style

Hrynkiewicz, Monika, Anna Iwaniak, Piotr Minkiewicz, Małgorzata Darewicz, and Wojciech Płonka. 2023. "Analysis of Structure–Activity Relationships of Food-Derived DPP IV-Inhibitory Di- and Tripeptides Using Interpretable Descriptors" Applied Sciences 13, no. 23: 12935. https://doi.org/10.3390/app132312935

APA Style

Hrynkiewicz, M., Iwaniak, A., Minkiewicz, P., Darewicz, M., & Płonka, W. (2023). Analysis of Structure–Activity Relationships of Food-Derived DPP IV-Inhibitory Di- and Tripeptides Using Interpretable Descriptors. Applied Sciences, 13(23), 12935. https://doi.org/10.3390/app132312935

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