Interscaffolding additivity. Association of P1 variants of eglin c and of turkey ovomucoid third domain with serine proteinases

Identification of Kazal Inhibitor Scaffolds with Identical Canonical Binding Loops and Their Effects on Binding Properties

Kazal inhibitors hold high potential as scaffolds for therapeutic molecules, taking advantage of the easily exchangeable canonical binding loop. Different Kazal inhibitor backbones have been suggested to be therapeutically useful, but the impact of different Kazal-like scaffolds on binding properties is still largely unknown. Here, we identified trypsin-targeting human serine protease inhibitor Kazal type 1 (SPINK1) homologues in different mammalian species that cluster in two P2-P1 combinations, implying the coevolution of these residues. We generated loop exchange variants of human SPINK1 for comparison with Kazal inhibitors from related species. Using comprehensive biophysical characterization of the inhibitor-enzyme interactions, we found not only affinity but also pH resistance to be highly backbone-dependent. Differences are mostly observed in complex stability, which varies by over one order of magnitude. We provide clear evidence for high backbone dependency within the Kazal family. Hence, when designing Kazal inhibitor-based therapeutic molecules, testing different backbones after optimizing the canonical binding loop can be beneficial and may result in increased affinity, complex stability, specificity, and pH resistance.

Directed Evolution-Driven Increase of Structural Plasticity Is a Prerequisite for Binding the Complement Lectin Pathway Blocking MASP-Inhibitor Peptides

MASP-1 and MASP-2 are key activator proteases of the complement lectin pathway. The first specific mannose-binding lectin-associated serine protease (MASP) inhibitors had been developed from the 14-amino-acid sunflower trypsin inhibitor (SFTI) peptide by phage display, yielding SFTI-based MASP inhibitors, SFMIs. Here, we present the crystal structure of the MASP-1/SFMI1 complex that we analyzed in comparison to other existing MASP-1/2 structures. Rigidified backbone structure has long been accepted as a structural prerequisite for peptide inhibitors of proteases. We found that a hydrophobic cluster organized around the P2 Thr residue is essential for the structural stability of wild-type SFTI. We also found that the same P2 Thr prevents binding of the rigid SFTI-like peptides to the substrate-binding cleft of both MASPs as the cleft is partially blocked by large gatekeeper enzyme loops. Directed evolution removed this obstacle by replacing the P2 Thr with a Ser, providing the SFMIs with high-degree structural plasticity, which proved to be essential for MASP inhibition. To gain more insight into the structural criteria for SFMI-based MASP-2 inhibition, we systematically modified MASP-2-specific SFMI2 by capping its two termini and by replacing its disulfide bridge with varying length thioether linkers. By doing so, we also aimed to generate a versatile scaffold that is resistant to reducing environment and has increased stability in exopeptidase-containing biological environments. We found that the reduction-resistant disulfide-substituted l-2,3-diaminopropionic acid (Dap) variant possessed near-native potency. As MASP-2 is involved in the life-threatening thrombosis in COVID-19 patients, our synthetic, selective MASP-2 inhibitors could be relevant coronavirus drug candidates.

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.

Novel MASP-2 inhibitors developed via directed evolution of human TFPI1 are potent lectin pathway inhibitors

The lectin pathway (LP) of the complement system is an important antimicrobial defense mechanism, but it also contributes significantly to ischemia reperfusion injury (IRI) associated with myocardial infarct, stroke, and several other clinical conditions. Mannan-binding lectin-associated serine proteinase 2 (MASP-2) is essential for LP activation, and therefore, it is a potential drug target. We have previously developed the first two generations of MASP-2 inhibitors by in vitro evolution of two unrelated canonical serine proteinase inhibitors. These inhibitors were selective LP inhibitors, but their nonhuman origin rendered them suboptimal lead molecules for drug development. Here, we present our third-generation MASP-2 inhibitors that were developed based on a human inhibitor scaffold. We subjected the second Kunitz domain of human tissue factor pathway inhibitor 1 (TFPI1 D2) to directed evolution using phage display to yield inhibitors against human and rat MASP-2. These novel TFPI1-based MASP-2 inhibitor (TFMI-2) variants are potent and selective LP inhibitors in both human and rat serum. Directed evolution of the first Kunitz domain of TFPI1 had already yielded the potent kallikrein inhibitor, Kalbitor® (ecallantide), which is an FDA-approved drug to treat acute attacks of hereditary angioedema. Like hereditary angioedema, acute IRI is also related to the uncontrolled activation of a specific plasma serine proteinase. Therefore, TFMI-2 variants are promising lead molecules for drug development against IRI.

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.

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.

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

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.

Recombinant protease inhibitors for herbivore pest control: a multitrophic perspective

Protease inhibitors are a promising complement to Bt toxins for the development of insect-resistant transgenic crops, but their limited specificity against proteolytic enzymes and the ubiquity of protease-dependent processes in living organisms raise questions about their eventual non-target effects in agroecosystems. After a brief overview of the main factors driving the impacts of insect-resistant transgenic crops on non-target organisms, the possible effects of protease inhibitors are discussed from a multitrophic perspective, taking into account not only the target herbivore proteases but also the proteases of other organisms found along the trophic chain, including the plant itself. Major progress has been achieved in recent years towards the design of highly potent broad-spectrum inhibitors and the field deployment of protease inhibitor-expressing transgenic plants resistant to major herbivore pests. A thorough assessment of the current literature suggests that, whereas the non-specific inhibitory effects of recombinant protease inhibitors in plant food webs could often be negligible and their 'unintended' pleiotropic effects in planta of potential agronomic value, the innocuity of these proteins might always remain an issue to be assessed empirically, on a case-by-case basis.

Determinants of affinity and proteolytic stability in interactions of Kunitz family protease inhibitors with mesotrypsin

