The effects of exosite occupancy on the substrate specificity of thrombin

Effect of Sulfotyrosine and Negatively Charged Amino Acid of Leech-Derived Peptides on Binding and Inhibitory Activity Against Thrombin

Hirudins, natural sulfo(glyco)proteins, are clinical anticoagulants that directly inhibit thrombin, a key coagulation factor. Their potent thrombin inhibition primarily results from antagonistic interactions with both the catalytic and non-catalytic sites of thrombin. Hirudins often feature sulfate moieties on Tyr residues in their anionic C-terminus region, enabling strong interactions with thrombin exosite-I and effectively blocking its engagement with fibrinogen. Although sulfotyrosines have been identified in various hirudin variants, the precise relationship between sulfotyrosine and the number of negatively charged amino acids within the anionic-rich C-terminus peptide domain for the binding of thrombin has remained elusive. By using Fmoc-SPPS, hirudin dodecapeptides homologous to the C-terminus of hirudin variants from various leech species were successfully synthesized, and the effect of sulfotyrosine and the number of negatively charged amino acids on hirudin-thrombin interactions was investigated. Our findings did not reveal any synergistic effect between an increasing number of sulfotyrosines or negatively charged amino acids and their inhibitory activity on thrombin or fibrinolysis in the assays, despite a higher binding level toward thrombin in the sulfated dodecapeptide Hnip_Hirudin was observed in SPR analysis.

Thrombin - A Molecular Dynamics Perspective

Thrombin is a crucial enzyme involved in blood coagulation, essential for maintaining circulatory system integrity and preventing excessive bleeding. However, thrombin is also implicated in pathological conditions such as thrombosis and cancer. Despite the application of various experimental techniques, including X-ray crystallography, NMR spectroscopy, and HDXMS, none of these methods can precisely detect thrombin's dynamics and conformational ensembles at high spatial and temporal resolution. Fortunately, molecular dynamics (MD) simulation, a computational technique that allows the investigation of molecular functions and dynamics in atomic detail, can be used to explore thrombin behavior. This review summarizes recent MD simulation studies on thrombin and its interactions with other biomolecules. Specifically, the 17 studies discussed here provide insights into thrombin's switch between 'slow' and 'fast' forms, active and inactive forms, the role of Na+ binding, the effects of light chain mutation, and thrombin's interactions with other biomolecules. The findings of these studies have significant implications for developing new therapies for thrombosis and cancer. By understanding thrombin's complex behavior, researchers can design more effective drugs and treatments that target thrombin.

Structural and functional analysis of the simultaneous binding of two duplex/quadruplex aptamers to human α-thrombin

The long-range communication between the two exosites of human α-thrombin (thrombin) tightly modulates the protein-effector interactions. Duplex/quadruplex aptamers represent an emerging class of very effective binders of thrombin. Among them, NU172 and HD22 aptamers are at the forefront of exosite I and II recognition, respectively. The present study investigates the simultaneous binding of these two aptamers by combining a structural and dynamics approach. The crystal structure of the ternary complex formed by the thrombin with NU172 and HD22_27mer provides a detailed view of the simultaneous binding of these aptamers to the protein, inspiring the design of novel bivalent thrombin inhibitors. The crystal structure represents the starting model for molecular dynamics studies, which point out the cooperation between the binding at the two exosites. In particular, the binding of an aptamer to its exosite reduces the intrinsic flexibility of the other exosite, that preferentially assumes conformations similar to those observed in the bound state, suggesting a predisposition to interact with the other aptamer. This behaviour is reflected in a significant increase of the anticoagulant activity of NU172 when the inactive HD22_27mer is bound to exosite II, providing a clear evidence of the synergic action of the two aptamers.

Molecular dynamics simulations of human α-thrombin in different structural contexts: evidence for an aptamer-guided cooperation between the two exosites

