Converting trypsin to chymotrypsin: the role of surface loops

The High-Affinity Chymotrypsin Inhibitor Eglin C Poorly Inhibits Human Chymotrypsin-Like Protease: Gln192 and Lys218 Are Key Determinants

Eglin C, a small protein from the medicinal leech, has been long considered a general high-affinity inhibitor of chymotrypsins and elastases. Here, we demonstrate that eglin C inhibits human chymotrypsin-like protease (CTRL) weaker by several orders of magnitude than other chymotrypsins. In order to identify the underlying structural aspects of this unique deviation, we performed comparative molecular dynamics simulations on experimental and AlphaFold model structures of bovine CTRA and human CTRL. Our results indicate that in CTRL, the primary determinants of the observed weak inhibition are amino-acid positions 192 and 218 (using conventional chymotrypsin numbering), which participate in shaping the S1 substrate-binding pocket and thereby affect the stability of the protease-inhibitor complexes.

Replacement of a single residue changes the primary specificity of thrombin

BACKGROUND

Thrombin prefers substrates carrying Arg at the site of cleavage (P1) because of the presence of D189 in the primary specificity (S1) pocket but can also cleave substrates carrying Phe at P1. The structural basis of this property is unknown.

OBJECTIVES

Solve the X-ray structure of thrombin bound to a ligand carrying Phe at P1 and investigate the effects of replacing D189.

METHODS

X-ray crystallography is used to solve the structure of thrombin bound to the irreversible inhibitor H-D-Phe-Pro-Phe-CH2Cl (PPPCK). Residue D189 is mutated to Ala, Lys, Phe, and Ser.

RESULTS

The X-ray structure of the thrombin-PPPCK complex is solved at 2.5 Å resolution and compared to the structure of thrombin bound to H-D-Phe-Pro-Arg-CH2Cl (PPACK). PPPCK binds to thrombin in a conformation similar to that of PPACK, but Phe at P1 makes no contacts with D189. Replacement of D189 with Ala, Lys, Phe, or Ser reverses both substrate preference and stability enhancement from Arg to Phe.

CONCLUSION

D189 in the S1 pocket confers thrombin "trypsin-like" specificity for Arg at P1. However, the S1 pocket is wide enough to also enable "chymotrypsin-like" specificity for Phe at P1. Consistent with these structural features, a single amino acid replacement (D189A) switches thrombin specificity from trypsin-like to chymotrypsin-like, converting the substrate preference from H-D-Phe-Pro-Arg-p-nitroanilide to H-D-Phe-Pro-Phe-p-nitroanilide and preferential stability enhancement from PPACK to PPPCK. The observation that thrombin specificity is controlled mainly by a single residue establishes a new paradigm in the field of trypsin-like proteases.

Comprehensive Understanding of Laboratory Evolution for Food Enzymes: From Design to Screening Innovations

In the field of food processing, enzymes play a pivotal role in improving product quality and flavor, and extending shelf life. However, the exposure of traditional food enzymes to high temperatures during processing often leads to a decrease in activity or even inactivation, limiting the effectiveness of their application under high-temperature conditions. Therefore, the modification of thermostability and activity of enzymes to adapt to extreme conditions through protein engineering has become a key way to improve the efficiency and economic benefits of industrial production. Directed evolution and semirational design strategies in the laboratory have proven to be broadly applicable frameworks for biochemical researchers in the food field, including those who are beginners. In this review, we systematically summarize semirational design strategies and high-throughput screening strategies, and introduce some intuitive computer simulation software to improve the thermostability and enzyme activity of food enzymes. The application of these strategies and techniques provides a comprehensive guide for the optimization of food enzymes. In addition, the latest hot topics of genetically engineered food enzymes in the field of application are discussed.

Leveraging High-Resolution Ion Mobility-Mass Spectrometry for Cyclic Peptide Soft Spot Identification

Cyclic peptides are an important class of molecules that gained significant attention in the field of drug discovery due to their unique pharmacological characteristics and enhanced proteolytic stability. Yet, gastrointestinal degradation remains a major hurdle in the discovery of orally bioavailable cyclic peptides. Soft spot identification (SSID) of the regions in the cyclic peptide sequence susceptible to amide hydrolysis by proteases is used in the discovery stage to guide medicinal chemistry design. SSID can be an arduous task, traditionally performed using liquid chromatography-tandem mass spectrometry (LC-MS/MS), often resulting in complex and time-consuming manual analysis, particularly when isomeric linear peptide metabolites chromatographically coelute. Here, we present an alternative orthogonal approach that entails a high-resolution ion mobility (HRIM) system based on Structures for Lossless Ion Manipulation (SLIM) technology interfaced with quadrupole time-of-flight (QTOF) mass spectrometry to address some of the challenges associated with SSID. Two strategies were used to resolve linear isomeric peptide metabolites: labeled and label-free, both utilizing the HRIM platform. The label-free strategy leverages negative polarity to ionize the isomers which achieves better separation of the gas phase ions in the ion mobility (IM) dimension as compared to positive polarity, which is a more conventional approach when studying proteins and peptides. The second approach uses an isotope-labeled dimethyl tag on the terminal amine group, acting as a "shift reagent" to influence the mobility of isomers in the positive mode. This method resulted in baseline separation for the isomers of interest and produced unique product ions in the fragmentation spectra for unambiguous soft spot identification. Both label-free and labeled strategies demonstrated the ability to solve the challenges associated with SSID for cyclic peptides.

Let's get biophysical - How to get your favorite protein's digits

In these days of information overload and high-throughput analysis, it is easy to lose focus on the study of individual proteins. It is our conjecture that such investigations are still crucially important and offer uniquely penetrative insights. We thus present a discussion of biophysical methods to allow readers to get to know their protein of interest better. Although this perspective is not written with the expert in mind, we hope that for interdisciplinary scientists, or researchers who do not routinely perform biophysical analyses, the content will be helpful and inspiring.

Remote loop evolution reveals a complex biological function for chitinase enzymes beyond the active site

Loops are small secondary structural elements that play a crucial role in the emergence of new enzyme functions. However, the evolutionary molecular mechanisms how proteins acquire these loop elements and obtain new function is poorly understood. To address this question, we study glycoside hydrolase family 19 (GH19) chitinase-an essential enzyme family for pathogen degradation in plants. By revealing the evolutionary history and loops appearance of GH19 chitinase, we discover that one loop which is remote from the catalytic site, is necessary to acquire the new antifungal activity. We demonstrate that this remote loop directly accesses the fungal cell wall, and surprisingly, it needs to adopt a defined structure supported by long-range intramolecular interactions to perform its function. Our findings prove that nature applies this strategy at the molecular level to achieve a complex biological function while maintaining the original activity in the catalytic pocket, suggesting an alternative way to design new enzyme function.

Thrombin has dual trypsin-like and chymotrypsin-like specificity

