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 Peptides with Potential Antimicrobial Activity: An Approach Based on Bioinformatic Studies Bartłomiej Dziuba 1,* and Marta Dziuba 2 1

2

University of Warmia and Mazury in Olsztyn, Chair of Industrial and Food Microbiology, Cieszynski Square 1, Olsztyn 10-957, Poland 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 enzymes. The obtained peptides were characterized by the following parameters: molecular weight, isoelectric point, composition and number of amino acid 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 peptide 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 Amino Acid

Ala Arg Asn Asp Cys Gln Glu Gly His Ile Leu Lys Met Phe Pro Ser Thr Trp Tyr Val

Average Amino Acid Content (%) a

b

Number of Peptides Containing Given Amino Acid

6.8 8.2 3.4 2.3 5.8 3.7 2.3 10.5 2.1 6.1 9.9 10.3 1.0 4.0 5.0 5.1 3.4 2.1 2.6 5.9

9.8 12.3 6.1 5.7 13.6 7.1 5.6 12.5 5.8 8.1 11.4 12.7 3.7 6.2 9.7 7.5 6.3 5.0 7.0 8.2

286 278 234 166 176 220 175 349 152 312 339 337 113 258 215 282 226 178 157 300

AMPs from Milk Proteins Min.-Max. Amino Acid Content (%) 0–33.3 0–33.3 0–25.0 0–25.0 0–40.0 0–25.0 0–25.0 0–63.1 0–18.4 0–40.0 0–58.3 0–62.5 0–7.4 0–37.5 0–53.2 0–23.1 0–33.3 0–38.5 0–33.3 0–37.3

Average Amino Acid Content (%) a

b

Number of Peptides Containing Given Amino Acid

6.4 7.0 2.0 2.6 3.0 7.0 4.4 3.3 1.2 6.5 8.7 10.2 1.0 2.6 8.0 4.0 5.3 3.2 5.6 8.0

10.6 11.9 9.7 14.0 10.4 11.4 11.3 9.3 6.1 10.2 12.3 15.5 3.6 7.3 11.8 9.1 10.8 8.0 14.4 12.5

36 35 12 11 17 36 23 21 12 38 42 19 16 21 40 26 29 24 23 38

Min.-Max Amino Acid Content (%) 0–25.0 0–25.0 0–25.0 0–25.0 0–25.0 0–33.3 0–25.0 0–20.0 0–14.3 0–33.3 0–37.5 0–33.3 0–7.1 0–16.7 0–28.6 0–21.4 0–33.3 0–16.7 0–33.3 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-casein 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 Value

Min.–Max. Value

Predominant Value (%)

Mean Value

Min.–Max. Value

Predominant Value (%)

Molecular mass (Da) pI Net charge Instability index Aliphatic index GRAVY Boman Index (kcal/mol)

3242.9 9.3 3.7 31.7 83.9 −0.2 1.5

393.5–14,350.8 3.4–13.3 −7.0–20 −50.4–166.2 0–227.5 −3.51–3.6 −2.6–6.8

2000–4000 (47) 9–10 (29) 0–5 (56) 0–60 (71) 40–120 (65) −1–0 (47) 1–2 (28)

1906.9 8.1 2.2 41.5 89.6 −0.4 1.3

393.5–6707.4 3.4–12.0 −7.0–10.1 −30.9–157.7 0.0–226.7 −2.5–2.2 −6.9–5.0

393–1000 (39) 10–11 (25) −2–0 (34) 0–20 (24) 60–100 (49) −1–0 (49) 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

Number of cases

180 160 140 120 100 80 60 40 20 0 -10

-5

0

5

10

15

20

25

Net charge

Many host-defense peptides that are produced naturally by the immune system 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, trypsin, alcalase, chymotrypsin, papain and pancreatic elastase, 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 “Enzyme 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 pancreatic elastase 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],[4144],[45-45],[46-53],[54-57],[58-62],[63-65],[66-71],[72-72],[73-76],[77-81],[82-86],[87-92],[93-93],[9495],[96-98],[99-99],[100-101],[102-106],[107-109],[110-111],[112-112],[113-116],[117-120],[121126],[127-127],[128-129],[130-136],[137-137],[138-138],[139-142],[143-143],[144-149],[150-156],[157158],[159-162],[163-163],[164-167],[168-169],[170-170],[171-176],[177-182],[183-186],[187-187],[188198],[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). Sequence

AMP Origin/

Net

Isoelectric

Molecular Boman Instability Aliphatic

Position

Charge

Point pH

Mass (Da)

Index

Index

Index

0.1

7.72

1004.07

4.63

53.06

12.50

−2.625

1.000

1.0

10.08

2046.41

1.52

20.89

113.89

0.028

−2.0

3.7

2523.89

−0.22

16.36

117.20

−2.0

3.79

999.9

5.3

119.70

3.0

