Structure of the complex formed by bovine trypsin and bovine pancreatic trypsin inhibitor. II. Crystallographic refinement at 1.9 A resolution

Advancements in Inactivation of Soybean Trypsin Inhibitors

Soybean Trypsin Inhibitors (STIs) in soy-based foods have negative effects on soybean protein digestion and pancreatic health of humans. The inactivation of STIs is a critical unit operation aimed at enhancing the nutritional properties of soy-based foods during processing. This paper reviews the structure of STIs and soybean proteins, as well as the mechanisms of digestion. Various technologies (physical, chemical, biological) have been used to inactivate STIs. Their parameter settings, operating procedures, advantages, and disadvantages are also described. Mechanisms of inactivation of STIs (Kunitz trypsin inhibitor (KTI) and Bowman-Birk inhibitor (BBI)) conformations under different treatments are clarified. In addition, emerging technologies, e.g., Ohmic Heating, Electron Beam Irradiation, Dielectric-Barrier Discharge, and probiotics, have demonstrated great potential to inactivate STIs. We advise that multiple emerging technologies should combine with other unit operating systems to maximize inactivation efficiency.

Molecular Insights into Substrate Binding of the Outer Membrane Enzyme OmpT

The enzyme OmpT of the outer membrane of Escherichia coli shows proteolytic activity and cleaves peptides and proteins. Using molecular dynamics simulations in a fully hydrated lipid bilayer on a time scale of hundreds of nanoseconds, we draw a detailed atomic picture of substrate recognition in the OmpT-holo enzyme complex. Hydrogen bonds and salt bridges are essential for maintaining the integrity of the active site and play a central role for OmpT in recognizing its substrate. Electrostatic interactions are critical at all stages from approaching the substrate to docking at the active site. Computational alanine scanning based on the Molecular Mechanics Generalized Born Surface Area (MM-GBSA) approach confirms the importance of multiple residues in the active site that form salt bridges. The substrate fluctuates along the axis of the β-barrel, which is associated with oscillations of the binding cleft formed by the residue pairs D210-H212 and D83-D85. Principal component analysis suggests that substrate and protein movements are correlated. We observe the transient presence of putative catalytic water molecules near the active site, which may be involved in the nucleophilic attack on the cleavable peptide bond of the substrate.

Interfering with Host Proteases in SARS-CoV-2 Entry as a Promising Therapeutic Strategy

Due to its fast international spread and substantial mortality, the coronavirus disease COVID-19 evolved to a global threat. Since currently, there is no causative drug against this viral infection available, science is striving for new drugs and approaches to treat the new disease. Studies have shown that the cell entry of coronaviruses into host cells takes place through the binding of the viral spike (S) protein to cell receptors. Priming of the S protein occurs via hydrolysis by different host proteases. The inhibition of these proteases could impair the processing of the S protein, thereby affecting the interaction with the host-cell receptors and preventing virus cell entry. Hence, inhibition of these proteases could be a promising strategy for treatment against SARS-CoV-2. In this review, we discuss the current state of the art of developing inhibitors against the entry proteases furin, the transmembrane serine protease type-II (TMPRSS2), trypsin, and cathepsin L.

A Molecular Lid Mechanism of K+ Channel Blocker Action Revealed by a Cone Peptide

Many venomous organisms carry in their arsenal short polypeptides that block K⁺ channels in a highly selective manner. These toxins may compete with the permeating ions directly via a "plug" mechanism or indirectly via a "pore-collapse" mechanism. An alternative "lid" mechanism was proposed but remained poorly defined. Here we study the Drosophila Shaker channel block by Conkunitzin-S1 and Conkunitzin-C3, two highly similar toxins derived from cone venom. Despite their similarity, the two peptides exhibited differences in their binding poses and biophysical assays, implying discrete action modes. We show that while Conkunitzin-S1 binds tightly to the channel turret and acts via a "pore-collapse" mechanism, Conkunitzin-C3 does not contact this region. Instead, Conk-C3 uses a non-conserved Arg to divert the permeant ions and trap them in off-axis cryptic sites above the SF, a mechanism we term a "molecular-lid". Our study provides an atomic description of the "lid" K⁺ blocking mode and offers valuable insights for the design of therapeutics based on venom peptides.

Effect of Concentrated Salts Solutions on the Stability of Immobilized Enzymes: Influence of Inactivation Conditions and Immobilization Protocol

This paper aims to investigate the effects of some salts (NaCl, (NH₄)₂SO₄ and Na₂SO₄) at pH 5.0, 7.0 and 9.0 on the stability of 13 different immobilized enzymes: five lipases, three proteases, two glycosidases, and one laccase, penicillin G acylase and catalase. The enzymes were immobilized to prevent their aggregation. Lipases were immobilized via interfacial activation on octyl agarose or on glutaraldehyde-amino agarose beads, proteases on glyoxyl agarose or glutaraldehyde-amino agarose beads. The use of high concentrations of salts usually has some effects on enzyme stability, but the intensity and nature of these effects depends on the inactivation pH, nature and concentration of the salt, enzyme and immobilization protocol. The same salt can be a stabilizing or a destabilizing agent for a specific enzyme depending on its concentration, inactivation pH and immobilization protocol. Using lipases, (NH₄)₂SO₄ generally permits the highest stabilities (although this is not a universal rule), but using the other enzymes this salt is in many instances a destabilizing agent. At pH 9.0, it is more likely to find a salt destabilizing effect than at pH 7.0. Results confirm the difficulty of foreseeing the effect of high concentrations of salts in a specific immobilized enzyme.