BACKGROUND

The residue at the site of activation of protein C is Arg in all species except the ray-finned fish, where it is Trp. This feature raises the question of whether thrombin is the physiological activator of protein C across vertebrates.

OBJECTIVES

To establish if thrombin can cleave at Trp residues.

METHODS

The activity of wild-type thrombin and mutant D189S was tested with a library of chromogenic substrates and toward wild-type protein C and mutants carrying substitutions at the site of cleavage.

RESULTS

Thrombin has trypsin-like and chymotrypsin-like specificity and cleaves substrates at Arg or Trp residues. Cleavage at Arg is preferred, but cleavage at Trp is significant and comparable with that of chymotrypsin. The D189S mutant of thrombin has broad specificity and cleaves at basic and aromatic residues without significant preference. Thrombin also cleaves natural substrates at Arg or Trp residues, showing activity toward protein C across vertebrates, including the ray-finned fish. The rate of activation of protein C in the ray-finned fish is affected by the sequence preceding Trp at the scissile bond.

CONCLUSION

The results provide a possible solution for the paradoxical presence of a Trp residue at the site of cleavage of protein C in ray-finned fish and support thrombin as the physiological activator of protein C in all vertebrates. The dual trypsin-like and chymotrypsin-like specificity of thrombin suggests that the spectrum of physiological substrates of this enzyme is broader currently assumed.

The Impact of Second-Shell Nucleotides on Ligand Specificity in Cyclic Dinucleotide Riboswitches

Ligand specificity is an essential requirement for all riboswitches. Some variant riboswitches utilize a common structural motif, yet through subtle sequence differences, they are able to selectively respond to different small molecule ligands and regulate downstream gene expression. These variants discriminate between structurally and chemically similar ligands. Crystal structures provide insight into how specificity is achieved. However, ligand specificity cannot always be explained solely by nucleotides in direct contact with the ligand. The cyclic dinucleotide variant family contains two classes, cyclic-di-GMP and cyclic-AMP-GMP riboswitches, that were distinguished based on the identity of a single nucleotide in contact with the ligand. Here we report a variant riboswitch with a mutation at a second ligand-contacting position that is promiscuous for both cyclic-di-GMP and cyclic-AMP-GMP despite a predicted preference for cyclic-AMP-GMP. A high-throughput mutational analysis, SMARTT, was used to quantitatively assess thousands of sites in the first- and second-shells of ligand contact for impacts on ligand specificity and promiscuity. In addition to nucleotides in direct ligand contact, nucleotides more distal from the binding site, within the J1/2 linker and the terminator helix, were identified that impact ligand specificity. These findings provide an example of how nucleotides outside the ligand binding pocket influence the riboswitch specificity. Moreover, these distal nucleotides could be used to predict promiscuous sequences.

Sequence - dynamics - function relationships in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) are crucial regulators of cellular signaling. Their activity is regulated by the motion of a conserved loop, the WPD-loop, from a catalytically inactive open to a catalytically active closed conformation. WPD-loop motion optimally positions a catalytically critical residue into the active site, and is directly linked to the turnover number of these enzymes. Crystal structures of chimeric PTPs constructed by grafting parts of the WPD-loop sequence of PTP1B onto the scaffold of YopH showed WPD-loops in a wide-open conformation never previously observed in either parent enzyme. This wide-open conformation has, however, been observed upon binding of small molecule inhibitors to other PTPs, suggesting the potential of targeting it for drug discovery efforts. Here, we have performed simulations of both enzymes and show that there are negligible energetic differences in the chemical step of catalysis, but significant differences in the dynamical properties of the WPD-loop. Detailed interaction network analysis provides insight into the molecular basis for this population shift to a wide-open conformation. Taken together, our study provides insight into the links between loop dynamics and chemistry in these YopH variants specifically, and how WPD-loop dynamic can be engineered through modification of the internal protein interaction network.

Engineering of TEV protease variants with redesigned substrate specificity

Due to their ability to catalytically cleave proteins and peptides, proteases present unique opportunities for the use in industrial, biotechnological, and therapeutic applications. Engineered proteases with redesigned substrate specificities have the potential to expand the scope of practical applications of this enzyme class. We here apply a combinatorial protease engineering-based screening method that links proteolytic activity to the solubility and correct folding of a fluorescent reporter protein to redesign the substrate specificity of tobacco etch virus (TEV) protease. The target substrate EKLVFQA differs at three out of seven positions from the TEV consensus substrate sequence. Flow cytometric sorting of a semi-rational TEV protease library, consisting of focused mutations of the substrate binding pockets as well as random mutations throughout the enzyme, led to the enrichment of a set of protease variants that recognize and cleave the novel target substrate.

Loop dynamics and the evolution of enzyme activity

In the early 2000s, Tawfik presented his 'New View' on enzyme evolution, highlighting the role of conformational plasticity in expanding the functional diversity of limited repertoires of sequences. This view is gaining increasing traction with increasing evidence of the importance of conformational dynamics in both natural and laboratory evolution of enzymes. The past years have seen several elegant examples of harnessing conformational (particularly loop) dynamics to successfully manipulate protein function. This Review revisits flexible loops as critical participants in regulating enzyme activity. We showcase several systems of particular interest: triosephosphate isomerase barrel proteins, protein tyrosine phosphatases and β-lactamases, while briefly discussing other systems in which loop dynamics are important for selectivity and turnover. We then discuss the implications for engineering, presenting examples of successful loop manipulation in either improving catalytic efficiency, or changing selectivity completely. Overall, it is becoming clearer that mimicking nature by manipulating the conformational dynamics of key protein loops is a powerful method of tailoring enzyme activity, without needing to target active-site residues.

Bioinspired Framework Catalysts: From Enzyme Immobilization to Biomimetic Catalysis

Enzymatic catalysis has fueled considerable interest from chemists due to its high efficiency and selectivity. However, the structural complexity and vulnerability hamper the application potentials of enzymes. Driven by the practical demand for chemical conversion, there is a long-sought quest for bioinspired catalysts reproducing and even surpassing the functions of natural enzymes. As nanoporous materials with high surface areas and crystallinity, metal-organic frameworks (MOFs) represent an exquisite case of how natural enzymes and their active sites are integrated into porous solids, affording bioinspired heterogeneous catalysts with superior stability and customizable structures. In this review, we comprehensively summarize the advances of bioinspired MOFs for catalysis, discuss the design principle of various MOF-based catalysts, such as MOF-enzyme composites and MOFs embedded with active sites, and explore the utility of these catalysts in different reactions. The advantages of MOFs as enzyme mimetics are also highlighted, including confinement, templating effects, and functionality, in comparison with homogeneous supramolecular catalysts. A perspective is provided to discuss potential solutions addressing current challenges in MOF catalysis.

