Binding of amino acid side-chains to S1 cavities of serine proteinases

Alchemical Calculation of Relative Free Energies for Charge-Changing Mutations at Protein-Protein Interfaces Considering Fixed and Variable Protonation States

The calculation of relative free energies (ΔΔG) for charge-changing mutations at protein-protein interfaces through alchemical methods remains challenging due to variations in the system's net charge during charging steps, the possibility of mutated and contacting ionizable residues occurring in various protonation states, and undersampling issues. In this study, we present a set of strategies, collectively termed TIRST/TIRST-H+, to address some of these challenges. Our approaches combine thermodynamic integration (TI) with the prediction of pKa shifts to calculate ΔΔG values. Moreover, special sets of restraints are employed to keep the alchemically transformed molecules separated. The accuracy of the devised approaches was assessed on a large and diverse data set comprising 164 point mutations of charged residues (Asp, Glu, Lys, and Arg) to Ala at the protein-protein interfaces of complexes with known three-dimensional structures. Mean absolute and root-mean-square errors ranging from 1.38 to 1.66 and 1.89 to 2.44 kcal/mol, respectively, and Pearson correlation coefficients of ∼0.6 were obtained when testing the approaches on the selected data set using the GPU-TI module of Amber18 suite and the ff14SB force field. Furthermore, the inclusion of variable protonation states for the mutated acid residues improved the accuracy of the predicted ΔΔG values. Therefore, our results validate the use of TIRST/TIRST-H+ in prospective studies aimed at evaluating the impact of charge-changing mutations to Ala on the stability of protein-protein complexes.

Spatom: a graph neural network for structure-based protein-protein interaction site prediction

Accurate identification of protein-protein interaction (PPI) sites remains a computational challenge. We propose Spatom, a novel framework for PPI site prediction. This framework first defines a weighted digraph for a protein structure to precisely characterize the spatial contacts of residues, then performs a weighted digraph convolution to aggregate both spatial local and global information and finally adds an improved graph attention layer to drive the predicted sites to form more continuous region(s). Spatom was tested on a diverse set of challenging protein-protein complexes and demonstrated the best performance among all the compared methods. Furthermore, when tested on multiple popular proteins in a case study, Spatom clearly identifies the interaction interfaces and captures the majority of hotspots. Spatom is expected to contribute to the understanding of protein interactions and drug designs targeting protein binding.

Non-Canonical Amino Acids in Analyses of Protease Structure and Function

All known organisms encode 20 canonical amino acids by base triplets in the genetic code. The cellular translational machinery produces proteins consisting mainly of these amino acids. Several hundred natural amino acids serve important functions in metabolism, as scaffold molecules, and in signal transduction. New side chains are generated mainly by post-translational modifications, while others have altered backbones, such as the β- or γ-amino acids, or they undergo stereochemical inversion, e.g., in the case of D-amino acids. In addition, the number of non-canonical amino acids has further increased by chemical syntheses. Since many of these non-canonical amino acids confer resistance to proteolytic degradation, they are potential protease inhibitors and tools for specificity profiling studies in substrate optimization and enzyme inhibition. Other applications include in vitro and in vivo studies of enzyme kinetics, molecular interactions and bioimaging, to name a few. Amino acids with bio-orthogonal labels are particularly attractive, enabling various cross-link and click reactions for structure-functional studies. Here, we cover the latest developments in protease research with non-canonical amino acids, which opens up a great potential, e.g., for novel prodrugs activated by proteases or for other pharmaceutical compounds, some of which have already reached the clinical trial stage.

Scaffold Hopping Transformations Using Auxiliary Restraints for Calculating Accurate Relative Binding Free Energies

In silico screening of drug-target interactions is a key part of the drug discovery process. Changes in the drug scaffold via contraction or expansion of rings, the breaking of rings, and the introduction of cyclic structures from acyclic structures are commonly applied by medicinal chemists to improve binding affinity and enhance favorable properties of candidate compounds. These processes, commonly referred to as scaffold hopping, are challenging to model computationally. Although relative binding free energy (RBFE) calculations have shown success in predicting binding affinity changes caused by perturbing R-groups attached to a common scaffold, applications of RBFE calculations to modeling scaffold hopping are relatively limited. Scaffold hopping inevitably involves breaking and forming bond interactions of quadratic functional forms, which is highly challenging. A novel method for handling ring opening/closure/contraction/expansion and linker contraction/expansion is presented here. To the best of our knowledge, RBFE calculations on linker contraction/expansion have not been previously reported. The method uses auxiliary restraints to hold the atoms at the ends of a bond in place during the breaking and forming of the bonds. The broad applicability of the method was demonstrated by examining perturbations involving small-molecule macrocycles and mutations of proline in proteins. High accuracy was obtained using the method for most of the perturbations studied. The rigor of the method was isolated from the force field by validating the method using relative and absolute hydration free energy calculations compared to standard simulation results. Unlike other methods that rely on λ-dependent functional forms for bond interactions, the method presented here can be employed using modern molecular dynamics software without modification of codes or force field functions.

A Cloud Computing Platform for Scalable Relative and Absolute Binding Free Energy Predictions: New Opportunities and Challenges for Drug Discovery

Free energy perturbation (FEP) has become widely used in drug discovery programs for binding affinity prediction between candidate compounds and their biological targets. However, limitations of FEP applications also exist, including, but not limited to, high cost, long waiting time, limited scalability, and breadth of application scenarios. To overcome these problems, we have developed XFEP, a scalable cloud computing platform for both relative and absolute free energy predictions using optimized simulation protocols. XFEP enables large-scale FEP calculations in a more efficient, scalable, and affordable way, for example, the evaluation of 5000 compounds can be performed in 1 week using 50-100 GPUs with a computing cost roughly equivalent to the cost for the synthesis of only one new compound. By combining these capabilities with artificial intelligence techniques for goal-directed molecule generation and evaluation, new opportunities can be explored for FEP applications in the drug discovery stages of hit identification, hit-to-lead, and lead optimization based not only on structure exploitation within the given chemical series but also including evaluation and comparison of completely unrelated molecules during structure exploration in a larger chemical space. XFEP provides the basis for scalable FEP applications to become more widely used in drug discovery projects and to speed up the drug discovery process from hit identification to preclinical candidate compound nomination.

