Elucidation of the mode of interaction of thermolysin with a proteinaceous metalloproteinase inhibitor, SMPI, based on a model complex structure and a structural dynamics analysis

Synthesis, experimental evaluation and molecular modelling of hydroxamate derivatives as zinc metalloproteinase inhibitors

Enzymes of the M4 family of zinc-metalloproteinases are virulence factors secreted from gram-positive or gram-negative bacteria, and putative drug targets in the treatment of bacterial infections. In order to have a therapeutic value such inhibitors should not interfere with endogenous zinc-metalloproteinases. In the present study we have synthesised a series of hydroxamate derivatives and validated the compounds as inhibitors of the M4 enzymes thermolysin and pseudolysin, and the endogenous metalloproteinases ADAM-17, MMP-2 and MMP-9 using experimental binding studies and molecular modelling. In general, the compounds are stronger inhibitors of the MMPs than of the M4 enzymes, however, an interesting exception is LM2. The compounds bound stronger to pseudolysin than to thermolysin, and the molecular modelling studies showed that occupation of the S2(') subpocket by an aromatic group is favourable for strong interactions with pseudolysin.

Inhibition of pseudolysin and thermolysin by hydroxamate-based MMP inhibitors

In the present study, we have investigated the inhibition of thermolysin and pseudolysin by a series of compounds previously identified as matrix metalloproteinase (MMP) inhibitors using experimental binding studies and theoretical calculations. The experimental studies showed that some of the compounds were able to inhibit thermolysin and pseudolysin in the low μM range. The studies revealed that, in general, the compounds bound in the order MMPs > pseudolysin > thermolysin, and the strongest pseudolysin and thermolysin binders were compounds 8-12. Furthermore, compounds 8 and 9 were unique in that they bound much stronger to the two bacterial enzymes than to the MMPs. The docking calculations suggested that the phenyl group of the strongest binders (compounds 8 and 9) occupy the S2(')-subpocket, while a second ring system occupy the S1-subpocket in both thermolysin and pseudolysin. When the compounds possess two ring systems, the largest and most electron rich ring system seems to occupy the S1-subpocket.

Prokaryote-derived protein inhibitors of peptidases: A sketchy occurrence and mostly unknown function

In metazoan organisms protein inhibitors of peptidases are important factors essential for regulation of proteolytic activity. In vertebrates genes encoding peptidase inhibitors constitute up to 1% of genes reflecting a need for tight and specific control of proteolysis especially in extracellular body fluids. In stark contrast unicellular organisms, both prokaryotic and eukaryotic consistently contain only few, if any, genes coding for putative peptidase inhibitors. This may seem perplexing in the light of the fact that these organisms produce large numbers of proteases of different catalytic classes with the genes constituting up to 6% of the total gene count with the average being about 3%. Apparently, however, a unicellular life-style is fully compatible with other mechanisms of regulation of proteolysis and does not require protein inhibitors to control their intracellular and extracellular proteolytic activity. So in prokaryotes occurrence of genes encoding different types of peptidase inhibitors is infrequent and often scattered among phylogenetically distinct orders or even phyla of microbiota. Genes encoding proteins homologous to alpha-2-macroglobulin (family I39), serine carboxypeptidase Y inhibitor (family I51), alpha-1-peptidase inhibitor (family I4) and ecotin (family I11) are the most frequently represented in Bacteria. Although several of these gene products were shown to possess inhibitory activity, with an exception of ecotin and staphostatins, the biological function of microbial inhibitors is unclear. In this review we present distribution of protein inhibitors from different families among prokaryotes, describe their mode of action and hypothesize on their role in microbial physiology and interactions with hosts and environment.

Molecular insight into pseudolysin inhibition using the MM-PBSA and LIE methods

Pseudolysin, the extracellullar elastase of Pseudomonas aeruginosa (EC: 3.4.24.26) plays an important role in the pathogenesis of P. aeruginosa infections. In the present study, molecular dynamics simulations and theoretical affinity predictions were used to gain molecular insight into pseudolysin inhibition. Four low molecular weight inhibitors were docked at their putative binding sites and molecular dynamics (MD) simulations were performed for 5.0 ns, and the free energy of binding was calculated by the linear interaction energy method. The number and the contact surface area of stabilizing hydrophobic, aromatic, and hydrogen bonding interactions appears to reflect the affinity differences between the inhibitors. The proteinaceous inhibitor, Streptomyces metalloproteinase inhibitor (SMPI) was docked in three different binding positions and MD simulations were performed for 3.0 ns. The MD trajectories were used for molecular mechanics-Poisson-Boltzmann surface area analysis of the three binding positions. Computational alanine scanning of the average pseudolysin-SMPI complexes after MD revealed residues at the pseudolysin-SMPI interface giving the main contribution to the free energy of binding. The calculations indicated that SMPI interacts with pseudolysin via the rigid active site loop, but that also contact sites outside this loop contribute significantly to the free energy of association.

Comparative sequence and structure analysis reveal features of cold adaptation of an enzyme in the thermolysin family

Knowledge about the structural features underlying cold adaptation is important for designing enzymes of different industrial relevance. Vibriolysin from Antarctic bacterium strain 643 (VAB) is at present the only enzyme of the thermolysin family from an organism that thrive in extremely cold climate. In this study comparative sequence-structure analysis and molecular dynamics (MD) simulations were used to reveal the molecular features of cold adaptation of VAB. Amino acid sequence analysis of 44 thermolysin enzymes showed that VAB compared to the other enzymes has: (1) fewer arginines, (2) a lower Arg/(Lys + Arg) ratio, (3) a lower fraction of large aliphatic side chains, expressed by the (Ile + Leu)/(Ile + Leu + Val) ratio, (4) more methionines, (5) more serines, and (6) more of the thermolabile amino acid asparagine. A model of the catalytic domain of VAB was constructed based on homology with pseudolysin. MD simulations for 3 ns of VAB, pseudolysin, and thermolysin supported the assumption that cold-adapted enzymes have a more flexible three-dimensional (3D) structure than their thermophilic and mesophilic counterparts, especially in some loop regions. The structural analysis indicated that VAB has fewer intramolecular cation-pi electron interactions and fewer hydrogen bonds than its mesophilic (pseudolysin) and thermophilic (thermolysin) counterparts. Lysine is the dominating cationic amino acids involved in salt bridges in VAB, while arginine is dominating in thermolysin and pseudolysin. VAB has a greater volume of inaccessible cavities than pseudolysin and thermolysin. The electrostatic potentials on the surface of the catalytic domain were also more negative for VAB than for thermolysin and pseudolysin. Thus, the MD simulations, the structural patterns, and the amino acid composition of VAB relative to other enzymes of the thermolysin family suggest that VAB possesses the biophysical properties generally following adaptation to cold climate.

The Protein-protein Interactions Between SMPI and Thermolysin Studied by Molecular Dynamics and MM/PBSA Calculations

Thermolysin is a zinc-metalloendopeptidase secreted by the gram-positive thermophilic bacterium Bacillus thermoproteolyticus. Thermolysin belongs to the gluzinicin family of enzymes, which is selectively inhibited by Steptomyces metalloproteinase inhibitor (SMPI). Very little is known about the interaction between SMPI and thermolysin. Knowledge about the protein-protein interactions is very important for designing new thermolysin inhibitors with possible industrial or pharmaceutical applications. In the present study, two binding modes between SMPI and thermolysin were studied by 2300 picoseconds (ps) of comparative molecular dynamics (MD) simulations and calculation of the free energy of binding using the molecular mechanics-Poisson-Boltmann surface area (MM/PBSA) method. One of the positions, the 'horizontal arrow head docking' (HAHD) was similar to the previously proposed binding mode by Tate et al. (Tate, S., Ohno, A., Seeram, S. S., Hiraga, K., Oda, K., and Kainosho, M. J. Mol. Biol. 282, 435-446 (1998)). The other position, the 'vertical arrow head docking' (VAHD) was obtained by a manual docking guided by the shape and charge distribution of SMPI and the binding pocket of thermolysin. The calculations showed that SMPI had stronger interactions with thermolysin in the VAHD than in the HAHD complex, and the VAHD complex was considered more realistic than the HAHD complex. SMPI interacted with thermolysin not only at the active site but had auxiliary binding sites contributing to proper interactions. The VAHD complex can be used for designing small molecule inhibitors mimicking the SMPI-thermolysin binding interfaces.

The many faces of protease-protein inhibitor interaction

Proteases and their natural protein inhibitors are among the most intensively studied protein-protein complexes. There are about 30 structurally distinct inhibitor families that are able to block serine, cysteine, metallo- and aspartyl proteases. The mechanisms of inhibition can be related to the catalytic mechanism of protease action or include a mechanism-unrelated steric blockage of the active site or its neighborhood. The structural elements that are responsible for the inhibition most often include the N- or the C-terminus or exposed loop(s) either separately or in combination of several such elements. During complex formation, no major conformational changes are usually observed, but sometimes structural transitions of the inhibitor and enzyme occur. In many cases, convergent evolution, with respect to the inhibitors' parts that are responsible for the inhibition, can be inferred from comparisons of their structures or sequences, strongly suggesting that there are only limited ways to inhibit proteases by proteins.

NMR structure of Streptomyces killer toxin-like protein, SKLP: further evidence for the wide distribution of single-domain betagamma-crystallin superfamily proteins

A protein isolated from the culture supernatant of the soil bacterium, Streptomyces sp. F-287, exhibits cytocidal effects for both budding and fission yeasts, and causes morphological changes of yeasts and filamentous fungi. This protein, which was the first killer toxin-like protein for yeasts identified in the Streptomyces microorganism, was named SKLP (Streptomyces killer toxin-like protein). Since the amino acid sequence of the protein, as determined by sequential Edman degradations, seemed to be unique, we determined the structure by NMR spectroscopy. Although the actual target of SKLP in yeasts has not been determined yet, the structure might give us a clue to characterize the targets. The solution structure of SKLP determined by NMR, however, turned out to be a single-domain crystallin-like protein, with two Greek key motifs and a short extra beta-strand at the N terminus. The final ensemble of 20 NMR structures overlaid onto their mean coordinate with rmsd values of 0.32(+/-0.06) A for the backbone atoms involved in the secondary structure elements. As a yeast killer toxin, WmKT, isolated from the yeast strain Williopsis mrakii also has a Greek key beta-barrel fold, we have made a detailed comparison of the structural features of SKLP with the other crystallin superfamily proteins. It is very interesting that SKLP has a unique electrostatic potential distribution on the molecular surface. Namely, one surface of the beta-barrel fold in SKLP has a large negatively charged region, with an isolated positive charge of the Arg62 side-chain at the center. The edge of this surface is surrounded by positively charged residues, including Arg31, Arg65 and Arg74. The salient features of the charge distribution on this surface and the cluster of Arg residues might be related to the target binding of SKLP.

Analysis of the internal motion of free and ligand-bound human lysozyme by use of 15N NMR relaxation measurement: a comparison with those of hen lysozyme

Human lysozyme has a structure similar to that of hen lysozyme and differs in amino acid sequence by 51 out of 129 residues with one insertion at the position between 47 and 48 in hen lysozyme. The backbone dynamics of free or (NAG)3-bound human lysozyme has been determined by measurements of 15N nuclear relaxation. The relaxation data were analyzed using the Lipari-Szabo formalism and were compared with those of hen lysozyme, which was already reported (Mine S et al.. 1999, J Mol Biol 286:1547-1565). In this paper, it was found that the backbone dynamics of free human and hen lysozymes showed very similar behavior except for some residues, indicating that the difference in amino acid sequence did not affect the behavior of entire backbone dynamics, but the folded pattern was the major determinant of the internal motion of lysozymes. On the other hand, it was also found that the number of residues in (NAG)3-bound human and hen lysozymes showed an increase or decrease in the order parameters at or near active sites on the binding of (NAG)3, indicating the increase in picosecond to nanosecond. These results suggested that the immobilization of residues upon binding (NAG)3 resulted in an entropy penalty and that this penalty was compensated by mobilizing other residues. However, compared with the internal motions between both ligand-bound human and hen lysozymes, differences in dynamic behavior between them were found at substrate binding sites, reflecting a subtle difference in the substrate-binding mode or efficiency of activity between them.

Structural basis of the endoproteinase-protein inhibitor interaction

Proteolytic enzymes are potentially hazardous to their protein environment, so that their activity must be carefully controlled. Living organisms use protein inhibitors as a major tool to regulate the proteolytic activity of proteinases. Most of the inhibitors for which 3D structures are available are directed towards serine proteinases, interacting with the active sites in a 'canonical' i.e. substrate-like manner via an exposed reactive site loop of conserved conformation. More recently, some non-canonically binding serine proteinase inhibitors directed against coagulation factors, in particular thrombin, a few cysteine proteinase inhibitors inhibitory towards papain-like proteinases, and three zinc endopeptidase inhibitors directed against metzincins and thermolysin have been characterised in the free and complexed state, displaying novel mechanisms of inhibition with their target proteinases. These different interaction modes are presented and briefly discussed with respect to the different strategies applied by nature.

Analysis of the relationship between enzyme activity and its internal motion using nuclear magnetic resonance: 15N relaxation studies of wild-type and mutant lysozyme

A mutant lysozyme where R14 and H15 are deleted together has higher activity and a similar binding ability to an inhibitor, trimer of N-acetylglucosamine ((NAG)3), compared with wild-type lysozyme. Since this has been attributed to intrinsic protein dynamic properties, we investigated the relationship between the activity and the internal motions of proteins. Backbone dynamics of the free and the complex forms with the (NAG)3 have been studied by measurement of the 15N T1 and T2 relaxation rates and NOE determinations at 600 MHz. Analysis of the data using the model-free formalism showed that the generalized order parameters (S2) were almost the same in wild-type and mutant lysozyme in unbound state, indicating that the mutation had little effect on the global internal motions. On the other hand, in the presence of (NAG)3, although some signals located around the active site were broadened or decreased in intensity because of strong perturbation by (NAG)3, there were several residues that showed increased or decreased backbone S2 in the complexed lysozymes. A comparison of the internal motions of the wild-type and mutant complexes showed a number of distinct dynamic differences between them. In particular, many residues located at or near active-site regions (turn 1, strand 2, turn 2 and long loop), displayed greater backbone dynamics reflecting the order parameter in mutant complex relative to mutant free. Furthermore, the Rex values at the loop C-D region, which was considered to be important for enzymatic activity, significantly increased. From these results, it was suggested that variations in the dynamics of these regions may play an important role in the enzyme activity.

Studies on the formation and stability of a complex between Streptomyces proteinaceous metalloprotease inhibitor and thermolysin

The effects of certain physicochemical parameters on the formation and stability of a complex between Streptomyces proteinaceous metalloprotease inhibitor (SMPI) and thermolysin were investigated. SMPI had its lowest Ki value at a pH of around 6.5 (similar to the pH dependence of the kcat/K(m) of thermolysin catalysis), reflecting the splitting mechanism of the SMPI inhibition of thermolysin. This Ki increased with an increase in pressure, and in (Ki-1) was almost linear with respect to pressure. The volume of the reaction (delta Vcomp), which is the volume change accompanying enzyme-inhibitor complex formation, was calculated as +8.1 +/- 0.3 mL.mol-1, which has a sign opposite to delta Vcomp for neutral peptide inhibitors and acyl-peptide substrates. The temperature dependence of Ki-1 gave the reaction enthalpy (delta Hcomp) and reaction entropy (delta Scomp) of the complex formation as 34.6 +/- 1.4 kJ.mol-1 and 298 +/- 5 J.mol-1.K-1, respectively. These positive reaction volumes and reaction entropies were related to the electrostatic interactions and ionic strength dependence of Ki which corresponded to the key ionic interaction during complex formation. Complex formation with SMPI stabilized thermolysin against pressure perturbation as observed by the changes in the Trp fluorescence of thermolysin with increasing pressure. Thermal stability, however, was affected very little by complex formation with SMPI. Phosphoramidon, Cbz-Phe-Gly-NH2 and Cbz-Phe also positively affected the pressure-tolerance of thermolysin, in the following order: Cbz-Gly-Phe-NH2 < Cbz-Phe << phosphoramidon. The third compound exhibited stabilizing effects comparable with those of SMPI, which suggests that the interaction between SMPI and thermolysin was localized to the reactive site.

NMR structure of the Streptomyces metalloproteinase inhibitor, SMPI, isolated from Streptomyces nigrescens TK-23: another example of an ancestral beta gamma-crystallin precursor structure

The Streptomyces metalloproteinase inhibitor, SMPI, isolated from Streptomyces nigrescens TK-23, is a proteinaceous metalloproteinase inhibitor, and consists of 102 amino acid residues with two disulfide bridges. SMPI specifically inhibits metalloproteinases such as thermolysin. In the present work, the solution structure of SMPI was determined on the basis of 1536 nuclear Overhauser enhancement derived distance restraints and 52 dihedral angle restraints obtained from three-bond spin coupling constants. The final ensemble of 20 NMR structures overlaid onto their mean coordinate with backbone (N, Calpha, C') r.m.s.d. values of 0. 45(+/-0.11) A and 0.57(+/-0.18) A for residues 6 to 99 and the entire 102 residues, respectively. SMPI is essentially composed of two beta-sheets, each consisting of four antiparallel beta-strands. The structure can be considered as two Greek key motifs with 2-fold internal symmetry, a Greek key beta-barrel. One unique structural feature found in SMPI is in its extension between the first and second strands of the second Greek key motif. Interestingly, this extended segment is known to be involved in the inhibitory activity of SMPI. In the absence of sequence similarity, the SMPI structure shows clear similarity to both domains of the eye lens crystallins, both domains of the calcium sensor protein-S, as well as the single-domain yeast killer toxin. The yeast killer toxin structure was thought to be a precursor of the two-domain beta gamma-crystallin proteins, because of its structural similarity to each domain of the beta gamma-crystallins. SMPI thus provides another example of a single-domain protein structure that corresponds to the ancestral fold from which the two-domain proteins in the beta gamma-crystallin superfamily are believed to have evolved.