10.45

1003.25

1.67

2.0

10.6

1119.33

0.0

6.94

−2.0

DDKHYQKA (Pancreatic elastase EC 3.4.21.36,

αs2-casein, gen. var. A-11P

Leukocyte elastase EC 3.4.21.37)

f(74–81)

GQRDLLFKDSALGFLRIP

Lactoferrin

(Prolyl oligopeptidase EC 3.4.21.26)

f(294–311)

ADALNLDGGYIYTAGKCGLVPVLAE

Lactoferrin

(V-8 protease EC 3.4.21.19)

f(389–413)

QEQNQEQP (Prolyl oligopeptidase EC 3.4.21.26),

κ-casein, gen. var. A

Thermolysin EC 3.4.24.27)

f(1–8)

KKYKVPQL (Pepsin 1.3 EC 3.4.23.1,

αs1-casein, gen. var. B-8P

Pancreatic elastase EC 3.4.21.71)

f(102–109)

AVAVVKKGSNF (Chymase EC 3.4.212.39,

Lactoferrin

Metridin EC 3.4.21.3)

f(94–104)

EMPFPK(Ficain EC 3.4.22.3,

β-casein, gen. var. A2-5P

Bromelain EC 3.4.22.4)

f(108–113)

EPEQSL

β-lactoglobulin gen. var. B

(Ficain EC 3.4.22.3)

f(112–117)

ITRINKKIEKFQS

β-casein, gen. var. A2-5P

(Leukocyte elastase EC 3.4.21.37)

f(23–35)

ITRINKKIEKF

β-casein, gen. var. A2-5P

(Proteinase P1 (lactocepin) EC 3.4.21.96)

f(23–33)

ALFGKNGKNCPDKFCLFK

Lactoferrin

(Proteinase P1 (lactocepin) EC 3.4.21.96)

f(616–633)

a

GRAVY SVM a RFC b

ANN c

DAC d

0.622

NAMP

0.507

0.659

0.4835

AMP

0.834

0.536

0.568

0.5215 NAMP

0.816

0.00

−3.263

0.976

0.5845

AMP

0.000

46.29

85.00

−1.262

0.952

0.509

NAMP

0.916

0.13

−14.91

97.27

0.591

0.943

0.483

AMP

0.912

747.91

1.17

145.77

0.00

−0.983

1.000

0.571

AMP

0.190

3.13

701.73

2.94

174.73

65.00

−1.517

1.000

0.4695 NAMP

3.0

10.83

1604.92

2.98

88.52

90.00

−0.915

0.492

0.45

AMP

0.739

3.0

10.83

1389.71

2.71

62.85

106.36

−0.691

0.694

0.46

AMP

0.879

2.9

9.71

2030.45

1.06

−5.79

48.89

−0.317

0.706

0.7655

AMP

0.987

, SVM: support vector machines; b, RF: random forest; c, ANN: artificial neural networks; and d, DA: discriminant analysis.

0.569

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, glycine, 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 valine, 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 proteases. 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 caseins 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. Sequence

AMP Origin/Position

Net Charge

Isoelectric Point (pH)

Molecular Instability GRAVY Mass (Da) Index

Aliphatic Index

RPKHPIK

f(1–7) of αS1-casein

3.1

11.17

875.0

−2.029

19.84

55.71

RPKHPIKHQGLPQEVL NENLLRF

Isracidin f(1–23) of αS1-casein

2.2

9.99

2764.2

−0.987

46.49

97.39

GLPQE

f(10–14) of αS1-casein

−1.0

4.00

542.5

−1.040

123.56

78.00

IKHQGLPQE

f(6–14) of αS1-casein

0.1

6.75

1049.1

−1.256

63.66

86.67

VLNENLLR

f(15–22) of αS1-casein

0.0

5.97

970.1

0.075

32.83

182.5

LRLKKYKVPQL

αS1-casein, gen. var. B-8P f(99–109)

4

10.46

1385.7

−0.636

46.18

132.73

SDIPNPIGSENSEK

f(180–193) of αS1-casein

−2

4.14

1486.5

−1.314

32.69

55.71

KTKLTEEEKNRLNFLK KISQRYQKFALPQYLK TVYQHQK

casocidin-I f(150–188) of αS2-casein gen. var. A

7.1

10.08

4869.6

−1.323

27.85

70.00

VYQHQKAMKPWIQPK TKVIPYVRY

f(183–206) αS2-casein gen. var. A

5.1

10.12

3002.6

−0.783

2.25

72.92

LKKISQ

f(165–170) of ovine αS2-casein

2.0

10.00

715.8

−0.633

12.13

130.00

VYQHQKAMKPWIQPK TKVIPYVRYL

f(183–207) of bovine αS2-casein

5.1

10.12

3115.7

−0.600

2.56

85.60

Half Life Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli

1h 2 min 2 min 1h 2 min 2 min 10 h >20 h >20 h 20 h 30 min >10 h 100 h >20 h >10 h 5.5 h 3 min 2 min 1.9 h >20 h >10 h 1.3 h 3 min 2 min 100 h >20 h >10 h 5.5 h 3 min 2 min 100 h >20 h >10 h

Boman Index (kcal/mol) 3.68 2.42 1.29 1.92 2.02 1.67 2.89 2.83 1.37 1.7 1.12

S2 Table S1. Cont. Sequence

AMP Origin/Position

Net Charge

Isoelectric Point (pH)

LKKISQRYQKFALPQY

f(164–179) of bovine αS2-casein

4.0

10.17

2011.4

−0.938

25.45

79.38

PYVRYL

f(203–208) of ovine αS2-casein

1.0

9.00

809.9

−0.117

−4.23

113.33

KTVYQHQKAMKPWIQ PKTKVIPYVRYL

f(181–207) of ovine αS2-casein

6.1

10.22

3345.0

−0.726

3.11

79.26

LKTVYQHQKAMKPWI QPKTKVIPYVRYL

f(180–207) of ovine αS2-casein

6.1

10.22

3458.22

−0.564

0.32

90.36

f(175–207) of ovine αS2-casein

6.1

10.1

4030.8

−0.503

11.18

91.52

f(172–207) of ovine αS2-casein

7.1

10.17

4434.3

−0.589

11.08

83.89

f(164–207) of ovine αS2-casein

10.1

10.35

5451.6

−0.700

9.69

86.36

LKKISQRYQKFALPQ YL

f(164–180) of ovine αS2-casein

4.0

10.17

2124.5

−0.659

24.54

97.65

YQEPVLGPVRGPFPI

β-casein, gen. var. A2-5P f(193–207)

0

6.00

1668.9

−0.033

80.71

90.67

YQEPVLGPVRGPFPIIV

β-casein, genetic variant A2-5P f(193–209)

0

6.00

1881.2

0.482

67.4

120.00

MAIPPKKNQDKTEIPTI NTIASGEPTSTPTTEAV ESTVATLEDSPEVIESPP EINTVQVTSTAV

Kappacin f(107–169) of κ-casein gen. var. A

−7.0

4.04

6707.4

−0.370

64.17

77.66

ALPQYLKTVYQHQKA MKPWIQPKTKVIPYVR YL QKFALPQYLKTVYQHQ KAMKPWIQPKTKVIPY VRYL LKKISQRYQKFALPQY LKTVYQHQKAMKPWI QPKTKVIPYVRYL

Molecular Instability GRAVY Mass (Da) Index

Aliphatic Index

Half Life Mamalian 5.5 h Yeast 3 min E.coli 2 min Mamalian >20 h Yeast >20 h E.coli ? Mamalian 1.3 h Yeast 3 min E.coli 3 min Mamalian 5.5 Yeast 3min E.coli 3min Mamalian 4.4 h Yeast >20 min E.coli >10 min Mamalian 0.8 h Yeast 10 min E.coli 10 h Mamalian 5.5 h Yeast 2 min E.coli 2 min Mamalian 5.5 h Yeast 2 min E.coli 2 min Mamalian 2.8 h Yeast 10 min E.coli 2 min Mamalian 2.8 h Yeast 10 min E.coli 2 min Mamalian 30 h Yeast 20 h E.coli

10 h

Boman Index (kcal/mol) 2.02 1.04 1.33 1.11 0.91 1.06 1.44 1.61 0.3 −0.25

1.55

S3 Table S1. Cont. Sequence

AMP Origin/Position

Net Charge

Isoelectric Point (pH)

Molecular Mass (Da)

GRAVY

Instability Index

Aliphatic Index

YYQQKPVA

f(42–49) of bovine κ-casein

1.0

8.50

996.1

−1.137

62.9

48.75

VQVTSTAV

f(162–169) of bovine κ-casein

0.0

5.49

803.9

1.087

−11.29

121.25

STVATL

f(141–146) of bovine κ-casein

0.0

5.24

590.6

1.267

8.33

130.00

FSDKIAK

f(18–24) of bovine κ-casein

1.0

8.59

807.9

−0.429

−15.69

70.00

IQY

f(28–30) of bovine κ-casein

0.0

6.09

422.79

−0.100

−18.47

130.00

YVL

f(30–32) of bovine κ-casein

0.0

6.09

393.51

2.233

6.67

226.67

EIPT

f(118–121) of bovine κ-casein

−1.0

3.85

458.54

−0.325

3.85

97.50

VESTVATL

f(139–146) of bovine κ-casein

−1.0

4.00

818.9

1.038

32.83

133.75

PAAVRSPAQILQ

f(64–75) of bovine κ-casein

1.0

10.18

1250.4

0.200

157.72

114.17

