Int. J. Mol. Sci. 2014, 15, 14531-14545; doi:10.3390/ijms150814531 OPEN ACCESS International Journal of Molecular Sciences ISSN 1422-0067 www.mdpi.com/journal/ijms Article New Milk Protein-Derived with Potential Antimicrobial Activity: An Approach Based on Bioinformatic Studies

Bartłomiej Dziuba 1,* and Marta Dziuba 2

1 University of Warmia and Mazury in Olsztyn, Chair of Industrial and Food Microbiology, Cieszynski Square 1, Olsztyn 10-957, Poland 2 University of Warmia and Mazury in Olsztyn, Chair of Food Biochemistry, Cieszynski Square 1, Olsztyn 10-957, Poland; E-Mail: [email protected]

* Author to whom correspondence should be addressed; E-Mail: [email protected]; Tel.: +48-8-9523-3786.

Received: 6 June 2014; in revised form: 14 July 2014 / Accepted: 16 July 2014 / Published: 20 August 2014

Abstract: New peptides with potential antimicrobial activity, encrypted in milk protein sequences, were searched for with the use of bioinformatic tools. The major milk proteins were hydrolyzed in silico by 28 . The obtained peptides were characterized by the following parameters: molecular weight, isoelectric point, composition and number of residues, net charge at pH 7.0, aliphatic index, instability index, Boman index, and GRAVY index, and compared with those calculated for known 416 antimicrobial peptides including 59 antimicrobial peptides (AMPs) from milk proteins listed in the BIOPEP database. A simple analysis of physico-chemical properties and the values of biological activity indicators were insufficient to select potentially antimicrobial peptides released in silico from milk proteins by proteolytic enzymes. The final selection was made based on the results of multidimensional statistical analysis such as support vector machines (SVM), random forest (RF), artificial neural networks (ANN) and discriminant analysis (DA) available in the Collection of Anti-Microbial Peptides (CAMP database). Eleven new peptides with potential antimicrobial activity were selected from all peptides released during in silico proteolysis of milk proteins.

Keywords: milk proteins; in silico proteolysis; antimicrobial peptides; bioinformatics

Int. J. Mol. Sci. 2014, 15 14532

1. Introduction

In evolutionary terms, milk is a model example of a molecular system containing substances with exceptional ability to prevent and inhibit microbial infections [1]. The role of molecules with preventive effects should be explored to promote the development of new antimicrobial treatments, new natural food preservatives or nutraceuticals [2]. The overall antimicrobial effectiveness of milk resulting from the synergistic activity of milk peptides and proteins other than immunoglobulins, such as lactoferrin, lactoperoxidase and lysozyme, is much higher than that of individual molecules [3]. Antimicrobial milk components may demonstrate antibiotic-like activity, and they could pose a natural alternative to antibiotics [4]. Moreover, the milk proteins’ sequences contain several motifs that can be released during enzymatic hydrolysis to increase antimicrobial potential of milk proteins. In the past 20 years, numerous proteins and peptides with antimicrobial properties have been isolated from various organisms and species, ranging from bacteria to humans, or released from food proteins by proteolysis [5–8]. Proteins with low molecular mass and peptides containing up to 80 amino acids, with highly varied structure, high specificity and activity, have been most thoroughly researched. Many of these peptides demonstrating potent antimicrobial activity are low-molecular mass peptides with up to 20 amino acid residues. Intensive works are being undertaken for isolating, purifying and describing peptides for commercial applications. More than 1000 linear peptides have been identified to date [9]. Their number continues to increase because peptides can be now isolated from sources other than living organisms, for instance from food proteins [10]. The new generation of native peptides is referred to collectively as antimicrobial peptides (AMPs). Milk proteins are a natural source of bioactive peptides with diverse physiological and antimicrobial properties. Their activity is revealed after enzymatic proteolysis or after fermentation involving proteolytic microorganisms. Subject to proteolytic conditions, a diverse combination of bioactive peptides deriving from milk proteins can be obtained, but only few peptides have been identified and characterized for their antimicrobial activity [1,10–16]. Moreover, several milk protein-derived peptides demonstrate more than one type of activity, and they are referred to as multifunctional peptides [16,17]. A preliminary analysis of antimicrobial linear peptides of various origin and milk proteins revealed the presence of several motifs characterized by high structural similarity and, consequently, similar physicochemical properties and biological activity indicators. The above observation was used to formulate a research hypothesis that in addition to the identified fragments, milk protein sequences can contain several new motifs with antimicrobial activity. The objective of this study was to search for new milk protein-derived peptides with antimicrobial potential applying computer simulated proteolysis of milk proteins [18] and prediction algorithms such as SVM (support vector machines), RF (random forest), ANN (artificial neural networks) and DA (discriminant analysis) available at the interface of the CAMP database [19].

2. Results and Discussion

2.1. Physico-Chemical Characteristic of AMPs

Comparison of theoretically calculated physicochemical properties and amino acid content of 416 antimicrobial peptides listed in BIOPEP database [20] with 59 antimicrobial peptides originating

Int. J. Mol. Sci. 2014, 15 14533 from milk proteins is summarized in Tables 1 and 2. Physicochemical properties of individual AMPs from milk proteins are presented in Table S1. More than 80% of all peptides contain Lys, Gly and Leu amino acids. The Ile, Val, Ala, Arg, Ser, Phe, Asn, Thr, Gln and Pro are present in at least 50% of peptides, while Asp, Cys, Glu, His, Met, Trp and Tyr are their minor components. The content of amino acids in the sequences of the AMPs was calculated in reference to all of the examined peptides. The results indicate that certain amino acids are more common in the sequences and can influence their biological activity. Amino acids such as Arg, His, Lys, Phe, Tyr, Trp, Leu, Pro could be predominant in the sequences of biologically active peptides, depending on their type of activity [21–23].

Table 1. Amino acid content of analyzed antimicrobial peptides.

AMPs in BIOPEP AMPs from Milk Proteins Average Number of Average Number of Amino Min.-Max. Min.-Max Amino Acid Peptides Amino Acid Peptides Acid Amino Acid Amino Acid Content (%) Containing Given Content (%) Containing Given Content (%) Content (%) a b Amino Acid a b Amino Acid Ala 6.8 9.8 286 0–33.3 6.4 10.6 36 0–25.0 Arg 8.2 12.3 278 0–33.3 7.0 11.9 35 0–25.0 Asn 3.4 6.1 234 0–25.0 2.0 9.7 12 0–25.0 Asp 2.3 5.7 166 0–25.0 2.6 14.0 11 0–25.0 Cys 5.8 13.6 176 0–40.0 3.0 10.4 17 0–25.0 Gln 3.7 7.1 220 0–25.0 7.0 11.4 36 0–33.3 Glu 2.3 5.6 175 0–25.0 4.4 11.3 23 0–25.0 Gly 10.5 12.5 349 0–63.1 3.3 9.3 21 0–20.0 His 2.1 5.8 152 0–18.4 1.2 6.1 12 0–14.3 Ile 6.1 8.1 312 0–40.0 6.5 10.2 38 0–33.3 Leu 9.9 11.4 339 0–58.3 8.7 12.3 42 0–37.5 Lys 10.3 12.7 337 0–62.5 10.2 15.5 19 0–33.3 Met 1.0 3.7 113 0–7.4 1.0 3.6 16 0–7.1 Phe 4.0 6.2 258 0–37.5 2.6 7.3 21 0–16.7 Pro 5.0 9.7 215 0–53.2 8.0 11.8 40 0–28.6 Ser 5.1 7.5 282 0–23.1 4.0 9.1 26 0–21.4 Thr 3.4 6.3 226 0–33.3 5.3 10.8 29 0–33.3 Trp 2.1 5.0 178 0–38.5 3.2 8.0 24 0–16.7 Tyr 2.6 7.0 157 0–33.3 5.6 14.4 23 0–33.3 Val 5.9 8.2 300 0–37.3 8.0 12.5 38 0–37.5 a: For all antimicrobial peptides; and b: For peptides that contain given amino acid.

The highest amino acid content in all analyzed antimicrobial peptides was recorded for Lys, Gly, Arg and Leu, whereas the lowest content for Met, Trp, His, Asp, Glu and Tyr. However, evaluation of the amino acids content in the sequences where a given amino acid was found suggest that, in addition to the four predominant amino acids, Cys, Ala and Pro may also play a certain role. The majority of antimicrobial peptides are cationic amphiphatic peptides, containing one or no acidic residues and a high number of cationic Arg, Lys or His residues. Hydrophobic amino acids such as Trp, Val, Leu and Ileu account for 30%–50% of the peptide sequence and play a key role in secondary structure formation and interactions with bacterial membrane. The presence of cationic regions and the

Int. J. Mol. Sci. 2014, 15 14534 hydrophobic character of peptides influence the antibacterial activity [24,25]. Isracidin, the first antimicrobial cationic peptide derived from αs1- after chymosin cleavage, is a 23-mer peptide that contains 6 basic residues, including 5 in the N-terminal region, and 7 hydrophobic residues, with the net charge 2.2. Isracidin exhibits activity against a wide range of Gram-positive and Gram-negative bacteria [26]. The analysis of amino acid content of antimicrobial peptides revealed many fragments with the predominance of one or several amino acids [27,28]. In comparison with all AMPs listed in the BIOPEP database, amino acids such as Asn, Gly, Lys, Phe and Ser were less frequently encountered in AMPs from milk proteins. On average, the evaluated peptides had a higher content of Lys, Leu, Val and Pro. However, if the content of selected amino acids was evaluated only in sequences where they appeared, Asp, Glu, Tyr and Ala were encountered more frequently in several cases. The noted results are similar to those given by Wang and Wang [8] for AMPs in the Antimicrobial Peptide Database (APD database), but somewhat different from those reported by Hammami et al. [29]. The main difference relates to the presence of Ala residues in AMPs sequences in the Bactibase. The specific properties of the peptide structure can be attributed to its amino acid composition. The minimum, mean and maximum values of selected physicochemical properties of the analyzed antimicrobial peptides are presented in Table 2.