Insight to the residue in P2 position prevents the peptide inhibitor from being hydrolyzed by serine proteases

Peptidic inhibitors of proteases are attracting increasing interest not only as drug candidates but also for studying the function and regulation mechanisms of these enzymes. Previously, we screened out a cyclic peptide inhibitor of human uPA [Formula: see text] and found that Ala substitution of P2 residue turns upain-1 to a substrate. To further investigate the effect of P2 residue on the peptide behavior transformation, we constructed upain-1-W3F, which has Phe replacement in the P2 position. We determined K_D and K_i of upain-1-W3F and found that upain-1-W3F might still exist as an inhibitor. Furthermore, the high-resolution crystal structure of upain-1-W3F·uPA reveals that upain-1-W3F indeed stays as an intact inhibitor bind to uPA. We thus propose that the P2 residue plays a nonnegligible role in the conversion of upain-1 to a substrate. These results also proposed a strategy to optimize the pharmacological properties of peptide-based drug candidates by hydrophobicity and steric hindrance.Abbreviations : uPA: urokinase-type plasminogen activator; SPD: serine protease domain; S1 pocket: specific substrate-binding pocket.

Surface loops of trypsin-like serine proteases as determinants of function

Trypsin and chymotrypsin-like serine proteases from family S1 (clan PA) constitute the largest protease group in humans and more generally in vertebrates. The prototypes chymotrypsin, trypsin and elastase represent simple digestive proteases in the gut, where they cleave nearly any protein. Multidomain trypsin-like proteases are key players in the tightly controlled blood coagulation and complement systems, as well as related proteases that are secreted from diverse immune cells. Some serine proteases are expressed in nearly all tissues and fluids of the human body, such as the human kallikreins and kallikrein-related peptidases with specialization for often unique substrates and accurate timing of activity. HtrA and membrane-anchored serine proteases fulfill important physiological tasks with emerging roles in cancer. The high diversity of all family members, which share the tandem β-barrel architecture of the chymotrypsin-fold in the catalytic domain, is conferred by the large differences of eight surface loops, surrounding the active site. The length of these loops alters with insertions and deletions, resulting in remarkably different three-dimensional arrangements. In addition, metal binding sites for Na+, Ca2+ and Zn2+ serve as regulatory elements, as do N-glycosylation sites. Depending on the individual tasks of the protease, the surface loops determine substrate specificity, control the turnover and allow regulation of activation, activity and degradation by other proteins, which are often serine proteases themselves. Most intriguingly, in some serine proteases, the surface loops interact as allosteric network, partially tuned by protein co-factors. Knowledge of these subtle and complicated molecular motions may allow nowadays for new and specific pharmaceutical or medical approaches.

Structural analysis of a Vibrio phospholipase reveals an unusual Ser-His-chloride catalytic triad

Phospholipases can disrupt host membranes and are important virulence factors in many pathogens. VvPlpA is a phospholipase A₂ secreted by Vibrio vulnificus and essential for virulence. Its homologs, termed thermolabile hemolysins (TLHs), are widely distributed in Vibrio bacteria, but no structural information for this virulence factor class is available. Herein, we report the crystal structure of VvPlpA to 1.4-Å resolution, revealing that VvPlpA contains an N-terminal domain of unknown function and a C-terminal phospholipase domain and that these two domains are packed closely together. The phospholipase domain adopts a typical SGNH hydrolase fold, containing the four conserved catalytic residues Ser, Gly, Asn, and His. Interestingly, the structure also disclosed that the phospholipase domain accommodates a chloride ion near the catalytic His residue. The chloride is five-coordinated in a distorted bipyramid geometry, accepting hydrogen bonds from a water molecule and the amino groups of surrounding residues. This chloride substitutes for the most common Asp/Glu residue and forms an unusual Ser-His-chloride catalytic triad in VvPlpA. The chloride may orient the catalytic His and stabilize the charge on its imidazole ring during catalysis. Indeed, VvPlpA activity depended on chloride concentration, confirming the important role of chloride in catalysis. The VvPlpA structure also revealed a large hydrophobic substrate-binding pocket that is capable of accommodating a long-chain acyl group. Our results provide the first structure of the TLH family and uncover an unusual Ser-His-chloride catalytic triad, expanding our knowledge on the biological role of chloride.

Analysis of Kunitz inhibitors from plants for comprehensive structural and functional insights

Legume Kunitz type trypsin inhibitor (KTI) family is one of the most versatile families of proteins. A typical KTI features a single peptide folded in β-trefoil manner, with the molecular weight about 20-22kDa and two disulphide bonds. The members are known to inhibit a wide range of serpins proteases at the same time many of them possess unique features. Copaifera langsdorffii Trypsin inhibitor (CTI) has a β-trefoil fold made up of two non-covalently bound polypeptide chains with only a single disulfide bridge. Delonix regia Trypsin inhibitor (DrTI) has one amino acid insertion between P1 and P2 of the reactive site distorting its conformation. Bauhinia bauhinioides Cruzipain inhibitor (BbCI) has a conservative β-trefoil fold but lacks disulfide bonds. Such subtle differences in structures make Kunitz inhibitors different from other inhibitor families. Most of the studies on these inhibitors are focused towards their proposed role in defense from insect pests and wounding but their exact physiological role in nature is still uncharted. Thus, it would be very interesting to closely analyze the structural details of these inhibitors in order to ascertain their biological role and other fascinating applications.

