Molecular dissection of Na+ binding to thrombin

To be, or not to be cleaved: Directed evolution of a canonical serine protease inhibitor against active and inactive protease pair identifies binding loop residue critical for prevention of proteolytic cleavage

Canonical serine protease inhibitor proteins occupy the substrate-binding groove of their target enzyme via a surface loop. Unlike true substrates, inhibitors are cleaved by the target protease extremely slowly. Here, we applied an unbiased directed evolution approach to investigate which loop residues hamper proteolytic cleavage while maintaining high-affinity binding. As a protease inhibitor model system, we used human chymotrypsin C (CTRC) and Schistocerca gregaria protease inhibitor 2 (SGPI-2). We created an SGPI-2 library displayed on M13 phage by randomizing the binding loop amino acid positions, with the exception of the structurally indispensable Cys residues. We selected binding phage clones against active CTRC and the inactive mutant Ser195Ala. All CTRC-selected binders inhibited CTRC activity and also bound to the inactive Ser195Ala mutant, but the Ser195Ala-selected clones proved to be either inhibitors or substrates of active CTRC. Substrate-like behavior of SGPI-2 variants was associated with the absence of the P2 Thr, the residue next to the specificity determinant P1 amino acid. The selected SGPI-2 variants containing a P2 Thr bound strongly to CTRC even if the other loop residues deviated from the optimal inhibitory consensus sequence. In the absence of a P2 Thr, however, SGPI-2 variants became substrates unless all other loop residues were optimal for binding. Structural modeling confirmed that P2 Thr is important for organizing a stabilizing H-bond network. The observations indicate that binding loops of canonical serine protease inhibitors evolved amino acids not only to support tight binding to the target enzyme but also to inhibit proteolytic cleavage.

Replacement of a single residue changes the primary specificity of thrombin

BACKGROUND

Thrombin prefers substrates carrying Arg at the site of cleavage (P1) because of the presence of D189 in the primary specificity (S1) pocket but can also cleave substrates carrying Phe at P1. The structural basis of this property is unknown.

OBJECTIVES

Solve the X-ray structure of thrombin bound to a ligand carrying Phe at P1 and investigate the effects of replacing D189.

METHODS

X-ray crystallography is used to solve the structure of thrombin bound to the irreversible inhibitor H-D-Phe-Pro-Phe-CH2Cl (PPPCK). Residue D189 is mutated to Ala, Lys, Phe, and Ser.

RESULTS

The X-ray structure of the thrombin-PPPCK complex is solved at 2.5 Å resolution and compared to the structure of thrombin bound to H-D-Phe-Pro-Arg-CH2Cl (PPACK). PPPCK binds to thrombin in a conformation similar to that of PPACK, but Phe at P1 makes no contacts with D189. Replacement of D189 with Ala, Lys, Phe, or Ser reverses both substrate preference and stability enhancement from Arg to Phe.

CONCLUSION

D189 in the S1 pocket confers thrombin "trypsin-like" specificity for Arg at P1. However, the S1 pocket is wide enough to also enable "chymotrypsin-like" specificity for Phe at P1. Consistent with these structural features, a single amino acid replacement (D189A) switches thrombin specificity from trypsin-like to chymotrypsin-like, converting the substrate preference from H-D-Phe-Pro-Arg-p-nitroanilide to H-D-Phe-Pro-Phe-p-nitroanilide and preferential stability enhancement from PPACK to PPPCK. The observation that thrombin specificity is controlled mainly by a single residue establishes a new paradigm in the field of trypsin-like proteases.

Effect of prothrombin Belgrade mutation, causing antithrombin resistance, on fibrin clot properties

INTRODUCTION

Prothrombin Belgrade mutation is the result of the c.1787G>A substitution in the prothrombin gene. It is located in the antithrombin and sodium binding site and leads to impaired inactivation of thrombin by antithrombin, resulting in antithrombin resistance and thrombotic disorders. However, it negatively affects sodium binding and may have hypocoagulant effects. Considering that prothrombin Belgrade mutation mechanism is still not fully elucidated and that sodium binding is important for thrombin affinity towards fibrinogen, our aim was to determine whether this mutation affects fibrin clot formation and lysis.

METHODS

Using HEK293T cell line, recombinant wild type and mutated prothrombin were generated by transient transfection. Samples that correspond to plasma of a non-carrier, heterozygous and homozygous carriers were reconstituted using prothrombin deficient plasma and recombinant proteins. Reconstituted samples were used in OHP assay (Overall Hemostasis Potential) to determine kinetic profiles of coagulation and fibrinolysis. Clot turbidity assay was performed to observe kinetics of clot formation and lysis more closely. Fibrin clots formed in reconstituted plasma samples were analyzed by confocal microscopy to determine density of fibrin network. Fibrin clots were additionally observed using electron microscopy to determine thickness of individual fibrin fibers.

RESULTS

No significant difference found in OHP, OCP, OFP, and fibrin network density between wild type, heterozygous, and homozygous carrier reconstituted plasma samples. There were significant differences between samples for slope and slope time parameters in kinetic profiles and fibrin fiber thickness.

CONCLUSIONS

Results indicate that prothrombin Belgrade mutation has no significant impact on fibrinolysis, however it may affect kinetics of clot formation and its architecture.

The Prothrombin-Prothrombinase Interaction

The hemostatic response to vascular injury entails a sequence of proteolytic events where several inactive zymogens of the trypsin family are converted to active proteases. The cascade starts with exposure of tissue factor from the damaged endothelium and culminates with conversion of prothrombin to thrombin in a reaction catalyzed by the prothrombinase complex composed of the enzyme factor Xa, cofactor Va, Ca2+, and phospholipids. This cofactor-dependent activation is paradigmatic of analogous reactions of the blood coagulation and complement cascades, which makes elucidation of its molecular mechanism of broad significance to the large class of trypsin-like zymogens to which prothrombin belongs. Because of its relevance as the most important reaction in the physiological response to vascular injury, as well as the main trigger of pathological thrombotic complications, the mechanism of prothrombin activation has been studied extensively. However, a molecular interpretation of this mechanism has become available only recently from important developments in structural biology. Here we review current knowledge on the prothrombin-prothrombinase interaction and outline future directions for the study of this key reaction of the coagulation cascade.

Molecular mechanism by which spider-driving peptide potentiates coagulation factors

Hemostasis is a crucial process that quickly forms clots at injury sites to prevent bleeding and infections. Dysfunctions in this process can lead to hemorrhagic disorders, such as hemophilia and thrombocytopenia purpura. While hemostatic agents are used in clinical treatments, there is still limited knowledge about potentiators targeting coagulation factors. Recently, LCTx-F2, a procoagulant spider-derived peptide, was discovered. This study employed various methods, including chromogenic substrate analysis and dynamic simulation, to investigate how LCTx-F2 enhances the activity of thrombin and FXIIa. Our findings revealed that LCTx-F2 binds to thrombin and FXIIa in a similar manner, with the N-terminal penetrating the active-site cleft of the enzymes and the intermediate section reinforcing the peptide-enzyme connection. Interestingly, the C-terminal remained at a considerable distance from the enzymes, as evidenced by the retention of affinity for both enzymes using truncated peptide T-F2. Furthermore, results indicated differences in the bonding relationship of critical residues between thrombin and FXIIa, with His13 facilitating binding to thrombin and Arg7 being required for binding to FXIIa. Overall, our study sheds light on the molecular mechanism by which LCTx-F2 potentiates coagulation factors, providing valuable insights that may assist in designing drugs targeting procoagulation factors.