Table 2. The physico-chemical characteristics of peptides collected in BIOPEP and peptides from milk proteins.

AMPs Collected in BIOPEP AMPs from Milk Proteins Index Mean Min.–Max. Predominant Mean Min.–Max. Predominant Value Value Value (%) Value Value Value (%) Molecular mass (Da) 3242.9 393.5–14,350.8 2000–4000 (47) 1906.9 393.5–6707.4 393–1000 (39) pI 9.3 3.4–13.3 9–10 (29) 8.1 3.4–12.0 10–11 (25) Net charge 3.7 −7.0–20 0–5 (56) 2.2 −7.0–10.1 −2–0 (34) Instability index 31.7 −50.4–166.2 0–60 (71) 41.5 −30.9–157.7 0–20 (24) Aliphatic index 83.9 0–227.5 40–120 (65) 89.6 0.0–226.7 60–100 (49) GRAVY −0.2 −3.51–3.6 −1–0 (47) −0.4 −2.5–2.2 −1–0 (49) Boman Index (kcal/mol) 1.5 −2.6–6.8 1–2 (28) 1.3 −6.9–5.0 1–3 (61)

The molecular mass of the evaluated peptides ranged from 393 to 14,500 Da with an average of approximately 3200 Da. Molecules with a molecular mass in the 2000–4000 Da range were predominant in approximately 47% of all peptides. AMPs from milk proteins had lower molecular mass of 2000 Da on average, and they were characterized by a predominance of molecules with molecular mass under 1000 Da (approximately 39%). Out of 416 analyzed peptides, 57 had pI < pH 7 and 359 had pI > pH 7, whereas in the group of milk protein AMPs, 25 had pI < pH 7 and 34 had pI > pH 7. The above findings correlate with peptide net charges in neutral pH and the content of various amino acids discussed earlier. The net charge for all AMPs in BIOPEP ranged from −7.0 to 20.0. Thirty peptides were negatively charged, 344 had a positive charge, and the rest were neutral (Figure 1). The average net charge was 3.7, whereas the average net charge of peptides from the APD database was 4.56 [8]. In the analyzed group of milk protein AMPs, 12 were negatively charged, 35 had a positive charge and 12 were neutral. The average net charge was 2.2.

Int. J. Mol. Sci. 2014, 15 14535

Figure 1. Net charge distribution for all antimicrobial peptides from BIOPEP database.

260

240

220

200

180

160

140

120

Number of cases of Number 100

80

60

40

20

0 -10-50 5 10152025 Net charge

Many host-defense peptides that are produced naturally by the have positively charged molecules, which can be attributed to the presence of basic amino acids in their sequence. The physical forces responsible for antibacterial activity are net positive charge, hydrophobicity and flexibility [22,30,31]. In addition to cationic AMPs, significant numbers of anionic peptides from eukaryotic organisms, including milk proteins, have been reported [1,27,32]. Anionic peptides are generally abundant in Asp and Glu, they have a net charge of −1 or −2, and similarly to cationic peptides, they can form amphiphilic structures, such as α-helix or β-sheet, that are crucial for their activity or exert various effects on sensitive microorganisms [1,25,33]. The value of the instability index calculated for 271 peptides (approximately 65%) and 30 milk protein peptides (51%) was below 40, and it was indicative of peptide stability [34]. The aliphatic index, which is positively correlated with thermostability, is defined as the relative volume of Ala, Val, Ile and Leu side chains. In eight analyzed fragments, the aliphatic index was zero, which implies an absence of the above amino acids in the analyzed peptides. The value of the aliphatic index was in the range of 40–120 for 65% of all peptides and 60–100 for 49% of milk protein AMPs. The value of the GRAVY index, a measure of peptide solubility, was negative for 232 (56%) AMPs in the BIOPEP database and 30 milk protein AMPs (51%), and it was indicative of their hydrophilic nature [35]. The value of the Boman index, which is a measure of peptide affinity to proteins and its ability to establish biological interactions, ranged from −2.6 to 6.8 for all analyzed peptides and from −6.9 to 5.0 for milk protein AMPs. The predominant value of the Boman index was 1–2 for all AMPs (28%) and 1–3 for milk protein AMPs (61%). The physicochemical properties of peptides and their biological activity indicators were used to select new potentially antimicrobial peptides released in silico from milk proteins. However,

Int. J. Mol. Sci. 2014, 15 14536 the main criterion of performed prediction was the analysis made with statistical tools available in CAMP database [36].

2.2. In Silico Proteolysis of Milk Proteins

Several strategies are applied to identify and produce biologically active peptides from milk proteins [37]. At first, peptides present in the protein sequence have to be released by enzymatic hydrolysis with digestive enzymes [27,38–40], by fermenting milk with proteolytic starter cultures [41–43] or by the action of enzymes derived from microorganisms [12,44]. Proteolytic enzymes, such as pepsin, , alcalase, , and pancreatic , are most commonly applied to obtain bioactive peptides [24,45]. In the next stage, peptide fractions are separated and purified, their antimicrobial properties are determined, and they are identified usually by various MS techniques. Due to the fact that sequences of milk proteins and the cleavage sites of enzymes are known as well as the properties of antimicrobial peptides are well defined it is possible to apply bioinformatic tools in the search for new AMPs. This strategy is based on computer simulation of proteolysis and application of multivariate statistical methods to determine potentially antimicrobial motifs released from the analyzed proteins. In silico proteolysis can be performed with the use of several programs, such as PMAP (www.proteolysis.org/proteases) [46] and PeptideCutter (www.expasy.org/tools/peptidecutter/) [47], or a dedicated tool in the BIOPEP database. Detailed description of BIOPEP database has been presented in our earlier work [18]. Briefly, the “Record Operation” menu in BIOPEP contains the “ Action” option, which can be used to design the proteolytic process. In this study, 28 enzymes from the BIOPEP database were used to simulate the proteolysis of major milk proteins. The simulation produced thousands of milk protein fragments. Antimicrobial activity predictions were based solely on peptide chains containing 5 to 30 amino acids.

The results example of in silico proteolysis of αs1-casein var. gen. B by is presented in Figure 2. Only 15 out of 58 fragments were further evaluated.

Figure 2. Results window of peptides released from bovine αs1-casein var. gen. B by pancreatic elastase (EC 3.4.21.36) generated by BIOEP database.

Report of enzyme action

Results of enzyme action RPKHPI - KHQG - L - PQEV - L - NENL - L - RFFV - A - PFPQV - FG - KEKV - NEL - SKDI - G - SESTEDQA - MEDI - KEMEA - ESI - SSSEEI - V - PNSV - EQKHI - QKEDV - PSERYL - G - YL - EQL - L - RL - KKYKV - PQL - EI - V - PNSA - EERL - HSMKQG - I - HA - QQKEPMI - G - V - NQEL - A - YFYPEL - FRQFYQL - DA - YPSG - A - WYYV - PL - G - TQYTDA - PSFSDI - PNPI - G - SENSEKTTMPL - W

Location of released peptides [1-6],[7-10],[11-11],[12-15],[16-16],[17-20],[21-21],[22-25],[26-26],[27-31],[32-33],[34-37],[38-40],[41- 44],[45-45],[46-53],[54-57],[58-62],[63-65],[66-71],[72-72],[73-76],[77-81],[82-86],[87-92],[93-93],[94- 95],[96-98],[99-99],[100-101],[102-106],[107-109],[110-111],[112-112],[113-116],[117-120],[121- 126],[127-127],[128-129],[130-136],[137-137],[138-138],[139-142],[143-143],[144-149],[150-156],[157- 158],[159-162],[163-163],[164-167],[168-169],[170-170],[171-176],[177-182],[183-186],[187-187],[188- 198],[199-199]

Int. J. Mol. Sci. 2014, 15 14537

Additionally, the “search for active fragments” command was used to determine released fragments whose structure corresponded to that of peptides with known antimicrobial activity. Those fragments were not further examined.

2.3. Prediction of Antimicrobial Activity of Peptides Released During in Silico Proteolysis of Milk Proteins