A clip domain serine protease involved in moulting in the silkworm, Bombyx mori: cloning, characterization, expression patterns and functional analysis

Clip domain serine proteases (CLIPs), characterized by one or more conserved clip domains, are essential components of extracellular signalling cascades in various biological processes, especially in innate immunity and the embryonic development of insects. Additionally, CLIPs may have additional non-immune functions in insect development. In the present study, the clip domain serine protease gene Bombyx mori serine protease 95 (BmSP95), which encodes a 527-residue protein, was cloned from the integument of B. mori. Bioinformatics analysis indicated that BmSP95 is a typical CLIP of the subfamily D and possesses a clip domain at the N terminus, a trypsin-like serine protease (tryp_spc) domain at the C terminus and a conserved proline-rich motif between these two domains. At the transcriptional level, BmSP95 is expressed in the integument during moulting and metamorphosis, and the expression pattern is consistent with the fluctuating 20-hydroxyecdysone (20E) titre in B. mori. At the translational level, BmSP95 protein is synthesized in the epidermal cells, secreted as a zymogen and activated in the moulting fluid. Immunofluorescence revealed that BmSP95 is distributed into the old endocuticle in the moulting stage. The expression of BmSP95 was upregulated by 20E. Moreover, expression of BmSP95 was downregulated by pathogen infection. RNA interference-mediated silencing of BmSP95 led to delayed moulting from pupa to moth. These results suggest that BmSP95 is involved in integument remodelling during moulting and metamorphosis.

Strong and widespread action of site-specific positive selection in the snake venom Kunitz/BPTI protein family

S1 family of serine peptidases is the largest family of peptidases. They are specifically inhibited by the Kunitz/BPTI inhibitors. Kunitz domain is characterized by the compact 3D structure with the most important inhibitory loops for the inhibition of S1 peptidases. In the present study we analysed the action of site-specific positive selection and its impact on the structurally and functionally important parts of the snake venom Kunitz/BPTI family of proteins. By using numerous models we demonstrated the presence of large numbers of site-specific positively selected sites that can reach between 30-50% of the Kunitz domain. The mapping of the positively selected sites on the 3D model of Kunitz/BPTI inhibitors has shown that these sites are located in the inhibitory loops 1 and 2, but also in the Kunitz scaffold. Amino acid replacements have been found exclusively on the surface, and the vast majority of replacements are causing the change of the charge. The consequence of these replacements is the change in the electrostatic potential on the surface of the Kunitz/BPTI proteins that may play an important role in the precise targeting of these inhibitors into the active site of S1 family of serine peptidases.

Effects of Acids, Bases, and Heteroatoms on Proximal Radial Distribution Functions for Proteins

The proximal distribution of water around proteins is a convenient method of quantifying solvation. We consider the effect of charged and sulfur-containing amino acid side-chain atoms on the proximal radial distribution function (pRDF) of water molecules around proteins using side-chain analogs. The pRDF represents the relative probability of finding any solvent molecule at a distance from the closest or surface perpendicular protein atom. We consider the near-neighbor distribution. Previously, pRDFs were shown to be universal descriptors of the water molecules around C, N, and O atom types across hundreds of globular proteins. Using averaged pRDFs, a solvent density around any globular protein can be reconstructed with controllable relative error. Solvent reconstruction using the additional information from charged amino acid side-chain atom types from both small models and protein averages reveals the effects of surface charge distribution on solvent density and improves the reconstruction errors relative to simulation. Solvent density reconstructions from the small-molecule models are as effective and less computationally demanding than reconstructions from full macromolecular models in reproducing preferred hydration sites and solvent density fluctuations.

Molecular medicine - To be or not to be

Molecular medicine is founded on the synergy between Chemistry, Physics, Biology and Medicine, with the ambitious goal of tackling diseases from a molecular perspective. This Review aims at retracing a personal outlook of the birth and development of molecular medicine, as well as at highlighting some of the most urgent challenges linked to aging and represented by incurable neurodegenerative diseases caused by protein misfolding. Furthermore, we emphasize the emerging role of the retromer dysfunctions and improper protein sorting in Alzheimer's disease and other important neurological disordered.

Insights into the serine protease mechanism based on structural observations of the conversion of a peptidyl serine protease inhibitor to a substrate

BACKGROUND

Serine proteases are one of the most studied group of enzymes. Despite the extensive mechanistic studies, some crucial details remain controversial, for example, how the cleaved product is released in the catalysis reaction. A cyclic peptidyl inhibitor (CSWRGLENHRMC, upain-1) of a serine protease, urokinase-type plasminogen activator (uPA), was found to become a slow substrate and cleaved slowly upon the replacement of single residue (W3A).

METHODS

By taking advantage of the unique property of this peptide, we report the high-resolution structures of uPA in complex with upain-1-W3A peptide at four different pH values by X-ray crystallography.

RESULTS