Super-high procoagulant activity of gecko thrombin: A gift from sky dragon

AIMS

Gecko, the "sky dragon" named by Traditional Chinese Medicine, undergoes rapid coagulation and scarless regeneration following tail amputation in the natural ecology, providing a perfect opportunity to develop the efficient and safe drug for blood clotting. Here, gecko thrombin (gthrombin) was recombinantly prepared and comparatively studied on its procoagulant activity.

METHODS

The 3D structure of gthrombin was constructed using the homology modeling method of I-TASSER. The active gthrombin was prepared by the expression of gecko prethrombin-2 in 293 T cells, followed by purification with Ni2+ -chelating column chromatography prior to activation by snake venom-derived Ecarin. The enzymatic activities of gthrombin were assayed by hydrolysis of synthetic substrate S-2238 and the fibrinogen clotting. The vulnerable nerve cells were used to evaluate the toxicity of gthrombin at molecular and cellular levels.

RESULTS

The active recombinant gthrombin showed super-high catalytic and fibrinogenolytic efficiency than those of human under different temperatures and pH conditions. In addition, gthrombin made nontoxic effects on the central nerve cells including neurons, contrary to those of mammalian counterparts, which contribute to neuronal damage, astrogliosis, and demyelination.

CONCLUSIONS

A super-high activity but safe procoagulant candidate drug was identified from reptiles, which provided a promising perspective for clinical application in rapid blood clotting.

Augmented Degradation of Factors VIII and IX in the Intermittent Movement State

The most common clinical presentation of hemophilia A and hemophilia B is bleeding in large joints and striated muscles. It is unclear why bleeding has a predilection to affect joints and muscles. As muscles and joints are involved in intermittent movement, we explored whether this phenomenon could be associated with an impact on factor VIII and IX levels. Purified proteins and a mouse model were assessed using coagulation assays, Western blot analysis and immuno-staining. Movement caused an increase in thrombin activity and a decrease in factor VIII and factor IX activity. The decrease in factor VIII activity was more significant in the presence of thrombin and during movement. Under movement condition, sodium ions appeared to enhance the activity of thrombin that resulted in decreased factor VIII activity. Unlike factor VIII, the reduction in factor IX levels in the movement condition was thrombin-independent. High factor VIII levels were found to protect factor IX from degradation and vice versa. In mice that were in movement, factor VIII and IX levels decreased in the microcirculation of the muscle tissue compared with other tissues and to the muscle tissue at rest. Movement had no effect on von Willebrand factor levels. Movement induces reduction in factor VIII and IX levels. It enables an increase in the binding of sodium ions to thrombin leading to enhanced thrombin activity and augmented degradation of factor VIII. These data suggest a potential mechanism underlying the tendency of hemophilia patients to bleed in muscles and joints.

Simulations suggest double sodium binding induces unexpected conformational changes in thrombin

Thrombin is a Na[Formula: see text]-activated serine protease existing in two forms targeted to procoagulant and anticoagulant activities, respectively. There is one Na[Formula: see text]-binding site that has been the focus of the study of the thrombin. However, molecular dynamics (MD) simulations suggest that there might be actually two Na[Formula: see text]-binding sites in thrombin and that Na[Formula: see text] ions can even bind to two sites simultaneously. In this study, we performed 12 independent 2-µs all-atom MD simulations for the wild-type (WT) thrombin and we studied the effects of the different Na[Formula: see text] binding modes on thrombin. From the root-mean-square fluctuations (RMSF) for the [Formula: see text]-carbons, we see that the atomic fluctuations mainly change in the 60s, 170s, and 220s loops, and the connection (residue 167 to 170). The correlation matrices for different binding modes suggest regions that may play an important role in thrombin's allosteric response and provide us a possible allosteric pathway for the sodium binding. Amorim-Hennig (AH) clustering tells us how the structure of the regions of interest changes on sodium binding. Principal component analysis (PCA) shows us how the different regions of thrombin change conformation together with sodium binding. Solvent-accessible surface area (SASA) exposes the conformational change in exosite I and catalytic triad. Finally, we argue that the double binding mode might be an inactive mode and that the kinetic scheme for the Na[Formula: see text] binding to thrombin might be a multiple-step mechanism rather than a 2-step mechanism.

The role of thrombin in haemostasis

Thrombin is a multifunctional serine protease generated in injured cells. The generation of thrombin in coagulation plays a central role in the functioning of haemostasis. The last enzyme in the coagulation cascade is thrombin, with the function of cleaving fibrinogen to fibrin, which forms the fibrin clot of a haemostatic plug. Although thrombin primarily converts fibrinogen to fibrin, it also has many other positive regulatory effects on coagulation. Thrombin has procoagulant, inflammatory, cellular proliferation and anticoagulant effects. In coagulation system, thrombin has two very distinct roles. Firstly, it acts as a procoagulant when it converts fibrinogen into an insoluble fibrin clot, activates factor (F) XIII, activates thrombin activatable fibrinolysis inhibitor (TAFI) and activates FV, FVIII and FXI. Thrombin also enhances platelet adhesion by inactivating a disintegrin and metalloprotease with thrombospondin type1 motif (ADAMTS13). However, when thrombin activates protein C, it acts as an anticoagulant. A natural anticoagulant pathway that supplies regulation of the blood coagulation system contains protein C, which is the key component. This is accomplished by the specific proteolytic inactivation of FV and FVIII. In this review, the multiple roles of thrombin in the haemostatic response to injury are studied in addition to the cofactors that determine thrombin activity and how thrombin activity is thought to be coordinated.

A systematic approach for evaluating the role of surface-exposed loops in trypsin-like serine proteases applied to the 170 loop in coagulation factor VIIa

Proteases play a major role in many vital physiological processes. Trypsin-like serine proteases (TLPs), in particular, are paramount in proteolytic cascade systems such as blood coagulation and complement activation. The structural topology of TLPs is highly conserved, with the trypsin fold comprising two β-barrels connected by a number of variable surface-exposed loops that provide a surprising capacity for functional diversity and substrate specificity. To expand our understanding of the roles these loops play in substrate and co-factor interactions, we employ a systematic methodology akin to the natural truncations and insertions observed through evolution of TLPs. The approach explores a larger deletion space than classical random or directed mutagenesis. Using FVIIa as a model system, deletions of 1–7 amino acids through the surface exposed 170 loop, a vital allosteric regulator, was introduced. All variants were extensively evaluated by established functional assays and computational loop modelling with Rosetta. The approach revealed detailed structural and functional insights recapitulation and expanding on the main findings in relation to 170 loop functions elucidated over several decades using more cumbersome crystallization and single deletion/mutation methodologies. The larger deletion space was key in capturing the most active variant, which unexpectedly had a six-amino acid truncation. This variant would have remained undiscovered if only 2–3 deletions were considered, supporting the usefulness of the methodology in general protease engineering approaches. Our findings shed further light on the complex role that surface-exposed loops play in TLP function and supports the important role of loop length in the regulation and fine-tunning of enzymatic function throughout evolution.