The antimicrobial activity of peptides can be analyzed with the use of dedicated tools in several databases. Various search algorithms are available in APD [8,48], Bactibase [29], PhytAMP [49], AntiBP [50], AMPA [51] and CAMP [19,52]. APD relies on the physicochemical properties of peptides, and if those parameters are within the APD-defined space for natural AMPs, the program will align their sequences with those present in the database. PhytAMP and Bactibase use the Hidden Markov model for the sequence alignment. Above two databases are mainly dedicated to plant AMPs and bacteriocins, respectively. AntiBP analyses introduced sequences using three statistical methods: Artificial Neural Networks, Support Vector Machines and Quantitive Matrices. Prediction models are dedicated to peptides longer than 15 amino acids. The prediction tool available in AMPA is a highly interesting option. The AMPA algorithm uses the propensity scale to generate an antimicrobial profile. The calculation is based on antimicrobial indices defined for individual residues. The AMPA algorithm is useful in analyses of whole proteins rather than their fragments. The prediction tools available in CAMP, based on machine learning algorithms such as SVM, RF, ANN and DA, were most suitable for our research needs. According to Thomas et al. [19], the results accuracy for different statistical models ranges from 87% to 93%. The discussed tools can be used to predict the antimicrobial activity of peptides regardless of their length. In the first stage, the antimicrobial potential of milk protein fragments obtained by proteolysis was evaluated with the use of four statistical models available in the CAMP database. The results were used to manually select peptides where a positive result was reported in at least two algorithms. Overlapping peptides were excluded based on sequences of known AMPs from milk proteins. The above procedure supported the identification of 60 fragments whose sequences and properties are presented in Table S2. In the following stage, peptides whose sequences partially overlapped the sequences of the identified AMPs as well as peptides with low water solubility, as predicted by the peptide property calculator [53] (http://www.innovagen.com), were eliminated. A total of 11 potentially antimicrobial peptides with the prediction scores higher than 0.45 for at least three statistical models (Table 3) were thus selected. The sequences of potentially antimicrobial AMPs were not located nearby the sequences of the known AMPs. The peptides identified in CAMP were additionally analyzed with the use of the tools available in the APD database.

Int. J. Mol. Sci. 2014, 15 14538

Table 3. The characteristics of potential AMPs from milk proteins released during in silico proteolysis and predicted by statistical models available in CAMP database (SVM, RF, ANN and DA).

AMP Origin/ Net Isoelectric Molecular Boman Instability Aliphatic Sequence GRAVY SVM a RFC b ANN c DAC d Position Charge Point pH Mass (Da) Index Index Index

DDKHYQKA (Pancreatic elastase EC 3.4.21.36, αs2-casein, gen. var. A-11P 0.1 7.72 1004.07 4.63 53.06 12.50 −2.625 1.000 0.622 NAMP 0.507 Leukocyte elastase EC 3.4.21.37) f(74–81) GQRDLLFKDSALGFLRIP Lactoferrin 1.0 10.08 2046.41 1.52 20.89 113.89 0.028 0.659 0.4835AMP 0.834 (Prolyl oligopeptidase EC 3.4.21.26) f(294–311) ADALNLDGGYIYTAGKCGLVPVLAE Lactoferrin −2.0 3.7 2523.89 −0.22 16.36 117.20 0.536 0.568 0.5215NAMP 0.816 (V-8 EC 3.4.21.19) f(389–413) QEQNQEQP (Prolyl oligopeptidase EC 3.4.21.26), κ-casein, gen. var. A −2.0 3.79 999.9 5.3 119.70 0.00 −3.263 0.976 0.5845 AMP 0.000 Thermolysin EC 3.4.24.27) f(1–8)

KKYKVPQL (Pepsin 1.3 EC 3.4.23.1, αs1-casein, gen. var. B-8P 3.0 10.45 1003.25 1.67 46.29 85.00 −1.262 0.952 0.509 NAMP 0.916 Pancreatic elastase EC 3.4.21.71) f(102–109) AVAVVKKGSNF ( EC 3.4.212.39, Lactoferrin 2.0 10.6 1119.33 0.13 −14.91 97.27 0.591 0.943 0.483 AMP 0.912 EC 3.4.21.3) f(94–104) EMPFPK( EC 3.4.22.3, β-casein, gen. var. A2-5P 0.0 6.94 747.91 1.17 145.77 0.00 −0.983 1.000 0.571 AMP 0.190 EC 3.4.22.4) f(108–113) EPEQSL β-lactoglobulin gen. var. B −2.0 3.13 701.73 2.94 174.73 65.00 −1.517 1.000 0.4695 NAMP 0.569 (Ficain EC 3.4.22.3) f(112–117) ITRINKKIEKFQS β-casein, gen. var. A2-5P 3.0 10.83 1604.92 2.98 88.52 90.00 −0.915 0.492 0.45 AMP 0.739 (Leukocyte elastase EC 3.4.21.37) f(23–35) ITRINKKIEKF β-casein, gen. var. A2-5P 3.0 10.83 1389.71 2.71 62.85 106.36 −0.691 0.694 0.46 AMP 0.879 (Proteinase P1 () EC 3.4.21.96) f(23–33) ALFGKNGKNCPDKFCLFK Lactoferrin 2.9 9.71 2030.45 1.06 −5.79 48.89 −0.317 0.706 0.7655 AMP 0.987 (Proteinase P1 (lactocepin) EC 3.4.21.96) f(616–633) a, SVM: support vector machines; b, RF: random forest; c, ANN: artificial neural networks; and d, DA: discriminant analysis.

Int. J. Mol. Sci. 2014, 15 14539

In the group of 28 enzymes applied in in silico proteolysis, 11 of them released fragments whose antimicrobial potential was predicted in the described procedure. Four potentially antimicrobial AMPs were released from lactoferrin by chymase, metridin, prolyl oligopeptidase, proteinase P1 and V-8 protease, three from β-casein by bromelain, ficain, leukocyte elastase, proteinase P1, and one each from the remaining casein fractions by leukocyte elastase, pancreatic elastase, pepsin, prolyl oligopeptidase and thermolysin, and one from β-lactoglobulin by ficain. Low-molecular mass peptides in the range of 390–2500 Da were the predominant group of milk protein-derived peptides with antimicrobial properties (73%). The above observations were confirmed by the results of extensive research into milk proteins, in particular with regard to isracidin (αs1-casein f(7–42)) [26]. Four smaller fragments with antimicrobial activity were derived from isracidin: caseicin A and B [12], fragment 1–7 and fragment 10–14 [43]. The selected in our work, potentially antimicrobial AMPs contained 5 to 25 amino acids, which corresponds to molecular mass of 701 to 2520 Da. The general classification of AMPs include cationic peptides which can be divided in three subclasses: linear peptides forming helical structures, cysteine-rich open-ended peptides containing single or several disulfide bridges and molecules rich in specific amino acids, such as proline, , histidine and tryptophan [54]. Peptides with a positive charge constitute the majority of antimicrobial peptides known as cationic antimicrobial peptides (CAPs). In the group of known AMPs derived from milk proteins (Table S1), 35 have a positive charge ranging from 0.1 to 10.1. The second group contains anionic peptides that are generally rich in glutamic and aspartic acid [54]. Those peptides account for approximately 30% of all known milk protein AMPs. Another important feature characterizing AMPs is their ability to form amphipathic structures [24]. This process involves hydrophobic amino acids such as tryptophan and , which account for 30%–50% of the peptide sequence [22]. An analysis of selected peptides with the use of APD tools revealed that five fragments can form amphipathic helices. Fragment GQRDLLFKDSALGFLRIP (Table 3) contains three cationic and two anionic residues, two arginines, one proline, and it may have five residues on the same hydrophobic surface. Peptides ITRINKKIEKFQS and ITRINKKIEKF have four positively charged and one negatively charged residue, and they may have three residues on the same hydrophobic surface. Peptide ALFGKNGKNCPDKFCLFK has four positively charged and one negatively charged residue and, additionally, two cysteine residues that could form a disulfide bond stabilizing the α helix or the beta structure. The negatively charged ADALNLDGGYIYTAGKCGLVPVLAE peptide has three negatively charged residues (two aspartic acids and one glutamic acid), one positively charged residue, one proline and, possibly, five residues on the same hydrophobic surface. The AVAVVKKGSNF fragment is cationic due to the presence of two lysines with a high hydrophobic ratio of 54%, but only two residues have been predicted on the same hydrophobic surface. The remaining five peptides were very short, ranging from six to eight residues. According to Laverty et al. [23] and Strøm et al. [55], ultra short antimicrobial peptides consist of approximately four or five amino acid residues with a selection of amino acids that fulfill minimum functionalities required for effective antimicrobial activity. Those functionalities include charged moieties, such as arginine and the lipophilic tryptophan, which form an antimicrobial pharmacophore with the correct balance between charge and lipophilicity. In view of the above, the prediction scores for the five shortest peptides—DDKHYQKA, QNQNQEP, KKYKVPQL, EMPFPK and EPEQSL—are difficult to explain. According to the APD database, those peptides are too short to form a helix, and they do not have residues on the same hydrophobic surface.

Int. J. Mol. Sci. 2014, 15 14540

Two of them are neutral, two are anionic and one is cationic. Fragment DDKHYQKA contains three positively charged and two negatively charged residues, but due to the presence of glutamine, the overall charge is 0.1. Fragment QEQNQEQP contains proline and two glutamic acids residues typical for anionic peptides. The hydrophobic ratio for both fragments is low at 12% and 0%, respectively. Fragment EMPFPK is a neutral peptide with one positively charged and one negatively charged residue, but it also contains two prolines, and its hydrophobic ratio is 33%. Peptide EPEQSL has two negatively charged residues, it contains one proline and has a hydrophobic ratio of 16%. Fragment KKYKVPQL contains three positively charged lysines, a proline residue, and its hydrophobic ratio is 25%. However, it is worth noting that some small peptides derived from milk proteins, such as GLPQE, EQLTK, STVATL and YVL, which also do not form amphipathic structures, exhibit antimicrobial activity. Several peptides are known to impact more than one physiological function.

