Allosteric modulation of BPTI interaction with human alpha- and zeta-thrombin

Allosteric Modulation of Thrombin by Thrombomodulin: Insights from Logistic Regression and Statistical Analysis of Molecular Dynamics Simulations

Thrombomodulin (TM), a transmembrane receptor integral to the anticoagulant pathway, governs thrombin's substrate specificity via interaction with thrombin's anion-binding exosite I. Despite its established role, the precise mechanisms underlying this regulatory function are yet to be fully unraveled. In this study, we deepen the understanding of these mechanisms through eight independent 1 μs all-atom simulations, analyzing thrombin both in its free form and when bound to TM fragments TM456 and TM56. Our investigations revealed distinct and significant conformational changes in thrombin mediated by the binding of TM56 and TM456. While TM56 predominantly influences motions within exosite I, TM456 orchestrates coordinated alterations across various loop regions, thereby unveiling a multifaceted modulatory role that extends beyond that of TM56. A highlight of our study is the identification of critical hydrogen bonds that undergo transformations during TM56 and TM456 binding, shedding light on the pivotal allosteric influence exerted by TM4 on thrombin's structural dynamics. This work offers a nuanced appreciation of TM's regulatory role in blood coagulation, paving the way for innovative approaches in the development of anticoagulant therapies and expanding the horizons in oncology therapeutics through a deeper understanding of molecular interactions in the coagulation pathway.

Serine protease dynamics revealed by NMR analysis of the thrombin-thrombomodulin complex

Serine proteases catalyze a multi-step covalent catalytic mechanism of peptide bond cleavage. It has long been assumed that serine proteases including thrombin carry-out catalysis without significant conformational rearrangement of their stable two-β-barrel structure. We present nuclear magnetic resonance (NMR) and hydrogen deuterium exchange mass spectrometry (HDX-MS) experiments on the thrombin-thrombomodulin (TM) complex. Thrombin promotes procoagulative fibrinogen cleavage when fibrinogen engages both the anion binding exosite 1 (ABE1) and the active site. It is thought that TM promotes cleavage of protein C by engaging ABE1 in a similar manner as fibrinogen. Thus, the thrombin-TM complex may represent the catalytically active, ABE1-engaged thrombin. Compared to apo- and active site inhibited-thrombin, we show that thrombin-TM has reduced μs-ms dynamics in the substrate binding (S1) pocket consistent with its known acceleration of protein C binding. Thrombin-TM has increased μs-ms dynamics in a β-strand connecting the TM binding site to the catalytic aspartate. Finally, thrombin-TM had doublet peaks indicative of dynamics that are slow on the NMR timescale in residues along the interface between the two β-barrels. Such dynamics may be responsible for facilitating the N-terminal product release and water molecule entry that are required for hydrolysis of the acyl-enzyme intermediate.

Thermodynamic compensation upon binding to exosite 1 and the active site of thrombin

Several lines of experimental evidence including amide exchange and NMR suggest that ligands binding to thrombin cause reduced backbone dynamics. Binding of the covalent inhibitor dPhe-Pro-Arg chloromethyl ketone to the active site serine, as well as noncovalent binding of a fragment of the regulatory protein, thrombomodulin, to exosite 1 on the back side of the thrombin molecule both cause reduced dynamics. However, the reduced dynamics do not appear to be accompanied by significant conformational changes. In addition, binding of ligands to the active site does not change the affinity of thrombomodulin fragments binding to exosite 1; however, the thermodynamic coupling between exosite 1 and the active site has not been fully explored. We present isothermal titration calorimetry experiments that probe changes in enthalpy and entropy upon formation of binary ligand complexes. The approach relies on stringent thrombin preparation methods and on the use of dansyl-l-arginine-(3-methyl-1,5-pantanediyl)amide and a DNA aptamer as ligands with ideal thermodynamic signatures for binding to the active site and to exosite 1. Using this approach, the binding thermodynamic signatures of each ligand alone as well as the binding signatures of each ligand when the other binding site was occupied were measured. Different exosite 1 ligands with widely varied thermodynamic signatures cause a similar reduction in ΔH and a concomitantly lower entropy cost upon DAPA binding at the active site. The results suggest a general phenomenon of enthalpy-entropy compensation consistent with reduction of dynamics/increased folding of thrombin upon ligand binding to either the active site or exosite 1.