Isolation and functional diversity of Bowman-Birk type serine proteinase inhibitors from Hyacinthus orientalis

Bowman-Birk inhibitors (BBIs) are plant-derived serine proteinase inhibitors. Endogenously, they function as defense molecules against pathogens and insects, but they also have been explored for applications in cancer treatment and inflammatory disorders. Here, we isolated 15 novel BBIs from the bulb of Hyacinthus orientalis (termed HOSPIs). These isoinhibitors consisted of two or three chains, respectively, that are linked by disulfides bonds based on proposed cleavage sites in the canonical BBI reactive site loop. They strongly inhibited trypsin (Ki = 0.22-167 nM) and α-chymotrypsin (Ki = 19-1200 nM). Notably, HOSPI-B4 contains a six-residue reactive loop, which appears to be the smallest such motif discovered in BBIs to date. HOSPI-A6 and -A7 contain an unusual reactive site, i.e. Leu-Met at the P1-P1' position and have strong inhibitory activity against trypsin, α-chymotrypsin, and elastase. Analysis of the cDNA encoding HOSPIs revealed that the precursors have HOSPI-like domains repeated at least twice with a defined linker sequence connecting individual domains. Lastly, mutational analysis of HOSPIs suggested that the linker sequence does not affect the inhibitory activity, and a Thr residue at the P2 site and a Pro at the P3' site are crucial for elastase inhibition. Using mammalian proteases as representative model system, we gain novel insight into the sequence diversity and proteolytic activity of plant BBI. These results may aid the rational design of BBI peptides with potent and distinct inhibitory activity against human, pathogen, or insect serine proteinases.

Bioengineering of an elastase inhibitor from Caesalpinia echinata (Brazil wood) seeds

Protease inhibitors have been widely used in several therapeutic applications such as in the treatment of bleeding disorders, hypertension, cancer and pulmonary diseases. In a previous work, we demonstrated that a Kunitz-type serine protease inhibitor isolated from the seeds of Caesalpinia echinata (CeEI) exhibits pharmacological potential in lung inflammatory diseases in which neutrophil elastase plays a crucial role. However, an important challenge in the use of natural products is to ensure a commercially viable production. In this work, we report the cloning, expression and purification of two recombinant CeEI isoinhibitors with 700 base pairs encoding two proteins with 181 amino acid residues (rCeEI-4 and rCeEI-5). After the expression, each yielding 22 mg/L of active protein, both isoinhibitors presented a molecular mass of about 23.0 kDa, evaluated by SDS-PAGE. The inhibition constants for human neutrophil elastase (HNE) were 0.67 nM (rCeEI-4) and 0.57 nM (rCeEI-5), i.e., similar to the native inhibitor (1.90 nM). Furthermore, rCeEI-4 was used as a template to design smaller functional peptides flanking the inhibitor reactive site: rCeEI-36, delimited between the amino acid residues N³⁶ and S⁸⁸ containing a disulfide bond in the reactive-site loop, and rCeEI-46, delimited between S⁴⁶ and L⁷⁵ without the disulfide bond. The yields were 18 mg/L (rCeEI-36) and 12 mg/L (rCeEI-46). Both peptides inhibit HNE in the nanomolar range (K_i 0.30 ± 0.01 and 8.80 ± 0.23, respectively). Considering their size and the inhibitory efficiency, these peptides may be considered in strategies for the development of drugs targeting pulmonary disorders where elastase is involved.

Assessment of software methods for estimating protein-protein relative binding affinities

A growing number of computational tools have been developed to accurately and rapidly predict the impact of amino acid mutations on protein-protein relative binding affinities. Such tools have many applications, for example, designing new drugs and studying evolutionary mechanisms. In the search for accuracy, many of these methods employ expensive yet rigorous molecular dynamics simulations. By contrast, non-rigorous methods use less exhaustive statistical mechanics, allowing for more efficient calculations. However, it is unclear if such methods retain enough accuracy to replace rigorous methods in binding affinity calculations. This trade-off between accuracy and computational expense makes it difficult to determine the best method for a particular system or study. Here, eight non-rigorous computational methods were assessed using eight antibody-antigen and eight non-antibody-antigen complexes for their ability to accurately predict relative binding affinities (ΔΔG) for 654 single mutations. In addition to assessing accuracy, we analyzed the CPU cost and performance for each method using a variety of physico-chemical structural features. This allowed us to posit scenarios in which each method may be best utilized. Most methods performed worse when applied to antibody-antigen complexes compared to non-antibody-antigen complexes. Rosetta-based JayZ and EasyE methods classified mutations as destabilizing (ΔΔG < -0.5 kcal/mol) with high (83-98%) accuracy and a relatively low computational cost for non-antibody-antigen complexes. Some of the most accurate results for antibody-antigen systems came from combining molecular dynamics with FoldX with a correlation coefficient (r) of 0.46, but this was also the most computationally expensive method. Overall, our results suggest these methods can be used to quickly and accurately predict stabilizing versus destabilizing mutations but are less accurate at predicting actual binding affinities. This study highlights the need for continued development of reliable, accessible, and reproducible methods for predicting binding affinities in antibody-antigen proteins and provides a recipe for using current methods.

PhosphoShield: Improving Trypsin Digestion of Phosphoproteins by Shielding the Negatively Charged Phosphate Moiety

Protein phosphorylation is a post-translational modification that is essential to cellular signaling, cellular function, and associated disease progression. Bottom-up proteomics based on enzymatic digestion is the most widely used approach for identifying and quantifying phosphoproteins in complex biological samples. Researchers have largely optimized the experimental conditions for trypsin digestion, and it is now a routine procedure. However, trypsin digestion is impaired by the presence of phosphorylated residues in the protein sequence. This impairment arises from the fact that there are commonly salt bridges between a negatively charged phosphate group and the side chain of protonated arginine or lysine. On average, 55% of all phosphopeptides have their phosphosites located less than three amino acid residues from a cleavage site. Salt bridges reduce the cleavage accessibility for trypsin by masking the basic site chain groups of arginine and lysine. Thus, there are frequent missed cleavages in the vicinity of phosphorylation sites, thereby lessening both the depth of proteome coverage and the quantification accuracy of phosphoproteomics. In this work, we propose a method termed PhosphoShield to mitigate salt bridge formation by adding a digallium complex that exhibits a high binding affinity to the phosphate group. We tested our method using quantitative mass spectrometry analysis of the phosphoproteome of human liver cancer cells (HepG2). PhosphoShield enhances the cleavage frequency of at least 17% of tryptic phosphopeptides having cleavage sites close to the phosphate group.