An important functional property of protein protease inhibitors is their stability to proteolysis. Mesotrypsin is a human trypsin that has been implicated in the proteolytic inactivation of several protein protease inhibitors. We have found that bovine pancreatic trypsin inhibitor (BPTI), a Kunitz protease inhibitor, inhibits mesotrypsin very weakly and is slowly proteolyzed, whereas, despite close sequence and structural homology, the Kunitz protease inhibitor domain of the amyloid precursor protein (APPI) binds to mesotrypsin 100 times more tightly and is cleaved 300 times more rapidly. To define features responsible for these differences, we have assessed the binding and cleavage by mesotrypsin of APPI and BPTI reciprocally mutated at two nonidentical residues that make direct contact with the enzyme. We find that Arg at P(1) (versus Lys) favors both tighter binding and more rapid cleavage, whereas Met (versus Arg) at P'(2) favors tighter binding but has minimal effect on cleavage. Surprisingly, we find that the APPI scaffold greatly enhances proteolytic cleavage rates, independently of the binding loop. We draw thermodynamic additivity cycles analyzing the interdependence of P(1) and P'(2) substitutions and scaffold differences, finding multiple instances in which the contributions of these features are nonadditive. We also report the crystal structure of the mesotrypsin·APPI complex, in which we find that the binding loop of APPI displays evidence of increased mobility compared with BPTI. Our data suggest that the enhanced vulnerability of APPI to mesotrypsin cleavage may derive from sequence differences in the scaffold that propagate increased flexibility and mobility to the binding loop.

Systematic mutational analysis of peptide inhibition of the p53-MDM2/MDMX interactions

Inhibition of the interaction between the tumor suppressor protein p53 and its negative regulators MDM2 and MDMX is of great interest in cancer biology and drug design. We previously reported a potent duodecimal peptide inhibitor, termed PMI (TSFAEYWNLLSP), of the p53-MDM2 and -MDMX interactions. PMI competes with p53 for MDM2 and MDMX binding at an affinity roughly 2 orders of magnitude higher than that of (17-28)p53 (ETFSDLWKLLPE) of the same length; both peptides adopt nearly identical alpha-helical conformations in the complexes, where the three highlighted hydrophobic residues Phe, Trp, and Leu dominate PMI or (17-28)p53 binding to MDM2 and MDMX. To elucidate the molecular determinants for PMI activity and specificity, we performed a systematic Ala scanning mutational analysis of PMI and (17-28)p53. The binding affinities for MDM2 and MDMX of a total of 35 peptides including 10 truncation analogs were quantified, affording a complete dissection of energetic contributions of individual residues of PMI and (17-28)p53 to MDM2 and MDMX association. Importantly, the N8A mutation turned PMI into the most potent dual-specific antagonist of MDM2 and MDMX reported to date, registering respective K(d) values of 490 pM and 2.4 nM. The co-crystal structure of N8A-PMI-(25-109)MDM2 was determined at 1.95 A, affirming that high-affinity peptide binding to MDM2/MDMX necessitates, in addition to optimized intermolecular interactions, enhanced helix stability or propensity contributed by non-contact residues. The powerful empirical binding data and crystal structures present a unique opportunity for computational studies of peptide inhibition of the p53-MDM2/MDMX interactions.

Insights into affinity and specificity in the complexes of alpha-lytic protease and its inhibitor proteins: binding free energy from molecular dynamics simulation

We report the binding free energy calculation and its decomposition for the complexes of alpha-lytic protease and its protein inhibitors using molecular dynamics simulation. Standard mechanism serine protease inhibitors eglin C and OMTKY3 are known to have strong binding affinity for many serine proteases. Their binding loops have significant similarities, including a common P1 Leu as the main anchor in the binding interface. However, recent experiments demonstrate that the two inhibitors have vastly different affinity towards alpha-lytic protease (ALP), a bacterial serine protease. OMTKY3 inhibits the enzyme much more weakly (by approximately 10(6) times) than eglin C. Moreover, a variant of OMTKY3 with five mutations, OMTKY3M, has been shown to inhibit 10(4) times more strongly than the wild-type inhibitor. The underlying mechanisms for the unusually large difference in binding affinities and the effect of mutation are not well understood. Here we use molecular dynamics simulation with molecular mechanics-Poisson Boltzmann/surface area method (MM-PB/SA) to investigate quantitatively the binding specificity. The calculated absolute binding free energies correctly differentiate the thermodynamic stabilities of these protein complexes, but the magnitudes of the binding affinities are systematically overestimated. Analysis of the binding free energy components provides insights into the molecular mechanism of binding specificity. The large DeltaDeltaG(bind) between eglin C and wild type OMTKY3 towards ALP is mainly attributable to the stronger nonpolar interactions in the ALP-eglin C complex, arising from a higher degree of structural complementarity. Here the electrostatic interaction contributes to a lesser extent. The enhanced inhibition in the penta-mutant OMTKY3M over its wild type is entirely due to an overall improvement in the solvent-mediated electrostatic interactions in the ALP-OMTKY3M complex. The results suggest that for these protein-complexes and similar enzyme-inhibitor systems (1) the binding is driven by nonpolar interactions, opposed by overall electrostatic and solute entropy contributions; (2) binding specificity can be tuned by improving the complementarity in electrostatics between two associating proteins. Binding free energy decomposition into contributions from individual protein residues provides additional detailed information on the structural determinants and subtle conformational changes responsible for the binding specificity.

Inhibition mechanism and the effects of structure on activity of male reproduction-related peptidase inhibitor Kazal-type (MRPINK) of Macrobrachium rosenbergii

In our previous reports, a Kazal family serine protease inhibitor, male reproduction-related peptidase inhibitor Kazal-type (MRPINK) has been identified from the prawn, Macrobrachium rosenbergii, and discovered having an inhibitory effect on the sperm gelatinolytic activity. MRPINK was predicated to inhibit chymotrypsin since it contains leucine and proline at P(1) positions of the two domains, respectively. In this report, recombinant MRPINK was as expected found to specifically inhibit chymotrypsin, but no inhibition was detected against trypsin or thrombin. By the analysis of kinetic tests, the inhibition mechanism of MRPINK was determined to be typical competitive model with K (i) of 354 nM. To elucidate the effects of structure on activity of MRPINK, the mutants (domain-1 only, domain-2 only, MRPINK(P88I), MRPINK(L37K), MRPINK(L37A), and MRPINK(L37G)) were prepared and their inhibitory activities assayed. The results showed that domain-2 was the key contributor to the inhibition of chymotrypsin (K (i) of 416 nM) and P(1) Pro was crucial for the activity. Nevertheless, whether the P(1) amino acid residue was Leu, or even if it was replaced by Lys, Ala, or Gly, domain-1 was ineffective to the activity.

Tailoring the specificity of a plant cystatin toward herbivorous insect digestive cysteine proteases by single mutations at positively selected amino acid sites