Human α-thrombin (thrombin) is a multifunctional enzyme that plays a pivotal role in the coagulation pathway. Thrombin activity can be effectively modulated by G-quadruplex-based oligonucleotide aptamers that specifically interact with the two positively charged regions (exosites I and II) on the protein surface. Although insightful atomic-level snapshots of the recognition between thrombin and aptamers have been recently achieved through crystallographic analyses, some dynamic aspects of this interaction have not been fully characterized. We here report molecular dynamics simulations of thrombin in different association states: ligand-free and binary/ternary complexes with the aptamers TBA (at exosite I) and HD22_27mer (at exosite II). The simulations carried out on the binary and ternary complexes formed by thrombin with these aptamers provide a dynamic view of the interactions that stabilize them in a crystal-free environment. Interestingly, the analysis of the dynamics of the exosites in different thrombin binding states clearly indicates that the HD22_27mer binding at the exosite II favours conformations of exosite I that are prone to the TBA binding. Similar effects are observed upon the binding of TBA to the exosite I. These observations provide an atomic-level picture of the exosite inter-communication in thrombin and explain the experimentally detected cooperativity of the TBA/HD22_27mer binding.

High throughput protease profiling comprehensively defines active site specificity for thrombin and ADAMTS13

We have combined random 6 amino acid substrate phage display with high throughput sequencing to comprehensively define the active site specificity of the serine protease thrombin and the metalloprotease ADAMTS13. The substrate motif for thrombin was determined by >6,700 cleaved peptides, and was highly concordant with previous studies. In contrast, ADAMTS13 cleaved only 96 peptides (out of >10⁷ sequences), with no apparent consensus motif. However, when the hexapeptide library was substituted into the P3-P3′ interval of VWF73, an exosite-engaging substrate of ADAMTS13, 1670 unique peptides were cleaved. ADAMTS13 exhibited a general preference for aliphatic amino acids throughout the P3-P3′ interval, except at P2 where Arg was tolerated. The cleaved peptides assembled into a motif dominated by P3 Leu, and bulky aliphatic residues at P1 and P1′. Overall, the P3-P2′ amino acid sequence of von Willebrand Factor appears optimally evolved for ADAMTS13 recognition. These data confirm the critical role of exosite engagement for substrates to gain access to the active site of ADAMTS13, and define the substrate recognition motif for ADAMTS13. Combining substrate phage display with high throughput sequencing is a powerful approach for comprehensively defining the active site specificity of proteases.

Dabigatran and Argatroban Diametrically Modulate Thrombin Exosite Function

Thrombin is a highly plastic molecule whose activity and specificity are regulated by exosites 1 and 2, positively-charged domains that flank the active site. Exosite binding by substrates and cofactors regulates thrombin activity by localizing thrombin, guiding substrates, and by inducing allosteric changes at the active site. Although inter-exosite and exosite-to-active-site allostery have been demonstrated, the impact of active site ligation on exosite function has not been examined. To address this gap, we used surface plasmon resonance to determine the effects of dabigatran and argatroban, active site-directed inhibitors, on thrombin binding to immobilized γ_A/γ_A-fibrin or glycoprotein Ibα peptide via exosite 1 and 2, respectively, and thrombin binding to γ_A/γ′-fibrin or factor Va, which is mediated by both exosites. Whereas dabigatran attenuated binding, argatroban increased thrombin binding to γ_A/γ_A- and γ_A/γ′-fibrin and to factor Va. The results with immobilized fibrin were confirmed by examining the binding of radiolabeled thrombin to fibrin clots. Thus, dabigatran modestly accelerated the dissociation of thrombin from γ_A/γ_A-fibrin clots, whereas argatroban attenuated dissociation. Dabigatran had no effect on thrombin binding to glycoprotein Ibα peptide, whereas argatroban promoted binding. These findings not only highlight functional effects of thrombin allostery, but also suggest that individual active site-directed thrombin inhibitors uniquely modulate exosite function, thereby identifying potential novel mechanisms of action.

How does association process affect fibrinogen hydrolysis by thrombin?