In the structures obtained at low pH (pH4.6 and 5.5), the cyclic peptide upain-1-W3A was found to be intact and remained in the active site of uPA. At 7.4, the scissile bond of the peptide was found cleaved, showing that the peptide became a uPA substrate. At pH9.0, the C-terminal part of the substrate was no longer visible, and only the P1 residue occupying the S1 pocket was identified.

CONCLUSIONS

The analysis of these structures provides explanations why the upain-1-W3A is a slow substrate. In addition, we clearly identified the cleaved fragments of the peptide at both sides of the scissile bond in the active site of the enzyme, showing a slow release of the cleaved peptide.

GENERAL SIGNIFICANCE

This work indicates that the quick release of the cleaved P' fragment after the first step of hydrolysis may not always be needed for the second hydrolysis.

Drosophila melanogaster clip-domain serine proteases: Structure, function and regulation

Mammalian chymotrypsin-like serine proteases (SPs) are one of the best-studied family of enzymes with roles in a wide range of physiological processes, including digestion, blood coagulation, fibrinolysis and humoral immunity. Extracellular SPs can form cascades, in which one protease activates the zymogen of the next protease in the chain, to amplify physiological or pathological signals. These extracellular SPs are generally multi-domain proteins, with pro-domains that are involved in protein-protein interactions critical for the sequential organization of the cascades, the control of their intensity and their proper localization. Far less is known about invertebrate SPs than their mammalian counterparts. In insect genomes, SPs and their proteolytically inactive homologs (SPHs) constitute large protein families. In addition to the chymotrypsin fold, many of these proteins contain additional structural domains, often with conserved mammalian orthologues. However, the largest group of arthropod SP regulatory modules is the clip domains family, which has only been identified in arthropods. The clip-domain SPs are extracellular and have roles in the immune response and embryonic development. The powerful reverse-genetics tools in Drosophila melanogaster have been essential to identify the functions of clip-SPs and their organization in sequential cascades. This review focuses on the current knowledge of Drosophila clip-SPs and presents, when necessary, data obtained in other insect models. We will first cover the biochemical and structural features of clip domain SPs and SPHs. Clip-SPs are implicated in three main biological processes: the control of the dorso-ventral patterning during embryonic development; the activation of the Toll-mediated response to microbial infections and the prophenoloxydase cascade, which triggers melanization. Finally, we review the regulation of SPs and SPHs, from specificity of activation to inhibition by endogenous or pathogen-encoded inhibitors.

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.

The ABC of protein kinase conformations

Due to their involvement in human diseases, protein kinases are an important therapeutic target class. Conformation is a key concept for understanding how functional activity, inhibition and sequence are linked. We assemble and annotate the mammalian structural kinome from the Protein Data Bank on the basis of a universal residue nomenclature. We identify a torsion angle around the Gly of the DFG-motif whose sharp distribution profile corresponds to three eclipsed conformations. This allows the definition a small set of clusters whose distribution shows a bias for the active conformation. A common rationale links the active and inactive state: stabilization of the active conformation, as well as inactivation by displacement of helix-αC or the DFG-motif is governed by the interaction between helix-αC and the DFG motif. In particular, the conformation of the DFG-motif is tightly correlated with the propensity of helix-αC displacement. Our analysis reveals detailed mechanisms for the displacement of helix-αC and the DFG and improves our understanding of the role of individual residues. By pooling conformations from the whole structural kinome, the energetic contributions of sequence and extrinsic factors can be estimated in free energy analyses. This article is part of a Special Issue entitled: Inhibitors of Protein Kinases.

Inhibition of outer membrane proteases of the omptin family by aprotinin

Bacterial proteases are important virulence factors that inactivate host defense proteins and contribute to tissue destruction and bacterial dissemination. Outer membrane proteases of the omptin family, exemplified by Escherichia coli OmpT, are found in some Gram-negative bacteria. Omptins cleave a variety of substrates at the host-pathogen interface, including plasminogen and antimicrobial peptides. Multiple omptin substrates relevant to infection have been identified; nonetheless, an effective omptin inhibitor remains to be found. Here, we purified native CroP, the OmpT ortholog in the murine pathogen Citrobacter rodentium. Purified CroP was found to readily cleave both a synthetic fluorescence resonance energy transfer substrate and the murine cathelicidin-related antimicrobial peptide. In contrast, CroP was found to poorly activate plasminogen into active plasmin. Although classical protease inhibitors were ineffective against CroP activity, we found that the serine protease inhibitor aprotinin displays inhibitory potency in the micromolar range. Aprotinin was shown to act as a competitive inhibitor of CroP activity and to interfere with the cleavage of the murine cathelicidin-related antimicrobial peptide. Importantly, aprotinin was able to inhibit not only CroP but also Yersinia pestis Pla and, to a lesser extent, E. coli OmpT. We propose a structural model of the aprotinin-omptin complex in which Lys15 of aprotinin forms salt bridges with conserved negatively charged residues of the omptin active site.

Computer models of the human immunoglobulins Binding sites and molecular interactions

In last month's issue(1) Richard Pumphrey showed how low resolution computer models could be used to investigate the shape and segmental flexibility of human immunoglobulin molecules. In this article, he extends the use of these models to the study of interactions between immunoglobulins and other molecules, such as their receptors, J chain, secretory component and complement fragments. Although precise studies of molecular interaction require computation in atomic detail, there is still much to be learnt from low resolution modelling, in which atomic details are glossed over but general principles become more obvious.