Does chemical cross-linking with NHS esters reflect the chemical equilibrium of protein-protein noncovalent interactions in solution?

Chemical cross-linking in combination with mass spectrometry has emerged as a powerful tool to study noncovalent protein complexes. Nevertheless, there are still many questions to answer. Does the amount of detected cross-linked complex correlate with the amount of protein complex in solution? In which concentration and affinity range is specific cross-linking possible? To answer these questions, we performed systematic cross-linking studies with two complexes, using the N-hydroxysuccinimidyl ester disuccinimidyl suberate (DSS): (1) NCoA-1 and mutants of the interacting peptide STAT6Y, covering a K(D) range of 30 nM to >25 μM, and (2) α-thrombin and basic pancreatic trypsin inhibitor (BPTI), a system that shows a buffer-dependent K(D) value between 100 and 320 μM. Samples were analyzed by matrix-assisted laser desorption/ionization mass spectrometry (MALDI-MS). For NCoA-1•STAT6Y, a good correlation between the amount of cross-linked species and the calculated fraction of complex present in solution was observed. Thus, chemical cross-linking in combination with MALDI-MS can be used to rank binding affinities. For the mid-affinity range up to about K(D) ≈ 25 μM, experiments with a nonbinding peptide and studies of the concentration dependence showed that only specific complexes undergo cross-linking with DSS. To study in which affinity range specific cross-linking can be applied, the weak α-thrombin•BPTI complex was investigated. We found that the detected complex is a nonspecifically cross-linked species. Consequently, based on the experimental approach used in this study, chemical cross-linking is not suitable for studying low-affinity complexes with K(D) >> 25 μM.

A segment of gamma ENaC mediates elastase activation of Na+ transport

The epithelial Na⁺ channel (ENaC) that mediates regulated Na⁺ reabsorption by epithelial cells in the kidney and lungs can be activated by endogenous proteases such as channel activating protease 1 and exogenous proteases such as trypsin and neutrophil elastase (NE). The mechanism by which exogenous proteases activate the channel is unknown. To test the hypothesis that residues on ENaC mediate protease-dependent channel activation wild-type and mutant ENaC were stably expressed in the FRT epithelial cell line using a tripromoter human ENaC construct, and protease-induced short-circuit current activation was measured in aprotinin-treated cells. The amiloride-sensitive short circuit current (I_Na) was stimulated by aldosterone (1.5-fold) and dexamethasone (8-fold). Dexamethasone-treated cells were used for all subsequent studies. The serum protease inhibitor aprotinin decreased baseline I_Na by approximately 50% and I_Na could be restored to baseline control values by the exogenous addition of trypsin, NE, and porcine pancreatic elastase (PE) but not by thrombin. All protease experiments were thus performed after exposure to aprotinin. Because NE recognition of substrates occurs with a preference for binding valines at the active site, several valines in the extracellular loops of α and γ ENaC were sequentially substituted with glycines. This scan yielded two valine residues in γ ENaC at positions 182 and 193 that resulted in inhibited responses to NE when simultaneously changed to other amino acids. The mutations resulted in decreased rates of activation and decreased activated steady-state current levels. There was an ∼20-fold difference in activation efficiency of NE against wild-type ENaC compared to a mutant with glycine substitutions at positions 182 and 193. However, the mutants remain susceptible to activation by trypsin and the related elastase, PE. Alanine is the preferred P₁ position residue for PE and substitution of alanine 190 in the γ subunit eliminated I_Na activation by PE. Further, substitution with a novel thrombin consensus sequence (LVPRG) beginning at residue 186 in the γ subunit (γ_Th) allowed for I_Na activation by thrombin, whereas wild-type ENaC was unresponsive. MALDI-TOF mass spectrometric evaluation of proteolytic digests of a 23-mer peptide encompassing the identified residues (T¹⁷⁶-S¹⁹⁸) showed that hydrolysis occurred between residues V193 and M194 for NE and between A190 and S191 for PE. In vitro translation studies demonstrated thrombin cleaved the γ_Th but not the wild-type γ subunit. These results demonstrate that γ subunit valines 182 and 193 are critical for channel activation by NE, alanine 190 is critical for channel activation by PE, and that channel activation can be achieved by inserting a novel thrombin consensus sequence. These results support the conclusion that protease binding and perhaps cleavage of the γ subunit results in ENaC activation.