Dissecting the Energetics of Intrinsically Disordered Proteins via a Hybrid Experimental and Computational Approach

Intrinsically disordered proteins (IDPs) play important roles in biology, but little is known about the energetics of their inter-residue interactions. Methods that have been successfully applied to analyze the energetics of globular proteins are not applicable to the fluctuating partially ordered ensembles populated by IDPs. A combined computational experimental strategy is introduced for analyzing the energetic role of individual residues in the free state of IDPs. The approach combines experimental measurements of the binding of wild-type and mutant IDPs to their partners with alchemical free energy calculations of the structured complexes. These data allow quantitative information to be deduced about the free state via a thermodynamic cycle. The approach is validated by the analysis of the effects of mutations upon the binding free energy of the ovomucoid inhibitor third binding domain to its partners and is applied to the C-terminal domain of the measles virus nucleoprotein, a 125-residue IDP involved in the RNA transcription and replication of measles virus. The analysis reveals significant inter-residue interactions in the unbound IDP and suggests a biological role for them. The work demonstrates that advances in force fields and computational hardware have now led to the point where it is possible to develop methods, which integrate experimental and computational techniques to reveal insights that cannot be studied using either technique alone.

Directed Evolution of Canonical Loops and Their Swapping between Unrelated Serine Proteinase Inhibitors Disprove the Interscaffolding Additivity Model

Reversible serine proteinase inhibitors comprise 18 unrelated families. Each family has a distinct representative structure but contains a surface loop that adopts the same, canonical conformation in the enzyme-inhibitor complex. The Laskowski mechanism universally applies for the action of all canonical inhibitors independent of their scaffold, but it has two nontrivial extrapolations. Intrascaffolding additivity states that all enzyme-contacting loop residues act independently of each other, while interscaffolding additivity claims that these residues act independently of the scaffold. These theories have great importance for engineering proteinase inhibitors but have not been comprehensively challenged. Therefore, we tested the interscaffolding additivity theory by hard-randomizing all enzyme-contacting canonical loop positions of a Kazal- and a Pacifastin-scaffold inhibitor, displaying the variants on M13 phage, and selecting the libraries on trypsin and chymotrypsin. Directed evolution delivered different patterns on both scaffolds against both enzymes, which contradicts interscaffolding additivity. To quantitatively assess the extent of non-additivity, we measured the affinities of the optimal binding loop variants and their binding loop-swapped versions. While optimal variants have picomolar affinities, swapping the evolved loops results in up to 200,000-fold affinity loss. To decipher the underlying causes, we characterized the stability, overall structure and dynamics of the inhibitors with differential scanning calorimetry, circular dichroism and NMR spectroscopy and molecular dynamic simulations. These studies revealed that the foreign loop destabilizes the lower-stability Pacifastin scaffold, while the higher-stability Kazal scaffold distorts the foreign loop. Our findings disprove interscaffolding additivity and show that loop and scaffold form one integrated unit that needs to be coevolved to provide high-affinity inhibition.

Potent, multi-target serine protease inhibition achieved by a simplified β-sheet motif

Engagement of an extended β-sheet is a common substrate/inhibitor interaction at the active site of serine proteases and is an important feature of Laskowski mechanism inhibitors that present a substrate-like loop to a target protease. This loop is cleaved but subsequently relegated forming a stable inhibitor/protease complex. Laskowski inhibitors are ubiquitous in nature and are used extensively in serine protease inhibitor design. However, most studies concentrate on introducing new sidechain interactions rather than the direct contributions of the substrate-like β-sheet to enzyme inhibition. Here we report the crystal structure of an simplified β-sheet inhibitory motif within the Sunflower Trypsin Inhibitor (SFTI) in complex with trypsin. We show that the intramolecular hydrogen bond network of this SFTI variant (SFTI-TCTR) engages the inhibitor sidechains that would normally interact with a target protease, giving mainchain interactions a more prominent role in complex formation. Despite having reduced sidechain interactions, this SFTI variant is remarkably potent and inhibits a diverse range of serine proteases. Crystal structural analysis and molecular modelling of SFTI-TCTR complexes again indicates an interface dominated by β–sheet interactions, highlighting the importance of this motif and the adaptability of SFTI as a scaffold for inhibitor design.

Catalytic Adaptation of Psychrophilic Elastase

The class I pancreatic elastase from Atlantic salmon is considered to be a cold-adapted enzyme in view of the cold habitat, the reduced thermostability of the enzyme, and the fact that it is faster than its mesophilic porcine counterpart at room temperature. However, no experimental characterization of its catalytic properties at lower temperatures has actually been reported. Here we use extensive computer simulations of its catalytic reaction, at different temperatures and with different peptide substrates, to compare its characteristics with those of porcine pancreatic elastase, with which it shares 67% sequence identity. We find that both enzymes have a preference for smaller aliphatic residues at the P1 position, while the reaction rate with phenylalanine at P1 is predicted to be substantially lower. With the former class of substrates, the calculated reaction rates for salmon enzyme are consistently higher than those of the porcine ortholog at all temperatures examined, and the difference is most pronounced at the lowest temperature. As observed for other cold-adapted enzymes, this is caused by redistribution of the activation free energy in terms of enthalpy and entropy and can be linked to differences in the mobility of surface-exposed loops in the two enzymes. Such mobility changes are found to be reflected by characteristic sequence conservation patterns in psychrophilic and mesophilic species. Hence, calculations of mutations in a single surface loop show that the temperature dependence of the catalytic reaction is altered in a predictable way.

[Prospects for the design of new therapeutically significant protease inhibitors based on knottins and sunflower seed trypsin inhibitor (SFTI 1)]