Plant cystatins, similar to other defense proteins, include hypervariable, positively selected amino acid sites presumably impacting their biological activity. Using 29 single mutants of the eighth domain of tomato (Solanum lycopersicum) multicystatin, SlCYS8, we assessed here the potential of site-directed mutagenesis at positively selected amino acid sites to generate cystatin variants with improved inhibitory potency and specificity toward herbivorous insect digestive cysteine (Cys) proteases. Compared to SlCYS8, several mutants (22 out of 29) exhibited either improved or lowered potency against different model Cys proteases, strongly suggesting the potential of positively selected amino acids as target sites to modulate the inhibitory specificity of the cystatin toward Cys proteases of agronomic significance. Accordingly, mutations at positively selected sites strongly influenced the inhibitory potency of SlCYS8 against digestive Cys proteases of the insect herbivore Colorado potato beetle (Leptinotarsa decemlineata). In particular, several variants exhibited improved potency against both cystatin-sensitive and cystatin-insensitive digestive Cys proteases of this insect. Of these, some variants also showed weaker activity against leaf Cys proteases of the host plant (potato [Solanum tuberosum]) and against a major digestive Cys protease of the two-spotted stinkbug Perillus bioculatus, an insect predator of Colorado potato beetle showing potential for biological control. Overall, these observations suggest the usefulness of site-directed mutagenesis at positively selected amino acid sites for the engineering of recombinant cystatins with both improved inhibitory potency toward the digestive proteases of target herbivores and weaker potency against nontarget Cys proteases in the host plant or the environment.

When the Surface Tells What Lies Beneath: Combinatorial Phage-display Mutagenesis Reveals Complex Networks of Surface–Core Interactions in the Pacifastin Protease Inhibitor Family

Pacifastin protease inhibitors are small cysteine-rich motifs of approximately 35 residues that were discovered in arthropods. The family is divided into two related groups on the basis of the composition of their minimalist inner core. In group I, the core is governed by a Lys10-Trp26 interaction, while in group II it is organized around Phe10. Group I inhibitors exhibit intriguing taxon specificity: potent arthropod-trypsin inhibitors from this group are almost inactive against vertebrate enzymes. The group I member SGPI-1 and the group II member SGPI-2 are extensively studied inhibitors. SGPI-1 is taxon-selective, while SGPI-2 is not. Individual mutations failed to explain the causes underlying this difference. We deciphered this phenomenon using comprehensive combinatorial mutagenesis and phage display. We produced a complete chimeric SGPI-1 / SGPI-2 inhibitor-phage library, in which the two sequences were shuffled at the highest possible resolution of individual residues. The library was selected for binding to bovine trypsin and crayfish trypsin. Sequence analysis of the selectants revealed that taxon specificity is due to an intra-molecular functional coupling between a surface loop and the Lys10-Trp26 core. Five SGPI-2 surface residues transplanted into SGPI-1 resulted in a variant that retained the "taxon-specific" core, but potently inhibited both vertebrate and arthropod enzymes. An additional rational point mutation resulted in a picomolar inhibitor of both trypsins. Our results challenge the generally accepted view that surface residues are the exclusive source of selectivity for canonical inhibitors. Moreover, we provide important insights into general principles underlying the structure-function properties of small disulfide-rich polypeptides, molecules that exist at the borderline between peptides and proteins.

Screened nonbonded interactions in native proteins manipulate optimal paths for robust residue communication

A protein structure is represented as a network of residues whereby edges are determined by intramolecular contacts. We introduce inhomogeneity into these networks by assigning each edge a weight that is determined by amino acid pair potentials. Two methodologies are utilized to calculate the average path lengths (APLs) between pairs: to minimize i), the maximum weight in the strong APL, and ii), the total weight in the weak APL. We systematically screen edges that have higher than a cutoff potential and calculate the shortest APLs in these reduced networks, while keeping chain connectivity. Therefore, perturbations introduced at a selected region of the residue network propagate to remote regions only along the nonscreened edges that retain their ability to disseminate the perturbation. The shortest APLs computed from the reduced homogeneous networks with only the strongest few nonbonded pairs closely reproduce the strong APLs from the weighted networks. The rate of change in the APL in the reduced residue network as compared to its randomly connected counterpart remains constant until a lower bound. Upon further link removal, this property shows an abrupt increase toward a random coil behavior. Under different perturbation scenarios, diverse optimal paths emerge for robust residue communication.

Modulating the proteinase inhibitory profile of a plant cystatin by single mutations at positively selected amino acid sites

Cysteine proteinase inhibitors of the cystatin superfamily have several important functions in plants, including the inhibition of exogenous cysteine proteinases during herbivory or infection. Here we used a maximum-likelihood approach to assess whether plant cystatins, like other proteins implicated in host-pest interactions, have been subject to positive selection during the course of their evolution. Several amino acid sites were identified as being positively selected in cystatins from either Poaceae (monocots) and Solanaceae (dicots). These hypervariable sites were located at strategic positions on the protein: on each side of the conserved glycine residues in the N-terminal trunk, within the first and second inhibitory loops entering the active site of target enzymes, and surrounding the larfav motif, a sequence of unknown function conserved among plant cystatins. Supporting the assumption that positively selected, hypervariable sites are indicative of amino acid sites implicated in functional diversity, mutants of the 8th cystatin unit of tomato multicystatin including alternative residues at positively selected sites in the N-terminal trunk exhibited highly variable affinities for the cysteine proteases papain, cathepsin B and cathepsin H. Overall, these observations support the hypothesis that plant cystatins have been under selective pressure to evolve in response to predatory challenges by herbivorous enemies. They also indicate the potential of site-directed mutagenesis at positively selected sites for the generation of cystatins with improved binding properties.

Predicting proteinase specificities from free energy calculations

The role of the primary binding residue (P1) in complexes between three different subtilases (subtilisin Carlsberg, thermitase and proteinase K) and their canonical protein inhibitor eglin c have been studied by free energy calculations. Based on the crystal structures of eglin c in complex with subtilisin Carlsberg and thermitase, and a homology model of the eglin c-proteinase K complex, a total of 57 mutants have been constructed and docked into their host proteins. The binding free energy was then calculated using molecular dynamics (MD) simulations combined with the linear interaction energy (LIE) method for all complexes differing only in the nature of the amino acid at the P1 position. LIE calculations for 19 different complexes for each subtilase were thus carried out excluding proline. The effects of substitutions at the P1 position on the binding free energies are found to be very large, and positively charged residues (Arg, Lys and His) are particularly deleterious for all three enzymes. The charged variants of the acidic side chains are found to bind more favorably as compared to their protonated states in all three subtilases. Furthermore, hydrophobic amino acids are accommodated most favorably at the S1-site in all three enzymes. Comparison of the three series of binding free energies shows only minor differences in the 19 computed relative binding free energies among these subtilases. This is further reflected in the correlation coefficient between the 23 relative binding free energies obtained, including the possible protonation states of ionizable side chains, but excluding the P1 Pro, for subtilisin Carlsberg versus thermitase (0.95), subtilisin versus proteinase K (0.94) and thermitase versus proteinase K (0.96).