Structure-activity relationships of dihydropyrimidone inhibitors against native and auto-processed human neutrophil elastase

BACKGROUND

Human neutrophil elastase (HNE) is a key driver of systemic and cardiopulmonary inflammation. Recent studies have established the existence of a pathologically active auto-processed form of HNE with reduced binding affinity against small molecule inhibitors.

METHOD

AutoDock Vina v1.2.0 and Cresset Forge v10 software were used to develop a 3D-QSAR model for a series of 47 DHPI inhibitors. Molecular Dynamics (MD) simulations were carried out using AMBER v18 to study the structure and dynamics of sc (single-chain HNE) and tcHNE (two-chain HNE). MMPBSA binding free energies of the previously reported clinical candidate BAY 85-8501 and the highly active BAY-8040 were calculated with sc and tcHNE.

RESULTS

The DHPI inhibitors occupy the S1 and S2 subsites of scHNE. The robust 3D-QSAR model showed acceptable predictive and descriptive capability with regression coefficient of r² = 0.995 and cross-validation regression coefficient q² = 0.579 for the training set. The key descriptors of shape, hydrophobics and electrostatics were mapped to the inhibitory activity. In auto-processed tcHNE, the S1 subsite undergoes widening and disruption. All the DHPI inhibitors docked with the broadened S1'-S2' subsites of tcHNE with lower AutoDock binding affinities. The MMPBSA binding free energy of BAY-8040 with tcHNE reduced in comparison with scHNE while the clinical candidate BAY 85-8501 dissociated during MD. Thus, BAY-8040 may have lower inhibitory activity against tcHNE whereas the clinical candidate BAY 85-8501 is likely to be inactive.

CONCLUSION

SAR insights gained from this study will aid the future development of inhibitors active against both forms of HNE.

Insights into the importance of WPD-loop sequence for activity and structure in protein tyrosine phosphatases

Protein tyrosine phosphatases (PTPs) possess a conserved mobile catalytic loop, the WPD-loop, which brings an aspartic acid into the active site where it acts as an acid/base catalyst. Prior experimental and computational studies, focused on the human enzyme PTP1B and the PTP from Yersinia pestis, YopH, suggested that loop conformational dynamics are important in regulating both catalysis and evolvability. We have generated a chimeric protein in which the WPD-loop of YopH is transposed into PTP1B, and eight chimeras that systematically restored the loop sequence back to native PTP1B. Of these, four chimeras were soluble and were subjected to detailed biochemical and structural characterization, and a computational analysis of their WPD-loop dynamics. The chimeras maintain backbone structural integrity, with somewhat slower rates than either wild-type parent, and show differences in the pH dependency of catalysis, and changes in the effect of Mg2+. The chimeric proteins' WPD-loops differ significantly in their relative stability and rigidity. The time required for interconversion, coupled with electrostatic effects revealed by simulations, likely accounts for the activity differences between chimeras, and relative to the native enzymes. Our results further the understanding of connections between enzyme activity and the dynamics of catalytically important groups, particularly the effects of non-catalytic residues on key conformational equilibria.

Rational Design of a Calcium-Independent Trypsin Variant

Trypsin is a long-known serine protease widely used in biochemical, analytical, biotechnological, or biocatalytic applications. The high biotechnological potential is based on its high catalytic activity, substrate specificity, and catalytic robustness in non-physiological reaction conditions. The latter is mainly due to its stable protein fold, to which six intramolecular disulfide bridges make a significant contribution. Although trypsin does not depend on cofactors, it essentially requires the binding of calcium ions to its calcium-binding site to obtain complete enzymatic activity and stability. This behavior is inevitably associated with a limitation of the enzyme’s applicability. To make trypsin intrinsically calcium-independent, we removed the native calcium-binding site and replaced it with another disulfide bridge. The resulting stabilized apo-trypsin (aTn) retains full catalytic activity as proven by enzyme kinetics. Studies using Ellmann’s reagent further prove that the two inserted cysteines at positions Glu70 and Glu80 are in their oxidized state, creating the desired functional disulfide bond. Furthermore, aTn is independent of calcium ions, possesses increased thermal and functional stability, and significantly reduced autolysis compared to wildtype trypsin. Finally, we confirmed our experimental data by solving the X-ray crystal structure of aTn.

Modulation of global stability, ligand binding and catalytic properties of trypsin by anions

Specific salts effect is well-known on stability and solubility of proteins, however, relatively limited knowledge is known regarding the effect on catalytic properties of enzymes. Here, we examined the effect of four sodium anions on thermal stability and catalytic properties of trypsin and binding of the fluorescent probe, p-aminobenzamidine (PAB), to the enzyme. We show that the specific anions effect on trypsin properties agrees with the localization of the anions in the Hofmeister series. Thermal stability of trypsin, T_m, the affinity of the fluorescent probe to the binding site, K_d, and the rate constant, k_cat, of trypsin-catalyzed hydrolysis of the substrate N-benzoyl-L-arginine ethyl ester (BAEE) increase with increasing kosmotropic character of anions in the order: perchlorate<bromide<chloride<sulfate, while the value of Michaelis constant, K_M, decreases. Correlations between the values of T_m, K_d for PAB, k_cat, and K_M for BAEE in the presence of 1 M studied salts suggest interrelation among these parameters of the enzyme. Global stabilization as well as increased rigidity of trypsin is accompanied by strengthening of interaction with fluorescent probe PAB and in accordance with decreasing values of K_M for the substrate BAEE. Strong correlations between parameters characterizing the trypsin properties with the charge densities of anions clearly indicate direct electrostatic interaction as a basis of the specific anion effect on the conformational and functional properties of the enzyme.

A systematic approach for evaluating the role of surface-exposed loops in trypsin-like serine proteases applied to the 170 loop in coagulation factor VIIa

Proteases play a major role in many vital physiological processes. Trypsin-like serine proteases (TLPs), in particular, are paramount in proteolytic cascade systems such as blood coagulation and complement activation. The structural topology of TLPs is highly conserved, with the trypsin fold comprising two β-barrels connected by a number of variable surface-exposed loops that provide a surprising capacity for functional diversity and substrate specificity. To expand our understanding of the roles these loops play in substrate and co-factor interactions, we employ a systematic methodology akin to the natural truncations and insertions observed through evolution of TLPs. The approach explores a larger deletion space than classical random or directed mutagenesis. Using FVIIa as a model system, deletions of 1–7 amino acids through the surface exposed 170 loop, a vital allosteric regulator, was introduced. All variants were extensively evaluated by established functional assays and computational loop modelling with Rosetta. The approach revealed detailed structural and functional insights recapitulation and expanding on the main findings in relation to 170 loop functions elucidated over several decades using more cumbersome crystallization and single deletion/mutation methodologies. The larger deletion space was key in capturing the most active variant, which unexpectedly had a six-amino acid truncation. This variant would have remained undiscovered if only 2–3 deletions were considered, supporting the usefulness of the methodology in general protease engineering approaches. Our findings shed further light on the complex role that surface-exposed loops play in TLP function and supports the important role of loop length in the regulation and fine-tunning of enzymatic function throughout evolution.