Small AMPs derived from bovine κ-casein YLV f(30–32), IQD f(28–30) and ovine αS2-casein LKKISQ f(165–170), PYVRYL f(203–208) also demonstrate antioxidative and ACE-inhibitory activity, respectively [56]. The mechanism of action of short peptides has not been fully explained and may be linked to a different mode of action [57]. Stability is yet another important feature of AMPs. Molecules with a stability index higher than 40 are regarded as unstable, and therefore they are characterized by lower bioavailability and shorter half-life. In the group of selected peptides (Table 3), six were predicted to be stable. Some processing techniques, such as chemical modification or incorporation of synthetic amino acids, can be applied to increase peptide stability and, consequently, lower susceptibility to hydrolysis by . In the peptide property calculator application [53] (http://www.innovagen.com), many peptides that had been predicted to be antimicrobial were determined as non-soluble due to their high hydrophobicity. Selected examples, characterized by high prediction scores in the applied statistical models, include LF f(149–175) AVAKFFSASCVPCIDRQAYPNLCQLCKG, LF f(147–153) GAVAKFF, LF f(60–65) DGGMVF, LF f(124–137) GWIIPMGILRPYLS, LF f(432–443) GYLAVAVVKKAN, LF f(666–686) VTAIANLKKCSTSPLLEACAF, LF f(126–135) IIPMGILRPY, LF f(120–132) GRSAGWIIPMGIL, LF f(517–525) LQGAVAKFFSASCVP, κ-casein, gen. var. A f(57–63) PYPYYAK. One of the drawbacks associated with the application of such AMPs is low solubility, which inhibits peptide delivery to target cells and tissues and peptide transport across membranes. The solubility problem may be addressed by gentle heating with sonification or through the application of DMSO (dimethyl sulfoxide), which is often used in pharmacology [58]. Delivery issues may be resolved through encapsulation in coordination polymers or particles and the use of micro- or nano-sized liquid marbles. Possible solutions to inhibited peptide transport across membranes include the development of soluble pro-drugs, the use of liquid dispersions and encapsulation in liposomes [44,45].

3. Experimental Section

3.1. Materials

Milk proteins as potential precursors of new antimicrobial peptides (AMPs) and their physicochemical properties were evaluated based on major milk proteins and AMPs listed in the

Int. J. Mol. Sci. 2014, 15 14541

BIOPEP database. The amino acid sequences of reference proteins for αs1-, αs2-, β-, κ-casein, β-lactoglobulin, α-lactalbumin, lactoferrin and 416 AMPs were analyzed.

3.2. Physicochemical Properties of AMPs

The physicochemical properties of antimicrobial peptides listed in the BIOPEP database were calculated with the use of software and algorithms freely available on the internet. Molecular mass, isoelectric point, the number and composition of amino acid residues, values of the instability index, aliphatic index and GRAVY index (grand average of hydropathicity) of antimicrobial peptides were computed in the ProtParam application [59]. The Boman index was calculated with APD2 algorithm: Antimicrobial Peptide Calculator and Predictor [60], and net charge was determined using Innovagen’s Peptide Property Calculator [53].

3.3. In Silico Proteolysis of Milk Proteins

The amino acid sequences of reference milk proteins were subjected to in silico proteolysis in the BIOPEP database [20]. The “Enzyme(s) action” application was used to determine the release of peptides from precursor proteins. The above supported the hydrolysis of proteins with up to three randomly-selected proteolytic enzymes out of the 28 available in the database. The option involving a single enzyme was applied for protein hydrolysis in silico in this experiment. Released peptides composed of 5 to 30 amino acid residues were submitted to statistical prediction of antimicrobial activity.

3.4. Prediction of Antimicrobial Activity of Peptides Released during in Silico Proteolysis of Milk Proteins

The antimicrobial activity of the released peptides was determined with the use of the Prediction Antimicrobial Peptides tool in the CAMP database [19,36]. Four multivariate statistical methods were used for prediction: Random Forest (RF), Support Vector Machines (SVM), Artificial Neural Network (ANN) and Discriminant Analysis (DA). The process of model development and evaluation has been described in detail by Waghu et al. [52]. Prediction results are presented with the relevant scores, excluding ANN, and peptides are classified as AMPs or non-AMPs. In this study, peptides were classified as AMPs if the resulting score was higher than 0.45 and if a positive recognition was obtained for at least three statistical methods. The selected peptides were processed in the APD database to describe their amino acid content and structure and to determine the presence of residues on the same hydrophobic surface of the molecules [48].

4. Conclusions

Research into bioactive peptides expands our knowledge of correlations between heath and diet. Increasing attention has been focused on AMPs due to their potential as a novel therapeutics. The results of the study indicate that milk proteins can be a source of new, potentially antimicrobial peptides. The bioinformatics applications that proved to be most useful for the needs of our research were the AMP prediction tools in the CAMP database that rely on machine learning algorithms such as SVM, RF, DA and ANN. Our results suggest that major milk proteins contain 60 potential AMPs that

Int. J. Mol. Sci. 2014, 15 14542 can be released by proteolytic enzymes. The highest prediction scores were reported for 11 potential AMPs in at least three statistical models. Those peptides that can be released by bromelain, chymase, ficain, leukocyte elastase, metridin, pancreatic elastase, pepsin, prolyl oligopeptidase, proteinase P1 thermolysin, V-8 protease were selected for chemical synthesis and further in vitro bioactivity analyses, which would verify the validity of the bioinformatics approach proposed in this study.

Acknowledgments

This research was supported by grant No. NN 312 484240 from the National Science Center in Poland.

Author Contributions

B.D. participated in the study design, conducted analysis and manuscript preparation; and M.D. participated in data analysis and manuscript preparation.

Conflicts of Interest

The authors declare no conflict of interest.

References

1. Floris, R.; Recio, I.; Berkhout, B.; Visser S. Antibacterial and antiviral effects of milk proteins and derivatives thereof. Curr. Pharm. Des. 2003, 9, 1257–1275. 2. Keymanesh, K.; Soltani, S.; Sardari, S. Application of antimicrobial peptides in agriculture and food industry. World J. Microbiol. Biotechnol. 2009, 25, 933–944. 3. Clare, D.A.; Swaisgood, H.E. Bioactive milk peptides: A prospectus. J. Dairy Sci. 2000, 83, 1187–1195. 4. Mils, S.; Ross, R.P.; Hill, C.; Fitzgerald, G.F.; Stanton, C. Milk intelligence: Mining milk for bioactive substances associated with human health. Int. Dairy J. 2011, 21, 377–401. 5. Dziuba, B. Antimicrobial peptides. In Bioactive Peptides and Food Proteins; Dziuba, J., Fornal, Ł., Eds.; WNT: Warsaw, Poland, 2009; pp. 35–46. 6. Giuliani, A.; Pirri, G.; Nicoletto, S.F. Antimicrobial peptides: An overview of a promising class of therapeutics. Cent. Eur. J. Biol. 2007, 2, 1–33. 7. Li-Chan, E.C.Y. Antimicrobial peptides. In Nutraceutical Proteins and Peptides in Health and Disease; Mine, Y., Shahidi, F., Eds.; Taylor & Francis Group, CRC: New York, NY, USA, 2006; pp. 99–136. 8. Wang, Z.; Wang, G. APD: The antimicrobial peptide database. Nucleic Acids Res. 2004, 32, D590–D592. 9. Hilpert, K.; Fjell, C.D.; Cherkasov, A. Short linear cationic antimicrobial peptides: Screening, optimizing, and prediction. In Peptide-Based Drug Design; Humana Press: Totowa, NJ, USA 2008; pp.127–159. 10. Pellegrini, A. Antimicrobial peptides from food proteins. Curr. Pharm. Des. 2003, 9, 1225–1238.

Int. J. Mol. Sci. 2014, 15 14543

11. Dashper, S.G.; Liu, S.W.; Reynolds, E.C. Antimicrobial peptides and their potential as oral therapeutic agents. Int. J. Pept. Res. Ther. 2007, 13, 505–516. 12. Hayes, M.; Ross, R.P.; Fitzgerald, G.F.; Hill, C.; Stanton, C. Casein-derived antimicrobial peptides generated by Lactobacillus acidophilus DPC6026. Appl. Environ. Microbiol. 2006, 72, 2260–2264. 13. López-Expósito, I.; Recio, I. Antibacterial activity of peptides and folding variants from milk proteins. Int. Dairy J. 2006, 16, 1294–1305. 14. Minervini, F.; Algaron, F.; Rizzello, C.G.; Fox, P.F.; Monnet, V.; Gobbetti, M. Angiotensin I-converting-enzyme-inhibitory and antibacterial peptides from Lactobacillus helveticus PR4 proteinase-hydrolyzed of milk from six species. Appl. Environ. Microbiol. 2003, 69, 5297–5305. 15. Recio, I.; Visser, S. Identification of two distinct antibacterial domains within the sequence of

bovine αs2-casein. Biochim. Biophys. Acta 1999, 1428, 314–326. 16. Tomita, M.; Wakabayashi, H.; Yamauchi, K.; Teraguchi, S.; Hayasawa, H. Bovine lactoferrin and lactoferricin derived from milk: Production and applications. Biochem. Cell Biol. 2002, 80, 109–112. 17. López-Expósito, I.; Quirós, A.; Amigo, L.; Recio I. Casein hydrolysates as a source of antimicrobial, antioxidant and antihypertensive peptides. Lait 2007, 87, 241–249. 18. Dziuba, M.; Dziuba B. In silico analysis of bioactive peptides. In Bioactive Proteins and Peptides as Functional Foods and Nutraceuticals; Mine, Y., Li-Chan, E.C.Y., Jiang, B., Eds.; Blackwell Publishing Ltd. and Institute of Food Technologists: London, UK, 2010; pp. 325–340. 19. Thomas, S.; Karnik, S.; Barai, R.S.; Jayaraman, V.K.; Idicula-Thomas, S. CAMP: A useful resource for research on antimicrobial peptides. Nucleic Acids Res. 2010, 38, D774–D780. 20. BIOPEP. Available online: http://www.uwm.edu.pl/biochemia/index.php/pl/biopep (accessed from 2 to 31 January 2014). 21. Epand, R.M.; Vogel, H.J. Diversity of antimicrobial peptides and their mechanisms of action. BBA Biomembr. 1999, 1462, 11–28. 22. Laverty, G.; Gorman, S.P.; Gilmore, B.F. The Potential of antimicrobial peptides as biocides. Int. J. Mol. Sci. 2011, 12, 6566–6596. 23. Zamiatnin, A.A.; Poronina, O.L. Computer biochemistry and molecular physiology of the natural peptide ligands: different functional groups. In Peptide Science—Present and Future. Proceedings of the 1st International Peptide Symposium; Shimonishi, Y., Ed.; Kluwer: Boston, MA, USA, 1999; pp. 212–214. 24. Akalin, A.S. Dairy-derived antimicrobial peptides: Action mechanisms, pharmaceutical uses and production proposals. Trends Food Sci. Technol. 2014, 36, 79–95. 25. Clare, D.A.; Catignani, G.L.; Swaisgood, H.E. Biodefense properties of milk: The role of antimicrobial proteins and peptides. Curr. Pharm. Des. 2003, 9, 1239–1255. 26. Lahov, E.; Regelson, W. Antibacterial and immunostimulating casein-derived substances from milk: Casecidin, isracidin peptides. Food Chem. Toxicol. 1996, 34, 131–145. 27. Benkerroum, N. Antimicrobial peptides generated from milk proteins: A survey and prospects for application in the industry—A review. Int. J. Dairy Technol. 2010, 63, 320–338. 28. McPhee, J.B.; Hancock, R.W. Function and therapeutic potential of host defense peptides. J. Peptide Sci. 2005, 11, 677–687.

Int. J. Mol. Sci. 2014, 15 14544

29. Hammami, R.; Zouhir, A.; Ben Hamida, J.; Fliss, I. BACTIBASE: A new web-accessible database for bacteriocin characterization. BMC Microbiol. 2007, 7, 89. 30. Steinstraesser, L.; Kraneburg, U.; Jacobson, F.; Al-Benna, S. Host defense peptides and their antimicrobial-immunomodulatory duality. Immunobiology 2011, 216, 322–333. 31. Zasloff, M. Antimicrobial peptides of multicellular organisms. Nature 2002, 415, 389–395. 32. Harris, F.; Dennison, S.R.; Phoenix, D.A. Anionic antimicrobial peptides from eukaryotic organisms. Curr. Protein Pept. Sci. 2009, 10, 585–606. 33. Wakabayashi, H.; Takase, M.; Tomita, M. Lactoferricin derived from milk protein lactoferrin. Curr. Pharm. Des. 2003, 9, 1277–1287. 34. Guruprasad, K.; Reddy, B.V.B.; Pandit, M.W. Correlation between stability of a protein and its dipeptide composition—A novel-approach for predicting in vivo stability of a protein from its primary sequence. Protein Eng. 1990, 4, 155–161. 35. Kyte, J.; Doolittle, R.F. A simple method for displaying the hydropathic character of a protein. J. Mol. Biol. 1982, 157, 105–132. 36. CAMP. Available online: http://www.camp.bicnirrh.res.in/ (accessed from 10 February to 30 April 2014). 37. Korhonen, H.; Pihlanto, A. Bioactive peptides: Production and functionality. Int. Dairy J. 2006, 16, 945–960. 38. López-Expósito, I.; Minervini, F.; Amigo, L.; Recio, I. Identification of antibacterial peptides from bovine kappa-casein. J. Food Prot. 2006, 69, 2992–2997. 39. McCann, K.B.; Shiell, B.J.; Michalski, W.P.; Lee, A.; Wan, J.; Roginski, H.; Coventry, M.J.

Isolation and characterisation of a novel antibacterial peptide from bovine αS1-casein. Int. Dairy J. 2006, 16, 316–323. 40. Pellegrini, A.; Dettling, C.; Thomas, U.; Hunziker, P. Isolation and identification of four bactericidal domains in the bovine β-lactoglobulin. Biochim. Biophys. Acta 2001, 1526, 131–140. 41. Lignitto, L.; Segato, S.; Balzan, S.; Cavatorta, V.; Oulahal, N.; Sforza, S.; Degraeve, P.; Galaverna, P.; Novelli, E. Preliminary investigation on the presence of peptides inhibiting the growth of Listeria innocua and Listeria monocytogenes in Asiago d’Allevo cheese. Dairy Sci. Technol. 2012, 92, 297–308. 42. Meira, S.M.M.; Daroit, D.J.; Helfer, V.E.; Corrêa, A.P.F.; Jéferson, S.; Carro, S.; Brandelli, A. Bioactive peptides in water-soluble extracts of ovine cheeses from Southern Brazil and Uruguay. Food Res. Int. 2012, 48, 322–329. 43. Rizzello, C.G.; Gobbetti, L.M.; Carbonara, T.; de Bari, M.D.; Zambonin, P.G. Antibacterial activities of peptides from the water-soluble extracts of Italian cheese varieties. J. Dairy Sci. 2005, 88, 2348–2360. 44. Agyei, D.; Danquah, M.K. Rethinking food-derived bioactive peptides for antimicrobial and immunomodulatory activities. Trends Food Sci. Technol. 2012, 23, 62–69. 45. Agyei, D.; Danquah, M.K. Industrial-scale manufacturing of pharmaceutical-grade bioactive peptides. Biotechnol. Adv. 2011, 29, 272–277. 46. PMAP. Available online: http://www.proteolysis.org/proteases (accessed on 7 January 2014). 47. PeptideCutter. Available online: http://www.expasy.org/tools/peptidecutter/ (accessed on 7 January 2014).

Int. J. Mol. Sci. 2014, 15 14545

48. Wang, G.S.; Li, X.; Wang, Z. APD2: The updated antimicrobial peptide database and its application in peptide design. Nucleic Acids Res. 2009, 37, D933–D937. 49. Hammami, R.; Ben Hamida, J.; Vergoten, G.; Fliss, I. PhytAMP: A database dedicated to antimicrobial plant peptides. Nucleic Acids Res. 2009, 37, D963–D968. 50. Lata, S.; Sharma, BK.; Raghava, G.P.S. Analysis and prediction of antibacterial peptides. BMC Bioinform. 2007, 8, 263. 51. Torrent, M.; di Tommaso, P.; Pulido, D.; Nogués, M.V.; Notredame, C.; Boix, E.; Andreu, D. AMPA: An automated web server for prediction of protein antimicrobial regions. Bioinformatics 2012, 28, 130–131. 52. Waghu, F.H.; Gopi, L.; Barai, R.S.; Ramteke, P.; Nizami, B.; Idicula-Thomas, S. CAMP: Collection of sequences and structures of antimicrobial peptides. Nucleic Acids Res. 2014, 42, D1154–D1158. 53. Peptide Property Calculator. Available online: http://www.innovagen.se/custom-peptide-synthesis/ peptide-property-calculator/peptide-property-calculator.asp (accessed from 7 January to 28 April 2014). 54. Marshall, S.H.; Arenas, G. Antimicrobial peptides: A natural alternative to chemical antibiotics and a potential for applied biotechnology. Elektron. J. Biotechnol. 2003, 6, 271–282. 55. Strøm, M.B.; Haug, B.E.; Skar, M.L.; Stensen, W.; Stiberg, T.; Svendsen, J.S. The pharmacophore of short cationic antibacterial peptides. J. Med. Chem. 2003, 46, 1567–1570. 56. López-Expósito, I.; Gómez-Ruiz, J.A.; Amigo, L.; Recio, I. Identification of antibacterial peptides

from ovine αS2-casein. Int. Dairy J. 2006, 16, 1072–1080. 57. Meng, S.; Xu, H.; Wang, F. Research advances of antimicrobial peptides and applications in food industry and agriculture. Curr. Protein Pept. Sci. 2010, 11, 264–273. 58. Vignes, R. Dimethyl Sulfoxide (DMSO)—A “new” clean, unique, superior solvent. In American Chemical Society Annual Meeting; Washington, DC, USA, 2000; pp. 1–20. 59. ProtParam. Available online: http://web.expasy.org/protparam/ (accessed from 7 January to 28 February 2014). 60. APD2: Antimicrobial Peptide Calculator and Predictor. Available online: http://aps.unmc.edu/AP/ prediction/prediction_main.php (accessed on 30 April 2014).

© 2014 by the authors; licensee MDPI, Basel, Switzerland. This article is an open access article distributed under the terms and conditions of the Creative Commons Attribution license (http://creativecommons.org/licenses/by/3.0/).

Supplementary Information

Table S1. The characteristics of known AMPs derived from milk proteins. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian 1 h RPKHPIK f(1–7) of αS1-casein 3.1 11.17 875.0 −2.029 19.84 55.71 Yeast 2 min 3.68 E.coli 2 min Mamalian 1 h RPKHPIKHQGLPQEVL Isracidin f(1–23) of α -casein 2.2 9.99 2764.2 −0.987 46.49 97.39 Yeast 2 min 2.42 NENLLRF S1 E.coli 2 min Mamalian 10 h GLPQE f(10–14) of αS1-casein −1.0 4.00 542.5 −1.040 123.56 78.00 Yeast >20 h 1.29 E.coli >20 h Mamalian 20 h IKHQGLPQE f(6–14) of αS1-casein 0.1 6.75 1049.1 −1.256 63.66 86.67 Yeast 30 min 1.92 E.coli >10 h Mamalian 100 h VLNENLLR f(15–22) of αS1-casein 0.0 5.97 970.1 0.075 32.83 182.5 Yeast >20 h 2.02 E.coli >10 h Mamalian 5.5 h LRLKKYKVPQL αS1-casein, gen. var. B-8P f(99–109) 4 10.46 1385.7 −0.636 46.18 132.73 Yeast 3 min 1.67 E.coli 2 min Mamalian 1.9 h SDIPNPIGSENSEK f(180–193) of αS1-casein −2 4.14 1486.5 −1.314 32.69 55.71 Yeast >20 h 2.89 E.coli >10 h KTKLTEEEKNRLNFLK Mamalian 1.3 h casocidin-I f(150–188) of α -casein gen. KISQRYQKFALPQYLK S2 7.1 10.08 4869.6 −1.323 27.85 70.00 Yeast 3 min 2.83 var. A TVYQHQK E.coli 2 min Mamalian 100 h VYQHQKAMKPWIQPK f(183–206) α -casein gen. var. A 5.1 10.12 3002.6 −0.783 2.25 72.92 Yeast >20 h 1.37 TKVIPYVRY S2 E.coli >10 h Mamalian 5.5 h LKKISQ f(165–170) of ovine αS2-casein 2.0 10.00 715.8 −0.633 12.13 130.00 Yeast 3 min 1.7 E.coli 2 min Mamalian 100 h VYQHQKAMKPWIQPK f(183–207) of bovine α -casein 5.1 10.12 3115.7 −0.600 2.56 85.60 Yeast >20 h 1.12 TKVIPYVRYL S2 E.coli >10 h

S2

Table S1. Cont. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian 5.5 h LKKISQRYQKFALPQY f(164–179) of bovine αS2-casein 4.0 10.17 2011.4 −0.938 25.45 79.38 Yeast 3 min 2.02 E.coli 2 min Mamalian >20 h PYVRYL f(203–208) of ovine αS2-casein 1.0 9.00 809.9 −0.117 −4.23 113.33 Yeast >20 h 1.04 E.coli ? Mamalian 1.3 h KTVYQHQKAMKPWIQ f(181–207) of ovine α -casein 6.1 10.22 3345.0 −0.726 3.11 79.26 Yeast 3 min 1.33 PKTKVIPYVRYL S2 E.coli 3 min Mamalian 5.5 LKTVYQHQKAMKPWI f(180–207) of ovine α -casein 6.1 10.22 3458.22 −0.564 0.32 90.36 Yeast 3min 1.11 QPKTKVIPYVRYL S2 E.coli 3min ALPQYLKTVYQHQKA Mamalian 4.4 h MKPWIQPKTKVIPYVR f(175–207) of ovine αS2-casein 6.1 10.1 4030.8 −0.503 11.18 91.52 Yeast >20 min 0.91 YL E.coli >10 min QKFALPQYLKTVYQHQ Mamalian 0.8 h KAMKPWIQPKTKVIPY f(172–207) of ovine αS2-casein 7.1 10.17 4434.3 −0.589 11.08 83.89 Yeast 10 min 1.06 VRYL E.coli 10 h LKKISQRYQKFALPQY Mamalian 5.5 h LKTVYQHQKAMKPWI f(164–207) of ovine αS2-casein 10.1 10.35 5451.6 −0.700 9.69 86.36 Yeast 2 min 1.44 QPKTKVIPYVRYL E.coli 2 min Mamalian 5.5 h LKKISQRYQKFALPQ f(164–180) of ovine α -casein 4.0 10.17 2124.5 −0.659 24.54 97.65 Yeast 2 min 1.61 YL S2 E.coli 2 min Mamalian 2.8 h YQEPVLGPVRGPFPI β-casein, gen. var. A2-5P f(193–207) 0 6.00 1668.9 −0.033 80.71 90.67 Yeast 10 min 0.3 E.coli 2 min Mamalian 2.8 h YQEPVLGPVRGPFPIIV β-casein, genetic variant A2-5P f(193–209) 0 6.00 1881.2 0.482 67.4 120.00 Yeast 10 min −0.25 E.coli 2 min MAIPPKKNQDKTEIPTI Mamalian 30 h NTIASGEPTSTPTTEAV Kappacin f(107–169) of κ-casein gen. Yeast 20 h −7.0 4.04 6707.4 −0.370 64.17 77.66 1.55 ESTVATLEDSPEVIESPP var. A E.coli 10 h EINTVQVTSTAV

S3

Table S1. Cont. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian 2.8 h YYQQKPVA f(42–49) of bovine κ-casein 1.0 8.50 996.1 −1.137 62.9 48.75 Yeast 10 min 1.38 E.coli 2 min Mamalian 100 h VQVTSTAV f(162–169) of bovine κ-casein 0.0 5.49 803.9 1.087 −11.29 121.25 Yeast >20 h 0.01 E.coli >10 h Mamalian 1.9 h STVATL f(141–146) of bovine κ-casein 0.0 5.24 590.6 1.267 8.33 130.00 Yeast >20 h −0.37 E.coli >10 h Mamalian 1.1 h FSDKIAK f(18–24) of bovine κ-casein 1.0 8.59 807.9 −0.429 −15.69 70.00 Yeast 3 min 1.93 E.coli 2 min Mamalian 20 h IQY f(28–30) of bovine κ-casein 0.0 6.09 422.79 −0.100 −18.47 130.00 Yeast >30 min 0.25 E.coli >10 h Mamalian 20 h YVL f(30–32) of bovine κ-casein 0.0 6.09 393.51 2.233 6.67 226.67 Yeast 30 min −2.59 E.coli >10 h Mamalian 1 h EIPT f(118–121) of bovine κ-casein −1.0 3.85 458.54 −0.325 3.85 97.50 Yeast 30 min 1.11 E.coli >10 h Mamalian 100 h VESTVATL f(139–146) of bovine κ-casein −1.0 4.00 818.9 1.038 32.83 133.75 Yeast >20 h 0.08 E.coli >10 h Mamalian >20 h PAAVRSPAQILQ f(64–75) of bovine κ-casein 1.0 10.18 1250.4 0.200 157.72 114.17 Yeast >20 h 0.84 E.coli ? Mamalian 4.4 h AVESTVATLEDSPEVIE f(136–156) of bovine κ-casein −6.0 3.40 2199.3 −0.243 134.70 88.10 Yeast >20 h 1.54 SPPE E.coli >10 h Mamalian 1h LTD1-f(1–5) of α-laktoalbumin, gen. EQLTK 0.0 6.10 617.7 −1.560 46.52 78.00 Yeast 30 min 3.11 var. B E.coli >10 h

S4

Table S1. Cont. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian >30 h LDT 2, α-chain f(17–31) of GYGGVSLPEWVCTTF −1.0 4.00 1615.8 0.453 44.95 64.67 Yeast >20 h −0.46 α- laktoalbumin, gen. var. B E.coli >10 h Mamalian 4.4 h LDT 2, β-chain f(109–114) of ALCSEK 0.0 6.04 649.7 −0.017 40.43 81.67 Yeast >20 h 1.29 α-laktoalbumin, . gen. varB E.coli >10 h Mamalian 1.2 h LDC, α-chain f(61–68) of CKDDQNPH −1.0 5.21 956.0 −2.525 5.15 0.00 Yeast >20 h 4.81 α-laktoalbumin, gen. var. B E.coli >10 h Mamalian 20 h LDC, β-chain chain f(75–80) of ISCDKF 0.0 5.83 711.8 0.267 80.62 65.00 Yeast 30 min 1.41 α-laktoalbumin, gen. var. B E.coli >10 h Mamalian 100 h VAGTWY f(15–20) β-laktoglobulin, gen. var. A 0.0 5.49 695.7 0.450 −30.87 65.00 Yeast >20 h −1.06 E.coli >10 h Mamalian 20 h IPAVFK f(78–83) β-laktoglobulin, gen. var. A 1.0 8.75 673.8 1.300 10.58 130.00 Yeast 30 min −1.36 E.coli >10 h Mamalian 4.4 h AASDISLLDAQSAPLR f(25–40) β-laktoglobulin, gen. var. A −1.0 4.21 1627.8 0.256 60.91 122.50 Yeast >20 h 1.32 E.coli >10 h Mamalian 100 h VLVLDTDYK f(92–100) β-laktoglobulin, gen.var. A −1.0 4.21 1065.2 0.344 −7.81 151.11 Yeast >20 h 0.86 E.coli >10 h Mamalian 1.1 h FKCRRWQWRMKKLG bovine lactoferricin B f(17–41) 7.9 11.84 3125.8 −0.576 77.92 50.80 Yeast 3 min 2.75 APSITCVRRAF of lactoferrin E.coli 2 min Mamalian 4.4 h APRKNVRW f(1–8) of lactoferrin 3.0 12.01 1026.2 −1.613 95.00 48.75 Yeast >20 h 4.23 E.coli >10 h Mamalian 1.1 h FKCRRWQWRMKKLG f(17–42) of lactoferrin 7.9 11.84 3196.9 −0.485 75.31 52.69 Yeast 3 min 2.57 APSITCVRRAFA E.coli 2 min

S5

Table S1. Cont. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian 1.1 h FKCRRWQWRMKKLGAPSITCVR f(17–43) of lactoferrin 7.9 11.84 3310.0 −0.326 72.89 65.19 Yeast 3 min 2.29 RAFAL E.coli 2 min Mamalian 4.4 h APRKNVRWCTISQPEW f(1–16) of lactoferrin 2.0 9.51 1971.2 −1.056 98.42 48.75 Yeast >20 h 2.72 E.coli >10 h Mamalian 1.2 h CIRA f(45–48) of lactoferrin 1.0 8.55 461.61 1.075 7.50 122.5 Yeast >20 h 3.69 E.coli >10 h Mamalian 1.1 h FKCRRWQWRMKKLGAPSITCVR f(17–47) of lactoferrin 7.9 11.37 3811.6 −0.316 78.95 69.35 Yeast 3 min 2.5 RAFALECIR E.coli 2 min Mamalian 4.4 h APRKNVRWCTI f(1–11) of lactoferrin 3.0 10.86 1343.6 −0.600 101.45 70.91 Yeast >20 h 2.74 E.coli >10 h Mamalian 1.2 h CRRWQWRMKKLGAPSITCV f(19–37) of lactofferin 4.9 10.93 2319.8 −0.468 77.14 61.58 Yeast >20 h 2.16 E.coli >10 h Mamalian 1.1 h FKCRRWQWRMKKLG f(17–30) of lactofferin 6.0 11.74 1923.3 −1.421 74.31 27.86 Yeast 3 min 3.55 E.coli 2 min Mamalian 5.5 h LECIRA f(43–48) of lactofferin 0.0 5.99 703.8 0.767 81.57 146.67 Yeast 3 min 1.46 E.coli 2 min Mamalian 4.4 h APRKNVRWCTISQPEWFKCRRW f(1–42) of lactofferin 9.9 11.63 5150.1 −0.702 84.35 51.19 Yeast >20 h 2.63 QWRMKKLGAPSITCVRRAFA E.coli >10 h Mamalian 1.1 h FKCRRWQWRMKKLGAPSITCVR f(17–48) of lactofferin 7.9 11.37 3883.7 −0.250 76.80 70.31 Yeast 3 min 2.36 RAFALECIRA E.coli 2 min Mamalian 2.8 h Lactoferrampin f(268–284) WKLLSKAQEKFGKNKSR 5.0 10.58 2048.4 −1.482 36.01 51.73 Yeast 3 min 2.97 of lactofferin E.coli 2 min

S6

Table S1. Cont. Boman Net Isoelectric Molecular Instability Aliphatic Sequence AMP Origin/Position GRAVY Half Life Index Charge Point (pH) Mass (Da) Index Index (kcal/mol) Mamalian 4.4 h Antiviral f(241–249) of ADRDQYELL −2.0 4.03 1122.2 −1.156 −25.68 97.78 Yeast >20 h 3.68 lactofferin E.coli >10 h TRVVWCAVG Lactoferrin (fr. 343–351) 1 9.01 990.1 1.156 −0.54 107.78 Mamalian 7.2 h −0.11 GRDPYKLRPV Lactoferrin (fr. 68–77) 2 9.99 1200.4 −1.33 58.29 68.00 Yeast >20 h 3.43 Mamalian 1 h Antiviral f(263–268) of EDLIWK −1.0 4.37 802.9 −0.583 40.43 130.00 Yeast 30 min 1.48 lactofferin E.coli >10 h

S7

Table S2. The characteristics of potential AMPs from milk proteins released during in silico proteolysis according to positive prediction scores of two or more statistical models available in CAMP database (SVM, RF, ANN and DA).

Boman AMP Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) Index GRAVY SVM RFC ANN DAC Origin/Position Charge Point pH Mass Da Index Index kcal/mol

DDKHYQKA (Pancreatic elastase EC 3.4.21.36, αs2-casein, gen. var. 0.1 7.72 1004.07 4.63 53.06 12.50 −2.625 1.000 0.622 NAMP 0.507 Leukocyte elastase EC 3.4.21.37) A-11P f(74–81)

αs2-casein, gen. var. FTKKTKL (Pancreatic elastase EC 3.4.21.36) 3.0 10.84 865.08 1.98 −3.56 55.71 −0.929 0.998 0.5215 AMP 0.940 A-11P f(147–153) DTDYKKYL (Pancreatic elastase EC 3.4.21.36, β-lactoglobulin, gen. Pepsin 1.3 EC 3.4.23.1, 0.0 6.75 1045.16 3.3 −10.04 48.75 −1.788 1.000 0.4615 NAMP 0.556 var. B f(96–103) Pancreatic elastase EC 3.4.21.71) KNTMEHV (Pancreatic elastase EC 3.4.21.36, αs2-casein, gen. var. Leukocyte elastase EC 3.4.21.37, Proteinase P1 0.1 6.75 857.9 2.83 −14.33 41.43 −1.243 0.928 0.609 NAMP 0.029 A-11P f(1–7) (lactocepin) EC 3.4.21.96)

αs2-casein, gen. var. NEEEYSI (Pancreatic elastase EC 3.4.21.36) −3 3.67 882.8 3.66 101.71 55.71 −1.657 1.000 0.5665 NAMP 0.152 A-11P f(48–54) α-lactalbumin, gen. QNNDSTEYG (Pancreatic elastase EC 3.4.21.36) −2 3,67 1026.9 4.39 42.26 0.00 −2.300 1.000 0.5875 NAMP 0.027 var. B f(43–51) TKESPQTHYYA Lactoferrin f(84–94) 0.1 6.41 1324.4 2.68 156.75 9.09 −1.700 0.999 0.524 NAMP 0.027 (Pancreatic elastase EC 3.4.21.36) Lactoferrin PNSKEKYYG (Pancreatic elastase EC 3.4.21.36) 1.0 9.47 1085.18 3.03 13.17 0.00 −2.244 1.000 0.252 NAMP 0.873 f(517–525) LQGAVAKFFSASCVP (Prolyl oligopeptidase EC Lactoferrin 1.0 8.97 1524.81 −0.58 54.95 84.67 0.980 0.7850.2385 AMP 0.716 3.4.21.26) Poor water solubility f(145–159) GQRDLLFKDSALGFLRIP Lactoferrin 1.0 10.08 2046.41 1.52 20.89 113.89 0.028 0.659 0.4835 AMP 0.834 (Prolyl oligopeptidase EC 3.4.21.26) f(294–311) QEQNQEQP (Prolyl oligopeptidase EC 3.4.21.26), κ-casein, gen. var. A −2 3.79 999.9 5.3 119.70 0.00 −3.263 0.976 0.5845 AMP 0.000 Thermolysin EC 3.4.24.27) f(1–8)