Energetics and Cooperativity of the Hydrogen Bonding and Anchor Interactions that Bind Peptides to MHC Class II Protein

The complexity of the interaction between major histocompatibility complex class II (MHC II) proteins and peptide ligands has been revealed through structural studies and crystallographic characterization. Peptides bind through side-chain "anchor" interactions with MHC II pockets and an extensive array of genetically conserved hydrogen bonds to the peptide backbone. Here we quantitatively investigate the kinetic hierarchy of these interactions. We present results detailing the impact of single side-chain mutations of peptide anchor residues on dissociation rates, utilizing two I-A(d)-restricted peptides, one of which has a known crystal structure, and 24 natural and non-natural amino acid mutant variants of these peptides. We find that the N-terminal P1, P4 and P6 anchor-pocket interactions can make significant contributions to binding stability. We also investigate the interactions of these peptides with four I-A(d) MHC II proteins, each mutated to disrupt conserved hydrogen bonds to the peptide backbone. These complexes exhibit kinetic behavior suggesting that binding energy is disproportionately invested near the peptide N terminus for backbone hydrogen bonds. We then evaluate the effects of simultaneously modifying both anchor and hydrogen bonding interactions. A quantitative analysis of 71 double mutant cycles reveals that there is little apparent cooperativity between anchor residue interactions and hydrogen bonds, even when they are directly adjacent (<5A).

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.

Emerging therapies for hepatitis C virus infection

Hepatitis C virus (HCV) has infected millions of people worldwide and has emerged as a global health crisis. The currently available therapy is interferon (IFN) either alone or in combination with ribavirin. However, the disappointing efficacy of IFN has led to the considerable need for improved treatments and a number of new therapies are under evaluation in clinical trials. These include pegylated IFNs, which have altered physiochemical characteristics allowing once-weekly dosing. Combination of pegylated IFN with ribavirin should further improve sustained response rates. However, not all patients are successfully treated with IFNs, particularly those infected with genotype 1 of the virus, and it is likely that potent, specific drugs will be required. The majority of new approaches currently trying to combat this viral disease are aimed at inhibition of viral targets. Most effort has been directed towards inhibition of the NS3 serine protease, and potent inhibitors have now been described. However, a clinical candidate is yet to emerge against this difficult target. Considerable work by leading researchers has provided crystal structures of the key replicative enzymes, NS3 protease, NS3 helicase, NS5B polymerase and full-length NS3 protease-helicase, and there is much hope that such structural information will bear fruit. More recently, inhibition of host targets, particularly inosine monophosphate dehydrogenase (IMPDH), has become of interest and there are on-going clinical trials with such inhibitors. Research aimed at novel treatments for HCV disease is gathering pace and very recent developments in cell-based assay systems can only hasten the discovery of improved therapies.

A physical reference state unifies the structure-derived potential of mean force for protein folding and binding

Extracting knowledge-based statistical potential from known structures of proteins is proved to be a simple, effective method to obtain an approximate free-energy function. However, the different compositions of amino acid residues at the core, the surface, and the binding interface of proteins prohibited the establishment of a unified statistical potential for folding and binding despite the fact that the physical basis of the interaction (water-mediated interaction between amino acids) is the same. Recently, a physical state of ideal gas, rather than a statistically averaged state, has been used as the reference state for extracting the net interaction energy between amino acid residues of monomeric proteins. Here, we find that this monomer-based potential is more accurate than an existing all-atom knowledge-based potential trained with interfacial structures of dimers in distinguishing native complex structures from docking decoys (100% success rate vs. 52% in 21 dimer/trimer decoy sets). It is also more accurate than a recently developed semiphysical empirical free-energy functional enhanced by an orientation-dependent hydrogen-bonding potential in distinguishing native state from Rosetta docking decoys (94% success rate vs. 74% in 31 antibody-antigen and other complexes based on Z score). In addition, the monomer potential achieved a 93% success rate in distinguishing true dimeric interfaces from artificial crystal interfaces. More importantly, without additional parameters, the potential provides an accurate prediction of binding free energy of protein-peptide and protein-protein complexes (a correlation coefficient of 0.87 and a root-mean-square deviation of 1.76 kcal/mol with 69 experimental data points). This work marks a significant step toward a unified knowledge-based potential that quantitatively captures the common physical principle underlying folding and binding. A Web server for academic users, established for the prediction of binding free energy and the energy evaluation of the protein-protein complexes, may be found at http://theory.med.buffalo.edu.

Selection of low-molecular-mass trypsin and chymotrypsin inhibitors based on the binding loop of CMTI-III using combinatorial chemistry methods

Using a combinatorial chemistry approach, a decapaptide library containing the N-terminal fragment of trypsin inhibitor CMTI-III was synthesized by the solid-phase method. The peptide library was screened for trypsin and chymotrypsin inhibitory activity applying the iterative method in solution. Two decapeptides were selected and resynthesized for each enzyme. The association equilibrium constants ((1.1+/-0.2)x10(8) and (7.3+/-1.6)x10(7)) determined for peptides with trypsin inhibitory activity indicate that they are 3-4-fold less active than the CMTI inhibitors. On the other hand, they are significantly more effective as compared with the starting sequence. Two peptides selected as chymotrypsin inhibitors displayed about 10 times higher activity (1.7+/-0.4)x10(7) and (1.1+/-0.2)x10(7), respectively) than those monosubstituted in position P(1) of the CMTI-III analogue. Considering low molecular weight of peptides selected and the lack of conformational constraints in their structures, the results are promising. They are good templates as starting sequences for further selection of small, peptidomimetic proteinase inhibitors.

Structure and energetics of protein-protein interactions: the role of conformational heterogeneity in OMTKY3 binding to serine proteases

Proteins with flexible binding surfaces can interact with numerous binding partners. However, this promiscuity is more difficult to understand in "rigid-body" proteins, whose binding results in little, or no, change in the position of backbone atoms. The binding of Kazal inhibitors to serine proteases is considered a classic case of rigid-body binding, although they bind to a wide range of proteases. We have studied the thermodynamics of binding of the Kazal serine protease inhibitor, turkey ovomucoid third domain (OMTKY3), to the serine protease subtilisin Carlsberg using isothermal titration calorimetry and have determined the crystal structure of the complex at very high resolution (1.1A). Comparison of the binding energetics and structure to other OMTKY3 interactions demonstrates that small changes in the position of side-chains can make significant contributions to the binding thermodynamics, including the enthalpy of binding. These effects emphasize that small, "rigid-body" proteins are still dynamic structures, and these dynamics make contributions to both the enthalpy and entropy of binding interactions.