Prospection of putative genes for digestive enzymes based on functional genome of the hepatopancreas of Amazon river prawn

Over recent years, Macrobrachium amazonicum has become a popular species for shrimp farming due to their fast growth, high survival rates, and marketability. Several studies have focused on the development of new technology for the culture of this species, but many aspects of their nutrition and physiology remain unknown. Thus, the goal of the present study was to obtain transcripts of putative genes encoding digestive enzymes, based on a library of the cDNA from the hepatopancreas of M. amazonicum, sequenced in the Ion TorrentTM platform. We identified fragments of nine genes related to digestive enzymes, acting over proteins, carbohydrates and lipids. Endo and exoproteases were also recorded in the hepatopancreas, indicating adaptation to the digestion of protein-rich foods. Nonetheless, the enzymes involved in the carbohydrate metabolism formed the largest functional group in M. amazonicum, including enzymes related to the digestion of starch, chitin, and cellulose. These findings indicate that the species has a genetic apparatus of a well-adapted omnivorous animal. This information may provide important insights for the selection of ingredients for the formulation of a more appropriate diet to the enzymatic repertoire of M. amazonicum.

Structural and biochemical analyses of selectivity determinants in chimeric Streptococcus Class A sortase enzymes

Sequence variation in related proteins is an important characteristic that modulates activity and selectivity. An example of a protein family with a large degree of sequence variation is that of bacterial sortases, which are cysteine transpeptidases on the surface of gram-positive bacteria. Class A sortases are responsible for attachment of diverse proteins to the cell wall to facilitate environmental adaption and interaction. These enzymes are also used in protein engineering applications for sortase-mediated ligations (SML) or sortagging of protein targets. We previously investigated SrtA from Streptococcus pneumoniae, identifying a number of putative β7-β8 loop-mediated interactions that affected in vitro enzyme function. We identified residues that contributed to the ability of S. pneumoniae SrtA to recognize several amino acids at the P1' position of the substrate motif, underlined in LPXTG, in contrast to the strict P1' Gly recognition of SrtA from Staphylococcus aureus. However, motivated by the lack of a structural model for the active, monomeric form of S. pneumoniae SrtA, here, we expanded our studies to other Streptococcus SrtA proteins. We solved the first monomeric structure of S. agalactiae SrtA which includes the C-terminus, and three others of β7-β8 loop chimeras from S. pyogenes and S. agalactiae SrtA. These structures and accompanying biochemical data support our previously identified β7-β8 loop-mediated interactions and provide additional insight into their role in Class A sortase substrate selectivity. A greater understanding of individual SrtA sequence and structural determinants of target selectivity may also facilitate the design or discovery of improved sortagging tools.

The active site region plays a critical role in Na+ binding to thrombin

The catalytic activity of thrombin and other enzymes of the blood coagulation and complement cascades is enhanced significantly by binding of Na+ to a site >15 Å away from the catalytic residue S195, buried within the 180 and 220 loops that also contribute to the primary specificity of the enzyme. Rapid kinetics support a binding mechanism of conformational selection where the Na+-binding site is in equilibrium between open (N) and closed (N∗) forms and the cation binds selectively to the N form. Allosteric transduction of this binding step produces enhanced catalytic activity. Molecular details on how Na+ gains access to this site and communicates allosterically with the active site remain poorly defined. In this study, we show that the rate of the N∗→N transition is strongly correlated with the analogous E∗→E transition that governs the interaction of synthetic and physiologic substrates with the active site. This correlation supports the active site as the likely point of entry for Na+ to its binding site. Mutagenesis and structural data rule out an alternative path through the pore defined by the 180 and 220 loops. We suggest that the active site communicates allosterically with the Na+ site through a network of H-bonded water molecules that embeds the primary specificity pocket. Perturbation of the mobility of S195 and its H-bonding capabilities alters interaction with this network and influences the kinetics of Na+ binding and allosteric transduction. These findings have general mechanistic relevance for Na+-activated proteases and allosteric enzymes.

Structure-Based Identification of Potent Lysine-Specific Demethylase 1 Inhibitor Peptides and Temporary Cyclization to Enhance Proteolytic Stability and Cell Growth-Inhibitory Activity

Peptides are attractive drug candidates, but their utility is greatly limited by their inherent susceptibility to proteolytic degradation and their inability to pass through the cell membrane. Here, we employ a strategy of temporary cyclization to develop a cell-active lysine-specific demethylase 1 (LSD1/KDM1A) inhibitor peptide. We first identified a highly potent LSD1-inhibitory linear peptide, with the assistance of X-ray crystal structure data of inhibitor peptide-bound LSD1·CoREST. The peptide was converted to a redox-activatable cyclic peptide incorporating cell-penetrating peptide (CPP), expecting selective activation under intracellular reducing conditions. The cyclic peptide moiety exhibited enhanced stability to protease and was converted to the linear, unmodified LSD1 inhibitor peptide under reducing conditions. The cyclic peptide with CPP inhibited the proliferation of human acute myeloid leukemia cells (HL-60) in the low micromolar concentration range.

The Picornavirus Precursor 3CD Has Different Conformational Dynamics Compared to 3Cpro and 3Dpol in Functionally Relevant Regions

Viruses have evolved numerous strategies to maximize the use of their limited genetic material, including proteolytic cleavage of polyproteins to yield products with different functions. The poliovirus polyprotein 3CD is involved in important protein-protein, protein-RNA and protein-lipid interactions in viral replication and infection. It is a precursor to the 3C protease and 3D RNA-dependent RNA polymerase, but has different protease specificity, is not an active polymerase, and participates in other interactions differently than its processed products. These functional differences are poorly explained by the known X-ray crystal structures. It has been proposed that functional differences might be due to differences in conformational dynamics between 3C, 3D and 3CD. To address this possibility, we conducted nuclear magnetic resonance spectroscopy experiments, including multiple quantum relaxation dispersion, chemical exchange saturation transfer and methyl spin-spin relaxation, to probe conformational dynamics across multiple timescales. Indeed, these studies identified differences in conformational dynamics in functionally important regions, including enzyme active sites, and RNA and lipid binding sites. Expansion of the conformational ensemble available to 3CD may allow it to perform additional functions not observed in 3C and 3D alone despite having nearly identical lowest-energy structures.

Origin and Expansion of the Serine Protease Repertoire in the Myelomonocyte Lineage