Thrombin, a key enzyme in the blood coagulation cascade, hydrolyzes fibrinogen into fibrin, which specifically associates into the fibers that build up a thrombus scaffold. The assembly of fibrin involves a set of stepwise reactions, for which a complete and detailed kinetic portrait is needed. Existing kinetic models focus on particular parts of the process, for example the mechanism of enzyme action itself or the kinetics of formation of fibrin assemblies. The current study considers a thorough model of the process from fibrinogen hydrolysis to the assembly of fibrin. Composing the model requires taking into account several reaction intermediates, stepwise removal of fibrinopeptides, and association of partially hydrolyzed fibrin, in particular desAA fibrin. The model is versatile enough to adopt new data both on fibrinogen hydrolysis and fibrin association. In addition, the model could be considered as an example of a kinetic description of other complex enzyme systems having several intermediates and feedbacks, such as the blood coagulation cascade and signal transduction.

Ligand binding to anion-binding exosites regulates conformational properties of thrombin

Thrombin participates in coagulation, anticoagulation, and initiation of platelet activation. To fulfill its diverse roles and maintain hemostasis, this serine protease is regulated via the extended active site region and anion-binding exosites (ABEs) I and II. For the current project, amide proton hydrogen-deuterium exchange coupled with MALDI-TOF mass spectrometry was used to characterize ligand binding to individual exosites and to investigate the presence of exosite-active site and exosite-exosite interactions. PAR3(44-56) and PAR1(49-62) were observed to bind to thrombin ABE I and then to exhibit long range effects over to ABE II. By contrast, Hirudin(54-65) focused more on ABE I and did not transmit influences over to ABE II. Although these three ligands were each directed to ABE I, they did not promote the same conformational consequences. D-Phe-Pro-Arg-chloromethyl ketone inhibition at the thrombin active site led to further local and long range consequences to thrombin-ABE I ligand complexes with the autolysis loop often most affected. When Hirudin(54-65) was bound to ABE I, it was still possible to bind GpIbα(269-286) or fibrinogen γ'(410-427) to ABE II. Each ligand exerted its predominant influences on thrombin and also allowed interexosite communication. The results obtained support the proposal that thrombin is a highly dynamic protein. The transmission of ligand-specific local and long range conformational events is proposed to help regulate this multifunctional enzyme.

Designing allosteric regulators of thrombin. Monosulfated benzofuran dimers selectively interact with Arg173 of exosite 2 to induce inhibition

Earlier, we reported on the design of sulfated benzofuran dimers (SBDs) as allosteric inhibitors of thrombin (Sidhu et al. J. Med. Chem.201154 5522-5531). To identify the site of binding of SBDs, we studied thrombin inhibition in the presence of exosite 1 and 2 ligands. Whereas hirudin peptide and heparin octasaccharide did not affect the IC(50) of thrombin inhibition by a high affinity SBD, the presence of full-length heparin reduced inhibition potency by 4-fold. The presence of γ' fibrinogen peptide, which recognizes Arg93, Arg97, Arg173, Arg175, and other residues, resulted in a loss of affinity that correlated with the ideal Dixon-Webb competitive profile. Replacement of several arginines and lysines of exosite 2 with alanine did not affect thrombin inhibition potency, except for Arg173, which displayed a 22-fold reduction in IC(50). Docking studies suggested a hydrophobic patch around Arg173 as a plausible site of SBD binding to thrombin. The absence of the Arg173-like residue in factor Xa supported the observed selectivity of inhibition by SBDs. Cellular toxicity studies indicated that SBDs are essentially nontoxic to cells at concentrations as high as 250 mg/kg. Overall, the work presents the localization of the SBD binding site, which could lead to allosteric modulators of thrombin that are completely different from all clinically used anticoagulants.

Role of electrostatic interactions in binding of thrombin to the fibrinogen γ' chain