Cleavage specificity analysis of six type II transmembrane serine proteases (TTSPs) using PICS with proteome-derived peptide libraries

Background

Type II transmembrane serine proteases (TTSPs) are a family of cell membrane tethered serine proteases with unclear roles as their cleavage site specificities and substrate degradomes have not been fully elucidated. Indeed just 52 cleavage sites are annotated in MEROPS, the database of proteases, their substrates and inhibitors.

Methodology/Principal Finding

To profile the active site specificities of the TTSPs, we applied Proteomic Identification of protease Cleavage Sites (PICS). Human proteome-derived database searchable peptide libraries were assayed with six human TTSPs (matriptase, matriptase-2, matriptase-3, HAT, DESC and hepsin) to simultaneously determine sequence preferences on the N-terminal non-prime (P) and C-terminal prime (P’) sides of the scissile bond. Prime-side cleavage products were isolated following biotinylation and identified by tandem mass spectrometry. The corresponding non-prime side sequences were derived from human proteome databases using bioinformatics. Sequencing of 2,405 individual cleaved peptides allowed for the development of the family consensus protease cleavage site specificity revealing a strong specificity for arginine in the P1 position and surprisingly a lysine in P1′ position. TTSP cleavage between R↓K was confirmed using synthetic peptides. By parsing through known substrates and known structures of TTSP catalytic domains, and by modeling the remainder, structural explanations for this strong specificity were derived.

Conclusions

Degradomics analysis of 2,405 cleavage sites revealed a similar and characteristic TTSP family specificity at the P1 and P1′ positions for arginine and lysine in unfolded peptides. The prime side is important for cleavage specificity, thus making these proteases unusual within the tryptic-enzyme class that generally has overriding non-prime side specificity.

Trypsin isozymes in the lobster Panulirus argus (Latreille, 1804): from molecules to physiology

Trypsin enzymes have been studied in a wide variety of animal taxa due to their central role in protein digestion as well as in other important physiological and biotechnological processes. Crustacean trypsins exhibit a high number of isoforms. However, while differences in properties of isoenzymes are known to play important roles in regulating different physiological processes, there is little information on this aspect for decapod trypsins. The aim of this review is to integrate recent findings at the molecular level on trypsin enzymes of the spiny lobster Panulirus argus, into higher levels of organization (biochemical, organism) and to interpret those findings in relation to the feeding ecology of these crustaceans. Trypsin in lobster is a polymorphic enzyme, showing isoforms that differ in their biochemical features and catalytic efficiencies. Molecular studies suggest that polymorphism in lobster trypsins may be non-neutral. Trypsin isoenzymes are differentially regulated by dietary proteins, and it seems that some isoenzymes have undergone adaptive evolution coupled with a divergence in expression rate to increase fitness. This review highlights important but poorly studied issues in crustaceans in general, such as the relation among trypsin polymorphism, phenotypic (digestive) flexibility, digestion efficiency, and feeding ecology.

Shape dependence of the radial distribution function of hydration water around proteins

The characterization of hydration water in proteins is important to understand their structure, function and folding properties. A calculation of the shape parameters reveals considerable asymmetry in the shapes of globular proteins. The present study suggests a generalized approach for the calculation of radial distribution of hydration water by accounting for the shape asymmetry in proteins. The surface and radial distribution function is analyzed for three groups of high resolution globular proteins and nonglobular proteins. This generalized approach depicts a considerable difference in the hydration water distribution pattern around aspherical proteins as compared to the earlier method and this difference is more pronounced for nonglobular proteins. The peaks for normalized RDF are found to be sharper compared to those of bulk water. The normalized RDF of hydration water exactly coincides with that of the bulk beyond 8.0 Å. The radial distribution of hydration water as a function of the water-protein distance matches with the experimentally observed distribution of hydration water around myoglobin. The results reveal that the distribution of hydration water is dependent on the shapes of proteins and hence a generalized approach should be used for the calculation of hydration water distribution around proteins, especially for nonglobular proteins.

Structure-kinetic relationship study of CDK8/CycC specific compounds

In contrast with the very well explored concept of structure-activity relationship, similar studies are missing for the dependency between binding kinetics and compound structure of a protein ligand complex, the structure-kinetic relationship. Here, we present a structure-kinetic relationship study of the cyclin-dependent kinase 8 (CDK8)/cyclin C (CycC) complex. The scaffold moiety of the compounds is anchored in the kinase deep pocket and extended with diverse functional groups toward the hinge region and the front pocket. These variations can cause the compounds to change from fast to slow binding kinetics, resulting in an improved residence time. The flip of the DFG motif ("DMG" in CDK8) to the inactive DFG-out conformation appears to have relatively little influence on the velocity of binding. Hydrogen bonding with the kinase hinge region contributes to the residence time but has less impact than hydrophobic complementarities within the kinase front pocket.

Structure of the recombinant BPTI/Kunitz-type inhibitor rShPI-1A from the marine invertebrate Stichodactyla helianthus

