Abstract
Thrombin is a serine protease that catalyzes a large number of different reactions including proteolytic cleave of fibrinogen to make the fibrin clot (procoagulant activity), of the protease activated receptors (for cell signaling) and of protein C generating activated protein C (anticoagulant activity). Thrombin has an effector binding site called the anion binding exosite 1 that is allosterically coupled to the active site. In this review, we survey results from thermodynamic characterization of the allosteric coupling as well as hydrogen-deuterium exchange mass spectrometry to reveal which parts of the thrombin structure are changed upon effector binding and/or mutagenesis, and finally NMR spectroscopy to characterize the different timescales of motions elicited by the effectors. We also relate the experimental work to computational network analysis of the thrombin-thrombomodulin complex.
1 Introduction
Thrombin is one of the best-studied serine proteases. The mechanism of serine proteases is multistep and involves covalent intermediates (Figure 1A). Depending on the substrate, reactions are slow compared to simpler enzymes (kcat for cleavage of fibrinogen by thrombin is 84 s-1 () whereas the kcat for adenylate kinase is 190 s-1 () and kcat for triosephosphate isomerase is 4,300 s-1 (). Allostery usually implies alteration in binding affinity by binding of an effector molecule at a second site. The case of thrombin is very interesting because it has a large number of different substrates and effectors change the activity of the enzyme towards one or another. Na + binding, for example, favors the cleavage of fibrinogen to form the fibrin clot (procoagulant activity) whereas thrombomodulin (TM) binding favors the proteolytic activation of protein C to form activated protein C (anticoagulant activity). Allostery also usually involves a structural change that can be detected by crystallography, however the structure of thrombin does not seem in many cases to be reflective of its activity.
FIGURE 1
Observations that thrombin may be an allosteric protein began to appear in the 1980s, even before the first crystal structures appeared. In 1981, Esmon discovered a cofactor, later named thrombomodulin (TM), that “switched” the activity of thrombin from pro-coagulant to anti-coagulant (). This “switch” behavior fascinated biophysicists who began to examine the conformational changes induced by TM once a smaller functional fragment of TM was isolated (). Using biophysical probes such as spin labels () and fluorophoric (), biophysicists demonstrated that TM-binding caused detectable alterations in the active site of thrombin. As with other allosteric systems, site-directed mutagenesis of key residues could partially shift the equilibrium of states to favor one or the other as demonstrated by mutation of a key TM-binding residue Glu192 to Gln (). Later structural studies revealed conformational changes caused by this mutation ().
Binding of adenosine nucleotides also modulated the activity of thrombin, and a full allosteric linkage analysis demonstrated these molecules were allosteric effectors (). In 1992, Di Cera’s group discovered that thrombin activity is modulated by Na + binding (), a phenomenon that had previously been reported for other blood coagulation serine proteases (). Here again, once structural studies emerged (), it was possible to find residues such as Trp 60D, far from the Na + binding site, which shifted the equilibrium of states (). Furthermore, despite the very small interaction surface expected for Na + binding, a large heat capacity change was observed upon Na + binding hinting that Na + binding may be a dynamic allosteric event. Most interestingly, the Na + binding and TM binding appeared to be allosterically coupled (). That meant that thrombin has at least 3 binding sites that might be coupled–the active site, the TM-binding site and the Na + -binding site.
2 Thermodynamic signatures of thrombin allostery
The crystal structure of the thrombin-TM complex revealed that TM only contacts thrombin via its fifth EGF-like domain despite the fact that the fourth EGF-like domain is required to “switch” the catalytic activity of thrombin () (Figure 1B). Our group attempted to measure the thrombin-TM binding by surface plasmon resonance and isothermal titration calorimetry with the goal of determining the coupling between Na + -binding and TM binding. By SPR, we could measure a Kd of 4 nM for the smallest active fragment of TM comprised of the 4th, 5th, and 6th EGF-like domains () (Figure 2A). The interaction was strongly electrostatically-driven as might be expected from the highly positively-charged TM binding site on thrombin and the fact that TM is highly negatively charged. Having demonstrated that we could detect binding, we attempted to measure the thermodynamics of binding by isothermal titration calorimetry (ITC). We could not measure any binding of TM to thrombin by ITC () (Figure 2C). As a control, we purchased a monoclonal antibody that was reported to bind in the TM binding site. By SPR, the antibody bound with the same binding affinity as TM, but the association and dissociation rates were much slower than the rates observed for TM binding (Figure 2B). The antibody binding was driven by a favorable ΔH as detected by ITC (Figure 2D) (). We surmised that the two interactions must have different thermodynamic driving forces with the thrombin-TM binding being entropically-driven. This was the first hint that thrombin may undergo dynamic allostery.
FIGURE 2
It is important to point out that the TM binding site on thrombin (Figure 1B) is called anion binding exosite 1 (ABE1). Fibrinogen interacts with thrombin at both the active site AND ABE1 (
Several groups have demonstrated the thermodynamic linkage between ABE1 and the thrombin active site using ITC. Binding of the lectin, bothrojaracin at ABE1 was shown to cause changes in the thrombin active site (
FIGURE 3

Isothermal titration calorimetry traces of various ligands binding to thrombin. (A) DAPA (70 µM in the syringe) binding to thrombin (6 µM in the cell). (B) DAPA (70 µM in the syringe) binding to the thrombin-aptamer complex (6 µM in the cell). (C) DAPA (70 µM in the syringe) binding to the thrombin-TM complex (6 µM in the cell). These data are from
3 Structure/function analyses of thrombin allostery
We revisited the amide exchange in thrombin later when mass spectrometry instrumentation had improved resulting in near 100% sequence coverage. Using the more up-to-date HDX-MS platform, we discovered an important aspect of TM-mediated allostery that was missed due to inadequate sequence coverage in the original thrombin HDX experiments. TM markedly stabilizes the active site conformation in which the new N-terminus is docked in the isoleucine binding pocket. All serine proteases follow the same activation mechanism in which they exist in a zymogen form until proteolytic cleavage at the activation peptide. As is the case for thrombin, many serine proteases contain cysteine disulfide bonds and the cleaved polypeptide remains attached to the protease via a disulfide bond. In thrombin, cleavage is after Arg, and the new N-terminal residue is an Ile. In all crystal structures of thrombin, one can observe the new N-terminal amino group interacting with Asp 194, which is adjacent to the catalytic Ser 195 (Figure 4A). We were surprised that in thrombin that did not have an inhibitor bound, the new N-terminal region was highly exchanging (Figure 4B). TM-binding markedly stabilizes this region, even though its binding site is distant (cf Figure 1B). In addition, TM binding stabilized the regions adjacent to and including the Na + -binding site (Figures 4C, D).
FIGURE 4

HDX-MS analysis of the thrombin-TM complex as compared to thrombin alone. (A) structure of thrombin showing the new N-terminus (blue) with Ile16 and the catalytic triad in sticks. (B) Deuterium incorporation into the new N-terminus of thrombin alone (blue) vs. in the thrombin-TM complex (black). (C) Deuterium incorporation into the 180 s loop of thrombin alone (blue) vs. in the thrombin-TM complex (black). (D) Deuterium incorporation into the Na + binding loop of thrombin alone (blue) vs. in the thrombin-TM complex (black).
Besides analyzing the thrombin-TM interaction, we also analyzed a number of thrombin mutants that had also been studied by other groups and for which crystal structures were available. Three mutants of particular interest will be discussed. The first is Glu 192 which when mutated to Gln was reported to mimic the catalytic switch induced by TM (
We used HDX-MS and accelerated molecular dynamics (AMD) to study the backbone dynamics of the W215A mutant thrombin and compared it to wild type thrombin (
FIGURE 5

Structure of WT thrombin (PDB 1PPB) highlighting 156–181 (residues 197–222; green). The side chains of Trp215 (blue) and the catalytic triad (black) are shown as sticks. Deuterium incorporation over 5 min into residues 156–180, 161–170, and 161–181 is shown for WT thrombin at 100 mM NaCl (gray) and 300 mM NaCl (black), and for the W215A mutant at 100 mM NaCl (cyan) and 300 mM NaCl (blue).
4 Computational network analysis of thrombin allostery
Accelerated molecular dynamics has the challenge that it is not trivial to determine how much to accelerate (i.e., how much to flatten the energy landscape) to explore realistic motions of proteins on longer timescales. We demonstrated that the AMD simulations could be calibrated by comparing experimental residual dipolar couplings with those back-calculated from the Boltzmann-reweighted AMD ensemble (
5 The multiple timescales of dynamic allostery in thrombin
The evidence from ITC and HDX-MS showed clearly that binding of TM to thrombin or mutation of specific amino acids could markedly alter the dynamics of thrombin far from the site of perturbation. However, the timescales of the motions and how they changed remained unknown. To answer this part of the question, we turned to NMR relaxation dispersion experiments. First, it is important to mention that the HDX-MS and NMR are highly complementary. Backbone resonances for many of the surface loops are simply missing from the NMR spectra indicating motions on the ns-µs timescale (
FIGURE 6

Relaxation dispersion results for (A) thrombin and (B) PPACK-inhibited thrombin. The spheres in each structure are those amide groups showing significant Rex. (C) Example dispersion curves and data points for apo-thrombin at 800 MHz (red) and 600 MHz (brown); and for PPACK-thrombin at 800 MHz (blue) and 600 MHz (black).
We also were able to use relaxation dispersion experiments to measure the motions in the thrombin-TM complex (
FIGURE 7

(A) Examples of thrombin-TM456 HSQC (red) doublet resonances compared to apo-thrombin (blue) resonances for residues at the N-terminus of the heavy chain. (B) Location of residues 16–23 (37-44seq) in the structure of thrombin (wheat) bound to TM456 (light blue) [PDB ID: 1DX5] (hot pink). The resonance corresponding to Glu 23CT (44seq) experienced similar μs-ms dynamics in thrombin-TM456 and apo-thrombin (purple sphere). (C) CPMG plot for the resonance corresponding to residue Glu 23 (44seq). The red and orange symbols are from spectra collected on thrombin-TM at 800 MHz and 600 MHz respectively, the blue and cyan curves are from apo-thrombin at 800 MHz and 600 MHz respectively, and the black and grey curves are from PPACK-thrombin at 800 MHz and 600 MHz respectively. These data are taken from
FIGURE 8

HDX-MS and NMR relaxation dispersion data showing the reduced dynamics in the C-terminal β-barrel of thrombin upon TM binding. (A) Structure of thrombin (wheat) bound to TM456 (light blue) [PDB ID: 1DX5]. The residues within the peptides identified in panel (A) are colored accordingly. Residues with missing HSQC resonances (black) and residues corresponding to resonances experiencing reduced μs-ms dynamics in thrombin-TM456 compared to apo-thrombin (blue spheres) are shown. (B) Deuterium uptake plots for the peptides spanning residues 161–180 (202-221seq), and 212–227 (260-275seq). (C) CPMG plots for resonances experiencing reduced μs-ms dynamics around the S1 pocket when TM is bound compared to apo-thrombin (symbols as in Figure 7). These data are taken from
6 A synthesis of the different observations of allostery
The picture that emerges from incorporating lessons from all these different experiments is that ABE1 and the active site are dynamically allosterically coupled. The dynamic allostery involves enthalpy-entropy compensation and from NMR relaxation dispersion experiments we see that whereas TM-binding causes increased dynamics in the N-terminal β-barrel, it causes decreased dynamics throughout the C-terminal β-barrel. In addition, TM causes slow motions that cause doublet signals near the active site. We did not detect the increased dynamics in the N-terminal β-barrel by HDX-MS, probably because the affected regions are β-strands interior to the protein. We did detect the decreased dynamics in the C-terminal β-barrel because the dynamic changes in this β-barrel largely involved loops in which the amide hydrogen bonds would be sensitive to slightly altered structures. The decreased dynamics in the C-terminal β-barrel are consistent with the structuring of the active site for substrate-binding and catalysis. Whereas Na + binding would be required to do this structuring in the absence of TM, Na + binding is not required with TM because TM itself is capable of structuring the active site. It is likely that each of these changes in dynamics plays a role in the overall catalytic mechanism of thrombin. TM-mediated structuring of the substrate binding site would increase the rate of association of protein C with thrombin as has been observed previously (
FIGURE 9

Schematic showing how the dynamics of thrombin are modulated by TM to facilitate cleavage of protein C. (A) TM reduces the dynamics and structures the C-terminal β-barrel facilitating substrate association. (B) TM drives the motions of the catalytic triad via µs-ms motions. (C) TM causes slow motions of the two β-barrels with respect to each other in order to facilitate release of the first product and subsequent entry of H2O into the active site.
Statements
Author contributions
EK is the sole author who wrote this review. The data in the figures is from previous publications and is re-presented here. All authors contributed to the article and approved the submitted version.
Funding
This work was funded by NIH R01HL127041.
Conflict of interest
The author declares that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.
Publisher’s note
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.
References
1
Baerga-OrtizA.BergqvistS.MandellJ.KomivesE. (2004). Two different proteins that compete for binding to thrombin have opposite kinetic and thermodynamic profiles. Protein Sci.13 (1), 166–176. 10.1110/ps.03120604
2
Baerga-OrtizA.RezaieA.KomivesE. (2000). Electrostatic dependence of the thrombin-thrombomodulin interaction. J. Mol. Biol.296 (2), 651–658. 10.1006/jmbi.1999.3447
3
BahA.CarrellC.ChenZ.GandhiP.Di CeraE. (2009). Stabilization of the E* form turns thrombin into an anticoagulant. J. Biol. Chem.284, 20034–20040. 10.1074/jbc.M109.012344
4
CantwellA.Di CeraE. (2000). Rational design of a potent anticoagulant thrombin. J. Biol. Chem.275, 39827–39830. 10.1074/jbc.C000751200
5
DangO.VindigniA.Di CeraE. (1995). An allosteric switch controls the procoagulant and anticoagulant activities of thrombin. Proc. Natl. Acad. Sci. U. S. A.92 (13), 5977–5981. 10.1073/pnas.92.13.5977
6
De CristofaroR.LandolfiR.Di CeraE. (1990). The linkage between adenosine nucleotide binding and amidase activity in human alpha-thrombin. Biophys. Chem.36, 77–84. 10.1016/0301-4622(90)85010-4
7
Di CeraE.DangQ.AyalaY. (1997). Molecular mechanisms of thrombin function. Cell Mol. Life Sci.53 (9), 701–730. 10.1007/s000180050091
8
Di CeraE.GuintoE.VindigniA.DangQ.AyalaY.WuyiM.et al (1995). The Na+ binding site of thrombin. J. Biol. Chem.270, 22089–22092. 10.1074/jbc.270.38.22089
9
EsmonC.OwenW. (1981). Identification of an endothelial cell cofactor for thrombin-catalyzed activation of protein C. Proc. Natl. Acad. Sci. U. S. A.78 (4), 2249–2252. 10.1073/pnas.78.4.2249
10
Fuentes-PriorP.IwanagaY.HuberR.PagilaR.RumennikG.SetoM.et al (2000). Structural basis for the anticoagulant activity of the thrombin-thrombomodulin complex. Nature404 (6777), 518–525. 10.1038/35006683
11
FuglestadB.GasperP.TonelliM.McCammonJ.MarkwickP.KomivesE. (2012). The dynamic structure of thrombin in solution. Biophys. J.103 (1), 79–88. 10.1016/j.bpj.2012.05.047
12
GandhiP.ChenZ.MathewsF.Di CeraE. (2008). Structural identification of the pathway of long-range communication in an allosteric enzyme. Proc. Natl. Acad. Sci. U. S. A.105, 1832–1837. 10.1073/pnas.0710894105
13
GasperP.FuglestadB.KomivesE.MarkwickP.McCammonJ. (2012). Allosteric networks in thrombin distinguish procoagulant vs. anticoagulant activities. Proc. Natl. Acad. Sci. U. S. A.109 (52), 21216–21222. 10.1073/pnas.1218414109
14
GibbsC.CoutreS.TsiangM.LiW.JainA.DunnK.et al (1995). Conversion of thrombin into an anticoagulant by protein engineering. Nature378, 413–416. 10.1038/378413a0
15
GruberA.CantwellA.Di CeraE.HansonS. (2002). The thrombin mutant W215A/E217A shows safe and potent anticoagulant and antithrombotic effects in vivo. J. Biol. Chem.277, 27581–27584. 10.1074/jbc.C200237200
16
GrütterM.PriestleJ.RahuelJ.GrossenbacherH.BodeW.HofsteengeJ.et al (1990). Crystal structure of the thrombin-hirudin complex: A novel mode of serine protease inhibition. EMBO J.9 (8), 2361–2365. 10.1002/j.1460-2075.1990.tb07410.x
17
GuintoE.Di CeraE. (1997). Critical role of W60d in thrombin allostery. Biophys. Chem.64 (1-3), 103–109. 10.1016/s0301-4622(96)02211-9
18
HandleyL.FuglestadB.StearnsK.TonelliM.FenwickR.MarkwickP.et al (2017). NMR reveals a dynamic allosteric pathway in thrombin. Sci. Rep.7, 39575. 10.1038/srep39575
19
HigginsD.LewisS.ShaferJ. (1983). Steady state kinetic parameters for the thrombin-catalyzed conversion of human fibrinogen to fibrin. J. Biol. Chem.258 (15), 9276–9282. 10.1016/s0021-9258(17)44663-1
20
KernsS.AgafonovR.ChoY.PontiggiaF.OttenR.PachovD.et al (2015). The energy landscape of adenylate kinase during catalysis. Nat. Struct. Mol. Biol.22, 124–131. 10.1038/nsmb.2941
21
Le BonniecB.EsmonC. (1991). Glu-192----Gln substitution in thrombin mimics the catalytic switch induced by thrombomodulin. Proc. Natl. Acad. Sci. U. S. A.88 (16), 7371–7375. 10.1073/pnas.88.16.7371
22
LechtenbergB.JohnsonD.FreundS.HuntingtonJ. (2010). NMR resonance assignments of thrombin reveal the conformational and dynamic effects of ligation. Proc. Natl. Acad. Sci. U. S. A.107, 14087–14092. 10.1073/pnas.1005255107
23
MandellJ.Baerga-OrtizA.AkashiS.TakioK.KomivesE. (2001). Solvent accessibility of the thrombin-thrombomodulin interface. J. Mol. Biol.306 (3), 575–589. 10.1006/jmbi.2000.4416
24
MandellJ.FalickA.KomivesE. (1998). Identification of protein-protein interfaces by decreased amide proton solvent accessibility. Proc. Natl. Acad. Sci. U. S. A.95 (25), 14705–14710. 10.1073/pnas.95.25.14705
25
MehD.SiebenlistK.MosessonM. (1996). Identification and characterization of the thrombin binding sites on fibrin. J. Biol. Chem.271 (38), 23121–23125. 10.1074/jbc.271.38.23121
26
MonteiroR.RapôsoJ.WisnerA.GuimarãesJ.BonC.ZingaliR. (1999). Allosteric changes of thrombin catalytic site induced by interaction of bothrojaracin with anion-binding exosites I and II. Biochem. Biophys. Res. Commun.262 (3), 819–822. 10.1006/bbrc.1999.1297
27
MusciG.BerlinerL.EsmonC. (1988). Evidence for multiple conformational changes in the active center of thrombin induced by complex formation with thrombomodulin: An analysis employing nitroxide spin-labels. Biochemistry27 (2), 769–773. 10.1021/bi00402a042
28
NiuW.ChenZ.Bush-PelcL.BahA.GandhiP.Di CeraE. (2009). Mutant N143P reveals how Na+ activates thrombin. J. Biol. Chem.284, 36175–36185. 10.1074/jbc.M109.069500
29
NiuW.ChenZ.GandhiP.VogtA.PozziN.PelcL.et al (2011). Crystallographic and kinetic evidence of allostery in a trypsin-like protease. Biochemistry50, 6301–6307. 10.1021/bi200878c
30
PeacockR.DavisJ.MarkwickP.KomivesE. (2018). Dynamic consequences of mutation of tryptophan 215 in thrombin. Biochemistry57 (18), 2694–2703. 10.1021/acs.biochem.8b00262
31
PeacockR.McGrannT.TonelliM.KomivesE. (2021). Serine protease dynamics revealed by NMR analysis of the thrombin-thrombomodulin complex. Sci. Rep.11 (1), 9354. 10.1038/s41598-021-88432-z
32
PinedaA.ChenZ.BahA.GarveyL.MathewsF.Di CeraE. (2006). Crystal structure of thrombin in a self-inhibited conformation. J. Biol. Chem.281 (43), 32922–32928. 10.1074/jbc.M605530200
33
StearnsD.KurosawaS.EsmonC. (1989). Microthrombomodulin. J. Biol. Chem.264 (6), 3352–3356. 10.1016/s0021-9258(18)94073-1
34
SteinerS.CastellinoF. (1985). Kinetic mechanism for stimulation by monovalent cations of the amidase activity of the plasma protease bovine activated protein C. Biochemistry24 (3), 609–617. 10.1021/bi00324a011
35
StrausD.RainesR.KawashimaE.KnowlesJ.GilbertW. (1985). Active site of triosephosphate isomerase: In vitro mutagenesis and characterization of an altered enzyme. Proc. Natl. Acad. Sci. U. S. A.82 (8), 2272–2276. 10.1073/pnas.82.8.2272
36
TreuheitN.BeachM.KomivesE. (2011). Thermodynamic compensation upon binding to exosite 1 and the active site of thrombin. Biochemistry50 (21), 4590–4596. 10.1021/bi2004069
37
van de LochtA.BodeW.HuberR.Le BonniecB.StoneS.EsmonC.et al (1997). The thrombin e192q-BPTI complex reveals gross structural rearrangements: Implications for the interaction with antithrombin and thrombomodulin. EMBO J.16 (11), 2977–2984. 10.1093/emboj/16.11.2977
38
WellsC.Di CeraE. (1992). Thrombin is a Na(+)-activated enzyme. Biochemistry31 (47), 11721–11730. 10.1021/bi00162a008
39
WuQ.TsiangM.Je.S. (1992). Localization of the single-stranded DNA binding site in the thrombin anion-binding exosite. J. Biol. Chem.267, 24408–24412. 10.1016/s0021-9258(18)35781-8
40
YeJ.EsmonN.EsmonC.JohnsonA. (1991). The active site of thrombin is altered upon binding to thrombomodulin. Two distinct structural changes are detected by fluorescence, but only one correlates with protein C activation. J. Biol. Chem.266 (34), 23016–23021. 10.1016/s0021-9258(18)54456-2
41
ZhangE.TulinskyA. (1997). The molecular environment of the Na+ binding site of thrombin. Biophys. Chem.63 (2-3), 185–200. 10.1016/s0301-4622(96)02227-2
Summary
Keywords
serine protease, thrombomodulin, isothermal titration calorimetry, HDX-MS (hydrogen/deuterium exchange-mass spectrometry), NMR relaxation dispersion
Citation
Komives EA (2023) Dynamic allostery in thrombin—a review. Front. Mol. Biosci. 10:1200465. doi: 10.3389/fmolb.2023.1200465
Received
05 April 2023
Accepted
19 June 2023
Published
29 June 2023
Volume
10 - 2023
Edited by
Emily Parker, Victoria University of Wellington, New Zealand
Reviewed by
Purnananda Guptasarma, Indian Institute of Science Education and Research Mohali, India
Daniel Deredge, University of Maryland, United States
Updates

Check for updates
Copyright
© 2023 Komives.
This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.
*Correspondence: Elizabeth A. Komives, ekomives@ucsd.edu
Disclaimer
All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.