Plant seed knottins, mainly from the Cucurbitacea family, and sunflower seed trypsin inhibitor (SFTI 1) are the most low-molecular canonical peptide inhibitors of serine proteases. High efficiency of inhibition of various serine proteases, structure rigidity together with the possibility of limited variations of amino acid sequences, high chemical stability, lack of toxic properties, opportunity of production by either chemical synthesis or use of heterologous expression systems make these inhibitors attractive templates for design of new compounds for regulation of therapeutically significant serine protease activities. Hence the design of such compounds represents a prospective research field. The review considers structural characteristics of these inhibitors, their properties, methods of preparation and design of new analogs. Examples of successful employment of natural serine protease inhibitors belonging to knottin family and SFTI 1 as templates for the design of highly specific inhibitors of certain proteases are given.

Polar Desolvation and Position 226 of Pancreatic and Neutrophil Elastases Are Crucial to their Affinity for the Kunitz-Type Inhibitors ShPI-1 and ShPI-1/K13L

The Kunitz-type protease inhibitor ShPI-1 inhibits human neutrophil elastase (HNE, K_i = 2.35·10⁻⁸ M) but does not interact with the porcine pancreatic elastase (PPE); whereas its P1 site variant, ShPI-1/K13L, inhibits both HNE and PPE (K_i = 1.3·10⁻⁹ M, and K_i = 1.2·10⁻⁸ M, respectively). By employing a combination of molecular modeling tools, e.g., structural alignment, molecular dynamics simulations and Molecular Mechanics Generalized-Born/Poisson-Boltzmann Surface Area free energy calculations, we showed that D226 of HNE plays a critical role in the interaction of this enzyme with ShPI-1 through the formation of a strong salt bridge and hydrogen bonds with K13 at the inhibitor’s P1 site, which compensate the unfavorable polar-desolvation penalty of the latter residue. Conversely, T226 of PPE is unable to establish strong interactions with K13, thereby precluding the insertion of K13 side-chain into the S1 subsite of this enzyme. An alternative conformation of K13 site-chain placed at the entrance of the S1 subsite of PPE, similar to that observed in the crystal structure of ShPI-1 in complex with chymotrypsin (PDB: 3T62), is also unfavorable due to the lack of stabilizing pair-wise interactions. In addition, our results suggest that the higher affinity of ShPI-1/K13L for both elastases mainly arises from the lower polar-desolvation penalty of L13 compared to that of K13, and not from stronger pair-wise interactions of the former residue with those of each enzyme. These results provide insights into the PPE and HNE inhibition and may contribute to the design of more potent and/or specific inhibitors toward one of these proteases.

Eliminating anti-nutritional plant food proteins: the case of seed protease inhibitors in pea

Several classes of seed proteins limit the utilisation of plant proteins in human and farm animal diets, while plant foods have much to offer to the sustainable intensification of food/feed production and to human health. Reduction or removal of these proteins could greatly enhance seed protein quality and various strategies have been used to try to achieve this with limited success. We investigated whether seed protease inhibitor mutations could be exploited to enhance seed quality, availing of induced mutant and natural Pisum germplasm collections to identify mutants, whilst acquiring an understanding of the impact of mutations on activity. A mutant (TILLING) resource developed in Pisum sativum L. (pea) and a large germplasm collection representing Pisum diversity were investigated as sources of mutations that reduce or abolish the activity of the major protease inhibitor (Bowman-Birk) class of seed protein. Of three missense mutations, predicted to affect activity of the mature trypsin / chymotrypsin inhibitor TI1 protein, a C77Y substitution in the mature mutant inhibitor abolished inhibitor activity, consistent with an absolute requirement for the disulphide bond C77-C92 for function in the native inhibitor. Two further classes of mutation (S85F, E109K) resulted in less dramatic changes to isoform or overall inhibitory activity. The alternative strategy to reduce anti-nutrients, by targeted screening of Pisum germplasm, successfully identified a single accession (Pisum elatius) as a double null mutant for the two closely linked genes encoding the TI1 and TI2 seed protease inhibitors. The P. elatius mutant has extremely low seed protease inhibitory activity and introgression of the mutation into cultivated germplasm has been achieved. The study provides new insights into structure-function relationships for protease inhibitors which impact on pea seed quality. The induced and natural germplasm variants identified provide immediate potential for either halving or abolishing the corresponding inhibitory activity, along with associated molecular markers for breeding programmes. The potential for making large changes to plant protein profiles for improved and sustainable food production through diversity is illustrated. The strategy employed here to reduce anti-nutritional proteins in seeds may be extended to allergens and other seed proteins with negative nutritional effects. Additionally, the novel variants described for pea will assist future studies of the biological role and health-related properties of so-called anti-nutrients.

Molecular Cloning and Functional Studies of Two Kazal-Type Serine Protease Inhibitors Specifically Expressed by Nasonia vitripennis Venom Apparatus

Two cDNA sequences of Kazal-type serine protease inhibitors (KSPIs) in Nasonia vitripennis, NvKSPI-1 and NvKSPI-2, were characterized and their open reading frames (ORFs) were 198 and 264 bp, respectively. Both NvKSPI-1 and NvKSPI-2 contained a typical Kazal-type domain. Real-time quantitative PCR (RT-qPCR) results revealed that NvKSPI-1 and NvKSPI-2 mRNAs were mostly detected specifically in the venom apparatus, while they were expressed at lower levels in the ovary and much lower levels in other tissues tested. In the venom apparatus, both NvKSPI-1 and NvKSPI-2 transcripts were highly expressed on the fourth day post eclosion and then declined gradually. The NvKSPI-1 and NvKSPI-2 genes were recombinantly expressed utilizing a pGEX-4T-2 vector, and the recombinant products fused with glutathione S-transferase were purified. Inhibition of recombinant GST-NvKSPI-1 and GST-NvKSPI-2 to three serine protease inhibitors (trypsin, chymotrypsin, and proteinase K) were tested and results showed that only NvKSPI-1 could inhibit the activity of trypsin. Meanwhile, we evaluated the influence of the recombinant GST-NvKSPI-1 and GST-NvKSPI-2 on the phenoloxidase (PO) activity and prophenoloxidase (PPO) activation of hemolymph from a host pupa, Musca domestica. Results showed PPO activation in host hemolymph was inhibited by both recombinant proteins; however, there was no significant inhibition on the PO activity. Our results suggested that NvKSPI-1 and NvKSPI-2 could inhibit PPO activation in host hemolymph and trypsin activity in vitro.