The BPTI/Kunitz-type inhibitor family includes several extremely potent serine protease inhibitors. To date, the inhibitory mechanisms have only been studied for mammalian inhibitors. Here, the first crystal structure of a BPTI/Kunitz-type inhibitor from a marine invertebrate (rShPI-1A) is reported to 2.5 Å resolution. Crystallization of recombinant rShPI-1A required the salt-induced dissociation of a trypsin complex that was previously formed to avoid intrinsic inhibitor aggregates in solution. The rShPI-1A structure is similar to the NMR structure of the molecule purified from the natural source, but allowed the assignment of disulfide-bridge chiralities and the detection of an internal stabilizing water network. A structural comparison with other BPTI/Kunitz-type canonical inhibitors revealed unusual ϕ angles at positions 17 and 30 to be a particular characteristic of the family. A significant clustering of ϕ and ψ angle values in the glycine-rich remote fragment near the secondary binding loop was additionally identified, but its impact on the specificity of rShPI-1A and similar molecules requires further study.

[Interaction of sea amemone Heteractis crispa Kunitz type polypeptides with pain vanilloid receptor TRPV1: in silico investigation]

Using methods of molecular biology we defined the structures of the 31 sea anemone Heteractis crispa genes encoding polypeptides which are structurally homologous to the Kunitz proteinase inhibitor family. Identified amino acid sequences have point residue substitutions, high degree of homology with sequences of known H. crispa Kunitz family members, and represent a combinatorial library of polypeptides. We generated their three-dimensional structures by homologous modeling methods. Analysis of their molecular electrostatic potential enabled us to divide given polypeptides into three clusters. One of them includes polypeptides APHC1, APHC2 and APHC3, which were earlier shown to possess a unique property of inhibiting of the pain vanilloid receptor TRPV1 in vitro and providing the analgesic effects in vivo in addition to their trypsin inhibitory activity. Molecular docking made possible establishing the spatial structure of the complexes, the nature of the polypeptides binding with TRPV1, as well as functionally important structural elements involved in the complex formation. Structural models have enabled us to propose a hypothesis contributing to understanding the APHC1-3 impact mechanism for the pain signals transduction by TRPV1: apparently, there is an increase of the receptor relaxation time resulted in binding of its two chains with the polypeptide molecule, which disrupt the functioning of the TRPV1 and leads to partial inhibition of signal transduction in electrophysiological experiments.

Structural insights into serine protease inhibition by a marine invertebrate BPTI Kunitz-type inhibitor

