Substrate-driven mapping of the degradome by comparison of sequence logos

Predicted salivary human protease activity in experimental gingivitis revealed by endoProteo-FASP approach

Gingivitis is a highly prevalent oral condition that can be studied in humans via the 21-d experimental gingivitis model, which allows for investigations into the induction and resolution of gingivitis. In this study, we used the autolysis of saliva as a source of peptides to predict the activity of human proteases in saliva during induction and resolution of inflammation. Healthy volunteers, with no remarkable oral or systemic conditions, were recruited into the study and stimulated saliva samples were collected at days 0, 21, and 35 of experimental gingivitis. Plaque and gingival indices were recorded to ensure clinical induction and resolution. Saliva was auto-digested at 37°C for 18 h before identification of peptides by mass spectrometry. Protease prediction was carried out using Proteasix in silico with the identified peptides. A comparison of day 0 to days 21 and 35 showed changes in predicted protease activity. Correlation network analysis revealed that at day 21 the proteases became less connected and showed a potential for a dysregulated system; by day 35 the connectivity was returning towards similar conditions at day 0. This study demonstrates that changes in predicted proteases are apparent even in saliva collected from donors experiencing inflammation around three teeth.

Revealing favorable and unfavorable residues in cooperative positions in protease cleavage sites

Proteases play critical roles in a wide variety of fundamental biological functions, and numerous protease inhibitors have been developed to treat various diseases including cancer. A wide range of experimental and computational methods have been developed to investigate the specificity and catalytic mechanisms of proteases. However, these methods only focused on the preferences of a single position around a cleavage site in a substrate, rarely on the compositionality of the subsites. We present new methods to quantify the specificity of proteases by considering the combinatorial patterns of amino acid residuals of cleavage sites in substrates. By incorporating the preference at positions, we modeled three types of favorable combinations of residues in cleavage sites. Moreover, by constructing a relationship weight matrix of residues between two positions, we can easily identify unfavorable combinations of residues at the positions. Applying these methods to a set of known cleavage sites of proteases, we revealed numerous favorable and unfavorable residues in cooperative positions in the protease cleavage sites. The results can help understand the specificity and catalytic mechanisms of proteases. To our knowledge, this is the first study that quantifies unfavorable combinations of amino acids between two sites. Furthermore, this method is not limited to the study of proteases and cleavage sites, and can be generalized to uncover the relationships of residues at meaningful sites in other proteins.

Electrostatic recognition in substrate binding to serine proteases

Serine proteases of the Chymotrypsin family are structurally very similar but have very different substrate preferences. This study investigates a set of 9 different proteases of this family comprising proteases that prefer substrates containing positively charged amino acids, negatively charged amino acids, and uncharged amino acids with varying degree of specificity. Here, we show that differences in electrostatic substrate preferences can be predicted reliably by electrostatic molecular interaction fields employing customized GRID probes. Thus, we are able to directly link protease structures to their electrostatic substrate preferences. Additionally, we present a new metric that measures similarities in substrate preferences focusing only on electrostatics. It efficiently compares these electrostatic substrate preferences between different proteases. This new metric can be interpreted as the electrostatic part of our previously developed substrate similarity metric. Consequently, we suggest, that substrate recognition in terms of electrostatics and shape complementarity are rather orthogonal aspects of substrate recognition. This is in line with a 2‐step mechanism of protein‐protein recognition suggested in the literature.

Secretome profiling reveals temperature-dependent growth of Aspergillus fumigatus

Aspergillus fumigatus is a ubiquitous opportunistic fungus. In this study, systematic analyses were carried out to study the temperature adaptability of A. fumigatus. A total of 241 glycoside hydrolases and 69 proteases in the secretome revealed the strong capability of A. fumigatus to degrade plant biomass and protein substrates. In total, 129 pathogenesis-related proteins detected in the secretome were strongly correlated with glycoside hydrolases and proteases. The variety and abundance of proteins remained at temperatures of 34°C-45°C. The percentage of endo-1,4-xylanase increased when the temperature was lowered to 20°C, while the percentage of cellobiohydrolase increased as temperature was increased, suggesting that the strain obtains carbon mainly by degrading xylan and cellulose, and the main types of proteases in the secretome were aminopeptidases and carboxypeptidases. Only half of the proteins were retained and their abundance declined to 9.7% at 55°C. The activities of the remaining β-glycosidases and proteases were merely 35% and 24%, respectively, when the secretome was treated at 60°C for 2 h. Therefore, temperatures >60°C restrict the growth of A. fumigatus.

Block-based characterization of protease specificity from substrate sequence profile

Background

The mechanism of action of proteases has been widely studied based on substrate specificity. Prior research has been focused on the amino acids at a single amino acid site, but rarely on combinations of amino acids around the cleavage bond.

Results

We propose a novel block-based approach to reveal the potential combinations of amino acids which may regulate the action of proteases. Using the entropies of eight blocks centered at a cleavage bond, we created a distance matrix for 61 proteases to compare their specificities. After quantitative analysis, we discovered a number of prominent blocks, each of which consists of successive amino acids near a cleavage bond, intuitively characterizing the site cooperation of the substrate sequences.