Calculation of binding free energies

Molecular dynamics simulations enable access to free energy differences governing the driving force underlying all biological processes. In the current chapter we describe alchemical methods allowing the calculation of relative free energy differences. We concentrate on the binding free energies that can be obtained using non-equilibrium approaches based on the Crooks Fluctuation Theorem. Together with the theoretical background, the chapter covers practical aspects of hybrid topology generation, simulation setup, and free energy estimation. An important aspect of the validation of a simulation setup is illustrated by means of calculating free energy differences along a full thermodynamic cycle. We provide a number of examples, including protein-ligand and protein-protein binding as well as ligand solvation free energy calculations.

Bowman-Birk inhibitors from legumes as colorectal chemopreventive agents

Aberrant functioning of serine proteases in inflammatory and carcinogenic processes within the gastrointestinal tract (GIT) has prompted scientists to investigate the potential of serine protease inhibitors, both natural and synthetic, as modulators of their proteolytic activities. Protease inhibitors of the Bowman-Birk type, a major protease inhibitor family in legume seeds, which inhibit potently and specifically trypsin- and chymotrypsin-like proteases, are currently being investigated as colorectal chemopreventive agents. Physiologically relevant amounts of Bowman-Birk inhibitors (BBI) can reach the large intestine in active form due to their extraordinary resistance to extreme conditions within the GIT. Studies in animal models have proven that dietary BBI from several legume sources, including soybean, pea, lentil and chickpea, can prevent or suppress carcinogenic and inflammatory processes within the GIT. Although the therapeutic targets and the action mechanism of BBI have not yet been elucidated, the emerging evidence suggests that BBI exert their preventive properties via protease inhibition; in this sense, serine proteases should be considered as primary targets in early stages of carcinogenesis. The validation of candidate serine proteases as therapeutic targets together with the identification, within the wide array of natural BBI variants, of the most potent and specific protease inhibitors, are necessary to better understand the potential of this protein family as colorectal chemopreventive agents.

N-tetradecylcarbamyl lipopeptides as novel agonists for Toll-like receptor 2

New analogues (UPam) of triacylated lipopeptide Pam3CysSK4, a popular agonist of Toll-like receptor 2 (TLR2), were designed making use of the cocrystal structure of a TLR2 heterodimer (with TLR1) with Pam3CysSK4. Twenty-two UPam derivatives that feature an N-tetradecylcarbamyl chain to mimic the native N-palmitoyl moiety and various small amino acids residues at the penultimate N-terminal position were prepared via solid-phase synthesis. In vitro evaluation of immunostimulatory properties revealed new potent TLR2 ligands.

Interaction of Subtilisin BPN′ and Recombinant Fungal Protease Inhibitor F from Silkworm with Substituted P1Site Residues

Fungal protease inhibitor F (FPI-F) from silkworm inhibits subtilisin and fungal proteases. FPI-F mutants P(1) residues of which, Thr(29), were replaced with Glu, Phe, Gly, Leu, Met, and Arg, were prepared. The inhibitory activities of mutated FPI-F against subtilisin and other mammalian proteases indicated that FPI-F might be a specific inhibitor toward subtilisin-type protease.

Using a Caesalpinia echinata Lam. protease inhibitor as a tool for studying the roles of neutrophil elastase, cathepsin G and proteinase 3 in pulmonary edema

Acute lung injury (ALI) is characterized by neutrophil infiltration and the release of proteases, mainly elastase (NE), cathepsin G (Cat G) and proteinase 3 (PR3), which can be controlled by specific endogenous inhibitors. However, inhibitors of these proteases have been isolated from different sources, including plants. For this study, CeEI, or Caesalpinia echinata elastase inhibitor, was purified from C. echinata (Brazil-wood) seeds after acetone fractionation, followed by ion exchange and reversed phase chromatographic steps. Characterization with SDS-PAGE, stability assays, amino acid sequencing and alignment with other protein sequences confirmed that CeEI is a member of the soybean Kunitz trypsin inhibitor family. Like other members of this family, CeEI is a 20 kDa monomeric protein; it is stable within a large pH and temperature range, with four cysteine residues forming two disulfide bridges, conserved amino acid residues and leucine-isoleucine residues in the reactive site. CeEI was able to inhibit NE and Cat G at a nanomolar range (with K(i)s of 1.9 and 3.6 nM, respectively) and inhibited PR3 within a micromolar range (K(i) 3.7 μM), leading to hydrolysis of specific synthetic substrates. In a lung edema model, CeEI reduced the lung weight and pulmonary artery pressure until 180 min after the injection of zymosan-activated polymorphonuclear neutrophils. In experiments performed in the presence of a Cat G and PR3, but not an NE inhibitor, lung edema was reduced only until 150 min and pulmonary artery pressure was similar to that of the control. These results confirm that NE action is crucial to edema establishment and progression. Additionally, CeEI appears to be a useful tool for studying the physiology of pulmonary edema and provides a template for molecular engineering and drug design for ALI therapy.

Additivity-based design of the strongest possible turkey ovomucoid third domain inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB)

We describe here successful designs of strong inhibitors for porcine pancreatic elastase (PPE) and Streptomyces griseus protease B (SGPB). For each enzyme two inhibitor variants were designed. In one, the reactive site residue (position 18) was retained and the best residues were substituted at contact positions 13, 14, and 15. In the other variant the best residues were substituted at all contact positions except the reactive site where a Gly was substituted. The four designed variants were: for PPE, T(13)E(14)Y(15)-OMTKY3 and T(13)E(14)Y(15)G(18)M(21)P(32)V(36)-OMTKY3, and for SGPB, S(13)D(14)Y(15)-OMTKY3 and S(13)D(14)Y(15)G(18)I(19)K(21)-OMTKY3. The free energies of association (ΔG(0)) of expressed variants have been measured with the proteases for which they were designed as well as with five other serine proteases and the results are discussed.