Activation and thermal stabilization of a recombinant γ-glutamyltranspeptidase from Bacillus licheniformis ATCC 27811 by monovalent cations

The regulation of enzyme activity through complexation with certain metal ions plays an important role in many biological processes. In addition to divalent metals, monovalent cations (MVCs) frequently function as promoters for efficient biocatalysis. Here, we examined the effect of MVCs on the enzymatic catalysis of a recombinant γ-glutamyltranspeptidase (BlrGGT) from Bacillus licheniformis ATCC 27,811 and the application of a metal-activated enzyme to L-theanine synthesis. The transpeptidase activity of BlrGGT was enhanced by Cs⁺ and Na⁺ over a broad range of concentrations with a maximum of 200 mM. The activation was essentially independent of the ionic radius, but K⁺ contributed the least to enhancing the catalytic efficiency. The secondary structure of BlrGGT remained mostly unchanged in the presence of different concentrations of MVCs, but there was a significant change in its tertiary structure under the same conditions. Compared with the control, the half-life (t_1/2) of the Cs⁺-enriched enzyme at 60 and 65 °C was shown to increase from 16.3 and 4.0 min to 74.5 and 14.3 min, respectively. The simultaneous addition of Cs⁺ and Mg²⁺ ions exerted a synergistic effect on the activation of BlrGGT. This was adequately reflected by an improvement in the conversion of substrates to L-theanine by 3.3-15.1% upon the addition of 200 mM MgCl₂ into a reaction mixture comprising the freshly desalted enzyme (25 μg/mL), 250 mM L-glutamine, 600 mM ethylamine, 200 mM each of the MVCs, and 50 mM borate buffer (pH 10.5). Taken together, our results provide interesting insights into the complexation of MVCs with BlrGGT and can therefore be potentially useful to the biocatalytic production of naturally occurring γ-glutamyl compounds. KEY POINTS: • The transpeptidase activity of B. licheniformis γ-glutamyltranspeptidase can be activated by monovalent cations. • The thermal stability of the enzyme was profoundly increased in the presence of 200 mM Cs⁺. • The simultaneous addition of Cs⁺and Mg²⁺ions to the reaction mixture improves L-theanine production.

Structural insights into the activation of blood coagulation factor XI zymogen by thrombin: A computational molecular dynamics study

Activation of human blood coagulation factor XI zymogen to factor XIa plays a significant role in the upstream coagulation pathway, in which factor XIa activates factor IX zymogen. The mechanistic details of the proteolytic activation of factor XI by the activating enzyme thrombin are not well-understood at atomic level. In this study, we employed a combination of molecular docking and microsecond time-scale molecular dynamics simulations to identify the key regions of interaction between fXI and thrombin. The activating complex between the substrate and enzyme was modeled to represent the initial acylation step of the serine-protease hydrolysis mechanism. The proposed solution structural complex, fIX:fIIa, obtained from 3 microseconds of MD refinement, suggests that the activation of factor XI is mediated by thrombin's anion binding exosite-II interactions with A3 and A4 domains. We predict that the two positively charged arginine residues (Arg409 and Arg413) in the exosite-2 region, the β- and γ-insertion loops of thrombin play an important structural role in the initial activating complex between fXI and thrombin.

The active site region plays a critical role in Na+ binding to thrombin

The catalytic activity of thrombin and other enzymes of the blood coagulation and complement cascades is enhanced significantly by binding of Na+ to a site >15 Å away from the catalytic residue S195, buried within the 180 and 220 loops that also contribute to the primary specificity of the enzyme. Rapid kinetics support a binding mechanism of conformational selection where the Na+-binding site is in equilibrium between open (N) and closed (N∗) forms and the cation binds selectively to the N form. Allosteric transduction of this binding step produces enhanced catalytic activity. Molecular details on how Na+ gains access to this site and communicates allosterically with the active site remain poorly defined. In this study, we show that the rate of the N∗→N transition is strongly correlated with the analogous E∗→E transition that governs the interaction of synthetic and physiologic substrates with the active site. This correlation supports the active site as the likely point of entry for Na+ to its binding site. Mutagenesis and structural data rule out an alternative path through the pore defined by the 180 and 220 loops. We suggest that the active site communicates allosterically with the Na+ site through a network of H-bonded water molecules that embeds the primary specificity pocket. Perturbation of the mobility of S195 and its H-bonding capabilities alters interaction with this network and influences the kinetics of Na+ binding and allosteric transduction. These findings have general mechanistic relevance for Na+-activated proteases and allosteric enzymes.

Probing activation-driven changes in coagulation factor IX by mass spectrometry

BACKGROUND

Activated factor IX (FIXa) is an inefficient enzyme that needs activated factor VIII (FVIII) for full activity. Recently, we identified a network of FVIII-driven changes in FIXa employing hydrogen-deuterium eXchange mass spectrometry (HDX-MS). Some changes also occurred in active-site inhibited FIXa, but others were not cofactor-driven, in particular those within the 220-loop (in chymotrypsin numbering).

OBJECTIVE

The aim of this work is to better understand the zymogen-to-enzyme transition in FIX, with specific focus on substrate-driven changes at the catalytic site.

METHODS

Footprinting mass spectrometry by HDX and Tandem-Mass Tags (TMT) labelling were used to explore changes occurring upon the conversion from FIX into FIXa. Mutagenesis and kinetic studies served to assess the role of the 220-loop.

RESULTS

HDX-MS displayed remarkably few differences between FIX and FIXa. In comparison with FIX, FIXa did exhibit decreased deuterium uptake at the N-terminus region. This was more prominent when the FIXa active site was occupied by an irreversible inhibitor. TMT-labelling showed that the N-terminus is largely protected from labelling, and that inhibitor binding increases protection to a minor extent. Occupation of the active site also reduced deuterium uptake within the 220-loop backbone. Mutagenesis within the 220-loop revealed that a putative H-bond network contributes to FIXa activity. TMT labeling of the N-terminus suggested that these 220-loop variants are more zymogen-like than wild-type FIXa.

CONCLUSION

In the absence of cofactor and substrate, FIXa is predominantly zymogen-like. Stabilization in its enzyme-like form involves, apart from FVIII-binding, also interplay between the 220-loop, N-terminus, and the substrate binding site.

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.

Proton Bridging in Catalysis by and Inhibition of Serine Proteases of the Blood Cascade System