Conclusion

This approach will help in the discovery of specific substrate sequences which may bridge between proteases and cleavage substrate as more substrate information becomes available.

Electronic supplementary material

The online version of this article (10.1186/s12859-017-1851-1) contains supplementary material, which is available to authorized users.

Molecular Connectivity Predefines Polypharmacology: Aliphatic Rings, Chirality, and sp3 Centers Enhance Target Selectivity

Dark chemical matter compounds are small molecules that have been recently identified as highly potent and selective hits. For this reason, they constitute a promising class of possible candidates in the process of drug discovery and raise the interest of the scientific community. To this purpose, Wassermann et al. (2015) have described the application of 2D descriptors to characterize dark chemical matter. However, their definition was based on the number of reported positive assays rather than the number of known targets. As there might be multiple assays for one single target, the number of assays does not fully describe target selectivity. Here, we propose an alternative classification of active molecules that is based on the number of known targets. We cluster molecules in four classes: black, gray, and white compounds are active on one, two to four, and more than four targets respectively, whilst inactive compounds are found to be inactive in the considered assays. In this study, black and inactive compounds are found to have not only higher solubility, but also a higher number of chiral centers, sp3 carbon atoms and aliphatic rings. On the contrary, white compounds contain a higher number of double bonds and fused aromatic rings. Therefore, the design of a screening compound library should consider these molecular properties in order to achieve target selectivity or polypharmacology. Furthermore, analysis of four main target classes (GPCRs, kinases, proteases, and ion channels) shows that GPCR ligands are more selective than the other classes, as the number of black compounds is higher in this target superfamily. On the other side, ligands that hit kinases, proteases, and ion channels bind to GPCRs more likely than to other target classes. Consequently, depending on the target protein family, appropriate screening libraries can be designed in order to minimize the likelihood of unwanted side effects early in the drug discovery process. Additionally, synergistic effects may be obtained by library design toward polypharmacology.

3D proteochemometrics: using three-dimensional information of proteins and ligands to address aspects of the selectivity of serine proteases

The high similarity between certain sub-pockets of serine proteases may lead to low selectivity of protease inhibitors. Therefore the application of proteochemometrics (PCM), which quantifies the relationship between protein/ligand descriptors and affinity for multiple ligands and targets simultaneously, is useful to understand and improve the selectivity profiles of potential inhibitors. In this study, protein field-based PCM that uses knowledge-based and WaterMap derived fields to describe proteins in combination with 2D (RDKit and MOE fingerprints) and 3D (4 point pharmacophoric fingerprints and GRIND) ligand descriptors was used to model the bioactivities of 24 homologous serine proteases and 5863 inhibitors in an integrated fashion. Of the multiple field-based PCM models generated based on different ligand descriptors, RDKit fingerprints showed the best performance in terms of external prediction with Rtest2 of 0.72 and RMSEP of 0.81. Further, visual interpretation of the models highlights sub-pocket specific regions that influence affinity and selectivity of serine protease inhibitors.

Interactome disassembly during apoptosis occurs independent of caspase cleavage

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.

Determinants of Macromolecular Specificity from Proteomics-Derived Peptide Substrate Data

BACKGROUND

Recent advances in proteomics methodologies allow for high throughput profiling of proteolytic cleavage events. The resulting substrate peptide distributions provide deep insights in the underlying macromolecular recognition events, as determinants of biomolecular specificity identified by proteomics approaches may be compared to structure-based analysis of corresponding protein-protein interfaces.

METHOD

Here, we present an overview of experimental and computational methodologies and tools applied in the area and provide an outlook beyond the protein class of proteases.

RESULTS AND CONCLUSION

We discuss here future potential, synergies and needs of the emerging overlap disciplines of proteomics and structure-based modelling.

Protease Inhibitors in View of Peptide Substrate Databases

Protease substrate profiling has nowadays almost become a routine task for experimentalists, and the knowledge on protease peptide substrates is easily accessible via the MEROPS database. We present a shape-based virtual screening workflow using vROCS that applies the information about the specificity of the proteases to find new small-molecule inhibitors. Peptide substrate sequences for three to four substrate positions of each substrate from the MEROPS database were used to build the training set. Two-dimensional substrate sequences were converted to three-dimensional conformations through mutation of a template peptide substrate. The vROCS query was built from single amino acid queries for each substrate position considering the relative frequencies of the amino acids. The peptide-substrate-based shape-based virtual screening approach gives good performance for the four proteases thrombin, factor Xa, factor VIIa, and caspase-3 with the DUD-E data set. The results show that the method works for protease targets with different specificity profiles as well as for targets with different active-site mechanisms. As no structure of the target and no information on small-molecule inhibitors are required to use our approach, the method has significant advantages in comparison with conventional structure- and ligand-based methods.

Identification of Protease Specificity by Combining Proteome-Derived Peptide Libraries and Quantitative Proteomics