AVESTVATLEDSPEVIE SPPE

f(136–156) of bovine κ-casein

−6.0

3.40

2199.3

−0.243

134.70

88.10

EQLTK

LTD1-f(1–5) of α-laktoalbumin, gen. var. B

0.0

6.10

617.7

−1.560

46.52

78.00

Half Life Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli

2.8 h 10 min 2 min 100 h >20 h >10 h 1.9 h >20 h >10 h 1.1 h 3 min 2 min 20 h >30 min >10 h 20 h 30 min >10 h 1h 30 min >10 h 100 h >20 h >10 h >20 h >20 h ? 4.4 h >20 h >10 h 1h 30 min >10 h

Boman Index (kcal/mol) 1.38 0.01 −0.37 1.93 0.25 −2.59 1.11 0.08 0.84 1.54 3.11

S4 Table S1. Cont. Sequence

AMP Origin/Position

Net Charge

Isoelectric Point (pH)

Molecular Mass (Da)

GRAVY

Instability Index

Aliphatic Index

GYGGVSLPEWVCTTF

LDT 2, α-chain f(17–31) of α- laktoalbumin, gen. var. B

−1.0

4.00

1615.8

0.453

44.95

64.67

ALCSEK

LDT 2, β-chain f(109–114) of α-laktoalbumin, . gen. varB

0.0

6.04

649.7

−0.017

40.43

81.67

CKDDQNPH

LDC, α-chain f(61–68) of α-laktoalbumin, gen. var. B

−1.0

5.21

956.0

−2.525

5.15

0.00

ISCDKF

LDC, β-chain chain f(75–80) of α-laktoalbumin, gen. var. B

0.0

5.83

711.8

0.267

80.62

65.00

VAGTWY

f(15–20) β-laktoglobulin, gen. var. A

0.0

5.49

695.7

0.450

−30.87

65.00

IPAVFK

f(78–83) β-laktoglobulin, gen. var. A

1.0

8.75

673.8

1.300

10.58

130.00

AASDISLLDAQSAPLR

f(25–40) β-laktoglobulin, gen. var. A

−1.0

4.21

1627.8

0.256

60.91

122.50

VLVLDTDYK

f(92–100) β-laktoglobulin, gen.var. A

−1.0

4.21

1065.2

0.344

−7.81

151.11

FKCRRWQWRMKKLG APSITCVRRAF

bovine lactoferricin B f(17–41) of lactoferrin

7.9

11.84

3125.8

−0.576

77.92

50.80

APRKNVRW

f(1–8) of lactoferrin

3.0

12.01

1026.2

−1.613

95.00

48.75

FKCRRWQWRMKKLG APSITCVRRAFA

f(17–42) of lactoferrin

7.9

11.84

3196.9

−0.485

75.31

52.69

Half Life Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli

>30 h >20 h >10 h 4.4 h >20 h >10 h 1.2 h >20 h >10 h 20 h 30 min >10 h 100 h >20 h >10 h 20 h 30 min >10 h 4.4 h >20 h >10 h 100 h >20 h >10 h 1.1 h 3 min 2 min 4.4 h >20 h >10 h 1.1 h 3 min 2 min

Boman Index (kcal/mol) −0.46 1.29 4.81 1.41 −1.06 −1.36 1.32 0.86 2.75 4.23 2.57

S5 Table S1. Cont. Net Isoelectric Molecular Charge Point (pH) Mass (Da)

GRAVY

Instability Index

Aliphatic Index

3310.0

−0.326

72.89

65.19

9.51

1971.2

−1.056

98.42

48.75

1.0

8.55

461.61

1.075

7.50

122.5

f(17–47) of lactoferrin

7.9

11.37

3811.6

−0.316

78.95

69.35

APRKNVRWCTI

f(1–11) of lactoferrin

3.0

10.86

1343.6

−0.600

101.45

70.91

CRRWQWRMKKLGAPSITCV

f(19–37) of lactofferin

4.9

10.93

2319.8

−0.468

77.14

61.58

FKCRRWQWRMKKLG

f(17–30) of lactofferin

6.0

11.74

1923.3

−1.421

74.31

27.86

LECIRA

f(43–48) of lactofferin

0.0

5.99

703.8

0.767

81.57

146.67

APRKNVRWCTISQPEWFKCRRW QWRMKKLGAPSITCVRRAFA

f(1–42) of lactofferin

9.9

11.63

5150.1

−0.702

84.35

51.19

FKCRRWQWRMKKLGAPSITCVR RAFALECIRA

f(17–48) of lactofferin

7.9

11.37

3883.7

−0.250

76.80

70.31

WKLLSKAQEKFGKNKSR

Lactoferrampin f(268–284) of lactofferin

5.0

10.58

2048.4

−1.482

36.01

51.73

Sequence

AMP Origin/Position

FKCRRWQWRMKKLGAPSITCVR RAFAL