S8

Table S2. Cont.

Boman AMP Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) Index GRAVY SVM RFC ANN DAC Origin/Position Charge Point pH Mass Da Index Index kcal/mol ADALNLDGGYIYTAGKCGLVPVLAE Lactoferrin f(389– −2.0 3.7 2523.89 −0.22 16.36 117.20 0.536 0.568 0.5215 NAMP 0.816 (V-8 proteaza EC 3.4.21.19) 413) GYLAVAVVKKAN Lactoferrin 2.0 10.18 1232.49 −0.46 −0.31 130.0 0.733 0.8930.4395 AMP 0.928 (Glutamyl II EC 3.4.21.82) f(432–443)

αS1-casein, gen. var. KKYKVPQL (Pepsin 1.3 EC 3.4.23.1, B-8P 3.0 10.45 1003.25 1.67 46.29 85.00 −1.262 0.952 0.509 NAMP 0.916 Pancreatic elastase EC 3.4.21.71) f(102–109) NVPGEIVESL (Pepsin 1.3 EC 3.4.23.1, β-casein gen. var. B G EC 3.4.21.20, Metridin EC 3.4.21.3, −2 3.79 1056.1 0.47 50.75 136.00 0.340 0.586 0.516 NAMP 0.001 f(7–17) Pancreatic elastase II EC 3.4.21.71) GRSAGWIIPMGIL (Pepsin 1.3 EC 3.4.23.1) Lactoferrin 1.0 11.04 1370.68 −0.82 52.79 127.69 0.923 0.1710.564 AMP 0.645 Poor water solubility f(120–132)