We present protease specificity profiling based on quantitative proteomics in combination with proteome-derived peptide libraries. Peptide libraries are generated by endoproteolytic digestion of proteomes without chemical modification of primary amines before exposure to a protease under investigation. After incubation with a test protease, treated and control libraries are differentially isotope-labeled using cost-effective reductive dimethylation. Upon analysis by liquid chromatography-tandem mass spectrometry, cleavage products of the test protease appear as semi-specific peptides that are enriched for the corresponding isotope label. We validate our workflow with two proteases with well-characterized specificity profiles: trypsin and caspase-3. We provide the first specificity profile of a protease encoded by a human endogenous retrovirus and for chlamydial protease-like activity factor (CPAF). For CPAF, we also highlight the structural basis of negative subsite cooperativity between subsites S1 and S2'. For A disintegrin and metalloproteinase with thrombospondin motifs (ADAMTS) -4, -5, and -15, we show a canonical preference profile, including glutamate in P1 and glycine in P3'. In total, we report nearly 4000 cleavage sites for seven proteases. Our protocol is fast, avoids enrichment or synthesis steps, and enables probing for lysine selectivity as well as subsite cooperativity. Due to its simplicity, we anticipate usability by most proteomic laboratories.

Polypharmacology: in silico methods of ligand design and development

How to design a ligand to bind multiple targets, rather than to a single target, is the focus of this review. Rational polypharmacology draws on knowledge that is both broad ranging and hierarchical. Computer-aided multitarget ligand design methods are described according to their nested knowledge level. Ligand-only and then receptor-ligand strategies are first described; followed by the metabolic network viewpoint. Subsequently strategies that view infectious diseases as multigenomic targets are discussed, and finally the disease level interpretation of medicinal therapy is considered. As yet there is no consensus on how best to proceed in designing a multitarget ligand. The current methodologies are bought together in an attempt to give a practical overview of how polypharmacology design might be best initiated.

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 PXD002265.

Quantitative Correlation of Conformational Binding Enthalpy with Substrate Specificity of Serine Proteases

Members of the same protease family show different substrate specificity, even if they share identical folds, depending on the physiological processes they are part of. Here, we investigate the key factors for subpocket and global specificity of factor Xa, elastase, and granzyme B which despite all being serine proteases and sharing the chymotrypsin-fold show distinct substrate specificity profiles. We determined subpocket interaction potentials with GRID for static X-ray structures and an in silico generated ensemble of conformations. Subpocket interaction potentials determined for static X-ray structures turned out to be insufficient to explain serine protease specificity for all subpockets. Therefore, we generated conformational ensembles using molecular dynamics simulations. We identified representative binding site conformations using distance-based hierarchical agglomerative clustering and determined subpocket interaction potentials for each representative conformation of the binding site. Considering the differences in subpocket interaction potentials for these representative conformations as well as their abundance allowed us to quantitatively explain subpocket specificity for the nonprime side for all three example proteases on a molecular level. The methods to identify key regions determining subpocket specificity introduced in this study are directly applicable to other serine proteases, and the results provide starting points for new strategies in rational drug design.

Involvement of Kallikrein-Related Peptidases in Normal and Pathologic Processes

Human kallikrein-related peptidases (KLKs) are a subgroup of serine proteases that participate in proteolytic pathways and control protein levels in normal physiology as well as in several pathological conditions. Their complex network of stimulatory and inhibitory interactions may induce inflammatory and immune responses and contribute to the neoplastic phenotype through the regulation of several cellular processes, such as proliferation, survival, migration, and invasion. This family of proteases, which includes one of the most useful cancer biomarkers, kallikrein-related peptidase 3 or PSA, also has a protective effect against cancer promoting apoptosis or counteracting angiogenesis and cell proliferation. Therefore, they represent attractive therapeutic targets and may have important applications in clinical oncology. Despite being intensively studied, many gaps in our knowledge on several molecular aspects of KLK functions still exist. This review aims to summarize recent data on their involvement in different processes related to health and disease, in particular those directly or indirectly linked to the neoplastic process.

Characterizing Protease Specificity: How Many Substrates Do We Need?

Calculation of cleavage entropies allows to quantify, map and compare protease substrate specificity by an information entropy based approach. The metric intrinsically depends on the number of experimentally determined substrates (data points). Thus a statistical analysis of its numerical stability is crucial to estimate the systematic error made by estimating specificity based on a limited number of substrates. In this contribution, we show the mathematical basis for estimating the uncertainty in cleavage entropies. Sets of cleavage entropies are calculated using experimental cleavage data and modeled extreme cases. By analyzing the underlying mathematics and applying statistical tools, a linear dependence of the metric in respect to 1/n was found. This allows us to extrapolate the values to an infinite number of samples and to estimate the errors. Analyzing the errors, a minimum number of 30 substrates was found to be necessary to characterize substrate specificity, in terms of amino acid variability, for a protease (S4-S4’) with an uncertainty of 5 percent. Therefore, we encourage experimental researchers in the protease field to record specificity profiles of novel proteases aiming to identify at least 30 peptide substrates of maximum sequence diversity. We expect a full characterization of protease specificity helpful to rationalize biological functions of proteases and to assist rational drug design.

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.