The deepest evolutionary branches of the trypsin/chymotrypsin family of serine proteases are represented by the digestive enzymes of the gastrointestinal tract and the multi-domain proteases of the blood coagulation and complement system. Similar to the very old digestive system, highly diverse cleavage specificities emerged in various cell lineages of the immune defense system during vertebrate evolution. The four neutrophil serine proteases (NSPs) expressed in the myelomonocyte lineage, neutrophil elastase, proteinase 3, cathepsin G, and neutrophil serine protease 4, collectively display a broad repertoire of (S1) specificities. The origin of NSPs can be traced back to a circulating liver-derived trypsin-like protease, the complement factor D ancestor, whose activity is tightly controlled by substrate-induced activation and TNFα-induced locally upregulated protein secretion. However, the present-day descendants are produced and converted to mature enzymes in precursor cells of the bone marrow and are safely sequestered in granules of circulating neutrophils. The potential site and duration of action of these cell-associated serine proteases are tightly controlled by the recruitment and activation of neutrophils, by stimulus-dependent regulated secretion of the granules, and by various soluble inhibitors in plasma, interstitial fluids, and in the inflammatory exudate. An extraordinary dynamic range and acceleration of immediate defense responses have been achieved by exploiting the high structural plasticity of the trypsin fold.

Characterization of the ExoU activation mechanism using EPR and integrative modeling

ExoU, a type III secreted phospholipase effector of Pseudomonas aeruginosa, serves as a prototype to model large, dynamic, membrane-associated proteins. ExoU is synergistically activated by interactions with membrane lipids and ubiquitin. To dissect the activation mechanism, structural homology was used to identify an unstructured loop of approximately 20 residues in the ExoU amino acid sequence. Mutational analyses indicate the importance of specific loop amino acid residues in mediating catalytic activity. Engineered disulfide cross-links show that loop movement is required for activation. Site directed spin labeling EPR and DEER (double electron-electron resonance) studies of apo and holo states demonstrate local conformational changes at specific sites within the loop and a conformational shift of the loop during activation. These data are consistent with the formation of a substrate-binding pocket providing access to the catalytic site. DEER distance distributions were used as constraints in RosettaDEER to construct ensemble models of the loop in both apo and holo states, significantly extending the range for modeling a conformationally dynamic loop.

Multiplexed Enzyme Activity-Based Probe Display via Hybridization

Emulsions offer the means to miniaturize and parallelize high-throughput screening but require a robust method to localize activity-based fluorescent probes in each droplet. Multiplexing probes in droplets is impractical, though highly desirable for identifying library members that possess very specific activity. Here, we present multiplexed probe immobilization on library beads for emulsion screening. During library bead preparation, we quantitated ∼10⁶ primers per bead by fluorescence in situ hybridization, however emulsion PCR yielded only ∼10³ gene copies per bead. We leveraged the unextended bead-bound primers to hybridize complementary probe-oligonucleotide heteroconjugates to the library beads. The probe-hybridized bead libraries were then used to program emulsion in vitro transcription/translation reactions and analyzed by FACS to perform multiplexed activity-based screening of trypsin and chymotrypsin mutant libraries for novel proteolytic specificity. The approach's modularity should permit a high degree of probe multiplexing and appears extensible to other enzyme classes and library types.

Remodeling enzyme active sites by stepwise loop insertion

The remolding active site loops via residue insertion/deletion as well as substitution is thought to play a key role in enzyme divergent evolution. However, enzyme engineering by residue insertion in active site loops often severely perturbs the protein structural integrity and causes protein misfolding and activity loss. We have designed a stepwise loop insertion strategy (StLois), in which a pair of randomized residues is introduced in a stepwise manner, efficiently collating mutational fitness effects. The strategy of StLois constitutes three key steps. First, the target regions should be identified through structural and functional analysis on the counterpart enzymes. Second, pair residues can be introduced in loop regions through insertion with NNK codon degeneracy. Third, the best hit used as a template for the next round mutagenesis. The residue insertion process can repeat as many times as necessary. By using the StLois method, we have evolved the substrate preference of a lactonase to phosphotriesterase. In this chapter, we describe the detailed StLois technique, which efficiently expands the residue in the loop region and remolds the architecture of enzyme active site for novel catalytic properties.

Engineered ChymotrypsiN for Mass Spectrometry-Based Detection of Protein Glycosylation

We have engineered the substrate specificity of chymotrypsin to cleave after Asn by high-throughput screening of large libraries created by comprehensive remodeling of the substrate binding pocket. The engineered variant (chymotrypsiN, ChyB-Asn) demonstrated an altered substrate specificity with an expanded preference for Asn-containing substrates. We confirmed that protein engineering did not compromise the stability of the enzyme by biophysical characterization. Comparison of wild-type ChyB and ChyB-Asn in profiling lysates of HEK293 cells demonstrated both qualitative and quantitative differences in the nature of the peptides and proteins identified by liquid chromatography and tandem mass spectrometry. ChyB-Asn enabled the identification of partially glycosylated Asn sites within a model glycoprotein and in the extracellular proteome of Jurkat T cells. ChymotrypsiN is a valuable addition to the toolkit of proteases to aid the mapping of N-linked glycosylation sites within proteins and proteomes.

Role of the I16-D194 ionic interaction in the trypsin fold

Activity in trypsin-like proteases is the result of proteolytic cleavage at R15 followed by an ionic interaction that ensues between the new N terminus of I16 and the side chain of the highly conserved D194. This mechanism of activation, first proposed by Huber and Bode, organizes the oxyanion hole and primary specificity pocket for substrate binding and catalysis. Using the clotting protease thrombin as a relevant model, we unravel contributions of the I16-D194 ionic interaction to Na⁺ binding, stability of the transition state and the allosteric E*-E equilibrium of the trypsin fold. The I16T mutation abolishes the I16-D194 interaction and compromises the architecture of the oxyanion hole. The D194A mutation also abrogates the I16-D194 interaction but, surprisingly, has no effect on the architecture of the oxyanion hole that remains intact through a new H-bond established between G43 and G193. In both mutants, loss of the I16-D194 ionic interaction compromises Na⁺ binding, reduces stability of the transition state, collapses the 215–217 segment into the primary specific pocket and abrogates the allosteric E*-E equilibrium in favor of a rigid conformation that binds ligand at the active site according to a simple lock-and-key mechanism. These findings refine the structural role of the I16-D194 ionic interaction in the Huber-Bode mechanism of activation and reveal a functional linkage with the allosteric properties of the trypsin fold like Na⁺ binding and the E*-E equilibrium.

Quantum Mechanics/Molecular Mechanics (QM/MM) Calculations Support a Concerted Reaction Mechanism for the Zika Virus NS2B/NS3 Serine Protease with Its Substrate

