Active site specificity profiling of the matrix metalloproteinase family: Proteomic identification of 4300 cleavage sites by nine MMPs explored with structural and synthetic peptide cleavage analyses

Secreted and membrane tethered matrix metalloproteinases (MMPs) are key homeostatic proteases regulating the extracellular signaling and structural matrix environment of cells and tissues. For drug targeting of proteases, selectivity for individual molecules is highly desired and can be met by high yield active site specificity profiling. Using the high throughput Proteomic Identification of protease Cleavage Sites (PICS) method to simultaneously profile both the prime and non-prime sides of the cleavage sites of nine human MMPs, we identified more than 4300 cleavages from P6 to P6′ in biologically diverse human peptide libraries. MMP specificity and kinetic efficiency were mainly guided by aliphatic and aromatic residues in P1′ (with a ~ 32–93% preference for leucine depending on the MMP), and basic and small residues in P2′ and P3′, respectively. A wide differential preference for the hallmark P3 proline was found between MMPs ranging from 15 to 46%, yet when combined in the same peptide with the universally preferred P1′ leucine, an unexpected negative cooperativity emerged. This was not observed in previous studies, probably due to the paucity of approaches that profile both the prime and non-prime sides together, and the masking of subsite cooperativity effects by global heat maps and iceLogos. These caveats make it critical to check for these biologically highly important effects by fixing all 20 amino acids one-by-one in the respective subsites and thorough assessing of the inferred specificity logo changes. Indeed an analysis of bona fide MEROPS physiological substrate cleavage data revealed that of the 37 natural substrates with either a P3-Pro or a P1′-Leu only 5 shared both features, confirming the PICS data. Upon probing with several new quenched-fluorescent peptides, rationally designed on our specificity data, the negative cooperativity was explained by reduced non-prime side flexibility constraining accommodation of the rigidifying P3 proline with leucine locked in S1′. Similar negative cooperativity between P3 proline and the novel preference for asparagine in P1 cements our conclusion that non-prime side flexibility greatly impacts MMP binding affinity and cleavage efficiency. Thus, unexpected sequence cooperativity consequences were revealed by PICS that uniquely encompasses both the non-prime and prime sides flanking the proteomic-pinpointed scissile bond

Modulation of homocysteine metabolism and redox homeostasis via regulation of cystathionine beta -synthase

Homocysteine is a thiol-containing amino acid that is generated by the hydrolysis of S-adenosylhomocysteine in methylation metabolism. Elevated levels of homocysteine (10--100 μM) are associated with cardiovascular diseases, neural tube defects, Alzheimer\u27s disease. The transsulfuration pathway represents an important route for homocysteine detoxification and provides a limiting reagent for glutathione synthesis, cysteine, thus connecting transmethylation and redox metabolism. The first step in the transsulfuration pathway is catalyzed by cystathionine β-synthase (CBS). Despite its significance, the in vivo regulation at the homocysteine junction is not well understood. In this study, we have characterized the impact of regulating CBS activity on homocysteine and redox homeostasis. First, we investigated a potential mechanistic link between the transsulfuration pathway and the methionine dependence of cancer cells, which prevents growth when methionine is replaced by homocysteine. We found that methionine restriction in vitro resulted in a reversible \u3e10-fold downregulation of CBS levels, which was caused by destabilization of CBS in the absence of a methionine metabolite, AdoMet. Lower CBS levels were accompanied by reduced cell viability under oxidative stress conditions. The AdoMet-related CBS downregulation was also observed in vivo, in a murine model for steatohepatitis and in human hepatocarcinoma. These findings provide a mechanistic link between the coordinate changes in methylation (low AdoMet) and redox (low GSH) metabolism seen in liver disease. Further, we employed a mathematical model to assess the effects of allosteric regulation of CBS by AdoMet. The model predicts an increase in homocysteine levels upon the loss of allosteric activation, thus explaining hyperhomocysteinemia observed in CBS D444N patients. Finally, we studied the basis of a sexual dimorphism in CBS regulation in mouse kidneys. We found that in mice, testosterone upregulates CBS mRNA, protein and activity levels resulting in ∼2-fold higher CBS in males than females. Lower CBS in females is correlated with an elevated plasma homocysteine, which is further amplified by a high methionine/low folate diet. These data suggest that the kidney transsulfuration pathway plays a major role in plasma homocysteine clearance