Additivity-based prediction of equilibrium constants for some protein-protein associations

For many protein families, such as serine proteinases or serine proteinase inhibitors, the family assignment predicts reactivity only in general terms. Both detailed specificity and quantitative reactivity are lacking. We believe that, for many such protein families, algorithms can be devised by defining the subset of n functionally important sequence positions, making the 19n possible single mutants and measuring their reactivity. Given the assumption that the contributions of the n positions are additive, the reactivities of the 20(n) variants can be predicted. This is illustrated by an almost complete algorithm for the Kazal family of protein inhibitors of serine proteinases.

Peptide mimics of the Bowman-Birk inhibitor reactive site loop

Bowman-Birk Inhibitors (BBIs) are small highly cross-linked proteins that typically display an almost symmetrical "double-headed" structure. Each "head" contains an independent proteinase binding domain. The realization that one BBI molecule could form a 1:1:1 complex with two enzymes led early workers to dissect this activity. Now, after three decades of research, it has been possible to isolate the antiproteinase activity as small ( approximately 11 residues), cyclic, synthetic peptides, which display most of the functional aspects of the protein. More recently, it has been found that these peptide fragments are not just a synthetic curiosity-a natural 14-residue cyclic peptide (SFTI-1), which too encapsulates the BBI inhibitory motif, is found to occur in sunflowers. This article reviews the properties of BBI-based peptides (including SFTI-1) and discusses the features that are important for inhibitory activity.

Analysis of serine proteinase-inhibitor interaction by alanine shaving

We analyzed the energetic importance of residues surrounding the hot spot (the P(1) position) of bovine pancreatic trypsin inhibitor (BPTI) in interaction with two proteinases, trypsin and chymotrypsin, by a procedure called molecular shaving. One to eight residues of the structural epitope, composed of two extended and exposed loops, were mutated to alanine(s). Although truncation of the side chains of residues surrounding the P(1) position to methyl groups caused a decrease in Delta G(den) values up to 6.4 kcal mole(-1), it did not influence the overall conformation of the inhibitor. We found that the replacement of up to six residues with alanines was fully additive at the level of protein stability. To analyze the influence of the structural epitope on the association energy, we determined association constants for BPTI variants and both enzymes and applied the additivity analysis. Shaving of two binding loops led to a progressive drop in the association energy, more pronounced for trypsin (decrease up to 9.6 kcal mole(-1)) than chymotrypsin (decrease up to 3.5 kcal mole(-1)). In the case of extensively mutated variants interacting with chymotrypsin, the association energies agreed very well with the values calculated from single mutational effects. However, when P(1)-neighboring residues were shaved to alanine(s), their contribution to the association energy was not fully removed because of the presence of methyl groups and main chain-main chain intermolecular hydrogen bonds. Moreover, the hot spot had a different contribution to the complex stability in the fully shaved BPTI variant compared with the wild type, which was caused by perturbations of the P(1)-S(1) electrostatic interaction.

Binding interactions between peptides and proteins of the class II major histocompatibility complex

The activation of helper T cells by peptides bound to proteins of the class II Major Histocompatibility Complex (MHC II) is pivotal to the initiation of an immune response. The primary functional requirement imposed on MHC II proteins is the ability to efficiently bind thousands of different peptides. Structurally, this is reflected in a unique architecture of binding interactions. The peptide is bound in an extended conformation within a groove on the membrane distal surface of the protein that is lined with several pockets that can accommodate peptide side-chains. Conserved MHC II protein residues also form hydrogen bonds along the length of the peptide main-chain. Here we review recent advances in the study of peptide-MHC II protein reactions that have led to an enhanced understanding of binding energetics. These results demonstrate that peptide-MHC II protein complexes achieve high affinity binding from the array of hydrogen bonds that are energetically segregated from the pocket interactions, which can then add to an intrinsic hydrogen bond-mediated affinity. Thus, MHC II proteins are unlike antibodies, which utilize cooperativity among binding interactions to achieve high affinity and specificity. The significance of these observations is discussed within the context of possible mechanisms for the HLA-DM protein that regulates peptide presentation in vivo and the design of non-peptide molecules that can bind MHC II proteins and act as vaccines or immune modulators.

Chromogenic substrates of bovine beta-trypsin: the influence of an amino acid residue in P1 position on their interaction with the enzyme

The Cucurbita maxima trypsin inhibitor CMTI-III molecule was used as a vehicle to design and synthesize a series of trypsin chromogenic substrates modified in position P1: Ac-Ala-Val-Abu-Pro-X-pNA, where X = Orn, Lys, Arg, Har, Arg(NO(2)), Cit, Hci, Phe(p-CN), Phe(p-NH(2)); pNA = p-nitroanilide. The most active compounds (as determined by specificity constant k(cat)/K(m)) were peptides with the Arg and Lys residues in the position discussed. Changes in the length and the decrease of the positive charge of the amino acid residue side chain in position P(1) resulted in the decrease or loss of the affinity towards bovine beta-trypsin. Among peptides containing amino acid residues with uncharged side chains in position P1, only one with p-cyano-l-Phe revealed activity. These results correspond well with trypsin inhibitory activity of CMTI-III analogues modified in the equivalent position, indicating the same type of interaction between position P1 of the substrate or inhibitor and S1 site specificity of trypsin.

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.

Recent advances in prevention and treatment of hepatitis C virus infections

Hepatitis C virus (HCV) is the leading cause of chronic hepatitis in humans. As members of the flavivirus family, HCVs are a group of small single-stranded, positive-sense RNA viruses. Upon translation of the genome, a polyprotein precursor is synthesized and further processed by both cellular and viral proteases to generate functional viral proteins. Treatment options are currently limited to the administration of alpha-interferon alone or in combination with ribavirin. Unfortunately, these approaches are characterized by relatively poor efficacy and an unfavorable side-effect profile. Therefore, intensive effort is directed at the discovery of novel molecules to treat this disease. These new approaches include the development of prophylactic and therapeutic vaccines, the identification of interferons with improved pharmacokinetic characteristics, and the discovery of novel drugs designed to inhibit the function of three major viral proteins: protease, helicase and polymerase. Finally, the HCV RNA genome itself, particularly the IRES element, is being actively exploited as an antiviral target using antisense molecules and catalytic ribozymes. This review summarizes the most recent findings in each of these areas. Although not intended to be comprehensive, it should serve as a first resource for those individuals who desire updated information in this rapidly changing field.