Zika virus (ZIKV) is mainly transmitted to humans by Aedes species mosquitoes and is associated with serious pathological disorders including microcephaly in newborns and Guillain-Barré syndrome in adults. Currently, there is no vaccine or anti-ZIKV drug available for preventing or controlling ZIKV infection. An attractive drug target for ZIKV treatment is a two-compartment (NS2B/NS3) serine protease that processes viral polyprotein during infection. Here, conventional molecular dynamics simulations of the ZIKV protease in complex with peptide substrate (TGKRS) sequence at the C-terminus of NS2B show that the substrate is in the active conformation for the cleavage reaction by ZIKV protease. Hybrid quantum mechanics/molecular mechanics (QM/MM) umbrella sampling simulations (PM6/ff14SB) of acylation results reveal that proton transfer from S135 to H51 and nucleophilic attack on the substrate by S135 are concerted. The rate-limiting step involves the formation of a tetrahedral intermediate. In addition, the single-point energy QM/MM calculations, precisely at the level of coupled cluster theory (LCCSD(T)/(aug)-cc-pVTZ), were performed to correct the potential energy profiles for the first step of the acylation process. The average computed activation barrier at this level of theory is 16.3 kcal mol^-1. Therefore, the computational approaches presented here are helpful for further designing of NS2B/NS3 inhibitors based on transition-state analogues.

A thermostable trypsin from freshwater fish Japanese dace (Tribolodon hakonensis): a comparison of the primary structures among fish trypsins

Trypsin from Japanese dace (Tribolodon hakonensis) (JD-T) living in freshwater (2-18 °C) was purified. JD-T represented typical fish trypsin characteristics regarding the effects of protease inhibitor, calcium-ion, and pH. For the effect of temperature, JD-T quite resembled to the trypsins from tropical-zone marine fish and freshwater fish (the catfish cultured in Thailand), i.e., the optimum temperature was 60 °C, and it was stable below 60 °C at pH 8.0 for 15 min incubation. From the data, it seemed that the trypsin from freshwater fish is thermostable in spite of the fact that their habitat temperatures are low. So, we determined the primary structure of JD-T to discuss its thermostability-structure relationship. JD-T possessed basic structural features of fish trypsin such as the catalytic triad, the Asp189 residue for substrate specificity, 12 Cys residues forming six disulfide-bridges, and the calcium-ion-binding loop. On the other hand, the contents of charged amino acid residues in whole JD-T molecule (16.2%) and N-terminal region (13.8%) were similar to those of tropical-zone marine fish and other freshwater fish trypsins. Then, JD-T conserved the five amino acid residues (Glu70, Asn72, Val75, Glu77, and Glu80) coordinate with calcium-ion, and the proportion of negatively charged amino acids to charged amino acids in the calcium-ion-binding region of JD-T (75.0%) was equivalent to those of tropical-zone marine fish and freshwater fish trypsins. Therefore, it was suggested that the high thermostability of JD-T are stemmed from these structural specificities.

Directed evolution of the 3C protease from coxsackievirus using a novel fluorescence-assisted intracellular method

Proteases are crucial for regulating biological processes in organisms through hydrolysis of peptide bonds. Recombinant proteases have moreover become important tools in biotechnological, and biomedical research and as therapeutics. We have developed a label-free high-throughput method for quantitative assessment of proteolytic activity in Escherichia coli. The screening method is based on co-expression of a protease of interest and a reporter complex. This reporter consists of an aggregation-prone peptide fused to a fluorescent protein via a linker that contains the corresponding substrate sequence. Cleavage of the substrate rescues the fluorescent protein from aggregation, resulting in increased fluorescence that correlates to proteolytic activity, which can be monitored using flow cytometry. In one round of flow-cytometric cell sorting, we isolated an efficiently cleaved tobacco etch virus (TEV) substrate from a 1:100 000 background of non-cleavable sequences, with around 6000-fold enrichment. We then engineered the 3C protease from coxsackievirus B3 (CVB3 3Cpro) towards improved proteolytic activity on the substrate LEVLFQ↓GP. We isolated highly proteolytic active variants from a randomly mutated CVB3 3Cpro library with up to 4-fold increase in activity. The method enables simultaneous measurement of proteolytic activity and protease expression levels and can therefore be applied for protease substrate profiling, as well as directed evolution of proteases.

Evolutionary and structural analyses uncover a role for solvent interactions in the diversification of cocoonases in butterflies

Multi-omic approaches promise to supply the power to detect genes underlying disease and fitness-related phenotypes. Optimal use of the resulting profusion of data requires detailed investigation of individual candidate genes, a challenging proposition. Here, we combine transcriptomic and genomic data with molecular modelling of candidate enzymes to characterize the evolutionary history and function of the serine protease cocoonase. Heliconius butterflies possess the unique ability to feed on pollen; recent work has identified cocoonase as a candidate gene in pollen digestion. Cocoonase was first described in moths, where it aids in eclosure from the cocoon and is present as a single copy gene. In heliconiine butterflies it is duplicated and highly expressed in the mouthparts of adults. At least six copies of cocoonase are present in Heliconius melpomene and copy number varies across H. melpomene sub-populations. Most cocoonase genes are under purifying selection, however branch-site analyses suggest cocoonase 3 genes may have evolved under episodic diversifying selection. Molecular modelling of cocoonase proteins and examination of their predicted structures revealed that the active site region of each type has a similar structure to trypsin, with the same predicted substrate specificity across types. Variation among heliconiine cocoonases instead lies in the outward-facing residues involved in solvent interaction. Thus, the neofunctionalization of cocoonase duplicates appears to have resulted from the need for these serine proteases to operate in diverse biochemical environments. We suggest that cocoonase may have played a buffering role in feeding during the diversification of Heliconius across the neotropics by enabling these butterflies to digest protein from a range of biochemical milieux.

Metal-Organic-Framework-Based Enzymatic Microfluidic Biosensor via Surface Patterning and Biomineralization

Recently, the biomineralization of enzyme in metal-organic-framework (enzyme-MOF) composite have shown a great potential to increase enzymes stability without compromising their activity; hence, it is desirable for its applications in biosensing devices. Nonetheless, most of the enzyme-MOF research has been focusing on enzyme encapsulation in particle form, which limits its synthesis flexibility for practical applications because of its requirement for postsynthesis immobilization onto solid support. Therefore, to develop a diagnostic device out of the biomineralized enzyme, surface patterning and integration of microfluidic system offers many advantages. In this work, mussel-inspired polydopamine/polyethyleneimine (PDA/PEI) coating is employed to pattern enzyme-MOF in microfluidic channels and exploit the wettability gradient for "pumpless transportation" effect. As a proof of concept, we combine a cascade reaction of glucose oxidase (GOx) and horseradish peroxidase (HRP) enzymes to detect glucose into a patterned zeolitic imidazole framework-8 (ZIF-8) thin film on a flexible polymeric substrate. The results show that the ZIF-8/GOx&HRP in situ composites on PDA/PEI patterns have good acid and thermal stability compared with samples without ZIF-8. ZIF-8/GOx&HRP in situ shows high selectivity toward glucose, linear sensitivity of 0.00303 Abs/μM, and the limit of detection of 8 μM glucose concentration. An unexpected benefit of this approach is the ability of the ZIF-8 thin-film structure to provide a diffusion limiting effect for substrate influx, thus, producing high range of linear response range (8 μM to 5 mM of glucose). This work provides insights into the spatial location of the enzymes in MOF thin films and the potential of such patterning techniques for MOF-based biosensors using other types of biological elements such as antibodies and aptamers.