Analysis of pathological defects in methionine metabolism using a simple mathematical model

AbstractDerangements in methionine metabolism are a hallmark of cancers and homocystinuria, an inborn error of metabolism. In this study, the metabolic consequences of the pathological changes associated with the key pathway enzymes, methionine adenosyl transferase (MAT), glycine N-methyl transferase (GNMT) and cystathionine β-synthase (CBS) as well as an activation of polyamine metabolism, were analyzed using a simple mathematical model describing methionine metabolism in liver. The model predicts that the mere loss of allosteric regulation of CBS by adenosylmethionine (AdoMet) leads to an increase in homocysteine concentration. This is consistent with the experimental data on the corresponding genetic defects, which specifically impair allosteric activation but not basal enzyme activity. Application of the characteristics of transformed hepatocytes to our model, i.e., substitution of the MATI/III isozyme by MATII, loss of GNMT activity and activation of polyamine biosynthesis, leads to the prediction of a significantly different dependence of methionine metabolism on methionine concentrations. The theoretical predictions were found to be in good agreement with experimental data obtained with the human hepatoma cell line, HepG2

Interactome disassembly during apoptosis occurs independent of caspase cleavage

Abstract Protein–protein interaction networks (interactomes) define the functionality of all biological systems. In apoptosis, proteolysis by caspases is thought to initiate disassembly of protein complexes and cell death. Here we used a quantitative proteomics approach, protein correlation profiling (PCP), to explore changes in cytoplasmic and mitochondrial interactomes in response to apoptosis initiation as a function of caspase activity. We measured the response to initiation of Fas‐mediated apoptosis in 17,991 interactions among 2,779 proteins, comprising the largest dynamic interactome to date. The majority of interactions were unaffected early in apoptosis, but multiple complexes containing known caspase targets were disassembled. Nonetheless, proteome‐wide analysis of proteolytic processing by terminal amine isotopic labeling of substrates (TAILS) revealed little correlation between proteolytic and interactome changes. Our findings show that, in apoptosis, significant interactome alterations occur before and independently of caspase activity. Thus, apoptosis initiation includes a tight program of interactome rearrangement, leading to disassembly of relatively few, select complexes. These early interactome alterations occur independently of cleavage of these protein by caspases

Deep profiling of the cleavage specificity and human substrates of snake venom metalloprotease HF3 by PICS using proteome derived peptide libraries and TAILS N-terminomics

Snakebite is a major medical concern in many parts of the world with metalloproteases playing important roles in the pathological effects of Viperidae venoms, including local tissue damage, hemorrhage, and coagulopathy. Hemorrhagic Factor 3 (HF3), a metalloprotease from Bothrops jararaca venom, induces local hemorrhage and targets extracellular matrix (ECM) components, including collagens and proteoglycans, and plasma proteins. However, the full substrate repertoire of this metalloprotease is unknown. We report positional proteomic studies identifying >2000 N-termini, including neo-N-termini of HF3 cleavage sites in mouse embryonic fibroblast secretome proteins. Terminal amine isotopic labeling of substrates (TAILS) analysis identified a preference for Leu at the P1′ position among candidate HF3 substrates including proteins of the ECM and focal adhesions and the cysteine protease inhibitor cystatin-C. Interestingly, 190 unique peptides matched to annotated cleavage sites in the TopFIND N-termini database, suggesting that these cleavages occurred at a site prone to cleavage or might have been generated by other proteases activated upon incubation with HF3, including caspases-3 and -7, cathepsins D and E, granzyme B, and MMPs 2 and 9. Using Proteomic identification of cleavage site specificity (PICS), a tryptic library derived from THP-1 monocytic cells was used as HF3 substrates for identifying protease cleavage sites and sequence preferences in peptides. A total of 799 unique cleavage sites were detected and, in accordance with TAILS analysis using native secreted protein substrates of MEF cells, revealed a clear preference for Leu at P1′. Taken together, these results greatly expand the known substrate degradome of HF3 and reveal potential new targets, which may serve as a basis to better elucidate the complex pathophysiology of snake envenomation

Transcriptomic and proteomic host response to Aspergillus fumigatus conidia in an air-liquid interface model of human bronchial epithelium.

