Rigidification of the autolysis loop enhances Na(+) binding to thrombin

Binding of Na⁺ to thrombin ensures high activity toward physiological substrates and optimizes the procoagulant and prothrombotic roles of the enzyme in vivo. Under physiological conditions of pH and temperature, the binding affinity of Na⁺ is weak due to large heat capacity and enthalpy changes associated with binding, and the K_d = 80 mM ensures only 64% saturation of the site at the concentration of Na⁺ in the blood (140 mM). Residues controlling Na⁺ binding and activation have been identified. Yet, attempts to improve the interaction of Na⁺ with thrombin and possibly increase catalytic activity under physiological conditions have so far been unsuccessful. Here we report how replacement of the flexible autolysis loop of human thrombin with the homologous rigid domain of the murine enzyme results in a drastic (up to 10-fold) increase in Na⁺ affinity and a significant improvement in the catalytic activity of the enzyme. Rigidification of the autolysis loop abolishes the heat capacity change associated with Na⁺ binding observed in the wild-type and also increases the stability of thrombin. These findings have general relevance to protein engineering studies of clotting proteases and trypsin-like enzymes.     

Probing the hirudin-thrombin interaction by incorporation of noncoded amino acids and molecular dynamics simulation

Thrombin is a primary target for the development of novel anticoagulants, since it plays two important and opposite roles in hemostasis: procoagulant and anticoagulant. All thrombin functions are influenced by Na+ binding, which triggers the transition of this enzyme from an anticoagulant (slow) form to a procoagulant (fast) form. In previous studies, we have conveniently produced by chemical synthesis analogues of the N-terminal fragment 1-47 of hirudin HM2 containing noncoded amino acids and displaying up to approximately 2700-fold more potent antithrombin activity, comparable to that of full-length hirudin. In the work presented here, we have exploited the versatility of chemical synthesis to probe the structural and energetic properties of the S3 site of thrombin through perturbations introduced in the structure of hirudin fragment 1-47. In particular, we have investigated the effects of systematic replacement of Tyr3 with noncoded amino acids retaining the aromatic nucleus of Tyr, as well as similar hydrophobic and steric properties, but possessing different electronic (e.g., p-fluoro-, p-iodo-, or p-nitro-Phe), charge (p-aminomethyl-Phe), or conformational (homo-Phe) properties. Our results indicate that the affinity of fragment 1-47 for thrombin is proportional to the desolvation free energy change upon complex formation, and is inversely related to the electric dipole moment of the amino acid side chain at position 3 of hirudin. In this study, we have also identified the key features that are responsible for the preferential binding of hirudin to the procoagulant (fast) form of thrombin. Strikingly, shaving at position 3, by Tyr --> Ala exchange, abolishes the differences in the affinity for thrombin allosteric forms, whereas a bulkier side chain (e.g., beta-naphthylalanine) improves binding preferentially to the fast form. These results provide strong, albeit indirect, evidence that the procoagulant (fast) form of thrombin is in a more open and accessible conformation with respect to the less forgiving structure it acquires in the slow form. This view is also supported by the results of molecular dynamics simulations conducted for 18 ns on free thrombin in full explicit water, showing that after approximately 5 ns thrombin undergoes a significant conformational transition, from a more open conformation (which we propose can be related to the fast form) to a more compact and closed one (which we propose can be related to the slow form). This transition mainly involves the Trp148 and Trp60D loop, the S3 site, and the fibrinogen binding site, whereas the S1 site, the Na+-binding site, and the catalytic pocket remain essentially unchanged. In particular, our data indicate that the S3 site of the enzyme is less accessible to water in the putative slow form. This structural picture provides a reasonable molecular explanation for the fact that physiological substrates related to the procoagulant activity of thrombin (fibrinogen, thrombin receptor 1, and factor XIII) orient a bulky side chain into the S3 site of the enzyme. Taken together, our results can have important implications for the design of novel thrombin inhibitors, of practical utility in the treatment of coagulative disorders.     

Identification of a 3rd Na+ Binding Site of the Glycine Transporter,GlyT2