αs2-casein, gen. var. FTKKTK (Thermolysin EC 3.4.24.27) 3.0 10.84 751.92 3.13 8.33 0.00 −1.717 1.000 0.5285 AMP 0.986 A-11P f(147–152)

αs2-casein, gen. var. ANEEEYS (Thermolysin EC 3.4.24.27) −3 3.67 840.8 4.11 101.71 14.29 −2.043 0.933 0.6535 NAMP 0.002 A-11P (fr 47–53) Lactoferrin AGDDQG (Thermolysin EC 3.4.24.27) −2 3.56 561.5 3.21 8.33 16.67 −1.583 0.624 0.551 NAMP 0.000 f(506–511)

ESTEVF ( EC 3.4.21.20, αs2-casein, gen. var. −2 3.80 710.7 2.09 72.53 48.33 −0.250 1.000 0.531 NAMP 0.005 (Pancreatic elastase II EC 3.4.21.71), A-8P f(152–147) PYPYYAK (Cathepsin G EC 3.4.21.20) κ-casein, gen. var. A 1.0 9.39 901.03 0.59 77.66 14.29 −1.314 0.999 0.3085 NAMP 0.719 Poor water solubility f(57–63)

S9

Table S2. Cont.

Boman AMP Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) Index GRAVY SVM RFC ANN DAC Origin/Position Charge Point pH Mass Da Index Index kcal/mol