Crystal structure of mutant carboxypeptidase T from Thermoactinomyces vulgaris with an implanted S1' subsite from pancreatic carboxypeptidase B

A site-directed mutagenesis method has been used to obtain the G215S/A251G/T257A/D260G/T262D mutant of carboxypeptidase T from Thermoactinomyces vulgaris (CPT), in which the amino-acid residues of the S1' subsite are substituted by the corresponding residues from pancreatic carboxypeptidase B (CPB). It was shown that the mutant enzyme retained the broad, mainly hydrophobic selectivity of wild-type CPT. The mutant containing the implanted CPB S1' subsite was crystallized and its three-dimensional structure was determined at 1.29 Å resolution by X-ray crystallography. A comparison of the three-dimensional structures of CPT, the G215S/A251G/T257A/D260G/T262D CPT mutant and CPB showed that the S1' subsite of CPT has not been distorted by the mutagenesis and adequately reproduces the structure of the CPB S1' subsite. The CPB-like mutant differs from CPB in substrate selectivity owing to differences between the two enzymes outside the S1' subsite. Moreover, the difference in substrate specificity between the enzymes was shown to be affected by residues other than those that directly contact the substrate.

Engineering potent mesotrypsin inhibitors based on the plant-derived cyclic peptide, sunflower trypsin inhibitor-1

Plants produce a diverse range of peptides and proteins that inhibit the activity of different serine proteases. The value of these inhibitors not only stems from their native role(s) in planta, but they are also regarded as promising templates for inhibitor engineering. Interest in this field has grown rapidly in recent years, particularly for therapeutic applications. The serine protease mesotrypsin has been implicated in several cancers, but is a challenging target for inhibitor engineering as a number of serine protease inhibitors that typically display broad-range activity show limited activity against mesotrypsin. In this study, we use a cyclic peptide isolated from sunflower seeds, sunflower trypsin inhibitor-1 (SFTI-1), as a scaffold for engineering potent mesotrypsin inhibitors. SFTI-1 comprises 14-amino acids and is a potent inhibitor of human cationic trypsin (K_i = 30 ± 0.8 pM) but shows 165,000-fold weaker activity against mesotrypsin (K_i = 4.96 ± 0.2 μM). Using an inhibitor library based on SFTI-1, we show that the inhibitor's P2' residue (Ile) is a key contributor to SFTI-1's limited activity against mesotrypsin. Substituting P2' Ile with chemically diverse amino acids, including non-canonical aromatic residues, produced new inhibitor variants that maintained a similar structure to SFTI-1 and showed marked improvements in activity (exceeding 100-fold). An assessment of the activity of the new inhibitors against closely-related trypsin paralogs revealed that the improved activity against mesotrypsin was accompanied by a loss in activity against off-target proteases, such that several engineered variants showed comparable activity against mesotrypsin and human cationic trypsin. Together, these findings identify potent mesotrypsin inhibitors that are suitable for further optimisation studies and demonstrate the potential gains in activity and selectivity that can be achieved by optimising the P2' residue, particularly for engineered SFTI-based inhibitors.

Proteolytic Cleavage-Mechanisms, Function, and "Omic" Approaches for a Near-Ubiquitous Posttranslational Modification

Proteases enzymatically hydrolyze peptide bonds in substrate proteins, resulting in a widespread, irreversible posttranslational modification of the protein's structure and biological function. Often regarded as a mere degradative mechanism in destruction of proteins or turnover in maintaining physiological homeostasis, recent research in the field of degradomics has led to the recognition of two main yet unexpected concepts. First, that targeted, limited proteolytic cleavage events by a wide repertoire of proteases are pivotal regulators of most, if not all, physiological and pathological processes. Second, an unexpected in vivo abundance of stable cleaved proteins revealed pervasive, functionally relevant protein processing in normal and diseased tissue-from 40 to 70% of proteins also occur in vivo as distinct stable proteoforms with undocumented N- or C-termini, meaning these proteoforms are stable functional cleavage products, most with unknown functional implications. In this Review, we discuss the structural biology aspects and mechanisms of catalysis by different protease classes. We also provide an overview of biological pathways that utilize specific proteolytic cleavage as a precision control mechanism in protein quality control, stability, localization, and maturation, as well as proteolytic cleavage as a mediator in signaling pathways. Lastly, we provide a comprehensive overview of analytical methods and approaches to study activity and substrates of proteolytic enzymes in relevant biological models, both historical and focusing on state of the art proteomics techniques in the field of degradomics research.

Engineering peptide ligase specificity by proteomic identification of ligation sites

Enzyme-catalyzed peptide ligation is a powerful tool for site-specific protein bioconjugation, but stringent enzyme–substrate specificity limits its utility. Here, we present an approach for comprehensive characterization of peptide ligase specificity for N termini using proteome-derived peptide libraries. We used this strategy to characterize the ligation efficiency for >25,000 enzyme–substrate pairs in the context of the engineered peptide ligase subtiligase and identified a family of 72 mutant subtiligases with activity toward N-terminal sequences that were previously recalcitrant to modification. We applied these mutants individually for site-specific bioconjugation of purified proteins including antibodies, and in algorithmically selected combinations for sequencing of the cellular N terminome with reduced sequence bias. We also developed a web application to enable algorithmic selection of the most efficient subtiligase variant(s) for bioconjugation to user-defined sequences. These studies provide a new toolbox of enzymes for site-specific protein modification and a general approach for rapidly defining and engineering peptide ligase specificity.

Phage-assisted continuous evolution of proteases with altered substrate specificity

Here we perform phage-assisted continuous evolution (PACE) of TEV protease, which canonically cleaves ENLYFQS, to cleave a very different target sequence, HPLVGHM, that is present in human IL-23. A protease emerging from ∼2500 generations of PACE contains 20 non-silent mutations, cleaves human IL-23 at the target peptide bond, and when pre-mixed with IL-23 in primary cultures of murine splenocytes inhibits IL-23-mediated immune signaling. We characterize the substrate specificity of this evolved enzyme, revealing shifted and broadened specificity changes at the six positions in which the target amino acid sequence differed. Mutational dissection and additional protease specificity profiling reveal the molecular basis of some of these changes. This work establishes the capability of changing the substrate specificity of a protease at many positions in a practical time scale and provides a foundation for the development of custom proteases that catalytically alter or destroy target proteins for biotechnological and therapeutic applications.Proteases are promising therapeutics to treat diseases such as hemophilia which are due to endogenous protease deficiency. Here the authors use phage-assisted continuous evolution to evolve a variant TEV protease with altered target peptide sequence specificities.