Inquiries into the participation of short hydrogen bonds in stabilizing transition states and intermediate states in the thrombin, factor Xa, plasmin and activated protein C-catalyzed reactions revealed that specific binding of effectors at Sn, n = 1-4 and S'n, n = 1-3 and at remote exosites elicit complex patterns of hydrogen bonding and involve water networks. The methods employed that yielded these discoveries include; (1) kinetics, especially partial or full kinetic deuterium solvent isotope effects with short cognate substrates and also with the natural substrates, (2) kinetic and structural probes, particularly low-field high-resolution nuclear magnetic resonance (1H NMR), of mechanism-based inhibitors and substrate-mimic peptide inhibitors. Short hydrogen bonds form at the transition states of the catalytic reactions at the active site of the enzymes as they do with mechanism-based covalent inhibitors of thrombin. The emergence of short hydrogen bonds at the binding interface of effectors and thrombin at remote exosites has recently gained recognition. Herein, I describe our contribution, a confirmation of this discovery, by low-field 1H NMR. The principal conclusion of this review is that proton sharing at distances below the sum of van der Waals radii of the hydrogen and both donor and acceptor atoms contribute to the remarkable catalytic prowess of serine proteases of the blood clotting system and other enzymes that employ acid-base catalysis. Proton bridges also play a role in tight binding in proteins and at exosites, i.e., allosteric sites, of enzymes.

Light Chain Mutation Effects on the Dynamics of Thrombin

Thrombin plays an important role in the process of hemostasis and blood coagulation. Studies in thrombin can help us find ways to treat cancer because thrombin is able to reduce the characteristic hypercoagulability of cancer. Thrombin is composed of two chains, the light chain and the heavy chain. The function of the heavy chain has been largely explored, while the function of the light chain was obscured until several disease-associated mutations in the light chain come to light. In this study, we want to explore the dynamic and conformation effects of mutations on the light chain further to determine possible associations between mutation, conformational changes, and disease. The study, which is a follow-up for our studies on apo thrombin and the mutant, ΔK9, mainly focuses on the mutants E8K and R4A. E8K is a disease-associated mutation, and R4A is used to study the role of Arg4, which is suggested experimentally to play a critical role for thrombin's catalytic activities. We performed five all-atom one microsecond-scale molecular dynamics (MD) simulations for both E8K and R4A, and quantified the changes in the conformational ensemble of the mutants. From the root-mean-square fluctuations (RMSF) for the α-carbons, we find that the atomic fluctuations change in the mutants in the 60s loop and γ loop. The correlation coefficients for the α-carbons indicate that the correlation relation for atom-pairs in the protein is also impacted. The clustering analysis and the principal component analysis (PCA) consistently tell us that the catalytic pocket and the regulatory loops are destabilized by the mutations. We also find that there are two binding modes for Na+ by clustering the vector difference between the Na+ ions and the 220s loop. After further analysis, we find that there is a relation between the Na+ binding and the rigidification of the γ loop, which may shed light on the mysterious role of the γ loop in thrombin.

Mechanisms of ligand binding

Many processes in chemistry and biology involve interactions of a ligand with its molecular target. Interest in the mechanism governing such interactions has dominated theoretical and experimental analysis for over a century. The interpretation of molecular recognition has evolved from a simple rigid body association of the ligand with its target to appreciation of the key role played by conformational transitions. Two conceptually distinct descriptions have had a profound impact on our understanding of mechanisms of ligand binding. The first description, referred to as induced fit, assumes that conformational changes follow the initial binding step to optimize the complex between the ligand and its target. The second description, referred to as conformational selection, assumes that the free target exists in multiple conformations in equilibrium and that the ligand selects the optimal one for binding. Both descriptions can be merged into more complex reaction schemes that better describe the functional repertoire of macromolecular systems. This review deals with basic mechanisms of ligand binding, with special emphasis on induced fit, conformational selection, and their mathematical foundations to provide rigorous context for the analysis and interpretation of experimental data. We show that conformational selection is a surprisingly versatile mechanism that includes induced fit as a mathematical special case and even captures kinetic properties of more complex reaction schemes. These features make conformational selection a dominant mechanism of molecular recognition in biology, consistent with the rich conformational landscape accessible to biological macromolecules being unraveled by structural biology.

Thrombin: An Approach to Developing a Higher-Order Reference Material and Reference Measurement Procedure for Substance Identity, Amount, and Biological Activities

Thrombin, the proteolytic enzyme that catalyzes the transformation of soluble fibrinogen to the polymerized fibrin clot, participates in multiple reactions in blood coagulation in addition to the clotting reaction. Although reference materials have existed for many years, structural characterization and measurement of biological activity have never been sufficient to permit claims of clear metrological traceability for the thrombin preparations. Our current state-of-the-art methods for protein characterization and determination of the catalytic properties of thrombin now make it practical to develop and characterize a metrologically acceptable reference material and reference measurement procedure for thrombin. Specifically, α-thrombin, the biologically produced protease formed during prothrombin activation, is readily available and has been extensively characterized. Dependences of thrombin proteolytic and peptide hydrolytic activities on a variety of substrates, pH, specific ions, and temperature are established, although variability remains for the kinetic parameters that describe thrombin enzymatic action. The roles of specific areas on the surface of the thrombin molecule (exosites) in substrate recognition and catalytic efficiency are described and characterized. It is opportune to develop reference materials of high metrological order and technical feasibility. In this article, we review the properties of α-thrombin important for its preparation and suggest an approach suitable for producing a reference material and a reference measurement procedure that is sensitive to thrombin’s catalytic competency on a variety of substrates.

19F NMR reveals the conformational properties of free thrombin and its zymogen precursor prethrombin-2

The conformational properties of trypsin-like proteases and their zymogen forms remain controversial because of a lack of sufficient information on their free forms. Specifically, it is unclear whether the free protease is zymogen-like and shifts to its mature form upon a ligand-induced fit or exists in multiple conformations in equilibrium from which the ligand selects the optimal fit via conformational selection. Here we report the results of ¹⁹F NMR measurements that reveal the conformational properties of a protease and its zymogen precursor in the free form. Using the trypsin-like, clotting protease thrombin as a relevant model system, we show that its conformation is quite different from that of its direct zymogen precursor prethrombin-2 and more similar to that of its fully active Na⁺-bound form. The results cast doubts on recent hypotheses that free thrombin is zymogen-like and transitions to protease-like forms upon ligand binding. Rather, they validate the scenario emerged from previous findings of X-ray crystallography and rapid kinetics supporting a pre-existing equilibrium between open (E) and closed (E*) forms of the active site. In this scenario, prethrombin-2 is more dynamic and exists predominantly in the E* form, whereas thrombin is more rigid and exists predominantly in the E form. Ligand binding to thrombin takes place exclusively in the E form without significant changes in the overall conformation. In summary, these results disclose the structural architecture of the free forms of thrombin and prethrombin-2, consistent with an E*-E equilibrium and providing no evidence that free thrombin is zymogen-like.

Role of the I16-D194 ionic interaction in the trypsin fold