The Na⁺/Cl^- dependent glycine transporters GlyT1 and GlyT2 regulate synaptic glycine concentrations. Glycine transport by GlyT2 is coupled to the co-transport of three Na⁺ ions, whereas transport by GlyT1 is coupled to the co-transport of only two Na⁺ ions. These differences in ion-flux coupling determine their respective concentrating capacities and have a direct bearing on their functional roles in synaptic transmission. The crystal structures of the closely related bacterial Na⁺-dependent leucine transporter, LeuT_Aa, and the Drosophila dopamine transporter, dDAT, have allowed prediction of two Na⁺ binding sites in GlyT2, but the physical location of the third Na⁺ site in GlyT2 is unknown. A bacterial betaine transporter, BetP, has also been crystallized and shows structural similarity to LeuT_Aa. Although betaine transport by BetP is coupled to the co-transport of two Na⁺ ions, the first Na⁺ site is not conserved between BetP and LeuT_Aa, the so called Na1'' site. We hypothesized that the third Na⁺ binding site (Na3 site) of GlyT2 corresponds to the BetP Na1'' binding site. To identify the Na3 binding site of GlyT2, we performed molecular dynamics (MD) simulations. Surprisingly, a Na⁺ placed at the location consistent with the Na1'' site of BetP spontaneously dissociated from its initial location and bound instead to a novel Na3 site. Using a combination of MD simulations of a comparative model of GlyT2 together with an analysis of the functional properties of wild type and mutant GlyTs we have identified an electrostatically favorable novel third Na⁺ binding site in GlyT2 formed by Trp263 and Met276 in TM3, Ala481 in TM6 and Glu648 in TM10.     

Crystal structure of the thrombin mutant D221A/D222K: the Asp222:Arg187 ion-pair stabilizes the fast form

The thrombin mutant D221A/D222K (ARK) does not bind Na+ and has interesting functional properties intermediate between those of the slow and fast forms of wild type. We solved the X-ray crystal structure of ARK bound at exosite I with a fragment of hirudin at 2.4-A resolution. The structure shows a slight collapse of the 186 and 220 loops into the Na+ binding site due to disruption of the Asp222:Arg187 ion-pair. The backbone O atoms of Arg221a and Lys224 are shifted into conformations that eliminate optimal interaction with Na+. A paucity of solvent molecules in the Na+ binding site is also noted, by analogy to what is seen in the structure of the slow form. These findings reinforce the crucial role of the Asp222:Arg187 ion-pair in stabilizing the fast form of thrombin.     

Free energy simulations of ligand binding to the aspartate transporter Glt(Ph)

Glutamate/Aspartate transporters cotransport three Na⁺ and one H⁺ ions with the substrate and countertransport one K⁺ ion. The binding sites for the substrate and two Na⁺ ions have been observed in the crystal structure of the archeal homolog Glt_Ph, while the binding site for the third Na⁺ ion has been proposed from computational studies and confirmed by experiments. Here we perform detailed free energy simulations of Glt_Ph, giving a comprehensive characterization of the substrate and ion binding sites, and calculating their binding free energies in various configurations. Our results show unequivocally that the substrate binds after the binding of two Na⁺ ions. They also shed light into Asp/Glu selectivity of Glt_Ph, which is not observed in eukaryotic glutamate transporters.     

Substrate binding and formation of an occluded state in the leucine transporter

Translocation through the extracellular vestibule and binding of leucine in the leucine transporter (LeuT) have been studied with molecular dynamics simulations. More than 0.1 μs of all-atom molecular dynamics simulations have been performed on different combinations of LeuT, bound substrate, and bound structural Na⁺ ions to describe molecular events involved in substrate binding and in the formation of the occluded state and to investigate the dynamics of this state. Three structural features are found to be directly involved in the initial steps of leucine transport: a Na⁺ ion directly coordinated to leucine (Na-1), two aromatic residues closing the binding site toward the extracellular vestibule (Tyr-108 and Phe-253), and a salt bridge in the extracellular vestibule (Arg-30 and Asp-404). These features account for observed differences between simulations of LeuT with and without bound substrate and for a possible pathway for leucine binding and thereby formation of the occluded LeuT binding site.     

Conformational changes in the selectivity filter of the open-state KcsA channel: an energy minimization study