Thrombomodulin Tightens the Thrombin Active Site Loops To Promote Protein C Activation

Thrombomodulin (TM) forms a 1:1 complex with thrombin. Whereas thrombin alone cleaves fibrinogen to make the fibrin clot, the thrombin-TM complex cleaves protein C to initiate the anticoagulant pathway. Crystallographic investigations of the complex between thrombin and TMEGF456 did not show any changes in the thrombin active site. Therefore, research has focused recently on how TM may provide a docking site for the protein C substrate. Previous work, however, showed that when the thrombin active site was occupied with substrate analogues labeled with fluorophores, the fluorophores responded differently to active (TMEGF1-6) versus inactive (TMEGF56) fragments of TM. To investigate this further, we have carried out amide H/(2)H exchange experiments on thrombin in the presence of active (TMEGF45) and inactive (TMEGF56) fragments of TM. Both on-exchange and off-exchange experiments show changes in the thrombin active site loops, some of which are observed only when the active TM fragment is bound. These results are consistent with the previously observed fluorescence changes and point to a mechanism by which TM changes the thrombin substrate specificity in favor of protein C cleavage.

Allosteric modulation of monomeric proteins*

Multimeric proteins (e.g. hemoglobin) are considered to be the prototypes of allosteric enzymes, whereas monomeric proteins (e.g. myoglobin) usually are assumed to be nonallosteric. However, the modulation of the functional properties of monomeric proteins by heterotropic allosteric effectors casts doubts on this assumption. Here, the allosteric properties of sperm whale myoglobin, human serum albumin, and human α-thrombin, generally considered as molecular models of monomeric proteins, are summarized.

A Natural Prothrombin Mutant Reveals an Unexpected Influence of A-chain Structure on the Activity of Human α-Thrombin

We have recently identified in two unrelated patients with bleeding tendency a homozygous mutation causing a deletion of one of the two contiguous Lys(9)/Lys(10) residues in the A-chain of alpha-thrombin (DeltaK9). We used in vitro expression analysis to clarify the role of the deletion of Lys(9) or Lys(10) in the thrombin function. The k(cat)/K(m) value of the hydrolysis by DeltaK9 of the synthetic substrate Phe-Pip-Arg-p-nitroanilide (where Pip represents l-pipecolyl) and fibrinopeptide A was 18- and 60-fold lower, respectively, compared with wild type (WT). Interaction with antithrombin was also reduced in the mutant, the association rate being about 20-fold lower than in the WT thrombin. The sensitivity to sodium ion of DeltaK9 was found significantly attenuated compared with the WT form. DeltaK9 has a very weak platelet-activating capacity, attributed to a severely defective PAR1 interaction, whereas the binding to the platelet glycoprotein Ibalpha was unaffected. Likewise, the interaction with protein C was severely impaired, whereas interaction with thrombomodulin had a normal K(d) value. At variance with these findings, both low affinity (basic pancreatic trypsin inhibitor) and high affinity (N-alpha-[2-naphthylsulfonyl-glycyl]-4-amidinophenylalanine-piperidide) thrombin inhibitors displayed a better binding to DeltaK9 than to the WT form, indicating a better accommodation of these inhibitors into the catalytic pocket of DeltaK9. A molecular dynamics simulation of the DeltaK9 thrombin in full explicit water solvent provided support to the role of the A-chain in affecting conformation and catalytic properties of the B-chain, especially in some insertion loops of the enzyme, such as the 60-loop, as well as in the geometry of the catalytic triad residues.

Prothrombinase assembly and S1 site occupation restore the catalytic activity of FXa impaired by mutation at the sodium-binding site