Design, docking, and evaluation of multiple libraries against multiple targets

We present a general approach to the design, docking, and virtual screening of multiple combinatorial libraries against a family of proteins. The method consists of three main stages: docking the scaffold, selecting the best substituents at each site of diversity, and comparing the resultant molecules within and between the libraries. The core "divide-and-conquer" algorithm for side-chain selection, developed from an earlier version (Sun et al., J Comp Aided Mol Design 1998;12:597-604), provides a way to explore large lists of substituents with linear rather than combinatorial time dependence. We have applied our method to three combinatorial libraries and three serine proteases: trypsin, chymotrypsin, and elastase. We show that the scaffold docking procedure, in conjunction with a novel vector-based orientation filter, reproduces crystallographic binding modes. In addition, the free-energy-based scoring procedure (Zou et al., J Am Chem Soc 1999;121:8033-8043) is able to reproduce experimental binding data for P1 mutants of macromolecular protease inhibitors. Finally, we show that our method discriminates between a peptide library and virtual libraries built on benzodiazepine and tetrahydroisoquinolinone scaffolds. Implications of the docking results for library design are explored.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Predicting the reactivity of proteins from their sequence alone: Kazal family of protein inhibitors of serine proteinases

An additivity-based sequence to reactivity algorithm for the interaction of members of the Kazal family of protein inhibitors with six selected serine proteinases is described. Ten consensus variable contact positions in the inhibitor were identified, and the 19 possible variants at each of these positions were expressed. The free energies of interaction of these variants and the wild type were measured. For an additive system, this data set allows for the calculation of all possible sequences, subject to some restrictions. The algorithm was extensively tested. It is exceptionally fast so that all possible sequences can be predicted. The strongest, the most specific possible, and the least specific inhibitors were designed, and an evolutionary problem was solved.

Identification of chymotrypsin inhibitors from a second-generation template assisted combinatorial peptide library

In an earlier study (McBride JD, Freeman N, Domingo GJ, Leatherbarrow RJ. Selection of chymotrypsin inhibitors from a conformationally-constrained combinatorial peptide library. J. Mol. Biol. 1996; 259: 819-827) we described a resin-bound cyclic peptide library, constructed based on the sequence of the anti-tryptic reactive site loop of Bowman Birk Inhibitor (BBI), a proteinase inhibitor protein. This library was used to identify re-directed chymotrypsin inhibitors with Ki values as low as 17 nM. We have now extended this work by constructing an enhanced library in which a further position, at the P4 site of the inhibitor, has been randomized. This new library has variation at three target locations (P4, P1 and P2) within the inhibitory loop region, producing 8,000 variants. Screening this library allowed selection of new inhibitor sequences with Ki values as low as 3.4 nM. The success of this approach is reflected by the fact that the inhibition constant given by the selected peptide sequence is slightly lower than that reported against chymotrypsin for the most studied full length BBI protein, Soybean BBI 2-IV.

Substitutions at the P(1) position in BPTI strongly affect the association energy with serine proteinases

The role of the S(1) subsite in trypsin, chymotrypsin and plasmin has been examined by measuring the association with seven different mutants of bovine pancreatic trypsin inhibitor (BPTI); the mutants contain Gly, Ala, Ser, Val, Leu, Arg, and Trp at the P(1) position of the reactive site. The effects of substitutions at the P(1) position on the association constants are very large, comprising seven orders of magnitude for trypsin and plasmin, and over five orders for chymotrypsin. All mutants showed a decrease of the association constant to the three proteinases in the same order: Ala>Gly>Ser>Arg>Val>Leu>Trp. Calorimetric and circular dichroism methods showed that none of the P1 substitutions, except the P1-Val mutant, lead to destabilisation of the binding loop conformation. The X-ray structure of the complex formed between bovine beta-trypsin and P(1)-Leu BPTI showed that the P(1)-Leu sterically conflicts with the side-chain of P(3)-Ile, which thereby is forced to rotate approximately 90 degrees. Ile18 (P(3)) in its new orientation, in turn interacts with the Tyr39 side-chain of trypsin. Introduction of a large side-chain at the P1' position apparently leads to a cascade of small alterations of the trypsin-BPTI interface that seem to destabilise the complex by it adopting a less optimized packing and by tilting the BPTI molecule up to 15 degrees compared to the native trypsin-BPTI complex.

Directed evolution study of temperature adaptation in a psychrophilic enzyme

We have used laboratory evolution methods to enhance the thermostability and activity of the psychrophilic protease subtilisin S41, with the goal of investigating the mechanisms by which this enzyme can adapt to different selection pressures. A combined strategy of random mutagenesis, saturation mutagenesis and in vitro recombination (DNA shuffling) was used to generate mutant libraries, which were screened to identify enzymes that acquired greater thermostability without sacrificing low-temperature activity. The half-life of seven-amino acid substitution variant 3-2G7 at 60 degrees C is approximately 500 times that of wild-type and far surpasses those of homologous mesophilic subtilisins. The dependence of half-life on calcium concentration indicates that enhanced calcium binding is largely responsible for the increased stability. The temperature optimum of the activity of 3-2G7 is shifted upward by approximately 10 degrees C. Unlike natural thermophilic enzymes, however, the activity of 3-2G7 at low temperatures was not compromised. The catalytic efficiency, k(cat)/K(M), was enhanced approximately threefold over a wide temperature range (10 to 60 degrees C). The activation energy for catalysis, determined by the temperature dependence of k(cat)/K(M) in the range 15 to 35 degrees C, is nearly identical to wild-type and close to half that of its highly similar mesophilic homolog, subtilisin SSII, indicating that the evolved S41 enzyme retained its psychrophilic character in spite of its dramatically increased thermostability. These results demonstrate that it is possible to increase activity at low temperatures and stability at high temperatures simultaneously. The fact that enzymes displaying both properties are not found in nature most likely reflects the effects of evolution, rather than any intrinsic physical-chemical limitations on proteins.

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.

What can the structures of enzyme-inhibitor complexes tell us about the structures of enzyme substrate complexes?