Potassium channels switch between closed and open conformations and selectively conduct K⁺ ions. There are at least two gates. The TM2 bundle at the intracellular site is the primary gate of KcsA, and rearrangements at the selectivity filter (SF) act as the second gate. The SF blocks ion flow via an inactivation process similar to C-type inactivation of voltage-gated K⁺ channels. We recently generated the open-state conformation of the KcsA channel. We found no major, possibly inactivating, structural changes in the SF associated with this massive inner-pore rearrangement, which suggests that the gates might act independently. Here we energy-minimize the open state of wild-type and mutant KcsA, validating in silico structures of energy-minimized SFs by comparison with crystallographic structures, and use these data to gain insight into how mutation, ion depletion, and K⁺ to Na⁺ substitution influence SF conformation. Both E71 or D80 protonations/mutations and the presence/absence of protein-buried water molecule(s) modify the H-bonding network stabilizing the P-loops, spawning numerous SF conformations. We find that the inactivated state corresponds to conformations with a partially unoccupied or an entirely empty SF. These structures, involving modifications in all four P-loops, are stabilized by H-bonds between amide H and carbonyl O atoms from adjacent P-loops, which block ion passage. The inner portions of the P-loops are more rigid than the outer parts. Changes are localized to the outer binding sites, with innermost site S4 persisting in the inactivated state. Strong binding by Na⁺ locally contracts the SF around Na⁺, releasing ligands that do not participate in Na⁺ coordination, and occluding the permeation pathway. K⁺ selectivity primarily appears to arise from the inability of the SF to completely dehydrate Na⁺ ions due to basic structural differences between liquid water and the “quasi-liquid” SF matrix.     

Effects of the Protonation State of the EEEE Motif of a Bacterial Na+-channel on Conduction and Pore Structure

A distinctive feature of prokaryotic Na⁺-channels is the presence of four glutamate residues in their selectivity filter. In this study, how the structure of the selectivity filter, and the free-energy profile of permeating Na⁺ ions are altered by the protonation state of Glu177 are analyzed. It was found that protonation of a single glutamate residue was enough to modify the conformation of the selectivity filter and its conduction properties. Molecular dynamics simulations revealed that Glu177 residues may adopt two conformations, with the side chain directed toward the extracellular entrance of the channel or the intracellular cavity. The likelihood of the inwardly directed arrangement increases when Glu177 residues are protonated. The presence of one glutamate residue with its chain directed toward the intracellular cavity increases the energy barrier for translocation of Na⁺ ions. These higher-energy barriers preclude Na⁺ ions to permeate the selectivity filter of prokaryotic Na⁺-channels when one or more Glu177 residues are protonated.     

Molecular insights on CALX-CBD12 interdomain dynamics from MD simulations,RDCs, and SAXS

Na⁺/Ca²⁺ exchangers (NCXs) are secondary active transporters that couple the translocation of Na⁺ with the transport of Ca²⁺ in the opposite direction. The exchanger is an essential Ca²⁺ extrusion mechanism in excitable cells. It consists of a transmembrane domain and a large intracellular loop that contains two Ca²⁺-binding domains, CBD1 and CBD2. The two CBDs are adjacent to each other and form a two-domain Ca²⁺ sensor called CBD12. Binding of intracellular Ca²⁺ to CBD12 activates the NCX but inhibits the NCX of Drosophila, CALX. NMR spectroscopy and SAXS studies showed that CALX and NCX CBD12 constructs display significant interdomain flexibility in the apo state but assume rigid interdomain arrangements in the Ca²⁺-bound state. However, detailed structure information on CBD12 in the apo state is missing. Structural characterization of proteins formed by two or more domains connected by flexible linkers is notoriously challenging and requires the combination of orthogonal information from multiple sources. As an attempt to characterize the conformational ensemble of CALX-CBD12 in the apo state, we applied molecular dynamics (MD) simulations, NMR (¹H-¹⁵N residual dipolar couplings), and small-angle x-ray scattering (SAXS) data in a combined strategy to select an ensemble of conformations in agreement with the experimental data. This joint approach demonstrated that CALX-CBD12 preferentially samples closed conformations, whereas the wide-open interdomain arrangement characteristic of the Ca²⁺-bound state is less frequently sampled. These results are consistent with the view that Ca²⁺ binding shifts the CBD12 conformational ensemble toward extended conformers, which could be a key step in the NCXs’ allosteric regulation mechanism. This strategy, combining MD with NMR and SAXS, provides a powerful approach to select ensembles of conformations that could be applied to other flexible multidomain systems.     

Coenzyme A binding to the aminoglycoside acetyltransferase (3)-IIIb increases conformational sampling of antibiotic binding site