Aspergillus fumigatus (A. fumigatus) is a wide-spread fungus that is a potent allergen in hypersensitive individuals but also an opportunistic pathogen in immunocompromised patients. It reproduces asexually by releasing airborne conidiospores (conidia). Upon inhalation, fungal conidia are capable of reaching the airway epithelial cells (AECs) in bronchial and alveolar tissues. Previous studies have predominantly used submerged monolayer cultures for studying this host-pathogen interaction; however, these cultures do not recapitulate the mucocililary differentiation phenotype of the in vivo epithelium in the respiratory tract. Thus, the aim of this study was to use well-differentiated primary human bronchial epithelial cells (HBECs) grown at the air-liquid interface (ALI) to determine their transcriptomic and proteomic responses following interaction with A. fumigatus conidia. We visualized conidial interaction with HBECs using confocal laser scanning microscopy (CLSM), and applied NanoString nCounter and shotgun proteomics to assess gene expression changes in the human cells upon interaction with A. fumigatus conidia. Western blot analysis was used to assess the expression of top three differentially expressed proteins, CALR, SET and NUCB2. CLSM showed that, unlike submerged monolayer cultures, well-differentiated ALI cultures of primary HBECs were estimated to internalize less than 1% of bound conidia. Nevertheless, transcriptomic and proteomic analyses revealed numerous differentially expressed host genes; these were enriched for pathways including apoptosis/autophagy, translation, unfolded protein response and cell cycle (up-regulated); complement and coagulation pathways, iron homeostasis, nonsense mediated decay and rRNA binding (down-regulated). CALR and SET were confirmed to be up-regulated in ALI cultures of primary HBECs upon exposure to A. fumigatus via western blot analysis. Therefore, using transcriptomics and proteomics approaches, ALI models recapitulating the bronchial epithelial barrier in the conductive zone of the respiratory tract can provide novel insights to the molecular response of bronchial epithelial cells upon exposure to A. fumigatus conidia

Active site specificity profiling datasets of matrix metalloproteinases (MMPs) 1, 2, 3, 7, 8, 9, 12, 13 and 14

The data described provide a comprehensive resource for the family-wide active site specificity portrayal of the human matrix metalloproteinase family. We used the high-throughput proteomic technique PICS (Proteomic Identification of protease Cleavage Sites) to comprehensively assay 9 different MMPs. We identified more than 4300 peptide cleavage sites, spanning both the prime and non-prime sides of the scissile peptide bond allowing detailed subsite cooperativity analysis. The proteomic cleavage data were expanded by kinetic analysis using a set of 6 quenched-fluorescent peptide substrates designed using these results. These datasets represent one of the largest specificity profiling efforts with subsequent structural follow up for any protease family and put the spotlight on the specificity similarities and differences of the MMP family. A detailed analysis of this data may be found in Eckhard et al. (2015) [1]. The raw mass spectrometry data and the corresponding metadata have been deposited in PRIDE/ProteomeXchange with the accession number http://www.ebi.ac.uk/pride/archive/projects/PXD002265

TAILS N-Terminomics and Proteomics Show Protein Degradation Dominates over Proteolytic Processing by Cathepsins in Pancreatic Tumors

Deregulated cathepsin proteolysis occurs across numerous cancers, but in vivo substrates mediating tumorigenesis remain ill-defined. Applying 8-plex iTRAQ terminal amine isotopic labeling of substrates (TAILS), a systems-level N-terminome degradomics approach, we identified cathepsin B, H, L, S, and Z in vivo substrates and cleavage sites with the use of six different cathepsin knockout genotypes in the Rip1-Tag2 mouse model of pancreatic neuroendocrine tumorigenesis. Among 1,935 proteins and 1,114 N termini identified by TAILS, stable proteolytic products were identified in wild-type tumors compared with one or more different cathepsin knockouts (17%–44% of 139 cleavages). This suggests a lack of compensation at the substrate level by other cathepsins. The majority of neo-N termini (56%–83%) for all cathepsins was consistent with protein degradation. We validated substrates, including the glycolytic enzyme pyruvate kinase M2 associated with the Warburg effect, the ER chaperone GRP78, and the oncoprotein prothymosin-alpha. Thus, the identification of cathepsin substrates in tumorigenesis improves the understanding of cathepsin functions in normal physiology and cancer