Comparison of periplasmic and intracellular expression of Arabidopsis thionin proproteins in E. coli

Thionins are antimicrobial plant peptides produced as preproproteins consisting of a signal peptide, the thionin domain, and a so-called acidic domain. Only thionin itself has been isolated from plants. To study the processing of the precursor, it has to be produced in a heterologous system. Since both domains contain several cysteines and, due to the known antimicrobial activity of the thionin, we tested the expression of all four Arabidopsis proproteins as fusion proteins. Periplasmic expression as fusion with maltose binding protein was not successful but cytoplasmic expression as His-tagged TRX fusion proteins with a TEV recognition sequence resulted in proteins of correct size. Use of the SHuffle strain C3030 further improved the expression. Fusion proteins inhibited growth of Escherichia coli. They could be cleaved by TEV protease, releasing authentic proproteins without any additional amino acid at the N-terminus.

SKEMPI: a Structural Kinetic and Energetic database of Mutant Protein Interactions and its use in empirical models

MOTIVATION

Empirical models for the prediction of how changes in sequence alter protein-protein binding kinetics and thermodynamics can garner insights into many aspects of molecular biology. However, such models require empirical training data and proper validation before they can be widely applied. Previous databases contained few stabilizing mutations and no discussion of their inherent biases or how this impacts model construction or validation.

RESULTS

We present SKEMPI, a database of 3047 binding free energy changes upon mutation assembled from the scientific literature, for protein-protein heterodimeric complexes with experimentally determined structures. This represents over four times more data than previously collected. Changes in 713 association and dissociation rates and 127 enthalpies and entropies were also recorded. The existence of biases towards specific mutations, residues, interfaces, proteins and protein families is discussed in the context of how the data can be used to construct predictive models. Finally, a cross-validation scheme is presented which is capable of estimating the efficacy of derived models on future data in which these biases are not present.

AVAILABILITY

The database is available online at http://life.bsc.es/pid/mutation_database/.

Secondary Binding Site of Trypsin: Revealed by Crystal Structure of Trypsin-Peptide Complex

Designed synthetic heterochiral peptides, when added to porcine trypsin, resulted in reduction of enzyme activity. The crystal structure of a complex formed between porcine trypsin and a heterochiral hepta peptide Boc-Pro-DAsp-Aib-Leu-Aib-Leu-Ala-NHMe has been determined at 1.9 A resolution. The hepta peptide does not bind at the active site, but is located in the interstitial region, and interacts with the calcium-binding loop (residues 60-80). The bound peptide interacts with the active site residue Ser195 through an acetate ion, and with Lys 60 mediated by water molecules. The structure, when compared with the other trypsin-peptide complexes, suggests that the flexibility of surface loops, concerted movement of the loops towards the active site, and the interaction of the bound peptide with Lys 60, may be responsible for the reduction in enzyme activity. This study provides a structural evidence for the earlier biochemical observation regarding the role of surface loops in the catalysis of the enzyme.

The P(2)' residue is a key determinant of mesotrypsin specificity: engineering a high-affinity inhibitor with anticancer activity

PRSS3/mesotrypsin is an atypical isoform of trypsin, the up-regulation of which has been implicated in promoting tumour progression. Mesotrypsin inhibitors could potentially provide valuable research tools and novel therapeutics, but small-molecule trypsin inhibitors have low affinity and little selectivity, whereas protein trypsin inhibitors bind poorly and are rapidly degraded by mesotrypsin. In the present study, we use mutagenesis of a mesotrypsin substrate, APPI (amyloid precursor protein Kunitz protease inhibitor domain), and of a poor mesotrypsin inhibitor, BPTI (bovine pancreatic trypsin inhibitor), to dissect mesotrypsin specificity at the key P(2)' position. We find that bulky and charged residues strongly disfavour binding, whereas acidic residues facilitate catalysis. Crystal structures of mesotrypsin complexes with BPTI variants provide structural insights into mesotrypsin specificity and inhibition. Through optimization of the P(1) and P(2)' residues of BPTI, we generate a stable high-affinity mesotrypsin inhibitor with an equilibrium binding constant K(i) of 5.9 nM, a >2000-fold improvement in affinity over native BPTI. Using this engineered inhibitor, we demonstrate the efficacy of pharmacological inhibition of mesotrypsin in assays of breast cancer cell malignant growth and pancreatic cancer cell invasion. Although further improvements in inhibitor selectivity will be important before clinical potential can be realized, the results of the present study support the feasibility of engineering protein protease inhibitors of mesotrypsin and highlight their therapeutic potential.

Tailoring elastase inhibition with synthetic peptides

Chronic wounds are the result of excessive amounts of tissue destructive proteases such as human neutrophil elastase (HNE). The high levels of this enzyme found on those types of wounds inactivate the endogenous inhibitor barrier thus, the search for new HNE inhibitors is required. This work presents two new HNE inhibitor peptides, which were synthesized based on the reactive-site loop of the Bowman-Birk inhibitor protein. The results obtained indicated that these new peptides are competitive inhibitors for HNE and, the inhibitory activity can be modulated by modifications introduced at the N- and C-terminal of the peptides. Furthermore, these peptides were also able to inhibit elastase from a human wound exudate while showing no cytotoxicity against human skin fibroblasts in vitro, greatly supporting their potential application in chronic wound treatment.

Cleavage of peptide bonds bearing ionizable amino acids at P(1) by serine proteases with hydrophobic S(1) pocket

Enzymatic hydrolysis of the synthetic substrate succinyl-Ala-Ala-Pro-Xxx-pNA (where Xxx=Leu, Asp or Lys) catalyzed by bovine chymotrypsin (CHYM) or Streptomyces griseus protease B (SGPB) has been studied at different pH values in the pH range 3-11. The pH optima for substrates having Leu, Asp, and Lys have been found to be 7.5-8.0, 5.5-6.0, and ∼10, respectively. At the normally reported pH optimum (pH 7-8) of CHYM and SGPB, the substrate with Leu at the reactive site is more than 25,000-fold more reactive than that with Asp. However, when fully protonated, Asp is nearly as good a substrate as Leu. The pK values of the side chains of Asp and Lys in the hydrophobic S(1) pocket of CHYM and SGPB have been calculated from pH-dependent hydrolysis data and have been found to be about 9 for Asp and 7.4 and 9.7 for Lys for CHYM and SGPB, respectively. The results presented in this communication suggest a possible application of CHYM like enzymes in cleaving peptide bonds contributed by acidic amino acids between pH 5 and 6.