Proteinases perform many beneficial functions that are essential to life, but they are also dangerous and must be controlled. Here we focus on one of the control mechanisms: the ubiquitous presence of protein proteinase inhibitors. We deal only with a subset of these: the standard mechanism, canonical protein inhibitors of serine proteinases. Each of the inhibitory domains of such inhibitors has one reactive site peptide bond, which serves all the cognate enzymes as a substrate. The reactive site peptide bond is in a combining loop which has an identical conformation in all inhibitors and in all enzyme-inhibitor complexes. There are at least 18 families of such inhibitors. They all share the conformation of the combining loops but each has its own global three-dimensional structure. Many three-dimensional structures of enzyme-inhibitor complexes were determined. They are frequently used to predict the conformation of substrates in very short-lived enzyme-substrate transition state complexes. Turkey ovomucoid third domain and eglin c have a Leu residue at P(1). In complexes with chymotrypsin, these P(1) Leu residues assume the same conformation. The relative free energies of binding of P(1) Leu (relative to either P(1) Gly or P(1) Ala) are within experimental error, the same for complexes of turkey ovomucoid third domain, eglin c, P(1) Leu variant of bovine pancreatic trypsin inhibitor and of a substrate with chymotrypsin. Therefore, the P(1) Leu conformation in transition state complexes is predictable. In contrast, the conformation of P(1) Lys(+) is strikingly different in the complexes of Lys(18) turkey ovomucoid third domain and of bovine pancreatic trypsin inhibitor with chymotrypsin. The relative free energies of binding are also quite different. Yet, the relative free energies of binding are nearly identical for Lys(+) in turkey ovomucoid third domain and in a substrate, thus allowing us to know the structure of the latter. Similar reasoning is applied to a few other systems.

Selection of human elastase inhibitors from a conformationally constrained combinatorial peptide library