Two loop segments (183-189 and 221-225) in the protease domain of factor Xa contribute to the formation of a Na(+)-binding site. Studies with factor Xa indicate that binding of a single Na(+) ion to this site influences its activity by altering the S1 specificity site, and substitution of Tyr(225) with Pro diminishes sensitivity to Na(+). Using full-length factor Xa(Y225P), the allosteric relationship between the Na(+) site and other structural determinants in factor Xa and prothrombinase was investigated. Direct binding and kinetic measurements with probes that target the S1 specificity pocket indicate that assembly of the mutant in prothrombinase corrected the impaired binding of these probes observed with free factor Xa(Y225P). This appears to result from the apparent allosteric linkage between the factor Va, S1, and Na(+)-binding sites, since binding of the cofactor to membrane-bound factor Xa(Y225P) enhances binding at the S1 site and vice versa. Additional studies revealed that the internal salt bridge (Ile(16)-Asp(194)) of factor Xa(Y225P) is partially destabilized, a process that is reversible upon occupation of the S1 site. The data establish that alterations at the factor Xa Na(+)-binding site shift the zymogen-protease equilibrium to a more zymogen-like state, and as a consequence binding of S1-directed probes and factor Va are adversely affected. Therefore, the zymogen-like characteristics of factor Xa(Y225P) have allowed for the apparent allosteric linkage between the S1, factor Va, and Na(+) sites to become evident and has provided insight into the structural transitions which accompany the conversion of factor X to factor Xa.

Regulation of blood coagulation

The protein C anticoagulant pathway converts the coagulation signal generated by thrombin into an anticoagulant response through the activation of protein C by the thrombin-thrombomodulin (TM) complex. The activated protein C (APC) thus formed interacts with protein S to inactivate two critical coagulation cofactors, factors Va and VIIIa, thereby dampening further thrombin generation. The proposed mechanisms by which TM switches the specificity of thrombin include conformational changes in thrombin, blocking access of normal substrates to thrombin and providing a binding site for protein C. The function of protein S appears to be to alter the cleavage site preferences of APC in factor Va, probably by changing the distance of the active site of APC relative to the membrane surface. The clinical relevance of this pathway is now established through the identification of deficient individuals with severe thrombotic complications and through the analysis of families with partial deficiencies in these components and an increased thrombotic tendency. One possible reason that even partial deficiencies are a thrombotic risk is that the function of the pathway can be down-regulated by inflammatory mediators. For instance, clinical studies have shown that the extent to which protein C levels decrease in patients with septic shock is predictive of a negative outcome. Initial clinical studies suggest that supplementation with protein C may be useful in the treatment of acute inflammatory diseases such as sepsis.

The role of Glu(192) in the allosteric control of the S(2)' and S(3)' subsites of thrombin

Thrombin is an allosteric protease controlled through exosites flanking the catalytic groove. Binding of a peptide derived from hirudin (Hir(52-65)) and/or of heparin to these opposing exosites alters catalysis. We have investigated the contribution of subsites S(2)' and S(3)' to this allosteric transition by comparing the hydrolysis of two sets of fluorescence-quenched substrates having all natural amino acids at positions P(2)' and P(3)'. Regardless of the amino acids, Hir(52-65) decreased, and heparin increased the k(cat)/K(m) value of hydrolysis by thrombin. Several lines of evidence have suggested that Glu(192) participates in this modulation. We have examined the role of Glu(192) by comparing the catalytic activity of thrombin and its E192Q mutant. Mutation substantially diminishes the selectivity of thrombin. The substrate with the "best" P(2)' residue was cleaved with a k(cat)/K(m) value only 49 times higher than the one having the "least favorable" P(2)' residue (versus 636-fold with thrombin). Mutant E192Q also lost the strong preference of thrombin for positively charged P(3)' residues and its strong aversion for negatively charged P(3)' residues. Furthermore, both Hir(52-65) and heparin increased the k(cat)/K(m) value of substrate hydrolysis. We conclude that Glu(192) is critical for the P(2)' and P(3)' specificities of thrombin and for the allostery mediated through exosite 1.