f(17–43) of lactoferrin

7.9

11.84

APRKNVRWCTISQPEW

f(1–16) of lactoferrin

2.0

CIRA

f(45–48) of lactoferrin

FKCRRWQWRMKKLGAPSITCVR RAFALECIR

Half Life Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli Mamalian Yeast E.coli

1.1 h 3 min 2 min 4.4 h >20 h >10 h 1.2 h >20 h >10 h 1.1 h 3 min 2 min 4.4 h >20 h >10 h 1.2 h >20 h >10 h 1.1 h 3 min 2 min 5.5 h 3 min 2 min 4.4 h >20 h >10 h 1.1 h 3 min 2 min 2.8 h 3 min 2 min

Boman Index (kcal/mol) 2.29 2.72 3.69 2.5 2.74 2.16 3.55 1.46 2.63 2.36 2.97

S6 Table S1. Cont. Net Isoelectric Molecular Charge Point (pH) Mass (Da)

GRAVY

Instability Index

Aliphatic Index

1122.2

−1.156

−25.68

97.78

9.01 9.99

990.1 1200.4

1.156 −1.33

−0.54 58.29

107.78 68.00

4.37

802.9

−0.583

40.43

130.00

Sequence

AMP Origin/Position

ADRDQYELL

Antiviral f(241–249) of lactofferin

−2.0

4.03

TRVVWCAVG GRDPYKLRPV

Lactoferrin (fr. 343–351) Lactoferrin (fr. 68–77)

1 2

EDLIWK

Antiviral f(263–268) of lactofferin

−1.0

Half Life Mamalian Yeast E.coli Mamalian Yeast Mamalian Yeast E.coli

4.4 h >20 h >10 h 7.2 h >20 h 1h 30 min >10 h

Boman Index (kcal/mol) 3.68 −0.11 3.43 1.48

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). Sequence (Enzyme)

AMP

Net

Isoelectric

Molecular

Origin/Position

Charge

Point pH

Mass Da

0.1

7.72

1004.07

3.0

10.84

0.0

DDKHYQKA (Pancreatic elastase EC 3.4.21.36,

αs2-casein, gen. var.

Leukocyte elastase EC 3.4.21.37)

A-11P f(74–81)

FTKKTKL (Pancreatic elastase EC 3.4.21.36) DTDYKKYL (Pancreatic elastase EC 3.4.21.36, Pepsin 1.3 EC 3.4.23.1, Pancreatic elastase EC 3.4.21.71) KNTMEHV (Pancreatic elastase EC 3.4.21.36, Leukocyte elastase EC 3.4.21.37, Proteinase P1 (lactocepin) EC 3.4.21.96) NEEEYSI (Pancreatic elastase EC 3.4.21.36) QNNDSTEYG (Pancreatic elastase EC 3.4.21.36) TKESPQTHYYA (Pancreatic elastase EC 3.4.21.36) PNSKEKYYG (Pancreatic elastase EC 3.4.21.36)

αs2-casein, gen. var. A-11P f(147–153) β-lactoglobulin, gen. var. B f(96–103) αs2-casein, gen. var. A-11P f(1–7) αs2-casein, gen. var. A-11P f(48–54) α-lactalbumin, gen. var. B f(43–51) Lactoferrin f(84–94) Lactoferrin f(517–525)

LQGAVAKFFSASCVP (Prolyl oligopeptidase EC

Lactoferrin

3.4.21.26) Poor water solubility

f(145–159)

GQRDLLFKDSALGFLRIP

Lactoferrin

(Prolyl oligopeptidase EC 3.4.21.26)

f(294–311)

QEQNQEQP (Prolyl oligopeptidase EC 3.4.21.26),

κ-casein, gen. var. A

Thermolysin EC 3.4.24.27)

f(1–8)

Boman Index

Instability Aliphatic

GRAVY SVM

RFC

ANN

DAC

Index

Index

4.63

53.06

12.50

−2.625

1.000 0.622 NAMP 0.507

865.08

1.98

−3.56

55.71

−0.929

0.998 0.5215

6.75

1045.16

3.3

−10.04

48.75

−1.788

1.000 0.4615 NAMP 0.556

0.1

6.75

857.9

2.83

−14.33

41.43

−1.243

0.928 0.609 NAMP 0.029

−3

3.67

882.8

3.66

101.71

55.71

−1.657

1.000 0.5665 NAMP 0.152

−2

3,67

1026.9

4.39

42.26

0.00

−2.300

1.000 0.5875 NAMP 0.027

0.1

6.41

1324.4

2.68

156.75

9.09

−1.700

0.999 0.524 NAMP 0.027

1.0

9.47

1085.18

3.03

13.17

0.00

−2.244

1.000 0.252 NAMP 0.873

1.0

8.97

1524.81

−0.58

54.95

84.67

0.980