Proteins isolated from marine invertebrates are frequently characterized by exceptional structural and functional properties. ShPI-1, a BPTI Kunitz-type inhibitor from the Caribbean Sea anemone Stichodactyla helianthus, displays activity not only against serine-, but also against cysteine-, and aspartate proteases. As an initial step to evaluate the molecular basis of its activities, we describe the crystallographic structure of ShPI-1 in complex with the serine protease bovine pancreatic trypsin at 1.7Å resolution. The overall structure and the important enzyme-inhibitor interactions of this first invertebrate BPTI-like Kunitz-type inhibitor:trypsin complex remained largely conserved compared to mammalian BPTI-Kunitz inhibitor complexes. However, a prominent stabilizing role within the interface was attributed to arginine at position P3. Binding free-energy calculations indicated a 10-fold decrease for the inhibitor affinity against trypsin, if the P3 residue of ShPI-1 is mutated to alanine. Together with the increased role of Arg(11) at P3 position, slightly reduced interactions at the prime side (Pn') of the primary binding loop and at the secondary binding loop of ShPI-1 were detected. In addition, the structure provides important information for site directed mutagenesis to further optimize the activity of rShPI-1A for biotechnological applications.

Atypical reactive center Kunitz-type inhibitor from the sea anemone Heteractis crispa

The primary structure of a new Kunitz-type protease inhibitor InhVJ from the sea anemone Heteractis crispa (Radianthus macrodactylus) was determined by protein sequencing and cDNA cloning. InhVJ amino acid sequence was shown to share high sequence identity (up to 98%) with the other known Kunitz-type sea anemones sequences. It was determined that the P1 Thr at the reactive site resulted in a decrease of the K(i) of InhVJ to trypsin and α-chymotrypsin (7.38 × 10(-8) M and 9.93 × 10(-7) M, respectively). By structure modeling the functional importance of amino acids at the reactive site as well as at the weak contact site were determined. The significant role of Glu45 for the orientation and stabilization of the InhVJ-trypsin complex was elucidated. We can suggest that there has been an adaptive evolution of the P1 residue at the inhibitor reactive site providing specialization or functional diversification of the paralogs. The appearance of a key so-called P1 Thr residue instead of Lys might lead to refinement of inhibitor specificity in the direction of subfamilies of serine proteases. The absence of Kv channel and TRPV1-receptor modulation activity was confirmed by electrophysiological screening tests.

Three genes expressing Kunitz domains in the epididymis are related to genes of WFDC-type protease inhibitors and semen coagulum proteins in spite of lacking similarity between their protein products

Background

We have previously identified a locus on human chromosome 20q13.1, encompassing related genes of postulated WFDC-type protease inhibitors and semen coagulum proteins. Three of the genes with WFDC motif also coded for the Kunitz-type protease inhibitor motif. In this report, we have reinvestigated the locus for homologous genes encoding Kunitz motif only. The identified genes have been analyzed with respect to structure, expression and function.

Results

We identified three novel genes; SPINT3, SPINT4 and SPINT5, and the structure of their transcripts were determined by sequencing of DNA generated by rapid amplification of cDNA ends. Each gene encodes a Kunitz domain preceded by a typical signal peptide sequence, which indicates that the proteins of 7.6, 8.7, and 9.7 kDa are secreted. Analysis of transcripts in 26 tissues showed that the genes predominantly are expressed in the epididymis. The recombinantly produced proteins could not inhibit the amidolytic activity of trypsin, chymotrypsin, plasmin, thrombin, coagulation factor Xa, elastase, urokinase and prostate specific antigen, whereas similarly made bovine pancreatic trypsin inhibitor (BPTI) had the same bioactivity as the protein isolated from bovine pancreas.

Conclusions

The similar organization, chromosomal location and site of expression, suggests that the novel genes are homologous with the genes of WFDC-type protease inhibitors and semen coagulum proteins, despite the lack of similarity in primary structure of their protein products. Their restricted expression to the epididymis suggests that they could be important for male reproduction. The recombinantly produced proteins are presumably bioactive, as demonstrated with similarly made BPTI, but may have a narrower spectrum of inhibition, as indicated by the lacking activity against eight proteases with differing specificity. Another possibility is that they have lost the protease inhibiting properties, which is typical of Kunitz domains, in favor of hitherto unknown functions.

Local order and mobility of water molecules around ambivalent helices

Water on a protein surface plays a key role in determining the structure and dynamics of proteins. Compared to the properties of bulk water, many aspects of the structure and dynamics of the water surrounding the proteins are less understood. It is interesting therefore to explore how the properties of the water within the solvation shell around the peptide molecule depend on its specific secondary structure. In this work we investigate the orientational order and residence times of the water molecules to characterize the structure, energetics, and dynamics of the hydration shell water around ambivalent peptides. Ambivalent sequences are identical sequences which display multiple secondary structures in different proteins. Molecular dynamics simulations of representative proteins containing variable helix, variable nonhelix, and conserved helix are also used to explore the local structure and mobility of water molecules in their vicinity. The results, for the first time, depict a different water distribution pattern around the conserved and variable helices. The water molecules surrounding the helical segments in variable helices are found to possess a less locally ordered structure compared to those around their corresponding nonhelical counterparts and conserved helices. The long conserved helices exhibit extremely high local residence times compared to the helical conformations of the variable helices, whereas the residence times of the nonhelical conformations of the variable helices are comparable to those of the short conserved helices. This differential pattern of the structure and dynamics of water molecules in the vicinity of conserved/variable helices may lend valuable insights for understanding the role of solvent effects in determining sequence ambivalency and help in improving the accuracy of water models used in the simulations of proteins.

Structure of hydration water in proteins: a comparison of molecular dynamics simulations and database analysis

Hydration layer water molecules play important structural and functional roles in proteins. Despite being a critical component in biomolecular systems, characterizing the properties of hydration water poses a challenge for both experiments and simulations. In this context we investigate the local structure of hydration water molecules as a function of the distance from the protein and water molecules respectively in 188 high resolution protein structures and compare it with those obtained from molecular dynamics simulations. Tetrahedral order parameter of water in proteins calculated from previous and present simulation studies show that the potential of bulk water overestimates the average tetrahedral order parameter compared to those calculated from crystal structures. Hydration waters are found to be more ordered at a distance between the first and second solvation shell from the protein surface. The values of the order parameter decrease sharply when the water molecules are located very near or far away from the protein surface. At small water-water distance, the values of order parameter of water are very low. The average order parameter records a maximum value at a distance equivalent to the first solvation layer with respect to the water-water radial distribution and asymptotically approaches a constant value at large distances. Results from present analysis will help to get a better insight into structure of hydration water around proteins. The analysis will also help to improve the accuracy of water models on the protein surface.

Crystal structures of a halophilic archaeal malate synthase from Haloferax volcanii and comparisons with isoforms A and G

BACKGROUND

Malate synthase, one of the two enzymes unique to the glyoxylate cycle, is found in all three domains of life, and is crucial to the utilization of two-carbon compounds for net biosynthetic pathways such as gluconeogenesis. In addition to the main isoforms A and G, so named because of their differential expression in E. coli grown on either acetate or glycolate respectively, a third distinct isoform has been identified. These three isoforms differ considerably in size and sequence conservation. The A isoform (MSA) comprises ~530 residues, the G isoform (MSG) is ~730 residues, and this third isoform (MSH-halophilic) is ~430 residues in length. Both isoforms A and G have been structurally characterized in detail, but no structures have been reported for the H isoform which has been found thus far only in members of the halophilic Archaea.

RESULTS

We have solved the structure of a malate synthase H (MSH) isoform member from Haloferax volcanii in complex with glyoxylate at 2.51 Å resolution, and also as a ternary complex with acetyl-coenzyme A and pyruvate at 1.95 Å. Like the A and G isoforms, MSH is based on a β8/α8 (TIM) barrel. Unlike previously solved malate synthase structures which are all monomeric, this enzyme is found in the native state as a trimer/hexamer equilibrium. Compared to isoforms A and G, MSH displays deletion of an N-terminal domain and a smaller deletion at the C-terminus. The MSH active site is closely superimposable with those of MSA and MSG, with the ternary complex indicating a nucleophilic attack on pyruvate by the enolate intermediate of acetyl-coenzyme A.

CONCLUSIONS

The reported structures of MSH from Haloferax volcanii allow a detailed analysis and comparison with previously solved structures of isoforms A and G. These structural comparisons provide insight into evolutionary relationships among these isoforms, and also indicate that despite the size and sequence variation, and the truncated C-terminal domain of the H isoform, the catalytic mechanism is conserved. Sequence analysis in light of the structure indicates that additional members of isoform H likely exist in the databases but have been misannotated.

Structural insight into serine protease Rv3671c that Protects M. tuberculosis from oxidative and acidic stress

Rv3671c, a putative serine protease, is crucial for persistence of Mycobacterium tuberculosis in the hostile environment of the phagosome. We show that Rv3671c is required for M. tuberculosis resistance to oxidative stress in addition to its role in protection from acidification. Structural and biochemical analyses demonstrate that the periplasmic domain of Rv3671c is a functional serine protease of the chymotrypsin family and, remarkably, that its activity increases on oxidation. High-resolution crystal structures of this protease in an active strained state and in an inactive relaxed state reveal that a solvent-exposed disulfide bond controls the protease activity by constraining two distant regions of Rv3671c and stabilizing it in the catalytically active conformation. In vitro biochemical studies confirm that activation of the protease in an oxidative environment is dependent on this reversible disulfide bond. These results suggest that the disulfide bond modulates activity of Rv3671c depending on the oxidative environment in vivo.

CAPRI targets T29-T42: proving ground for new docking procedures

The critical assessment of protein interactions (CAPRI) experiment provides a unique opportunity for unbiased assessment of docking procedures. The recent CAPRI targets T29-T42 entailed docking of bound, unbound, and modeled structures, presenting a wide range of prediction difficulty. We submitted accurate predictions for targets T40, T41, and T42, a good prediction for T32 and acceptable predictions for T29 and T34. The accuracy of our docking results generally matched the prediction difficulty; hence, docking of modeled proteins produced less accurate results. However, there were interesting exceptions: an accurate prediction was submitted for the dimer of modeled tetratricopeptide repeat (T42) and only an acceptable prediction for the bound/unbound case T29. The ensembles of docking models produced in the scans included an acceptable or better prediction for every target. We show here that our recently developed postscan reevaluation procedure, which tests propensity and solvation measures of the whole interface and the interface core, successfully distinguished these predictions from false docking models. For enzyme-inhibitor targets, we show that the distance of the interface from the enzyme's centroid ranked high native like docking models. Also, for one case we demonstrate that docking of an ensemble of conformers produced by normal modes analysis can improve the accuracy of the prediction.

Protease inhibitors as models for the study of oxidative folding

The correct balance between proteases and their natural protein inhibitors is of great importance in living systems. Protease inhibitors usually comprise small folds that are crosslinked by a high number of disulfide bonds, making them perfect models for the study of oxidative folding. To date, the oxidative folding of numerous protease inhibitors has been analyzed, revealing a great diversity of folding pathways that differ mainly in the heterogeneity and native disulfide-bond content of their intermediates. The two extremes of this diversity are represented by bovine pancreatic trypsin inhibitor and hirudin, which fold, respectively, via few native intermediates and heterogeneous scrambled isomers. Other proteins, such as leech carboxypeptidase inhibitor, share characteristics of both models displaying mixed folding pathways. The study of the oxidative folding of two-domain inhibitors, such as secretory leukocyte protease inhibitor, tick carboxypeptidase inhibitor, and Ascaris carboxypeptidase inhibitor, has provided some clues about how two-domain protease inhibitors may fold, that is, either by folding each domain autonomously or with one domain assisting in the folding of the other. Finally, the recent determination of the structures of the major intermediates of protease inhibitors has shed light on the molecular mechanisms guiding the oxidative folding of small disulfide-rich proteins.

MDockPP: A hierarchical approach for protein-protein docking and its application to CAPRI rounds 15-19

A hierarchical approach has been developed for protein-protein docking. In the first step, a Fast Fourier Transform (FFT)-based docking algorithm is used to globally sample all putative binding modes, in which the protein is represented by a reduced model, that is, each side chain on the protein surface is represented by its center of mass. Compared to conventional FFT docking with all-atom models, the FFT docking method with a reduced model is expected to generate more hits because it allows larger side-chain flexibility. Next, the filtered binding modes (normally several thousands) are refined by an iteratively derived knowledge-based scoring function ITScorePP and by considering backbone/loop flexibility using an ensemble docking algorithm. The distance-dependent potentials of ITScorePP were extracted by a physics-based iterative method, which circumvents the long-standing reference state problem in the knowledge-based approaches. With this hierarchical protocol, we have participated in the CAPRI experiments for Rounds 15-19 of 11 targets (T32-T42). In the predictor experiments, we achieved correct binding modes for six targets: three are with high accuracy (T40 for both distinct binding modes, T41, and T42), two are with medium accuracy (T34 and T37), and one is acceptable (T32). In the scorer experiments, of the seven target complexes that contain at least one acceptable mode submitted by the CAPRI predictor groups, we obtained correct binding modes for four targets: three are with high accuracy (T37, T40, and T41) and one is with medium accuracy (T34), suggesting good accuracy and robustness of ITScorePP.