Activity in trypsin-like proteases is the result of proteolytic cleavage at R15 followed by an ionic interaction that ensues between the new N terminus of I16 and the side chain of the highly conserved D194. This mechanism of activation, first proposed by Huber and Bode, organizes the oxyanion hole and primary specificity pocket for substrate binding and catalysis. Using the clotting protease thrombin as a relevant model, we unravel contributions of the I16-D194 ionic interaction to Na⁺ binding, stability of the transition state and the allosteric E*-E equilibrium of the trypsin fold. The I16T mutation abolishes the I16-D194 interaction and compromises the architecture of the oxyanion hole. The D194A mutation also abrogates the I16-D194 interaction but, surprisingly, has no effect on the architecture of the oxyanion hole that remains intact through a new H-bond established between G43 and G193. In both mutants, loss of the I16-D194 ionic interaction compromises Na⁺ binding, reduces stability of the transition state, collapses the 215–217 segment into the primary specific pocket and abrogates the allosteric E*-E equilibrium in favor of a rigid conformation that binds ligand at the active site according to a simple lock-and-key mechanism. These findings refine the structural role of the I16-D194 ionic interaction in the Huber-Bode mechanism of activation and reveal a functional linkage with the allosteric properties of the trypsin fold like Na⁺ binding and the E*-E equilibrium.

Surface loops of trypsin-like serine proteases as determinants of function

Trypsin and chymotrypsin-like serine proteases from family S1 (clan PA) constitute the largest protease group in humans and more generally in vertebrates. The prototypes chymotrypsin, trypsin and elastase represent simple digestive proteases in the gut, where they cleave nearly any protein. Multidomain trypsin-like proteases are key players in the tightly controlled blood coagulation and complement systems, as well as related proteases that are secreted from diverse immune cells. Some serine proteases are expressed in nearly all tissues and fluids of the human body, such as the human kallikreins and kallikrein-related peptidases with specialization for often unique substrates and accurate timing of activity. HtrA and membrane-anchored serine proteases fulfill important physiological tasks with emerging roles in cancer. The high diversity of all family members, which share the tandem β-barrel architecture of the chymotrypsin-fold in the catalytic domain, is conferred by the large differences of eight surface loops, surrounding the active site. The length of these loops alters with insertions and deletions, resulting in remarkably different three-dimensional arrangements. In addition, metal binding sites for Na+, Ca2+ and Zn2+ serve as regulatory elements, as do N-glycosylation sites. Depending on the individual tasks of the protease, the surface loops determine substrate specificity, control the turnover and allow regulation of activation, activity and degradation by other proteins, which are often serine proteases themselves. Most intriguingly, in some serine proteases, the surface loops interact as allosteric network, partially tuned by protein co-factors. Knowledge of these subtle and complicated molecular motions may allow nowadays for new and specific pharmaceutical or medical approaches.

Residues W215, E217 and E192 control the allosteric E*-E equilibrium of thrombin

A pre-existing, allosteric equilibrium between closed (E*) and open (E) conformations of the active site influences the level of activity in the trypsin fold and defines ligand binding according to the mechanism of conformational selection. Using the clotting protease thrombin as a model system, we investigate the molecular determinants of the E*-E equilibrium through rapid kinetics and X-ray structural biology. The equilibrium is controlled by three residues positioned around the active site. W215 on the 215-217 segment defining the west wall of the active site controls the rate of transition from E to E* through hydrophobic interaction with F227. E192 on the opposite 190-193 segment defining the east wall of the active site controls the rate of transition from E* to E through electrostatic repulsion of E217. The side chain of E217 acts as a lever that moves the entire 215-217 segment in the E*-E equilibrium. Removal of this side chain converts binding to the active site to a simple lock-and-key mechanism and freezes the conformation in a state intermediate between E* and E. These findings reveal a simple framework to understand the molecular basis of a key allosteric property of the trypsin fold.

Live-Cell Imaging of Physiologically Relevant Metal Ions Using Genetically Encoded FRET-Based Probes

Essential biochemical reactions and processes within living organisms are coupled to subcellular fluctuations of metal ions. Disturbances in cellular metal ion homeostasis are frequently associated with pathological alterations, including neurotoxicity causing neurodegeneration, as well as metabolic disorders or cancer. Considering these important aspects of the cellular metal ion homeostasis in health and disease, measurements of subcellular ion signals are of broad scientific interest. The investigation of the cellular ion homeostasis using classical biochemical methods is quite difficult, often even not feasible or requires large cell numbers. Here, we report of genetically encoded fluorescent probes that enable the visualization of metal ion dynamics within individual living cells and their organelles with high temporal and spatial resolution. Generally, these probes consist of specific ion binding domains fused to fluorescent protein(s), altering their fluorescent properties upon ion binding. This review focuses on the functionality and potential of these genetically encoded fluorescent tools which enable monitoring (sub)cellular concentrations of alkali metals such as K+, alkaline earth metals including Mg2+ and Ca2+, and transition metals including Cu+/Cu2+ and Zn2+. Moreover, we discuss possible approaches for the development and application of novel metal ion biosensors for Fe2+/Fe3+, Mn2+ and Na+.

Sodium-site in serine protease domain of human coagulation factor IXa: evidence from the crystal structure and molecular dynamics simulations study

Essentials Consensus sequence and biochemical data suggest a Na+ -site in the factor (F) IXa protease domain. X-ray structure of the FIXa EGF2/protease domain at 1.37 Å reveals a Na+ -site not observed earlier. Molecular dynamics simulations data support that Na+  ± Ca2+ promote FIXa protease domain stability. Sulfate ions found in the protease domain mimic heparin sulfate binding mode in FIXa. SUMMARY: Background Activated coagulation factor IX (FIXa) consists of a γ-carboxyglutamic acid domain, two epidermal growth factor-like (EGF) domains, and a C-terminal protease domain. Consensus sequence and biochemical data support the existence of a Na+ -site in the FIXa protease domain. However, soaking experiments or crystals grown in high concentration of ammonium sulfate did not reveal a Na+ -site in wild-type or mutant FIXa EGF2/protease domain structure. Objective Determine the structure of the FIXa EGF2/protease domain in the presence of Na+ ; perform molecular dynamics (MD) simulations to explore the role of Na+ in stabilizing FIXa structure. Methods Crystallography, MD simulations, and modeling heparin binding to FIXa. Results Crystal structure at 1.37-Å resolution revealed that Na+ is coordinated to carbonyl groups of residues 184A, 185, 221A, and 224 in the FIXa protease domain. The Na+ -site in FIXa is similar to that of FXa and is linked to the Asp189 S1-site. In MD simulations, Na+ reduced fluctuations in residues 217-225 (Na+ -loop) and 70-80 (Ca2+ -loop), whereas Ca2+ reduced fluctuations only in residues of the Ca2+ -loop. Ca2+ and Na+ together reduced fluctuations in residues of the Ca2+ -loop and Na+ -loop (residues 70-80, 183-194, and 217-225). Moreover, we observed four sulfate ions that make salt bridges with FIXa protease domain Arg/Lys residues, which have been implicated in heparin binding. Based upon locations of the sulfate ions, we modeled heparin binding to FIXa, which is similar to the heparin binding in thrombin. Conclusions The FIXa Na+ -site in association with Ca2+ contributes to stabilization of the FIXa protease domain. The heparin binding mode in FIXa is similar to that in thrombin.