0.785 0.2385

AMP

0.716

1.0

10.08

2046.41

1.52

20.89

113.89

0.028

0.659 0.4835

AMP

0.834

−2

3.79

999.9

5.3

119.70

0.00

−3.263

0.976 0.5845

AMP

0.000

kcal/mol

AMP

0.940

S8 Table S2. Cont. Sequence (Enzyme)

AMP

Net

Origin/Position

Charge

Point pH

Mass Da

−2.0

3.7

2523.89

2.0

10.18

3.0

ADALNLDGGYIYTAGKCGLVPVLAE

Lactoferrin f(389–

(V-8 proteaza EC 3.4.21.19)

413)

GYLAVAVVKKAN

Lactoferrin

(Glutamyl endopeptidase II EC 3.4.21.82)

f(432–443)

KKYKVPQL (Pepsin 1.3 EC 3.4.23.1, Pancreatic elastase EC 3.4.21.71) NVPGEIVESL (Pepsin 1.3 EC 3.4.23.1, Cathepsin G EC 3.4.21.20, Metridin EC 3.4.21.3, Pancreatic elastase II EC 3.4.21.71)

B-8P

β-casein gen. var. B f(7–17)

Poor water solubility

f(120–132)

AGDDQG (Thermolysin EC 3.4.24.27)

Instability

Aliphatic

Index

Index

−0.22

16.36

117.20

0.536

0.568 0.5215 NAMP 0.816

1232.49

−0.46

−0.31

130.0

0.733

0.893 0.4395

10.45

1003.25

1.67

46.29

85.00

−1.262

0.952 0.509 NAMP 0.916

−2

3.79

1056.1

0.47

50.75

136.00

0.340

0.586 0.516 NAMP 0.001

1.0

11.04

1370.68

−0.82

52.79

127.69

0.923

0.171 0.564

AMP

0.645

3.0

10.84

751.92

3.13

8.33

0.00

−1.717

1.000 0.5285

AMP

0.986

−3

3.67

840.8

4.11

101.71

14.29

−2.043

0.933 0.6535 NAMP 0.002

−2

3.56

561.5

3.21

8.33

16.67

−1.583

0.624 0.551 NAMP 0.000

−2

3.80

710.7

2.09

72.53

48.33

−0.250

1.000 0.531 NAMP 0.005

1.0

9.39

901.03

0.59

77.66

14.29

−1.314

0.999 0.3085 NAMP 0.719

Index kcal/mol

GRAVY SVM

RFC

ANN

AMP

DAC

0.928

f(102–109)

Lactoferrin

ANEEEYS (Thermolysin EC 3.4.24.27)

Boman

αS1-casein, gen. var.

GRSAGWIIPMGIL (Pepsin 1.3 EC 3.4.23.1)

FTKKTK (Thermolysin EC 3.4.24.27)

Isoelectric Molecular

αs2-casein, gen. var. A-11P f(147–152) αs2-casein, gen. var. A-11P (fr 47–53) Lactoferrin f(506–511)

ESTEVF (Cathepsin G EC 3.4.21.20,

αs2-casein, gen. var.

(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

Poor water solubility

f(57–63)

S9 Table S2. Cont. Sequence (Enzyme)

AMP

Net

Origin/Position

Charge

Point pH

Mass Da

−1

3.80

636.6

2.0

10.6

1.0

TDAPSF (Chymase EC 3.4.21.39),

αS1-casein, gen. var.

Metridin EC 3.4.21.3)

B-8P f(174–179)

AVAVVKKGSNF (Chymase EC 3.4.212.39,

Lactoferrin

Metridin EC 3.4.21.3)

f(94–104)

IIPMGILRPY (Chymase EC 3.4.21.39)

Lactoferrin

Poor water solubility

f(126–135)

VTAIANLKKCSTSPLLEACAF

Lactoferrin

(Chymase EC 3.4.21.39) Poor water solubility

f(666–686)

PQRDMPIQAF (Chymase EC 3.4.21.39,

β-casein gen. var. B

Metridin EC 3.4.21.3)

f(181–190)

EAGRDPY (Chymase EC 3.4.21.39),

Lactoferrin

Metridin EC 3.4.21.3)

f(66–72)

GTKESPQTHY (Chymase EC 3.4.21.39),

Lactoferrin

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. Poor water solubility DMPIQA (Papain 3.4.22.2, Bromelain EC 3.4.22.4) QPTTMA (Papain 3.4.22.2, Ficain EC 3.4.22.3, Leukocyte elastase EC 3.4.21.37, Bromelain EC 3.4.22.4) Poor water solubility

B-8P f(204–209) β-casein gen. var. B f(184–189) κ-casein, gen. var. A f(91–96)

Isoelectric Molecular

Boman

Instability

Aliphatic

Index

Index

1.65

72.53

16.67

−0.333