TDAPSF (Chymase EC 3.4.21.39), αS1-casein, gen. var. −1 3.80 636.6 1.65 72.53 16.67 −0.333 0.901 0.5085 NAMP 0.011 Metridin EC 3.4.21.3) B-8P f(174–179) AVAVVKKGSNF (Chymase EC 3.4.212.39, Lactoferrin 2.0 10.6 1119.33 0.13 −14.91 97.27 0.591 0.9430.483 AMP 0.912 Metridin EC 3.4.21.3) f(94–104) IIPMGILRPY (Chymase EC 3.4.21.39) Lactoferrin 1.0 9.84 1172.5 −0.79 47.87 156.0 0.980 0.3710.595 AMP 0.837 Poor water solubility f(126–135) VTAIANLKKCSTSPLLEACAF Lactoferrin 0.9 8.29 2180.62 0.00 43.29 107.14 0.748 0.685 0.917 AMP 0.852 (Chymase EC 3.4.21.39) Poor water solubility f(666–686) PQRDMPIQAF (Chymase EC 3.4.21.39, β-casein gen. var. B 0 6.27 1202.3 2.26 72.20 49.00 −0.720 0.715 0.697 NAMP 0.018 Metridin EC 3.4.21.3) f(181–190) EAGRDPY (Chymase EC 3.4.21.39), Lactoferrin −1 4.37 806.8 3.97 8.57 14.29 −1.857 0.848 0.5785 NAMP 0.000 Metridin EC 3.4.21.3) f(66–72) GTKESPQTHY (Chymase EC 3.4.21.39), Lactoferrin 0.1 6.75 1147.2 3.03 126.91 0.00 −1.96 0.718 0.625 NAMP 0.007 Metridin EC 3.4.21.3) f(83–92)