Na+-binding modes involved in thrombin's allosteric response as revealed by molecular dynamics simulations, correlation networks and Markov modeling

The monovalent sodium ion (Na+) is a critical modulator of thrombin. However, the mechanism of thrombin's activation by Na+ has been widely debated for more than twenty years. Details of the linkage between thrombin and Na+ remain vague due to limited temporal and spatial resolution in experiments. In this work, we combine microsecond scale atomic-detailed molecular dynamics simulations with correlation network analyses and hidden Markov modeling to probe the detailed thermodynamic and kinetic picture of Na+-binding events and their resulting allosteric responses in thrombin. We reveal that ASP189 and ALA190 comprise a stable Na+-binding site (referred as "inner" Na+-binding site) along with the previously known one (referred as "outer" Na+-binding site). The corresponding newly identified Na+-binding mode introduces significant allosteric responses in thrombin's regulatory regions by stabilizing selected torsion angles of residues responsive to Na+-binding. Our Markov model indicates that the bound Na+ prefers to transfer between the two Na+-binding sites when an unbinding event takes place. These results suggest a testable hypothesis of a substrate-driven Na+ migration (ΔG ∼ 1.7 kcal mol-1) from the "inner" Na+-binding site to the "outer" one during thrombin's catalytic activities. The binding of a Na+ ion at the "inner" Na+-binding site should be inferred as a prerequisite for thrombin's efficient recognition to the substrate, which opens a new angle for our understanding of Na+-binding's allosteric activation on thrombin and sheds light on detailed processes in thrombin's activation.

Dynamic Consequences of Mutation of Tryptophan 215 in Thrombin

Thrombin normally cleaves fibrinogen to promote coagulation; however, binding of thrombomodulin to thrombin switches the specificity of thrombin toward protein C, triggering the anticoagulation pathway. The W215A thrombin mutant was reported to have decreased activity toward fibrinogen without significant loss of activity toward protein C. To understand how mutation of Trp215 may alter thrombin specificity, hydrogen-deuterium exchange experiments (HDXMS), accelerated molecular dynamics (AMD) simulations, and activity assays were carried out to compare the dynamics of Trp215 mutants with those of wild type (WT) thrombin. Variation in NaCl concentration had no detectable effect on the sodium-binding (220s_CT) loop, but appeared to affect other surface loops. Trp215 mutants showed significant increases in amide exchange in the 170s_CT loop consistent with a loss of H-bonding in this loop identified by the AMD simulations. The W215A thrombin showed increased amide exchange in the 220s_CT loop and in the N-terminus of the heavy chain. The AMD simulations showed that a transient conformation of the W215A thrombin has a distorted catalytic triad. HDXMS experiments revealed that mutation of Phe227, which engages in a π-stacking interaction with Trp215, also caused significantly increased amide exchange in the 170s_CT loop. Activity assays showed that only the F227V mutant had wild type catalytic activity, whereas all other mutants showed markedly lower activity. Taken together, the results explain the reduced pro-coagulant activity of the W215A mutant and demonstrate the allosteric connection between Trp215, the sodium-binding loop, and the active site.

Missense mutations in the gene encoding prothrombin corresponding to Arg596 cause antithrombin resistance and thrombomodulin resistance

Antithrombin (AT) and thrombomodulin (TM) play important roles in the process of natural anticoagulation in vivo. Recently, we reported that the prothrombin Yukuhashi mutation (p.Arg596Leu) was associated with AT and TM resistance-related thrombophilia. To assess the AT and TM resistances associated with other missense mutations by single base substitution in the Arg596 codon, we generated recombinant variants (596Gln, 596Trp, 596Gly, and 596Pro) and investigated the effects on AT and TM anticoagulant functions. All variants except 596Pro were secreted in amounts comparable to that of the wild-type but exhibited variable procoagulant activities. After a 30-minute inactivation by AT, the relative residual activity of wild-type thrombin decreased to 15 ± 4.0%, in contrast to values of all variants were maintained at above 80%. The thrombin–AT complex formation, as determined by enzyme-linked immunosorbent assay, was reduced with all tested variants in the presence and absence of heparin. In the presence of soluble TM (sTM), the relative fibrinogen clotting activity of wild-type thrombin decreased to 16 ± 0.12%, whereas that of tested variants was 37%–56%. In a surface plasmon resonance assay, missense Arg596 mutations reduced thrombin–TM affinity to an extent similar to the reduction of fibrinogen clotting inhibition. In the presence of sTM or cultured endothelial-like cells, APC generation was enhanced differently by variant thrombins in a thrombin–TM affinity- dependent manner. These data indicate that prothrombin Arg596 missense mutations lead to AT and TM resistance in the variant thrombins and suggest that prothrombin Arg596 is important for AT- and TM- mediated anticoagulation.

Interplay between conformational selection and zymogen activation

Trypsin-like proteases are synthesized as zymogens and activated through a mechanism that folds the active site for efficient binding and catalysis. Ligand binding to the active site is therefore a valuable source of information on the changes that accompany zymogen activation. Using the physiologically relevant transition of the clotting zymogen prothrombin to the mature protease thrombin, we show that the mechanism of ligand recognition follows selection within a pre-existing ensemble of conformations with the active site accessible (E) or inaccessible (E*) to binding. Prothrombin exists mainly in the E* conformational ensemble and conversion to thrombin produces two dominant changes: a progressive shift toward the E conformational ensemble triggered by removal of the auxiliary domains upon cleavage at R271 and a drastic drop of the rate of ligand dissociation from the active site triggered by cleavage at R320. Together, these effects produce a significant (700-fold) increase in binding affinity. Limited proteolysis reveals how the E*-E equilibrium shifts during prothrombin activation and influences exposure of the sites of cleavage at R271 and R320. These new findings on the molecular underpinnings of prothrombin activation are relevant to other zymogens with modular assembly involved in blood coagulation, complement and fibrinolysis.

Structure of prothrombin in the closed form reveals new details on the mechanism of activation

The clotting factor prothrombin exists in equilibrium between closed and open conformations, but the physiological role of these forms remains unclear. As for other allosteric proteins, elucidation of the linkage between molecular transitions and function is facilitated by reagents stabilized in each of the alternative conformations. The open form of prothrombin has been characterized structurally, but little is known about the architecture of the closed form that predominates in solution under physiological conditions. Using X-ray crystallography and single-molecule FRET, we characterize a prothrombin construct locked in the closed conformation through an engineered disulfide bond. The construct: (i) provides structural validation of the intramolecular collapse of kringle-1 onto the protease domain reported recently; (ii) documents the critical role of the linker connecting kringle-1 to kringle-2 in stabilizing the closed form; and (iii) reveals novel mechanisms to shift the equilibrium toward the open conformation. Together with functional studies, our findings define the role of closed and open conformations in the conversion of prothrombin to thrombin and establish a molecular framework for prothrombin activation that rationalizes existing phenotypes associated with prothrombin mutations and points to new strategies for therapeutic intervention.

