No competing interest
This is an open-access article distributed under the terms of the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.
Matrix metalloproteinases (MMPs), the class of enzymes involved in the degradation of extra cellular matrix and surrounding cells are known to be expressed during cancer cell invasion, arthritis and metastasis and MMP-9, in particular seems to be a key protease associated with tumor progression. Therefore development of inhibitors of MMP-9 is a challenging task which can have therapeutic benefit for patients suffering from various cancers.β-lactoglobulin (β-LG), a core member of lipocalin (Lcn)family is an abundant whey protein in milk of ruminants and has attracted considerable attention as a rich source of bioactive peptides involved in diverse biological functions. Many of the lipocalins are shown to possess protease inhibitory activity including Lcn-2 and β-LG being a major lipocalin, its tryptic peptides were verified for the inhibitory effect on MMP-9.Hence, to investigate the effect of β-LG peptides on the protease activity, we have used β-LG purified from buffalo colostrum.The tryptic peptides of β-LG were separated and analyzed by LC-MS/MS. Based on in silico analysis the sequence information for singly charged peptide ions, m/z-573 IIAEK (P1), m/z-673 GLDIQK (P2), m/z-916 IDALNENK(P3) and m/z-933 IIVTQTMK (P4) were deduced and the interaction of these peptides with MMP-9 was verified by validated target of Matrix metalloproteinase-9 (PDB: 1GKC) using MolDock. A good correlation was observed in binding affinity of hypocholesterolemic peptides (IIAEK & GLDIQK) along with other two peptides (IDALNENK & IIVTQTMK) of β-LG from the buffalo colostrum implicating its utilization in functional foods.Keywords: β-lactoglobulin, LC-MS/MS, Collision induced dissociation, In silico analysis, Matrix metalloproteinases, Molecular docking.
Milk contains a large number of bioactive peptides with various biological attributes. These peptides either exist naturally in milk or released after milk protein proteolysis. The natural bioactive peptides found in milk include epidermal growth factor (EGF), transforming growth factor (TGF), nerve growth factor (NGF), insulin, and insulin like growth factors I and II (IGF-I and IGF-II) and the concentrations of these peptides are generally high in colostrum, often higher than those found in blood circulation. Due to a limited protease activity in the gastrointestinal tract of neonates and the existence of protease inhibitors in the milk, these milk borne bioactive peptides are likely to survive the gastrointestinal digestion in suckling neonates. It is speculated that colostrum-borne bioactive peptides may play a role in regulating postnatal gut development in suckling neonates. Research in the field of bioactive peptides has intensified during the past two decades and has been extensively reviewed and an increasing number of in vitro and in vivo studies revealed that peptides released after enzymatic digestion can play very important role in human health [1-4] and modulate several regulatory processes [5]. Thus the peptides derived from milk proteins have been shown to improve immunity in the newborn with varied bioactivities like antihypertensive, antimicrobial, antioxidative, antithrombotic, immunomodulatory, mineral-carrying and opioid. Hence there is increasing commercial interest in the production of bioactive peptides with the purpose of using them as active ingredients in functional foods. Presently peptides derived from natural sources are of great demand with multifunctional abilities.β-lactoglobulin (β-LG), a small soluble secretary protein is the major whey protein in ruminants belong to lipocalin family of proteins. Lipocalins are small secretary proteins consisting of a highly conserved eight stranded anti-parallel hydrophobic barrel bordered on one side by β-helix. More than 20 lipocalins have been identified so far involved in diverse biological functions varying from transport of small hydrophobic molecules to tumor progression and metastasis. The interest concerning lipocalins in cancer has so far been focused to the variations in concentration and the modification of lipocalin expression in distinct cancer forms [6].In addition, lipocalins have been assigned a role in cell regulation although influence of extracellular lipocalins on intracellular cell regulation events however is not fully understood. Some of the lipocalins also have protease inhibitory properties and possess the ability to interact with tumor specific proteases, revealing another pathway for lipocalins to interact with cancer cells. Lcn-1/ Tear lipocalin [7], Von Ebnerâs gland protein [8], Mouse urinary protein [9], Triabin[10] are reported to have cysteine and papain, trypsin, serine protease inhibitory properties while, Lcn-2 is a ligand for matrix metallo protenase-9 (MMP-9) upregulated in many of the solid tumors [11]. MMP-9, a Zn²± dependent endopeptidase is markedly increased during metastasis and after myocardial infarction [12]. Therefore inhibition of MMP-9 is considered to be a novel approach controlling cancer and hypertension.Development of low molecular weight protease inhibitors with minimal cytotoxicity may have effective therapeutic potential as a new class of anticancer drugs [13]. Hence histidine containing tetrapeptide sublibraries from different natural and non natural amino acids, from which new inhibitors of gelatinases (MMP-9 and MMP-2) have been generated as potential anticancer drugs [14]. MMPs catalytic domain containing N terminal HExxHxxGxxH Zn²+ binding consensus sequence and C terminal hemopexin-like domain are known to play key role in determining the substrate specificities of the various MMPs [15]. A conserved cysteine residue in the prodomain ligates the Zn²+ in the active site of the proteases, causing inactivation of the catalytic domain. Hence cysteine containing peptides mimik this to inactivate the enzyme activity and also the presence of free sulfhydryl or hydroxamate groups were found capable of inhibiting the activity by chelating Zn²+ in the active site of MMPs [16].The bioactivities reported for milk β-LG derived peptides exhibits antihypertensive, antimicrobial, antioxidant, opioid and immunostimulatory effects [17-19]. The peptides of β-LG released after the proteolysis with trypsin, chymotrypsin, thermolysin and protenase K were shown to possess angiotensin I converting enzyme (ACE) inhibitory property, a key enzyme that controls hypertension [20,21]. Further, the tryptic peptides of β-LG inhibited cholesterol absorption as evaluated by in vivo studies and was shown to be greater than β-sitosterol [12]. Hence, β-LG is incorporated as active ingredient in Biozate, a pharmaceutical preparation [1]. Although structure function relationship is well ascribed for the peptides having ACE inhibitory activity, no such correlation exists for vast majority of functional peptides of β-LG. Incidentally; several lipocalins are reported to have protease inhibitory activity and are known to interact with tumor specific proteases [7-10]. Hence, in the present study we have made an attempt to deduce the sequence information of some of the tryptic peptides of β-LG purified from buffalo colostrum by in-silico analysis and assign their functional attributes as inhibitors of MMP-9 by docking analysis.
The β-LG isolated and purified from buffalo colostrum [22], was used for the analysis. Isolation and purification of BLG-colBLG-col from buffalo (B. bubalis) colostrum was isolated by ammonium sulfate fractionation and purified by gel permeation chromatography on a Sephadex G-100 as described earlier [22]. Trypsin digestion of BLG-colThe enzyme digestion of BLG-col was performed using trypsin in 50mM NH4HCO3 (pH 8). Enzyme to BLG-col ratio was maintained at 1:50 and the reaction mixture was incubated at 37°C for 12h [21]. The digested sample was fractionated on C18 column (Vydac, 4.6Ã250mm, 5µm) using a linear gradient of acetonitrile (2%/min) containing 0.1% trifluroacetic acid. Fractions were collected and concentrated in a Speed Vac system and reconstituted in 5µL of 0.2% v/v TFA solution. LC-MS/MS analysisThe LC-MS/MS analysis of the tryptic peptides of BLG-col were obtained on a HCT ULTRA ETD II mass spectrometer (Bruker Daltonics, Bremen, Germany), equipped with two Octapole followed by an ion trap and a separate hexapole next to negative chemical ionization chamber (nCI). Fluoranthene was used as the electron transfer reagent and methane was used for the chemical ionization. The reaction time of electron transfer was typically set at 100ms. Helium was used as the collision gas for collision induced dissociation (CID) experiments.Database searchThe raw data files of MS/MS in .mgf file format were extracted using data analysis 4.0 software and further searched in Mascot MS/MS ion search database http://www.matrixscience.com/cgi/search_form.pl?FORMVER=2&SEARCH=MIS using MSDB (MSDB 20060831; 3239079 sequences; 1079594700 residues), non-redundant NCBI (NCBInr 20090820; 9511482 sequences; 3251602805 residues) and SWISS-PROT (SwissProt 57.6; 495880 sequences; 174780353 residues) databases. Molecular Docking analysisThe model for MMP-9 used in this study was derived from the co-ordinates of the structure labeled 1GKC in the protein data bank, which represents the human MMP-9 bound to NFH at 2A° resolution [23]. The structure of the peptide was constructed using automated ArgusLab 4.0.1 web server (http://www.arguslab.com). The possibility of binding, precise location of binding site and the mode of binding of the ligand was carried out using an automated docking software Molegro Virtual Docker 2008, version 3.2.1. The possible interacting conformations and orientations were analyzed by clustering methods embedded in MolDock [24].All docking studies were carried out using the construct MMP-9 molecular model complexed with inhibitor NFH (N-formylhydroxylamine). The molecular structure of MMP-9 active domain was prepared from the coordinate sets of 1GKC chain A as the template [23]. All water molecules were excluded where as, zinc and chloride atoms were retained in the active site fixed to their crystal position through out the docking process. The binding site of the enzyme was computed within a cavity volume 77.312 A° and surface 271.36 A° such that the binding site of MMP-9 was well sampled with grid resolution of 0.3A°. The value of population size and maximum interactions 100 and 10,000 respectively were used for each run and 5 best poses were retained for ligand NFH docked into the template structure as a positive control and the RMSD was obtained.
The LC-MS spectra of the tryptic peptides of colostrum β-LG clearly indicated the presence of tryptic peptides of masses m/z-573 [M+H]+, m/z-673 [M+H]+,m/z-916 [M+H]+ and m/z- 933 [M+H]+ in analogy to the PMF data reported earlier for BLG-col [22] These ions were numbered accordingly as P1-P4 respectively and analyzed for sequence information. LC-ESI-MS/MS analysis using ion trap mass spectrometer with CID capability facilitated in fragmentation analysis of singly charged ion pairs m/z-573 (P1, Illustration 1A), 673 m/z-( P2, Illustration 1B), m/z-916 (P3, Illustration 1C) and m/z-933 (P4, Illustration 1D). The direct submission of .mgf file for MS/MS of m/z- 573, 673, 916 and 933 to Mascot MS/MS ion search yielded peptide sequences IIAEK (P1), GLDIQK (P2), IDALNENK (P3) and IIVTQTMK (P4) respectively (Supplimentary file 1-4). Notably, the sequences of P1 and P2 completely matched to hypocholesterolaemic peptides of bovine milk β-LG and these sequences were deposited in the peptidome database with an accession number PSE 141. Similarly IIVTQTMK corresponded to N- terminal sequence of buffalo colostrum [22] and milk β-LG [25] while, IDALNENK matched to the milk β-LG. The in-silico analysis of CID-MS/MS based fragmented peptides thus revealed complimentary data to PMF data of colostrum β-LG, besides the sequence tags were in compliance to f1-8 (P4), f9-14 (P2), f71-75 (P1) and f84-91(P3) of the β-LG of buffalo milk and cDNA sequences deposited in the database [26].In the present study, peptide sequences IIAEK (m/z-573), GLDIQK (m/z-673), IDALNENK (m/z-916) and IIVTQTMK (m/z-933) from BLG-col were confirmed after in silico analysis. In order to understand the interaction of these peptides with MMP-9, we performed docking calculations. The catalytic centre of MMP-9 is composed of the active site Zn²+, co-ordinated by histidine (401, 405 and 411) and glutamic acid (402) residues. The cavity information of the enzyme was obtained as cavity volume 77.312 A° and surface 271.36 A°. We observed that hydrogen bond, electrostatic interaction and hydrophobic interaction were the main stabilizing forces in the enzyme-ligand interaction. The known MMP-9 inhibitor NFH interacted with His411, His405, Glu402, Leu188, Pro421, Gly186, Leu187, Leu188 and Tyr423 while, the interactions of the peptides IIAEK (P1), GLDIQK (P2), IDALNENK (P3) and IIVTQTMK (P4) to the active pocket amino acids of MMP-9 were as shown in (Illustration 2,3,4 &5). Analogous to the well known inhibitor, all the peptides were found to interact with the active site amino acid residues based on hydrogen bond parameters. To evaluate the effectiveness of docked tryptic peptides in the Molegro docking programme, we first docked NFH and compared the resulted geometry with the corresponding crystal structure of MMP-9 complex (PDB: 1GKC). The best return poses in presence of Zn(II) with the docking score showed same interactions to be considered for high potency of the peptide and interestingly binding energy of P3 was found much higher (-74.69 kJ/mol) than NFH (-120.745 kJ/mol). Similarly, the binding energies of P1(-137.536 kJ/mol), P2(-130.475 kJ/mol) and P4 (-123.404 kJ/mol) were found in close match to the effect seen with respect with the inhibitor. These validated results of the docking analysis demonstrate the inhibitory effect of BLG-col peptides on MMP-9, hence a newer role can be assigned for BLG-col.
Currently whey proteins have attracted attention as biologically active ingredients with important nutritional and functional properties. The peptides of β-LG are known to exhibit well-defined biological functions including ACE inhibitory activity. Further, as ACE inhibitory peptides also known to possess MMP-9 inhibition property, the present study would be instrumental to verify similar such effects for ACE inhibitory peptides of BLG-col. Since naturally derived peptides would be expected not to have side effects otherwise associated with synthetically produced drugs, our study explores structural specificity of the milk derived peptides and provides valuable insights in designing highly specific inhibitors of MMP-9 from food origin.
We thank University Grant Commission, New Delhi, India for sanction of grant (F31-294/2005-06) to undertake this investigation and Rohit A. Chougule thank ICMR for award of SRF.
1. Korhonen H, Pihlanto A. Bioactive peptides: Production and functionality. Int. Dairy J. 2006;16: 945-60.2. Bouhallab S, Bougle D. Biopeptides of milk: caseinophosphopeptides and mineral bioavailability. Reprod. Nutr. Develop.2004;44: 493-98.3. Tirelli A, Denoni I, Resmini P. Bioactive peptides in milk products. Ital. J. Food Sci.1997;9: 91-98.4. Xu RJ. Bioactive peptides in milk and their biological and health implications. Food Rev. Int. 1998; 14:1-16.5. Jorgensen ALW, Madsen HRJ, Stagsted J. Colostrum and bioactive, Colostral peptides differentially modulate the innate immune response of intestinal epithelial cells. J. Pep. Sci. 2010;16: 21-30. 6. Bratt, T. Lipocalins and cancer. Biochim Biophys Acta.2000;1482: 318-26.7. Redl B. Human tear lipocalin. Biochim Biophys Acta.2000;1482: 241-248.8. Wojnar P, van't Hof W, Merschak P, Lechner M, Redl B. The N-terminal part of recombinant human tear lipocalin/von Ebner's gland protein confers cysteine proteinase inhibition depending on the presence of the entire cystatin-like sequence motifs. J. Biol. Chem. 2001;382: 1515-20.9. Fuentes-Prior P, Noeske-Jungblut C, Donner P, Schleuning WD, Huber R. And Bode W. Proc. Natl. Acad. Sci. 1997; 94:11845-50.10. Mahadevan NR, Rodvold J, Almanza G, Perez AF, Wheeler MC and Zanetti M ER. Stress drives Lipocalin 2 upregulation in prostate cancer cells in an NF-kappaB-dependent manner. BMC cancer 2011;11: 229.11. Kai H, Lekda H, Yasukawa H, Kai M, Seki Y, Kuwahara F, Ueno T,. Sugi K, Imaizumi T. Peripheral blood levels of matrix metalloproteases-2 and -9 are elevated in patients with acute coronary syndromes. J. Am. Coll. Cardiol.1998;32: 368-72.12. Nagaoka S, Futamura Y, Miwa K, Awano T, Yamauchi K, Kanamaru Y, Tadashi K, Kuwata T. Identification of novel hypocholosterolemic peptides derived from bovine milk β-lactoglobulin. Biochem. Biophy. Res. Commun. 2000; 281:11-17.13. Ferry G, Boutin JA, Atassi G, Fauchère JL, Tucker GC. Selection of a histidine-containing inhibitor of gelatinases through deconvolution of combinatorial tetrapeptide libraries. Molecular Diversity 1996; 2:135-46.14. Brown, ED. Matrix metalloproteinase inhibitors: A new class of anticancer agent. Curr. Opin. Invest. Drugs 1993; 2: 617-626.15. Borkakoti N. Matrix metalloproteases: variations on a theme. Prog Biophys Mol Biol. 1998; 70(1):73-94.16. Jensen SA, Andersen P, Vrhovski B, Weiss AS. Rational design of tropoelastin peptide-based inhibitors of metalloproteinases. Arch Biochem. Biophys. 2003; 409(2):335-40.17. Biziulevicius GA, Kislukhina OV, Kazlauskaite J, Zukaite V. Food protein enzymatic hydrolsates possess both antimicrobial and immunostimulatory activities: a âcause and effectâ theory of bifunctionality FEMS Immunol. Med. Micro.2006;46: 131-38.18. Hernandez-Ledesma B, Davalos A, Bartolome B, Amigo LJ. Preparation of antioxidant enzymatic hydrolyzates from β-lactoglobulin and β-lactoglobulin. Agri. Food Chem. 2005;53: 588-93.19. Antila BN, Paakari I., Jarvinen A, Mattila MJ, Laukkanen M, pihlanto-Leppala A, Mantsala P, Hellman J. Opioid peptides derived from in-vitro proteolysis of bovine whey proteins. Int. Dairy J. 1991;1: 215-29.20. Pihlanto-Leppala A, Rokka T, Korhonen H. Angiotensin I converting enzyme inhibitory peptides derived form bovine milk proteins. Int. Dairy J. 1998;8: 325-31.21. Pihlanto-Leppala A, Korhonen H, Piilola K, Tupasela T, Coronen H. Angiotensin I-converting enzyme inhibitory properties of whey protein digests: concentration and characterization of active peptides J. Dairy Res. 2000;67: 53-64.22. Chougule RA, Aparna HS. Characterization of β-lactoglobulin from buffalo (Bubalus bubalis) colostrum and its possible interaction with erythrocyte    lipocalins-interacting membrane receptor. J. Biochem. 2011;150: 279-88.23. Rowsell S, Hawtin P, Minshull CA, Sarah HJ, Brockbank MV, Barratt DG, Slater AM, McPheat WL, Waterson D, Henney AM and Pauptit RA. Crystal structure of human MMP9 in complex with a reverse hydrox- amate inhibitor. J. Mol. Biol. 2002;319:173-81.24. Thomsen R, Chirstensen MHJ. A new tequnique for high-accuracy molecular docking. J. Med. Chem. 2006;49: 3315-21.25. Kolde HJ, Liberatori J, Braunitzer G . The amino acid sequence of the water buffalo beta-lactoglobulin. Milchwissenschaft. 1981;36:83-86.26. Das P, Jain S, Nayak S, Apparao KBC, Totey SM, Gerg LC. Molecular cloning and sequence analysis of the cDNA encoding beta-lactoglobulin in Bubalus bubalis. DNA Seq. 1999; 10:105-8.27. Yamamoto D, Takai S, Miyazaki M. Inhibitory profiles of captopril on matrix metalloproteinase-9 activity. Eur. J. Pharmcol. 2008;588: 277-79.