Trp-26 imparts functional versatility to human alpha-defensin HNP1

We performed a comprehensive alanine scan of human alpha-defensin HNP1 and tested the ability of the resulting analogs to kill Staphylococcus aureus, inhibit anthrax lethal factor, and bind human immunodeficiency virus-1 gp120. By far, the most deleterious mutation for all of these functions was W26A. The activities lost by W26A-HNP1 were restored progressively by replacing W26 with non-coded, straight-chain aliphatic amino acids of increasing chain length. The hydrophobicity of residue 26 also correlated with the ability of the analogs to bind immobilized wild type HNP1 and to undergo further self-association. Thus, the hydrophobicity of residue 26 is not only a key determinant of the direct interactions of HNP1 with target molecules, but it also governs the ability of this peptide to form dimers and more complex quaternary structures at micromolar concentrations. Although all defensin peptides are cationic, their amphipathicity is at least as important as their positive charge in enabling them to participate in innate host defense.

Engineering substrate preference in subtilisin: structural and kinetic analysis of a specificity mutant

Bacillus subtilisin has been a popular model protein for engineering altered substrate specificity. Although some studies have succeeded in increasing the specificity of subtilisin, they also demonstrate that high specificity is difficult to achieve solely by engineering selective substrate binding. In this paper, we analyze the structure and transient state kinetic behavior of Sbt160, a subtilisin engineered to strongly prefer substrates with phenylalanine or tyrosine at the P4 position. As in previous studies, we measure improvements in substrate affinity and overall specificity. Structural analysis of an inactive version of Sbt160 in complex with its cognate substrate reveals improved interactions at the S4 subsite with a P4 tyrosine. Comparison of transient state kinetic behavior against an optimal sequence (DFKAM) and a similar, but suboptimal, sequence (DVRAF) reveals the kinetic and thermodynamic basis for increased specificity, as well as the limitations of this approach. While highly selective substrate binding is achieved in Sbt160, several factors cause sequence specificity to fall short of that observed with natural processing subtilisins. First, for substrate sequences which are nearly optimal, the acylation reaction becomes faster than substrate dissociation. As a result, the level of discrimination among these substrates diminishes due to the coupling between substrate binding and the first chemical step (acylation). Second, although Sbt160 has 24-fold higher substrate affinity for the optimal substrate DFKAM than for DVRAF, the increased substrate binding energy is not translated into improved transition state stabilization of the acylation reaction. Finally, as interactions at subsites become stronger, the rate-determining step in peptide hydrolysis changes from acylation to product release. Thus, the release of the product becomes sluggish and leads to a low k(cat) for the reaction. This also leads to strong product inhibition of substrate turnover as the reaction progresses. The structural and kinetic analysis reveals that differences in the binding modes at subsites for substrates, transition states, and products are subtle and difficult to manipulate via straightforward protein engineering. These findings suggest several new strategies for engineering highly sequence selective enzymes.

High sequence variability among hemocyte-specific Kazal-type proteinase inhibitors in decapod crustaceans

Crustacean hemocytes were found to produce a large number of transcripts coding for Kazal-type proteinase inhibitors (KPIs). A detailed study performed with the crayfish Pacifastacus leniusculus and the shrimp Penaeus monodon revealed the presence of at least 26 and 20 different Kazal domains from the hemocyte KPIs, respectively. Comparisons with KPIs from other taxa indicate that the sequences of these domains evolve rapidly. A few conserved positions, e.g. six invariant cysteines were present in all domain sequences whereas the position of P1 amino acid, a determinant for substrate specificity, varied highly. A study with a single crayfish animal suggested that even at the individual level considerable sequence variability among hemocyte KPIs produced exist. Expression analysis of four crayfish KPI transcripts in hematopoietic tissue cells and different hemocyte types suggest that some of these KPIs are likely to be involved in hematopoiesis or hemocyte release as they were produced in particular hemocyte types or maturation stages only.

Kazal-type serine proteinase inhibitors from the black tiger shrimp Penaeus monodon and the inhibitory activities of SPIPm4 and 5

Serine proteinase inhibitors (SPIs) play important roles in physiological and immunological processes involving proteinases in all multicellular organisms. In black tiger shrimp Penaeus monodon, nine different Kazal-type SPIs, namely SPIPm1-9, were identified from the cDNA libraries of hemocyte, hepatopancreas, hematopoietic tissue, ovary and lymphoid organ. They are multi-domain SPIs containing 2-7 and possibly more Kazal domains. Two interesting cDNA clones, SPIPm4 and SPIPm5 coding for two-domain Kazal-type SPIs, were identified from the heat-treated hemocyte cDNA libraries. The SPIPm4 and SPIPm5 consist of open reading frames of 387 and 399 bp coding for polypeptides of 128 and 132 amino acids with putative signal peptides of 21 and 19 amino acid residues and mature SPIs of 107 and 113 amino acid residues, respectively. Recombinant expression in an Escherichia coli expression system yielded recombinant proteins, rSPIPm4 and rSPIPm5, with molecular masses of 12.862 and 13.433 kDa, respectively. The inhibitory activities of SPIPm4 and SPIPm5 were tested against trypsin, chymotrypsin, subtilisin and elastase. The SPIPm4 exhibited potent inhibitory activity against subtilisin and weakly against chymotrypsin whereas the SPIPm5 strongly inhibited subtilisin and elastase. The inhibition was a competitive type with inhibition constants (K(i)) of 14.95 nM for SPIPm4 against subtilisin, 4.19 and 59.64 nM, respectively, for SPIPm5 against subtilisin and elastase. They had no bacteriostatic effect against Gram-positive bacteria: Bacillus subtilis, Bacillus megaterium, Staphylococcus aureus, and Gram-negative bacteria: Vibrio harveyi 639, E. coli JM109. Gene expression study revealed that the SPIPm5 gene was up-regulated in response to heat treatment suggesting the involvement of SPIs in stress responses.