NMR spectroscopy experiments and molecular dynamics simulations were performed to describe the dynamic properties of the aminoglycoside acetyltransferase (3)-IIIb (AAC) in its apo and coenzyme A (CoASH) bound forms. The (15)N-(1)H HSQC spectra indicate a partial structural change and coupling of the CoASH binding site with another region in the protein upon the CoASH titration into the apo enzyme. Molecular dynamics simulations indicate a significant structural and dynamic variation of the long loop in the antibiotic binding domain in the form of a relatively slow (250 ns), concerted opening motion in the CoASH-enzyme complex and that binding of the CoASH increases the structural flexibility of the loop, leading to an interchange between several similar equally populated conformations.     

Interactions of Local Anesthetics with Voltage-gated Na<Superscript>+</Superscript> Channels

Voltage-gated Na⁺ channels are dynamic transmembrane proteins responsible for the rising phase of the action potential in excitable membranes. Local anesthetics (LAs) and structurally related antiarrhythmic and anticonvulsant compounds target specific sites in voltage-gated Na⁺ channels to block Na⁺ currents, thus reducing excitability in neuronal, cardiac, or central nervous tissue. A high-affinity LA block is produced by binding to open and inactivated states of Na⁺ channels rather than to resting states and suggests a binding site that converts from a low- to a high-affinity conformation during gating. Recent findings using site-directed mutagenesis suggest that multiple S6 segments together form an LA binding site within the Na⁺ channel. While the selectivity filter may form the more extracellular-located part of this binding site, the role of the fast inactivation gate in LA binding has not yet been resolved. The receptor of the neurotoxin batrachotoxin (BTX) is adjacent to or even overlaps with the LA binding site. The close proximity of the LA and BTX binding sites to residues critical for inactivation, together with gating transitions through S6 segments, might explain the strong impact of LAs and BTX on inactivation of voltage-gated Na⁺ channels and might help elucidate the mechanisms underlying voltage- and frequency-dependent LA block.     

Interaction of phospholipids (Lysophosphatidylethanolamines) with water and sodium cation

We have performed detailed ab initio SCF calculations on the intermolecular interaction energies for one Na⁺ ion and one water molecule with two molecular fragments, one exemplifying a phospholipid (PL) head (PLHD) and the other, a phospholipid tail (PLTL). A 6-12-1 atom-atom pair potential for the interaction of a Na⁺ ion and water with a lysophosphatidyl-ethanolamine (LPEA) was derived from these results by a fitting procedure. This fitted potential was used to obtain isoenergy maps that provide energy profiles of the Na⁺ ion and the water around the phospholipids. The interaction of the Na⁺ ion with PL, as well as the interaction of water with the PL, can be visualized from these maps, which, as expected, show regions of hydrophilicity and hydrophobicity for the water and indicate a very strong binding site for the Na⁺ ion on the phosphate. It appears to be a stationary site that would limit the Na⁺ ion mobility. This binding site is located near the double-bonded oxygen atom of the phosphate group; its binding energy for Na⁺ is 67 kcal/mol. On the other hand the NH⁺ group of PLHD ahows strong electrostatic repulsion of Na⁺ while interacting with water with a binding energy of 13 kcal/mol. This potential energy well region for water is separated from another of similar depth near the phosphate by a barrier and both regions are expected to act as binding sites for water.     

Ion-Controlled Conformational Dynamics in the Outward-Open Transition from an Occluded State of LeuT

Neurotransmitter:sodium symporter (NSS) proteins are secondary Na⁺-driven active transporters that terminate neurotransmission by substrate uptake. Despite the availability of high-resolution crystal structures of a bacterial homolog of NSSs—Leucine Transporter (LeuT)—and extensive computational and experimental structure-function studies, unanswered questions remain regarding the transport mechanisms. We used microsecond atomistic molecular-dynamics (MD) simulations and free-energy computations to reveal ion-controlled conformational dynamics of LeuT in relation to binding affinity and selectivity of the more extracellularly positioned Na⁺ binding site (Na1 site). In the course of MD simulations starting from the occluded state with bound Na⁺, but in the absence of substrate, we find a spontaneous transition of the extracellular vestibule of LeuT into an outward-open conformation. The outward opening is enhanced by the absence of Na1 and modulated by the protonation state of the Na1-associated Glu-290. Consistently, the Na⁺ affinity for the Na1 site is inversely correlated with the extent of outward-open character and is lower than in the occluded state with bound substrate; however, the Na1 site retains its selectivity for Na⁺ over K⁺ in such conformational transitions. To the best of our knowledge, our findings shed new light on the Na⁺-driven transport cycle and on the symmetry in structural rearrangements for outward- and inward-open transitions.     