0.901 0.5085 NAMP 0.011

1119.33

0.13

−14.91

97.27

0.591

0.943 0.483

AMP

0.912

9.84

1172.5

−0.79

47.87

156.0

0.980

0.371 0.595

AMP

0.837

0.9

8.29

2180.62

0.00

43.29

107.14

0.748

0.685 0.917

AMP

0.852

0

6.27

1202.3

2.26

72.20

49.00

−0.720

0.715 0.697 NAMP 0.018

−1

4.37

806.8

3.97

8.57

14.29

−1.857

0.848 0.5785 NAMP 0.000

0.1

6.75

1147.2

3.03

126.91

0.00

−1.96

0.718 0.625 NAMP 0.007

0

5.19

747.9

−0.74

121.03

65.0

0.300

0.976 0.5275 NAMP 0.046

−1

3.80

673.7

0.86

81.57

81.67

−0.067

0.768 0.525 NAMP 0.005

0

5.52

647.7

1.08

61.00

16.67

−0.467

1.000 0.6035 NAMP 0.015

Index kcal/mol

GRAVY SVM

RFC

ANN

DAC

S10 Table S2. Cont. Sequence (Enzyme)

EERLHSMK (Bromelain EC 3.4.22.4)

AMP

Net

Origin/Position

Charge

Point pH

Mass Da

0.1

6.86

1029.1

−3

3.57

0

αS1-casein, gen. var. B-8P f(117–124)

SESTEDQA (Ficain EC 3.4.22.3,

αS1-casein, gen. var.

Bromelain EC 3.4.22.4)

B-8P f(46–53)

PQRDMPIQA (Ficain EC 3.4.22.3),

β-casein, gen. var.

Bromelain EC 3.4.22.4)

A-5P f(181–189)

EMPFPK (Ficain EC 3.4.22.3,

β-casein, gen. var.

Bromelain EC 3.4.22.4)

A2-5P f(108–113)

EPEQSL (Ficain EC 3.4.22.3) ESPQTHY (Ficain EC 3.4.22.3, Bromelain EC 3.4.22.4)

β-lactoglobulin, gen. var. B f(112–117) Lactoferrin f(86–92)