Mechanistic insights into thrombin's switch between "slow" and "fast" forms

Thrombin is a multifunctional enzyme that plays an important role in blood coagulation, cell growth, and metastasis. Depending upon the binding of sodium ions, thrombin presents significantly different enzymatic activities. In the environment with sodium ions, thrombin is highly active in cleaving the coagulated substrates and this is referred to as the "fast" form; in the environment without sodium ions, thrombin turns catalytically less active and is in the "slow" form. Although many experimental studies over the last two decades have attempted to reveal the structural and kinetic differences between these two forms, it remains vague and disputed how the functional switch between the "fast" and "slow" forms is mediated by Na+ cations. In this work, we employ microsecond-scale all-atom molecular dynamics simulations to investigate the differences in the structural ensembles in sodium-bound/unbound and potassium-bound/unbound thrombin. Our calculations indicate that the regulatory regions, including the 60s, γ loops, and exosite I and II, are primarily affected by both the bound and unbound cations. Conformational free energy surfaces, estimated from principal component analysis, further reveal the existence of multiple conformational states. The binding of a cation introduces changes in the distribution of these states. Through comparisons with potassium-binding, the binding of sodium ions appears to shift the population toward conformational states that might be catalytically favorable. Our study of thrombin in the presence of sodium/potassium ions suggests Na+-mediated generalized allostery is the mechanism of thrombin's functional switch between the "fast" and "slow" forms.

Induced Fit Is a Special Case of Conformational Selection

Conformational selection (CS) and induced fit (IF) are two widely used interpretations of binding of a ligand to biological macromolecules. Both mechanisms envision a two-step reaction in which a conformational transition either precedes (CS) or follows (IF) the binding step. Under pseudo-first-order conditions where the ligand is in excess compared to the macromolecule, both mechanisms produce two relaxations. A fast one eventually increases linearly with ligand concentration and reflects the binding interaction. A slow one saturates to a constant value after decreasing or increasing hyperbolically with ligand concentration. This relaxation is the one most often accessible to experimental measurements and is potentially diagnostic of the mechanism involved. A relaxation that decreases unequivocally identifies CS, but a hyperbolic increase is compatible with both CS and IF. The potential ambiguity between the two mechanisms is more than qualitative. Here we show that the entire kinetic repertoire of IF is nothing but a mathematical special case of CS as revealed by a simple transformation of the rate constants, which emphasizes the need for independent support of either mechanism from additional experimental evidence. We discuss a simple strategy for distinguishing between IF and CS under the most common conditions encountered in practice, i.e., when the ligand is in excess compared to the macromolecule and a single relaxation is accessible to experimental measurements.

Role of glycine 221 in catalytic activity of hyaluronan-binding protein 2

HABP2 (hyaluronan-binding protein 2) is a Ca²⁺-dependent serine protease with putative roles in blood coagulation and fibrinolysis. A G221E substitution, known as the Marburg I polymorphism, reportedly affects HABP2 function and has been associated with increased risk for cardiovascular disease. However, the importance of Gly-221 for HABP2 activity is unclear. Here, we used G221E, G221A, and G221S mutants to assess the role of Gly-221 in HABP2 catalysis. The G221E variant failed to activate the single-chain urokinase-type plasminogen activator, and the G221A and G221S variants displayed moderately reduced single-chain urokinase-type plasminogen activator activation. Activity toward the peptide substrate S-2288 was markedly decreased in all HABP2 variants, with G221E being the most defective and G221A being the least defective. In the absence of Ca²⁺, S-2288 cleavage by wild-type HABP2 was Na⁺-dependent, with K_m decreasing from 3.0 to 0.6 mm upon titration from 0 to 0.3 m Na⁺ In the presence of 5 mm Ca²⁺, K_m was further reduced to 0.05 mm, but without an appreciable contribution of Na⁺ At physiological concentrations of Na⁺ and Ca²⁺, the three HABP2 variants, and particularly G221E, displayed a major K_m increase for S-2288. Chemical footprinting revealed that Ile-16 is significantly less protected from chemical modification in G221E than in wild-type HABP2, suggesting impaired insertion of the N terminus into the G221E protease domain, with a concomitant impact on catalytic activity. Homology modeling suggested that the Glu-221 side chain could sterically hinder insertion of the N terminus into the HABP2 protease domain, helping to explain the detrimental effects of Glu-221 substitution on HABP2 activity.

Clinical and biochemical characterization of the prothrombin Belgrade mutation in a large Serbian pedigree: new insights into the antithrombin resistance mechanism

Essentials Prothrombin Belgrade mutation leads to antithrombin resistance. Clinical and biochemical phenotypes in a large family with this mutation were investigated. In carriers, we detected decreased factor II activity and increased endogenous thrombin potential. Prothrombin Belgrade mutation represents a strong prothrombotic risk factor.

SUMMARY

Background The recently reported c.1787G>A mutation in the prothrombin gene leads to Arg596Gln replacement in the protein molecule (prothrombin Belgrade). This substitution impairs binding of antithrombin to thrombin and results in inherited thrombophilia, known as antithrombin resistance. Objectives We aimed to elucidate the clinical and biochemical characteristics of thrombophilia associated with antithrombin resistance in a large Serbian family with the prothrombin Belgrade mutation. Patients and methods Nineteen family members were investigated, among whom 10 were carriers of the c.1787G>A mutation. In all subjects the clinical phenotype was determined and laboratory investigations of hemostatic parameters were performed. Results Six out of the 10 mutation carriers developed thromboembolic events, mainly deep venous and mesenteric vein thrombosis. The median age of the first thrombotic event was 26.5 (12-41) years, whereas the incidence rate of first thrombosis was 2.2% per year. In all mutation carriers prothrombin activity was significantly decreased in comparison with non-carriers, clearly distinguishing each group. However, the presence of the mutation did not affect the prothrombin antigen level in plasma. The endogenous thrombin potential was significantly increased in all carriers in comparison with non-carriers, indicating the presence of blood hypercoagulability. Interestingly, levels of D-dimer and the F1+2 fragment were similar in both groups. Conclusions Although rare, the prothrombin Belgrade mutation represents strong thrombophilia with early onset of thrombosis in the investigated family. According to our results, decreased prothrombin activity may be a simple screening test for detection of this mutation in thrombotic patients.

Na+ Binding Is Ineffective in Forming a Primary Substrate Pocket of Thrombin