Conserved water molecules in the specificity pocket of serine proteases and the molecular mechanism of Na+ binding

Conservation of clusters of buried water molecules is a structural motif present throughout the serine protease family. Frequently, these clusters are shaped as water channels forming extensive hydrogen-bonding networks linked to the protein backbone. The most conspicuous example is the water channel present in the specificity pocket of trypsin and thrombin. In thrombin, other vitamin K-dependent proteases, and some complement factors, Na⁺ binds in this water channel and enhances allosterically the catalytic activity of the enzyme, whereas digestive and fibrinolytic proteases are devoid of such regulation. A comparative analysis of proteases with and without Na⁺ binding capability reveals the role of the water channel in maintaining the structural organization of the specificity pocket and in Na⁺ coordination. This enables the formulation of a molecular mechanism for Na⁺ binding in thrombin and leads to the identification of the structural changes necessary to engineer a functional Na⁺ site and enhanced catalytic activity in trypsin and other proteases. Proteins 30:34–42, 1998. © 1998 Wiley-Liss, Inc.     

Microseconds Simulations Reveal a New Sodium-binding Site and the Mechanism of Sodium-coupled Substrate Uptake by LeuT

The bacterial sodium-coupled leucine/alanine transporter LeuT is broadly used as a model system for studying the transport mechanism of neurotransmitters because of its structural and functional homology to mammalian transporters such as serotonin, dopamine, or norepinephrine transporters, and because of the resolution of its structure in different states. Although the binding sites (S1 for substrate, and Na1 and Na2 for two co-transported sodium ions) have been resolved, we still lack a mechanistic understanding of coupled Na⁺- and substrate-binding events. We present here results from extensive (>20 μs) unbiased molecular dynamics simulations generated using the latest computing technology. Simulations show that sodium binds initially the Na1 site, but not Na2, and, consistently, sodium unbinding/escape to the extracellular (EC) region first takes place at Na2, succeeded by Na1. Na2 diffusion back to the EC medium requires prior dissociation of substrate from S1. Significantly, Na⁺ binding (and unbinding) consistently involves a transient binding to a newly discovered site, Na1″, near S1, as an intermediate state. A robust sequence of substrate uptake events coupled to sodium bindings and translocations between those sites assisted by hydration emerges from the simulations: (i) bindings of a first Na⁺ to Na1″, translocation to Na1, a second Na⁺ to vacated Na1″ and then to Na2, and substrate to S1; (ii) rotation of Phe²⁵³ aromatic group to seclude the substrate from the EC region; and (iii) concerted tilting of TM1b and TM6a toward TM3 and TM8 to close the EC vestibule.     

Threshold Occupancy and Specific Cation Binding Modes in the Hammerhead Ribozyme Active Site are Required for Active Conformation

The relationship between formation of active in-line attack conformations and monovalent (Na⁺) and divalent (Mg²⁺) metal ion binding in hammerhead ribozyme (HHR) has been explored with molecular dynamics simulations. To stabilize repulsions between negatively charged groups, different requirements of the threshold occupancy of metal ions were observed in the reactant and activated precursor states both in the presence and in the absence of a Mg²⁺ in the active site. Specific bridging coordination patterns of the ions are correlated with the formation of active in-line attack conformations and can be accommodated in both cases. Furthermore, simulation results suggest that the HHR folds to form an electronegative recruiting pocket that attracts high local concentrations of positive charge. The present simulations help to reconcile experiments that probe the metal ion sensitivity of HHR catalysis and support the supposition that Mg²⁺, in addition to stabilizing active conformations, plays a specific chemical role in catalysis.     

Ion Selectivity of the Cytoplasmic Binding Sites of the Na,K-ATPase: I. Sodium Binding is Associated with a Conformational Rearrangement

To investigate Na⁺ binding to the ion-binding sites presented on the cytoplasmic side of the Na,K-ATPase, equilibrium Na⁺-titration experiments were performed using two fluorescent dyes, RH421 and FITC, to detect protein-specific actions. Fluorescence changes upon addition of Na⁺ in the presence of various Mg²⁺ concentrations were similar and could be fitted with a Hill function. The half-saturating concentrations and Hill coefficients determined were almost identical. As RH421 responds to binding of a Na⁺ ion to the third neutral site whereas FITC monitors conformational changes in the ATP-binding site or its environment, this result implies that electrogenic binding of the third Na⁺ ion is the trigger for a structural rearrangement of the ATP-binding moiety. This enables enzyme phosphorylation, which is accompanied by a fast occlusion of the Na⁺ ions and followed by the conformational transition E₁/E₂ of the protein. The coordinated action both at the ion and the nucleotide binding sites allows for the first time a detailed formulation of the mechanism of enzyme phosphorylation that occurs only when three Na⁺ ions are bound. Received: 8 October 1998/Revised: 29 December 1998     