TTMPLW (Papain 3.4.22.2, Bromelain EC 3.4.22.4) αS1-casein, gen. var. 0 5.19 747.9 −0.74 121.03 65.0 0.300 0.976 0.5275 NAMP 0.046 Poor water solubility B-8P f(204–209) β-casein gen. var. B DMPIQA (Papain 3.4.22.2, Bromelain EC 3.4.22.4) −1 3.80 673.7 0.86 81.57 81.67 −0.067 0.768 0.525 NAMP 0.005 f(184–189) QPTTMA (Papain 3.4.22.2, Ficain EC 3.4.22.3, κ-casein, gen. var. A Leukocyte elastase EC 3.4.21.37, 0 5.52 647.7 1.08 61.00 16.67 −0.467 1.000 0.6035 NAMP 0.015 f(91–96) Bromelain EC 3.4.22.4) Poor water solubility

S10

Table S2. Cont.

Boman AMP Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) Index GRAVY SVM RFC ANN DAC Origin/Position Charge Point pH Mass Da Index Index kcal/mol

αS1-casein, gen. var. EERLHSMK (Bromelain EC 3.4.22.4) 0.1 6.86 1029.1 4.35 49.50 48.75 −1.712 0.717 0.5785 NAMP 0.001 B-8P f(117–124)

SESTEDQA (Ficain EC 3.4.22.3, αS1-casein, gen. var. −3 3.57 865.8 4.43 105.05 12.50 −1.812 0.998 0.637 NAMP 0.004 Bromelain EC 3.4.22.4) B-8P f(46–53) PQRDMPIQA (Ficain EC 3.4.22.3), β-casein, gen. var. 0 6.27 1055.2 2.84 79.11 54.44 −1.111 0.827 0.6745 NAMP 0.010 Bromelain EC 3.4.22.4) A-5P f(181–189) EMPFPK (Ficain EC 3.4.22.3, β-casein, gen. var. 0 6.94 747.91 1.17 145.77 0.00 −0.983 1.000 0.571 AMP 0.190 Bromelain EC 3.4.22.4) A2-5P f(108–113) β-lactoglobulin, gen. EPEQSL (Ficain EC 3.4.22.3) −2.0 3.13 701.73 2.94 174.73 65.00 −1.517 1.000 0.4695 NAMP 0.569 var. B f(112–117) ESPQTHY (Ficain EC 3.4.22.3, Lactoferrin f(86–92) −0.9 5.24 860.8 3.3 189.14 0.00 −2.086 0.551 0.658 NAMP 0.004 Bromelain EC 3.4.22.4)

PQLEIV (Leukocyte elastase EC 3.4.21.37, αS1-casein, gen. var. −1 4.00 697.8 −0.25 58.38 178.33 0.650 0.945 0.553 NAMP 0.020 Prolidase L.lactis s.cremoris H61 EC 3.4.13.9) B-8P f(107–112)

EKTTMPLW (Leukocyte elastase EC 3.4.21.37), αS1-casein, gen. var. 0 6.10 1005.2 0.98 93.28 48.75 −0.700 1.000 0.644 NAMP 0.021 Proteinase P1 (lactocepin EC 3.4.21.37) B-8P f(192–199)

αs2-casein, gen. var. NEEEYS (Leukocyte elastase EC 3.4.21.37) −3 3.67 769.7 5.1 117.0 0.00 −2.683 1.000 0.6795 NAMP 0.006 A-11P f(48–53)

S11

Table S2. Cont.

Boman AMP Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) Index GRAVY SVM RFC ANN DAC Origin/Position Charge Point pH Mass Da Index Index kcal/mol

QETYKQEKNMA αs2-casein, gen. var. 0.0 6.99 1369.52 3.72 76.85 9.09 −2.145 0.697 0.3575 NAMP 0.856 (Leukocyte elastase EC 3.4.21.37) A-11P f(17–27)

FTKKTKLTEEEKNRLNFLKKIS αs2-casein, gen. var. 4.0 10.53 2696.19 3.06 29.53 70.91 −1.218 0.531 0.497 AMP 0.465 (Leukocyte elastase EC 3.4.21.37) A-11P f(147–168) ITRINKKIEKFQS β-casein, gen. var. 3.0 10.83 1604.92 2.98 88.52 90.00 −0.915 0.492 0.454 AMP 0.739 (Leukocyte elastase EC 3.4.21.37) A2-5P f(23–35) LTLTDV (Leukocyte elastase EC 3.4.21.37) β-casein, gen. var. −1 3.80 660.7 0.00 8.33 178.33 1.15 0.551 0.5625 NAMP 0.002 Poor water solubility A2-5P f(125–130) GWIIPMGILRPYLS (Leukocyte elastase EC Lactoferrin 1.0 9.84 1616.01 −0.9 45.86 139.29 0.821 0.5010.751 AMP 0.619 3.4.21.37) Poor water solubility f(124–137) PPTVMF (Metridin EC 3.4.21.3) β-casein, gen. var. 0 5.96 690.8 −1.13 40.43 48.33 0.833 1.000 0.644 NAMP 0.053 Poor water solubility A2-5P f(152–172) Lactoferrin RPTEGY (Metridin EC 3.4.21.3) 0 6.00 721.7 3.91 58.33 0.00 −2.0 1.000 0.4985 NAMP 0.025 f(428–433) DGGMVF (Metridin EC 3.4.21.3) Lactoferrin −1.0 3.1 624.72 −0.42 28.90 48.33 0.767 0.9980.644 AMP 0.001 Poor water solubility f(60–65) β-casein, gen. var. APKHKEM (Pancreatic elastase II EC 3.4.21.71) 1.1 8.64 840.0 2.62 69.91 14.29 −1.771 1.000 0.4915 AMP 0.020 A2-5P f(103–109) KILDKVGINYWLAHKALCSEKLDQWLCEKL α-lactalbumin, gen. 1.0 8.26 3559.26 0.8 12.53 120.33 −0.113 0.691 0.5285 AMP 0.968 ( EC 3.4.22.17) var. B f(94–123) INYWLAHKALCSEKLDQWLCEKL α-lactalbumin, gen. 0.0 7.0 2805.32 0.82 23.38 110.43 −0.183 0.554 0.556 AMP 0.763 (Glycyl endopeptidase EC 3.4.22.25) var. B f(101–123)

S12

Table S2. Cont.

Boman Net Isoelectric Molecular Instability Aliphatic Sequence (Enzyme) AMP Origin/Position Index GRAVY SVM RFC ANN DAC Charge Point pH Mass Da Index Index kcal/mol APSITCVRRAFALECIRAIAEKKADAV TLDG (Glycyl endopeptidase Lactoferrin f(31–61) 0.9 8.29 3289.88 1.4 46.85 104.19 0.335 0.5830.490 AMP 0.797 EC 3.4.22.25) AVAKFFSASCVPCIDRQAYPNLCQLC KG (Glycyl endopeptidase EC 3.4.22.25) Lactoferrin f(149–175) 1.8 8.3 3033.62 0.61 25.65 76.79 0.375 0.9030.9695 AMP 0.984 Poor water solubility

LEIVPN (Proteinase P1 (lactocepin) αS1-casein, gen. var. B-8P −1 4.0 683.6 −0.07 58.38 178.33 0.650 0.556 0.535 NAMP 0.026 EC 3.4.21.96) Poor water solubility f(109–114)

RDMPI Proteinase P1 (lactocepin) αS1-casein, gen. var. B-8P 0 5.84 630.7 3.27 95.88 78.00 −0.640 1.000 0.627 NAMP 0.015 EC 3.4.21.96 f(102-109) ITRINKKIEKF (Proteinase P1 β-casein, gen. var. A2-5P 3.0 10.83 1389.71 2.71 62.85 106.36 −0.691 0.694 0.460 AMP 0.879 (lactocepin) EC 3.4.21.96) f(23–33) FMAIPPKKN (Proteinase P1 (lactocepin) κ-casein, gen. var. A 2.0 10.6 1045.32 0.63 40.80 54.44 −0.389 0.890 0.381 AMP 0.684 EC 3.4.21.96) f(105–113) GAVAKFF (Proteinase P1 (lactocepin) Lactoferrin f(147–153) 1.0 10.1 738.89 −1.28 −3.56 70.0 1.30 0.6310.4245 AMP 0.597 EC 3.4.21.96) Poor water solubility ALFGKNGKNCPDKFCLFK (Proteinase Lactoferrin f(616–633) 2.9 9.71 2030.45 1.06 −5.79 48.89 −0.317 0.706 0.7655 AMP 0.987 P1 (lactocepin) EC 3.4.21.96) PKHKEMPF (Prolidase L.lactis β-casein, gen. var. A2-5P 1.1 9.88 1013.23 2.15 117.35 0.00 −1.625 0.999 0.4065 NAMP 0.552 s.cremoris H61 EC 3.4.13.9) f(104–111) Fragments of known AMPs are marked in red.

Copyright of International Journal of Molecular Sciences is the property of MDPI Publishing and its content may not be copied or emailed to multiple sites or posted to a listserv without the copyright holder's express written permission. However, users may print, download, or email articles for individual use.