Functional analysis of TcCTLP-5C2, a chymotrypsin-like serine protease needed for molting in Tribolium castaneum

In a previous study, we have characterized a gene family encoding chymotrypsin-like proteases from the red flour beetle, Tribolium castaneum (TcCTLPs). We identified 14 TcCTLP genes that were predominantly expressed in the midgut, where they presumably function in digestion. Two genes (TcCTLP-6C and TcCTLP-5C₂), however, additionally showed considerable expression in the carcass, and RNAi studies demonstrated that they are required for molting (Broehan et al., 2010; Insect Biochem. Mol. Biol 40, 274-83). Thus, the enzyme has distinct functions in different physiological environments. To study molecular adaptations that facilitate enzyme function in different environments, we performed an in-depth analysis of the molecular and enzymatic properties of TcCTLP-5C_2. We expressed different mutated versions of TcCTLP-5C₂ in form of factor Xa activatable pro-enzymes in insect cells using a baculoviral expression system, and purified the recombinant proteins by affinity chromatography. By measuring and comparing the enzyme activities, we obtained information about the significance of single amino acid residues in motifs that determine substrate specificity and pH tolerance. Further, we showed that TcCTLP-5C₂ is modified by N-glycosylation at amino acid position N137, which lies opposite to the catalytic cleft. Comparison of the enzymatic properties of non-glycosylated and glycosylated TcCTLP-5C₂ versions showed that N-glycosylation decreases V_max (maximum velocity) and k_cat (turnover) while leaving the K_m (specificity) unchanged. Thus, we provide evidence that N-glycosylation alters catalytic behavior by allosteric effects presumably due to altered structural dynamics as observed for chemically glycosylated enzymes.

Stepwise Loop Insertion Strategy for Active Site Remodeling to Generate Novel Enzyme Functions

The remodeling of active sites to generate novel biocatalysts is an attractive and challenging task. We developed a stepwise loop insertion strategy (StLois), in which randomized residue pairs are inserted into active site loops. The phosphotriesterase-like lactonase from Geobacillus kaustophilus (GkaP-PLL) was used to investigate StLois's potential for changing enzyme function. By inserting six residues into active site loop 7, the best variant ML7-B6 demonstrated a 16-fold further increase in catalytic efficiency toward ethyl-paraoxon compared with its initial template, that is a 609-fold higher, >10⁷ fold substrate specificity shift relative to that of wild-type lactonase. The remodeled variants displayed 760-fold greater organophosphate hydrolysis activity toward the organophosphates parathion, diazinon, and chlorpyrifos. Structure and docking computations support the source of notably inverted enzyme specificity. Considering the fundamental importance of active site loops, the strategy has potential for the rapid generation of novel enzyme functions by loop remodeling.

An investigation into the applicability of the semiempirical method PM7 for modeling the catalytic mechanism in the enzyme chymotrypsin

The catalytic cycle for the serine protease α-chymotrypsin was investigated in an attempt to determine the suitability of using the semiempirical method PM7 in the program MOPAC for investigating enzyme-catalyzed reactions. All six classical intermediates were modeled using standard methods, and were characterized as stable minima on the potential energy surface. Using a modified saddle point optimization method, five transition states were located and verified both by vibrational and by intrinsic reaction coordinate analysis. Some individual features, such as the hydrogen bonds in the oxyanion hole, the nature of various electrostatic interactions, and the role of Met192, were examined. This involved designing and running computational experiments to model mutations that would allow features of interest, in particular the energies involved, to be isolated. Three features within the enzyme were examined in detail: the reaction site itself, where covalent bonds were made and broken, the electrostatic effects of the buried aspartate anion, a passive but essential component of the catalytic triad, and the oxyanion hole, where hydrogen bonds help stabilize charged intermediates. With one minor exception, all phenomena investigated agreed with previously-reported descriptions. This result, along with the fact that all the techniques used were relatively straightforward, leads to the recommendation that PM7 and related methods, such as PM6-D3H4, are appropriate for modeling similar enzyme-catalyzed reactions. Graphical abstract Fifth of six transition states, showing water splitting into hydroxyl anion and a proton, to form the second tetrahedral intermediate and histidinium ion. Atoms of the water molecule involved in the hydrolysis are indicated by halos.

EnCOUNTer: a parsing tool to uncover the mature N-terminus of organelle-targeted proteins in complex samples

Background

Characterization of mature protein N-termini by large scale proteomics is challenging. This is especially true for proteins undergoing cleavage of transit peptides when they are targeted to specific organelles, such as mitochondria or chloroplast. Protein neo-N-termini can be located up to 100–150 amino acids downstream from the initiator methionine and are not easily predictable. Although some bioinformatics tools are available, they usually require extensive manual validation to identify the exact N-terminal position. The situation becomes even more complex when post-translational modifications take place at the neo-N-terminus. Although N-terminal acetylation occurs mostly in the cytosol, it is also observed in some organelles such as chloroplast. To date, no bioinformatics tool is available to define mature protein starting positions, the associated N-terminus acetylation status and/or yield for each proteoform. In this context, we have developed the EnCOUNTer tool (i) to score all characterized peptides using discriminating parameters to identify bona fide mature protein N-termini and (ii) to determine the N-terminus acetylation yield of the most reliable ones.

Results

Based on large scale proteomics analyses using the SILProNAQ methodology, tandem mass spectrometry favoured the characterization of thousands of peptides. Data processing using the EnCOUNTer tool provided an efficient and rapid way to extract the most reliable mature protein N-termini. Selected peptides were subjected to N-terminus acetylation yield determination. In an A. thaliana cell lysate, 1232 distinct proteotypic N-termini were characterized of which 648 were located at the predicted protein N-terminus (position 1/2) and 584 were located further downstream (starting at position > 2). A large number of these N-termini were associated with various well-defined maturation processes occurring on organelle-targeted proteins (mitochondria, chloroplast and peroxisome), secreted proteins or membrane-targeted proteins. It was also possible to highlight some protein alternative starts, splicing variants or erroneous protein sequence predictions.

Conclusions

The EnCOUNTer tool provides a unique way to extract accurately the most relevant mature proteins N-terminal peptides from large scale experimental datasets. Such data processing allows the identification of the exact N-terminus position and the associated acetylation yield.

Electronic supplementary material

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

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.

Influence of cholinium-based ionic liquids on the structural stability and activity of α-chymotrypsin

Unanticipated high thermal stability and sustained activity of CT was found in the presence of [Ch][Ac], [Ch][Cl] and [Ch][Dhp], while [Ch][Cit] and [Ch][OH] act as strong destabilizers for the CT structure.