Ion-Releasing State of a Secondary Membrane Transporter

The crystal structure of Na⁺-coupled galactose symporter (vSGLT) reports the transporter in its substrate-bound state, with a Na⁺ ion modeled in a binding site corresponding to that of a homologous protein, leucine transporter (LeuT). In repeated molecular dynamics simulations, however, we find the Na⁺ ion instable, invariably and spontaneously diffusing out of the transporter through a pathway lined by D189, which appears to facilitate the diffusion of the ion toward the cytoplasm. Further analysis of the trajectories and close structural examination, in particular, comparison of the Na⁺-binding sites of vSGLT and LeuT, strongly indicates that the crystal structure of vSGLT actually represents an ion-releasing state of the transporter. The observed dynamics of the Na⁺ ion, in contrast to the substrate, also suggests that the cytoplasmic release of the Na⁺ ion precedes that of the substrate, thus shedding light on a key step in the transport cycle of this secondary transporter.     

The Melibiose Transporter of Escherichia coli: CRITICAL CONTRIBUTION OF LYS-377 TO THE STRUCTURAL ORGANIZATION OF THE INTERACTING SUBSTRATE BINDING SITES*

We examine the role of Lys-377, the only charged residue in helix XI, on the functional mechanism of the Na⁺-sugar melibiose symporter from Escherichia coli. Intrinsic fluorescence, FRET, and Fourier transform infrared difference spectroscopy reveal that replacement of Lys-377 with either Cys, Val, Arg, or Asp disables both Na⁺ and melibiose binding. On the other hand, molecular dynamics simulations extending up to 200–330 ns reveal that Lys-377 (helix XI) interacts with the anionic side chains of two of the three putative ligands for cation binding (Asp-55 and Asp-59 in helix II). When Asp-59 is protonated during the simulations, Lys-377 preferentially interacts with Asp-55. Interestingly, when a Na⁺ ion is positioned in the Asp-55-Asp-59 environment, Asp-124 in helix IV (a residue essential for melibiose binding) reorients and approximates the Asp-55-Asp-59 pair, and all three acidic side chains act as Na⁺ ligands. Under these conditions, the side chain of Lys-377 interacts with the carboxylic moiety of these three Asp residues. These data highlight the crucial role of the Lys-377 residue in the spatial organization of the Na⁺ binding site. Finally, the analysis of the second-site revertants of K377C reveals that mutation of Ile-22 (in helix I) preserves Na⁺ binding, whereas that of melibiose is largely abolished according to spectroscopic measurements. This amino acid is located in the border of the sugar-binding site and might participate in sugar binding through apolar interactions.     

Ligand Binding Shuttles Thrombin along a Continuum of Zymogen- and Proteinase-like States

The critical and multiple roles of thrombin in blood coagulation are regulated by ligands and cofactors. Zymogen activation imparts proteolytic activity to thrombin and also affects the binding of ligands to its two principal exosites. We have used the activation peptide fragment 1.2 (F12), a ligand for anion-binding exosite 2, to probe the zymogenicity of thrombin by isothermal titration calorimetry. We show that F12 binding is sensitive to subtle aspects of proteinase formation beyond simply reporting on zymogen cleavage. Large thermodynamic differences in F12 binding distinguish between a series of thrombin species poised along the transition of zymogen to proteinase. Active-site ligands transitioned a zymogen-like state to a proteinase-like state. Conversely, removal of Na⁺ converted proteinase-like thrombin to a more zymogen-like form. Thrombin mutants, with deformed x-ray structures, previously considered to be emblematic of specific regulated states of the enzyme, are instead shown to be variously zymogen-like and can be made proteinase-like by active-site ligation. Thermodynamic linkage between anion-binding exosite 2, the Na⁺-binding site, and the active site arises from interconversions of thrombin between a continuum of zymogen- and proteinase-like states. These interconversions, reciprocally regulated by different ligands, cast new light on the problem of thrombin allostery and provide a thermodynamic framework to explain the regulation of thrombin by different ligands.