A resin-bound cyclic peptide library was constructed based on the sequence of the reactive-site loop of Bowman-Birk inhibitor, a proteinase inhibitor protein. The constrained loop sequence, which incorporates the minimal proteinase-binding motif, was retained throughout the library, but selected residues known to be important for inhibitor specificity were randomised. The approach was used to create a 'one bead, one peptide' library with 8000 variants resulting from randomization at three target locations in the sequence (P4, P1 and P2'). This library allows us to examine the degree to which variations in this proteinase-binding motif can redirect activity, as well as providing information about the binding specificity of a proteinase target. Screening this library for binding to human leucocyte elastase identified sequences with a strong consensus, and on resynthesis all were found to act as inhibitors, with Ki values as low as 65 nM. Human leucocyte elastase is known to have a substrate preference for small alkyl chains at the P1 locus, with valine being preferred. However, alanine and not the expected valine was found in 21 out of 23 identified sequences. The remaining two sequences had threonine at P1, a finding that would be hard to predict based on substrate specificity alone. Further analysis of resynthesized peptides demonstrated that valine substitution results in an analogue that is hydrolysed far more rapidly than ones having library-selected P1 residues. Testing of the human leucocyte elastase-selected sequences as inhibitors of porcine pancreatic elastase demonstrates a significant difference in the specificity of the P4 locus between these two proteinases.

Probing intermolecular backbone H-bonding in serine proteinase-protein inhibitor complexes

BACKGROUND

Intermolecular backbone H-bonding (N-H.O=C) is a common occurrence at the interface of protein-protein complexes. For instance, the amide NH groups of most residues in the binding loop of eglin c, a potent serine proteinase inhibitor from the leech Hirudo medicinalis, are H-bonded to the carbonyl groups of residues in the target enzyme molecules such as chymotrypsin, elastase and subtilisins. We sought to understand the energetic significance of these highly conserved backbone-backbone H-bonds in the enzyme-inhibitor complexes.

RESULTS

We synthesized an array of backbone-engineered ester analogs of eglin c using native chemical ligation to yield five inhibitor proteins each containing a single backbone ester bond from P3 to P2' (i.e. -CONH-to -COO-). The structure at the ligation site (P6-P5) is essentially unaltered as shown by a high-resolution analysis of the subtilisin-BPN'-eglin c complex. The free-energy changes (DeltaDeltaGNH-->O) associated with the binding of ester analogs at P3, P1 and P2' with bovine alpha-chymotrypsin, subtilisin Carlsberg and porcine pancreatic elastase range from 0-4.5 kcal/mol. Most markedly, the NH-->O substitution at P2 not only stabilizes the inhibitor but also enhances binding to the enzymes by as much as 500-fold.

CONCLUSIONS

Backbone H-bond contributions are context dependent in the enzyme-eglin c complexes. The interplay of rigidity and adaptability of the binding loop of eglin c seems to play a prominent role in defining the binding action.

Interscaffolding additivity: binding of P1 variants of bovine pancreatic trypsin inhibitor to four serine proteases

Different families of protein inhibitors of serine proteases share similar conformation of the enzyme-binding loop, while their scaffoldings are completely different. In the enzyme-inhibitor complex, the P1position of the loop makes numerous contacts within the S1pocket and significantly influences the energy of the interaction. Here, we determine the association energies (DeltaGavalues) for the interaction of coded P1variants of bovine pancreatic trypsin inhibitor (BPTI) with bovine beta-trypsin (BT), anionic salmon trypsin (AST), bovine alpha-chymotrypsin (BCHYM), and human neutrophil elastase (HNE). The respective DeltaGaranges are 15, 13, 9, and 8 kcal mol-1(1 cal=4.18 J). Next, through interscaffolding additivity cycles, we compare our set of DeltaGavalues determined for BCHYM and HNE with similar data sets available in the literature for three other inhibitor families. The analysis of the cycles shows that 27 to 83 % of cycles fulfil the criteria of additvity. In one particular case (comparison of associations of P1variants of BPTI and OMTKY3 with BCHYM) there is a structural basis for strongly non-additive behaviour. We argue that the interscaffolding additvity depends on sequential and conformational similarities of sites where the mutation(s) are introduced and on the particular substitution. In the second interscaffolding analysis, we compare binding of the same P1mutants to BT and AST. The high correlation coefficient shows that both trypsins recognize with comparable strength the non-cognate side-chains. However, the cognate Arg and Lys side-chains are recognized significantly more strongly by AST.

The crystal structures of the complexes between bovine beta-trypsin and ten P1 variants of BPTI

The high-resolution X-ray structures have been determined for ten complexes formed between bovine beta-trypsin and P1 variants (Gly, Asp, Glu, Gln, Thr, Met, Lys, His, Phe, Trp) of bovine pancreatic trypsin inhibitor (BPTI). All the complexes were crystallised from the same conditions. The structures of the P1 variants Asp, Glu, Gln and Thr, are reported here for the first time in complex with any serine proteinase. The resolution of the structures ranged from 1.75 to 2.05 A and the R-factors were about 19-20 %. The association constants of the mutants ranged from 1.5x10(4) to 1.7x10(13) M-1. All the structures could be fitted into well-defined electron density, and all had very similar global conformations. All the P1 mutant side-chains could be accomodated at the primary binding site, but relative to the P1 Lys, there were small local changes within the P1-S1 interaction site. These comprised: (1) changes in the number and dynamics of water molecules inside the pocket; (2) multiple conformations and non-optimal dihedral angles for some of the P1 side-chains, Ser190 and Gln192; and (3) changes in temperature factors of the pocket walls as well as the introduced P1 side-chain. Binding of the cognate P1 Lys is characterised by almost optimal dihedral angles, hydrogen bonding distances and angles, in addition to considerably lower temperature factors. Thus, the trypsin S1 pocket seems to be designed particularly for lysine binding.

Comparison of protein-protein interactions in serine protease-inhibitor and antibody-antigen complexes: implications for the protein docking problem

The protein-protein interaction energy of 12 nonhomologous serine protease-inhibitor and 15 antibody-antigen complexes is calculated using a molecular mechanics formalism and dissected in terms of the main-chain vs. side-chain contribution, nonrotameric side-chain contributions, and amino acid residue type involvement in the interface interaction. There are major differences in the interactions of the two types of protein-protein complex. Protease-inhibitor complexes interact predominantly through a main-chain-main-chain mechanism while antibody-antigen complexes interact predominantly through a side-chain-side-chain or a side-chain-main-chain mechanism. However, there is no simple correlation between the main-chain-main-chain interaction energy and the percentage of main-chain surface area buried on binding. The interaction energy is equally effected by the presence of nonrotameric side-chain conformations, which constitute approximately 20% of the interaction energy. The ability to reproduce the interface interaction energy of the crystal structure if original side-chain conformations are removed from the calculation is much greater in the protease-inhibitor complexes than the antibody-antigen complexes. The success of a rotameric model for protein-protein docking appears dependent on the extent of the main-chain-main-chain contribution to binding. Analysis of (1) residue type and (2) residue pair interactions at the interface show that antibody-antigen interactions are very restricted with over 70% of the antibody energy attributable to just six residue types (Tyr > Asp > Asn > Ser > Glu > Trp) in agreement with previous studies on residue propensity. However, it is found here that 50% of the antigen energy is attributable to just four residue types (Arg = Lys > Asn > Asp). On average just 12 residue pair interactions (6%) contribute over 40% of the favorable interaction energy in the antibody-antigen complexes, with charge-charge and charge/polar-tyrosine interactions being prominent. In contrast protease inhibitors use a diverse set of residue types and residue pair interactions.

Phage display selection of P1 mutants of BPTI directed against five different serine proteinases

The P1 position of protein inhibitors and oligopeptide substrates determines, to a large extent, association energy with many serine proteinases. To test the agreement of phage display selection with the existing thermodynamic data, a small library of all 20 P1 mutants of basic pancreatic trypsin inhibitor (BPTI) was created, fused to protein III, and displayed on the surface of M13 phage. The wild type of displayed inhibitor monovalently and strongly inhibited trypsin with an association constant of Ka = 3 x 10(11) M(-1). The library was applied to select BPTI variants active against five serine proteinases of different specificity (bovine trypsin and chymotrypsin, human leukocyte and porcine pancreatic elastases, human azurocidin). The results of enrichment with four proteinases agreed well with the available thermodynamic data. In the case of azurocidin, the phage display selection allowed determination of the P1 specificity of this protein with the following frequencies for selected P1 variants: 43% Lys, 36% Leu, 7% Met, 7% Thr, 7% Gln.

Interaction of second and third domains of Japanese quail ovomucoid with ten mammalian trypsins

Second and third domains were prepared from Japanese quail ovomucoid and association equilibrium constants, Kas, were measured at 25 degreesC and pH 8 for these domains with trypsins from ten mammalian species: cat, cow, dog, guinea pig, hog, horse, man, rabbit, rat, and sheep. The values ranged from 108 M-1 to 1010 M-1 for the second domain-trypsin associations and from 106 M-1 to 108 M-1 for the third domain-trypsin associations. Changes in Ka values for the interactions between the trypsins and each domain are attributed to slight changes in surface conformation caused by the residue changes in the inhibitor-binding region other than the S1 pocket of the trypsin species. The representative of such residue changes is assumed to be the one observed at residue 217 of trypsin molecule. Concerning each trypsin, the Ka value with the second domain was always higher than that with the third domain. However, the ratios between the two equilibrium constants varied from 3 to 60 depending upon trypsin species. This means that amino acid changes in enzyme-contact residues other than the P1 site of the Kazal-type inhibitor can make it possible to recognize even a slight difference in inhibitor-binding surface among the enzymes with the same S1 pocket and highly similar overall three-dimensional structure.

Specificity of human cathepsin G

A series of tetrapeptide p-nitroanilide substrates of the general formula: suc-Ala-Ala-Pro-Aaa-p-nitroanilide was used to map the S1 binding pocket of human cathepsin G. Based on the kcat/Km parameter, the following order of preference was found: Lys=Phe>Arg=Leu>Met>Nle=Nva>Ala>Asp. Thus, the enzyme exhibits clear dual and equal trypsin- and chymotrypsin-like specificities. Particularly deleterious were beta-branched side chains of Ile and Val. The P1 substrate preferences found for cathepsin G are distinctly different from many other serine proteinases, including fiddler crab collagenase and chymotrypsin. The kcat/Km values obtained for P1 Lys, Phe, Arg and Leu substrates correlate well with those determined for analogous P1 mutants of basic pancreatic trypsin inhibitor (BPTI) obtained through recombinant techniques. To characterise the subsite specificity of the enzyme, a series of Cucurbita maxima trypsin inhibitor I (CMTI I) mutants were used comprising P2-P3' and P12' positions. All the mutants obtained were inhibitors of cathepsin G with association constants in the range: 105-109 M-1. Some of the mutations destabilised complex formation. In particular, Met8-->Arg substitution at P3', which increased association constant for chymotrypsin 46-fold, led to a 7-fold decrease of binding with cathepsin G. In addition, mutation of Ile6 at position P1' either to Val or Asp was deleterious for cathepsin G. In two cases (Ala18-->Gly (P12') and Pro4-->Thr (P2)), about a 10-fold increase in association constants was observed.