Elastase digests: new ammunition for shotgun membrane proteomics

Despite many advances in membrane proteomics during the last decade the fundamental problem of accessing the transmembrane regions itself has only been addressed to some extent. The present study establishes a method for the nano-LC-based analysis of complex membrane proteomes on the basis of a methanolic porcine pancreatic elastase digest to increase transmembrane coverage. Halobacterium salinarium purple and Corynebacterium glutamicum membranes were successfully analyzed by using the new protocol. We demonstrated that elastase digests yield a large proportion of transmembrane peptides, facilitating membrane protein identification. The potential for characterization of a membrane protein through full sequence coverage using elastase is there but is restricted to the higher abundance protein components. Compatibility of the work flow with the two most common mass spectrometric ionization techniques, ESI and MALDI, was shown. Currently better results are obtained using ESI mainly because of the low response of MALDI for strictly neutral peptides. New findings concerning elastase specificity in complex protein mixtures reveal a new prospect beyond the application in shotgun experiments. Furthermore peptide mass fingerprinting with less specific enzymes might be done in the near future but requires an adaptation of current search algorithms to the new proteases.

Domain inhibitory and bacteriostatic activities of the five-domain Kazal-type serine proteinase inhibitor from black tiger shrimp Penaeus monodon

Serine proteinase inhibitors (SPIs) in multi-cellular organisms are important modulators of proteinase activities in various biological processes. A five-domain Kazal-type SPI SPIPm2 from the black tiger shrimp Penaeus monodon is presumably involved in innate immune response. The SPIPm2 with the domain P1 residues T, A, E, K and E was isolated from the hemocyte cDNA libraries and found to strongly inhibit subtilisin and elastase, and weakly inhibit trypsin. To unravel further the inhibitory activity of each domain, we subcloned, over-expressed and purified each individual SPI domain. Their inhibitory specificities against trypsin, subtilisin and elastase were determined. Domain 1 was found to be inactive. Domains 2, 3 and 5 inhibited subtilisin. Domain 2 inhibited also elastase. Domain 4 weakly inhibited subtilisin and trypsin. The intact SPIPm2 inhibitor was found to possess bacteriostatic activity against the Bacillus subtilis but not the Bacillus megaterium, Staphylococcus aureus, Vibrio harveyi 639 and Escherichia coli JM109. Domains 2, 4 and 5 contributed to this bacteriostatic activity.

Neisseria gonorrhoeae-induced human defensins 5 and 6 increase HIV infectivity: role in enhanced transmission

Sexually transmitted infections (STIs) increase the likelihood of HIV transmission. Defensins are part of the innate mucosal immune response to STIs and therefore we investigated their role in HIV infection. We found that human defensins 5 and 6 (HD5 and HD6) promoted HIV infection, and this effect was primarily during viral entry. Enhancement was seen with primary viral isolates in primary CD4(+) T cells and the effect was more pronounced with R5 virus compared with X4 virus. HD5 and HD6 promoted HIV reporter viruses pseudotyped with vesicular stomatitis virus and murine leukemia virus envelopes, indicating that defensin-mediated enhancement was not dependent on CD4 and coreceptors. Enhancement of HIV by HD5 and HD6 was influenced by the structure of the peptides, as loss of the intramolecular cysteine bonds was associated with loss of the HIV-enhancing effect. Pro-HD5, the precursor and intracellular form of HD5, also exhibited HIV-enhancing effect. Using a cervicovaginal tissue culture system, we found that expression of HD5 and HD6 was induced in response to Neisseria gonorrhoeae (GC, for gonococcus) infection and that conditioned medium from GC-exposed cervicovaginal epithelial cells with elevated levels of HD5 also enhanced HIV infection. Introduction of small interfering RNAs for HD5 or HD6 abolished the HIV-enhancing effect mediated by GC. Thus, the induction of these defensins in the mucosa in the setting of GC infection could facilitate HIV infection. Furthermore, this study demonstrates the complexity of defensins as innate immune mediators in HIV transmission and warrants further investigation of the mechanism by which defensins modulate HIV infection.

The photochemistry properties on interaction silver with tryptophan

UV-vis and fluorescence spectra for interaction silver(I) ion with tryptophan (Trp) have been studied. The influence of pH of media, multicomponent concentration, including amino acids, silver(I) ion, formaldehyde, sodium dodecyl sulfate (SDS), and trihydroxymethylaminomethane (Tris) as well as temperature, illumination, time, etc. on reaction were investigated, and the mechanism of reaction has been explored. The results shown that the optimum condition for the interaction between tryptophan and silver(I) ion were found.

1.2A-resolution crystal structures reveal the second tetrahedral intermediates of streptogrisin B (SGPB)

Streptogrisin B (SGPB) has served as one of the models for studying the catalytic activities of serine peptidases. Here we report its native crystal structures at pH 4.2 at a resolution of 1.18A, and at pH 7.3 at a resolution of 1.23A. Unexpectedly, outstanding electron density peaks occurred in the active site and the substrate-binding region of SGPB in the computed maps at both pHs. The densities at pH 4.2 were assigned as a tetrapeptide, Asp-Ala-Ile-Tyr, whereas those at pH 7.3 were assigned as a tyrosine molecule and a leucine molecule existing at equal occupancies in both of the SGPB molecules in the asymmetric unit. Refinement with relaxed geometric restraints resulted in molecular structures representing mixtures of the second tetrahedral intermediates and the enzyme-product complexes of SGPB existing in a pH-dependent equilibrium. Structural comparisons with the complexes of SGPB with turkey ovomucoid third domain (OMTKY3) and its variants have shown that, upon the formation of the tetrahedral intermediate, residues Glu192A to Gly193 of SGPB move towards the alpha-carboxylate O of residue P1 of the bound species, and adjustments in the side-chain conformational angles of His57 and Ser195 of SGPB favor the progression of the catalytic mechanism of SGPB.