Thrombin is a serine protease involved in the blood coagulation reaction, and it shows maximum enzymatic activity in the presence of Na⁺. It has been supposed that Na⁺ binding promotes conversion from the inactive form, with a collapsed primary substrate pocket (S1 pocket), to the active form, with a properly formed S1 pocket. However, the evidence supporting this activation mechanism was derived from the X-ray crystallographic structures solved under nonphysiological conditions and using thrombin mutants; thus, it still remains elusive whether the activation mechanism is actually attributed to Na⁺ binding. To address the problem, we employed all-atom molecular dynamics simulations for both active and inactive forms of thrombin in the presence and absence of Na⁺ binding and examined the effect of Na⁺ binding on S1-pocket formation. In contrast to the conventional supposition, we revealed that Na⁺ binding does not prevent S1-pocket collapse virtually, but rather, the bound Na⁺ can move to the S1 pocket, thus blocking substrate access directly. Additionally, it was clarified that Na⁺ binding does not promote S1-pocket formation. According to these insights, we concluded that Na⁺ binding is irrelevant to the interconversion between the inactive and active forms of thrombin.

Molecular Mechanisms of Enzyme Activation by Monovalent Cations

Regulation of enzymes through metal ion complexation is widespread in biology and underscores a physiological need for stability and high catalytic activity that likely predated proteins in the RNA world. In addition to divalent metals such as Ca²⁺, Mg²⁺, and Zn²⁺, monovalent cations often function as efficient and selective promoters of catalysis. Advances in structural biology unravel a rich repertoire of molecular mechanisms for enzyme activation by Na⁺ and K⁺ Strategies range from short-range effects mediated by direct participation in substrate binding, to more distributed effects that propagate long-range to catalytic residues. This review addresses general considerations and examples.

Elucidating factors important for monovalent cation selectivity in enzymes: E. coli β-galactosidase as a model

Many enzymes require a specific monovalent cation (M(+)), that is either Na(+) or K(+), for optimal activity. While high selectivity M(+) sites in transport proteins have been extensively studied, enzyme M(+) binding sites generally have lower selectivity and are less characterized. Here we study the M(+) binding site of the model enzyme E. coli β-galactosidase, which is about 10 fold selective for Na(+) over K(+). Combining data from X-ray crystallography and computational models, we find the electrostatic environment predominates in defining the Na(+) selectivity. In this lower selectivity site rather subtle influences on the electrostatic environment become significant, including the induced polarization effects of the M(+) on the coordinating ligands and the effect of second coordination shell residues on the charge distribution of the primary ligands. This work expands the knowledge of ion selectivity in proteins to denote novel mechanisms important for the selectivity of M(+) sites in enzymes.

Molecular Basis of Enhanced Activity in Factor VIIa-Trypsin Variants Conveys Insights into Tissue Factor-mediated Allosteric Regulation of Factor VIIa Activity

The complex of coagulation factor VIIa (FVIIa), a trypsin-like serine protease, and membrane-bound tissue factor (TF) initiates blood coagulation upon vascular injury. Binding of TF to FVIIa promotes allosteric conformational changes in the FVIIa protease domain and improves its catalytic properties. Extensive studies have revealed two putative pathways for this allosteric communication. Here we provide further details of this allosteric communication by investigating FVIIa loop swap variants containing the 170 loop of trypsin that display TF-independent enhanced activity. Using x-ray crystallography, we show that the introduced 170 loop from trypsin directly interacts with the FVIIa active site, stabilizing segment 215-217 and activation loop 3, leading to enhanced activity. Molecular dynamics simulations and novel fluorescence quenching studies support that segment 215-217 conformation is pivotal to the enhanced activity of the FVIIa variants. We speculate that the allosteric regulation of FVIIa activity by TF binding follows a similar path in conjunction with protease domain N terminus insertion, suggesting a more complete molecular basis of TF-mediated allosteric enhancement of FVIIa activity.

Kinetic dissection of the pre-existing conformational equilibrium in the trypsin fold

Structural biology has recently documented the conformational plasticity of the trypsin fold for both the protease and zymogen in terms of a pre-existing equilibrium between closed (E*) and open (E) forms of the active site region. How such plasticity is manifested in solution and affects ligand recognition by the protease and zymogen is poorly understood in quantitative terms. Here we dissect the E*-E equilibrium with stopped-flow kinetics in the presence of excess ligand or macromolecule. Using the clotting protease thrombin and its zymogen precursor prethrombin-2 as relevant models we resolve the relative distribution of the E* and E forms and the underlying kinetic rates for their interconversion. In the case of thrombin, the E* and E forms are distributed in a 1:4 ratio and interconvert on a time scale of 45 ms. In the case of prethrombin-2, the equilibrium is shifted strongly (10:1 ratio) in favor of the closed E* form and unfolds over a faster time scale of 4.5 ms. The distribution of E* and E forms observed for thrombin and prethrombin-2 indicates that zymogen activation is linked to a significant shift in the pre-existing equilibrium between closed and open conformations that facilitates ligand binding to the active site. These findings broaden our mechanistic understanding of how conformational transitions control ligand recognition by thrombin and its zymogen precursor prethrombin-2 and have direct relevance to other members of the trypsin fold.

Do the crystallographic forms of prethrombin-2 revert to a single form in solution?

It has been earlier established (Pozzi et al. Biochemistry 50 (2011) 10195-10202) that prethrombin-2 crystallizes into two similar but distinct forms: a collapsed form and an alternative form. We employed long molecular dynamics (MD) simulations for these two forms to obtain solvent-equilibrated forms. We find that, at 200ns, the simulated solution collapsed form is quite similar to the X-ray crystal collapsed form, while the simulated solution alternative form deviates from the X-ray crystal alternative form as well as from the solution collapsed form. A detailed structural analysis suggests that the fluctuation of the 140s-loop, in cross-talk with the 220s-loop, may alter the conformation of the W215-E217 segment near the nascent thrombin active site. A rationale is provided for the manner in which interactions of prethrombin-2 with FVa may affect the equilibrium between the two forms of prethrombin-2.

Toward understanding allosteric activation of thrombin: a conjecture for important roles of unbound Na(+) molecules around thrombin

We shed light on important roles of unbound Na(+) molecules in enzymatic activation of thrombin. Molecular mechanism of Na(+)-activation of thrombin has been discussed in the context of allostery. However, the recent challenge to redesign K(+)-activated thrombin revealed that the allosteric interaction is insufficient to explain the mechanism. Under these circumstances, we have examined the roles of unbound Na(+) molecule in maximization of thrombin-substrate association reaction rate. We performed all-atomic molecular dynamics (MD) simulations of thrombin in the presence of three different cations; Li(+), Na(+), and Cs(+). Although these cations are commonly observed in the vicinity of the S1-pocket of thrombin, smaller cations are distributed more densely and extensively than larger ones. This suggests the two observation rules: (i) thrombin surrounded by Na(+) is at an advantage in the initial step of association reaction, namely, the formation of an encounter complex ensemble, and (ii) the presence of Na(+) molecules does not necessarily have an advantage in the final step of association reaction, namely, the formation of the stereospecific complex. In conclusion, we propose a conjecture that unbound Na(+) molecules also affect the maximization of rate constant of thrombin-substrate association reaction through optimally forming an encounter complex ensemble.