PQLEIV (Leukocyte elastase EC 3.4.21.37,

αS1-casein, gen. var.

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.

Proteinase P1 (lactocepin EC 3.4.21.37)

B-8P f(192–199)

NEEEYS (Leukocyte elastase EC 3.4.21.37)

αs2-casein, gen. var. A-11P f(48–53)

Isoelectric Molecular

Boman

Instability

Aliphatic

Index

Index

4.35

49.50

48.75

−1.712

0.717 0.5785 NAMP 0.001

865.8

4.43

105.05

12.50

−1.812

0.998 0.637 NAMP 0.004

6.27

1055.2

2.84

79.11

54.44

−1.111

0.827 0.6745 NAMP 0.010

0

6.94

747.91

1.17

145.77

0.00

−0.983

1.000 0.571

−2.0

3.13

701.73

2.94

174.73

65.00

−1.517

1.000 0.4695 NAMP 0.569

−0.9

5.24

860.8

3.3

189.14

0.00

−2.086

0.551 0.658 NAMP 0.004

−1

4.00

697.8

−0.25

58.38

178.33

0.650

0.945 0.553 NAMP 0.020

0

6.10

1005.2

0.98

93.28

48.75

−0.700

1.000 0.644 NAMP 0.021

−3

3.67

769.7

5.1

117.0

0.00

−2.683

1.000 0.6795 NAMP 0.006

Index kcal/mol

GRAVY SVM

RFC

ANN

AMP

DAC

0.190

S11 Table S2. Cont. Sequence (Enzyme)

AMP

Net

Origin/Position

Charge

Point pH

Mass Da

0.0

6.99

1369.52

4.0

10.53

3.0

QETYKQEKNMA

αs2-casein, gen. var.

(Leukocyte elastase EC 3.4.21.37)

A-11P f(17–27)

FTKKTKLTEEEKNRLNFLKKIS

αs2-casein, gen. var.

(Leukocyte elastase EC 3.4.21.37)

A-11P f(147–168)

ITRINKKIEKFQS

β-casein, gen. var.

(Leukocyte elastase EC 3.4.21.37)

A2-5P f(23–35)

LTLTDV (Leukocyte elastase EC 3.4.21.37)

β-casein, gen. var.

Poor water solubility

A2-5P f(125–130)

GWIIPMGILRPYLS (Leukocyte elastase EC

Lactoferrin

3.4.21.37) Poor water solubility

f(124–137)

PPTVMF (Metridin EC 3.4.21.3)

β-casein, gen. var.

Poor water solubility

A2-5P f(152–172)

RPTEGY (Metridin EC 3.4.21.3)

Lactoferrin f(428–433)

DGGMVF (Metridin EC 3.4.21.3)

Lactoferrin

Poor water solubility

f(60–65)

APKHKEM (Pancreatic elastase II EC 3.4.21.71)

β-casein, gen. var. A2-5P f(103–109)

KILDKVGINYWLAHKALCSEKLDQWLCEKL

α-lactalbumin, gen.

(Calpain EC 3.4.22.17)

var. B f(94–123)

INYWLAHKALCSEKLDQWLCEKL

α-lactalbumin, gen.

(Glycyl endopeptidase EC 3.4.22.25)

var. B f(101–123)

Isoelectric Molecular

Boman

Instability

Aliphatic

Index

Index

3.72

76.85

9.09

−2.145

0.697 0.3575 NAMP 0.856

2696.19

3.06

29.53

70.91

−1.218

0.531 0.497

AMP

0.465

10.83

1604.92

2.98

88.52

90.00

−0.915

0.492 0.454

AMP

0.739

−1

3.80

660.7

0.00

8.33

178.33

1.15

0.551 0.5625 NAMP 0.002

1.0

9.84

1616.01

−0.9

45.86

139.29

0.821

0.501 0.751

0

5.96

690.8

−1.13

40.43

48.33

0.833

1.000 0.644 NAMP 0.053

0

6.00

721.7

3.91

58.33

0.00

−2.0

1.000 0.4985 NAMP 0.025

−1.0

3.1

624.72

−0.42

28.90

48.33

0.767

0.998 0.644

AMP

0.001

1.1

8.64

840.0

2.62

69.91

14.29

−1.771

1.000 0.4915

AMP

0.020

1.0

8.26

3559.26

0.8

12.53

120.33

−0.113

0.691 0.5285

AMP

0.968

0.0

7.0

2805.32

0.82

23.38

110.43

−0.183

0.554 0.556

AMP

0.763

Index kcal/mol

GRAVY SVM

RFC

ANN

AMP

DAC

0.619

S12 Table S2. Cont. Sequence (Enzyme)

Net

Isoelectric

Molecular

Charge

Point pH

Mass Da

Lactoferrin f(31–61)

0.9

8.29

3289.88

Lactoferrin f(149–175)

1.8

8.3

−1

Boman

Instability

Aliphatic

Index

Index

1.4

46.85

104.19

3033.62

0.61

25.65

4.0

683.6

−0.07

0

5.84

630.7

3.0

10.83

2.0

Lactoferrin f(147–153) Lactoferrin f(616–633)

AMP Origin/Position

Index

RFC

ANN

DAC

0.335

0.583 0.490

AMP

0.797

76.79

0.375

0.903 0.9695

AMP

0.984

58.38

178.33

0.650

0.556 0.535 NAMP 0.026

3.27

95.88

78.00

−0.640

1.000 0.627 NAMP 0.015

1389.71

2.71

62.85

106.36

−0.691

0.694 0.460

AMP

0.879

10.6

1045.32

0.63

40.80

54.44

−0.389

0.890 0.381

AMP

0.684

1.0

10.1

738.89

−1.28

−3.56

70.0

1.30

0.631 0.4245

AMP

0.597

2.9

9.71

2030.45

1.06

−5.79

48.89

−0.317

0.706 0.7655

AMP

0.987

1.1

9.88

1013.23

2.15

117.35

0.00

−1.625

0.999 0.4065 NAMP 0.552

kcal/mol

GRAVY SVM

APSITCVRRAFALECIRAIAEKKADAV TLDG (Glycyl endopeptidase EC 3.4.22.25) AVAKFFSASCVPCIDRQAYPNLCQLC KG (Glycyl endopeptidase EC 3.4.22.25) Poor water solubility LEIVPN (Proteinase P1 (lactocepin)

αS1-casein, gen. var. B-8P

EC 3.4.21.96) Poor water solubility

f(109–114)

RDMPI Proteinase P1 (lactocepin)

αS1-casein, gen. var. B-8P

EC 3.4.21.96

f(102-109)

ITRINKKIEKF (Proteinase P1

β-casein, gen. var. A2-5P

(lactocepin) EC 3.4.21.96)

f(23–33)

FMAIPPKKN (Proteinase P1 (lactocepin)

κ-casein, gen. var. A

EC 3.4.21.96)

f(105–113)

GAVAKFF (Proteinase P1 (lactocepin) EC 3.4.21.96) Poor water solubility ALFGKNGKNCPDKFCLFK (Proteinase P1 (lactocepin) EC 3.4.21.96) PKHKEMPF (Prolidase L.lactis

β-casein, gen. var. A2-5P

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.

New milk protein-derived peptides with potential antimicrobial activity: an approach based on bioinformatic studies.

New peptides with potential antimicrobial activity, encrypted in milk protein sequences, were searched for with the use of bioinformatic tools. The ma...
339KB Sizes 2 Downloads 6 Views