Thrombin binds to the highly anionic fibrinogen γ' chain through anion-binding exosite II. This binding profoundly alters thrombin's ability to cleave substrates, including fibrinogen, factor VIII, and PAR1. However, it is unknown whether this interaction is due mainly to general electrostatic complementarity between the γ' chain and exosite II or if there are critical charged γ' chain residues involved. We therefore systematically determined the contribution of negatively charged amino acids in the γ' chain, both individually and collectively, to thrombin binding affinity. Surface plasmon resonance binding experiments were performed using immobilized γ' chain peptides with charged-to-uncharged amino acid substitutions, i.e., Asp to Asn, Glu to Gln, and pTyr to Tyr. Individually, the substitution of uncharged for charged amino acids resulted in only minor changes in binding affinity, with a maximal change in K(d) from 0.440 to 0.705 μM for the Asp419Asn substitution. However, substitution of all three charged amino acids in a conserved β-turn that is predicted to contact thrombin, pTyr418Tyr, Asp419Asn, and pTyr422Tyr, resulted in the loss of measurable binding, as did substitution of all the flanking charged amino acids. In addition, the binding of the γ' chain to thrombin was weakened in a dose-dependent manner with increasing NaCl concentration, resulting in a net loss of three or four ion pairs between thrombin and the γ' chain. Therefore, although each of the individual charges in the γ' chain contributes only incrementally to the overall binding affinity, the ensemble of the combined charges plays a profound role in the thrombin-γ' chain interactions.

The extended cleavage specificity of human thrombin

Thrombin is one of the most extensively studied of all proteases. Its central role in the coagulation cascade as well as several other areas has been thoroughly documented. Despite this, its consensus cleavage site has never been determined in detail. Here we have determined its extended substrate recognition profile using phage-display technology. The consensus recognition sequence was identified as, P2-Pro, P1-Arg, P1'-Ser/Ala/Gly/Thr, P2'-not acidic and P3'-Arg. Our analysis also identifies an important role for a P3'-arginine in thrombin substrates lacking a P2-proline. In order to study kinetics of this cooperative or additive effect we developed a system for insertion of various pre-selected cleavable sequences in a linker region between two thioredoxin molecules. Using this system we show that mutations of P2-Pro and P3'-Arg lead to an approximate 20-fold and 14-fold reduction, respectively in the rate of cleavage. Mutating both Pro and Arg results in a drop in cleavage of 200-400 times, which highlights the importance of these two positions for maximal substrate cleavage. Interestingly, no natural substrates display the obtained consensus sequence but represent sequences that show only 1-30% of the optimal cleavage rate for thrombin. This clearly indicates that maximal cleavage, excluding the help of exosite interactions, is not always desired, which may instead cause problems with dysregulated coagulation. It is likely exosite cooperativity has a central role in determining the specificity and rate of cleavage of many of these in vivo substrates. Major effects on cleavage efficiency were also observed for residues as far away as 4 amino acids from the cleavage site. Insertion of an aspartic acid in position P4 resulted in a drop in cleavage by a factor of almost 20 times.

Bioinformatic approaches for predicting substrates of proteases

Proteases have central roles in "life and death" processes due to their important ability to catalytically hydrolyze protein substrates, usually altering the function and/or activity of the target in the process. Knowledge of the substrate specificity of a protease should, in theory, dramatically improve the ability to predict target protein substrates. However, experimental identification and characterization of protease substrates is often difficult and time-consuming. Thus solving the "substrate identification" problem is fundamental to both understanding protease biology and the development of therapeutics that target specific protease-regulated pathways. In this context, bioinformatic prediction of protease substrates may provide useful and experimentally testable information about novel potential cleavage sites in candidate substrates. In this article, we provide an overview of recent advances in developing bioinformatic approaches for predicting protease substrate cleavage sites and identifying novel putative substrates. We discuss the advantages and drawbacks of the current methods and detail how more accurate models can be built by deriving multiple sequence and structural features of substrates. We also provide some suggestions about how future studies might further improve the accuracy of protease substrate specificity prediction.

Predicting serpin/protease interactions

Proteases are tightly regulated by specific inhibitors, such as serpins, which are able to undergo considerable and irreversible conformational changes in order to trap their targets. There has been a considerable effort to investigate serpin structure and functions in the past few decades; however, the specific interactions between proteases and serpins remain elusive. In this chapter, we describe detailed experimental protocols to determine and characterize the extended substrate specificity of proteases based on a substrate phage display technique. We also describe how to employ a bioinformatics system to analyze the substrate specificity data obtained from this technique and predict the potential inhibitory serpin partners of a protease (in this case, the immune protease, granzyme B) in a step-by-step manner. The method described here could also be applied to other proteases for more generalized substrate specificity analysis and substrate discovery.