The Na+ binding site of thrombin

Dynamic allostery in thrombin-a review

Thrombin is a serine protease that catalyzes a large number of different reactions including proteolytic cleave of fibrinogen to make the fibrin clot (procoagulant activity), of the protease activated receptors (for cell signaling) and of protein C generating activated protein C (anticoagulant activity). Thrombin has an effector binding site called the anion binding exosite 1 that is allosterically coupled to the active site. In this review, we survey results from thermodynamic characterization of the allosteric coupling as well as hydrogen-deuterium exchange mass spectrometry to reveal which parts of the thrombin structure are changed upon effector binding and/or mutagenesis, and finally NMR spectroscopy to characterize the different timescales of motions elicited by the effectors. We also relate the experimental work to computational network analysis of the thrombin-thrombomodulin complex.

Reaction within the coulomb-cage; science in retrospect

The Coulomb-cage is defined as the space where the electrostatic interaction between two bodies is more intensive than the thermal energy (k_BT). For small molecule, the Coulomb-cage is a small sphere, extending only few water molecules towards the bulk and its radius is sensitive to the ionic strength of the solution. For charged proteins or membranal structures, the Coulomb-cage can engulf large fraction of the surface and provides a preferred pathway for ion propagation along the surface. Similarly, electrostatic potential at the inner space of a channel can form preferential trajectories passage for ions. The dynamics of ions inside the Coulomb-cage of ions was formulated by the studies of proton-anion recombination of excited photoacids. In the present article, we recount the study of intra- Coulomb-cage reaction taking place on the surface of macro-molecular bodies like micelles, membranes, proteins and intra-protein cavities. The study progressed stepwise, tracing the dynamics of a proton ejected from a photo-acid molecule located at defined sites (on membrane, inter-membrane space, active site of enzyme, inside Large Pore Channels etc.). Accumulation of experimental observations encouraged us to study of the reaction mechanism by molecular dynamics simulations of ions within the Coulomb-cage of proteins surface or inside large pores. The intra-Coulomb-cage proton transfer events follows closely the fine structure of the electrostatic field inside the cage and reflects the shape of nearby dielectric boundaries, the temporal ordering of the solvent molecules and the structural fluctuations of the charged side chains. The article sums some 40 years of research, which in retrospect clarifies the intra-Coulomb-cage reaction mechanism.

Comparative sequence analysis of vitamin K-dependent coagulation factors

BACKGROUND

Prothrombin, protein C, and factors VII, IX, and X are vitamin K (VK)-dependent coagulation proteins that play an important role in the initiation, amplification, and subsequent attenuation of the coagulation response. Blood coagulation evolved in the common vertebrate ancestor as a specialization of the complement system and immune response, which in turn bear close evolutionary ties with developmental enzyme cascades. There is currently no comprehensive analysis of the evolutionary changes experienced by these coagulation proteins during the radiation of vertebrates and little is known about conservation of residues that are important for zymogen activation and catalysis.

OBJECTIVES

To characterize the conservation level of functionally important residues among VK-dependent coagulation proteins from different vertebrate lineages.

METHODS

The conservation level of residues important for zymogen activation and catalysis was analyzed in >1600 primary sequences of VK-dependent proteins.

RESULTS

Functionally important residues are most conserved in prothrombin and least conserved in protein C. Some of the most profound functional modifications in protein C occurred in the ancestor of bony fish when the basic residue in the activation site was replaced by an aromatic residue. Furthermore, during the radiation of placental mammals from marsupials, protein C acquired a cysteine-rich insert that introduced an additional disulfide in the EGF1 domain and evolved a proprotein convertase cleavage site in the activation peptide linker that also became significantly elongated.

CONCLUSIONS

Sequence variabilities at functionally important residues may lead to interspecies differences in the zymogen activation and catalytic properties of orthologous VK-dependent proteins.

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 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.

Antithrombin Resistance Rescues Clotting Defect of Homozygous Prothrombin-Y510N Dysprothrombinemia

A patient with hematuria in our clinic was diagnosed with urolithiasis. Analysis of the patient's plasma clotting time indicated that both activated partial thromboplastin time (52.6 seconds) and prothrombin time (19.4 seconds) are prolonged and prothrombin activity is reduced to 12.4% of normal, though the patient exhibited no abnormal bleeding phenotype and a prothrombin antigen level of 87.9%. Genetic analysis revealed the patient is homozygous for prothrombin Y510N mutation. We expressed and characterized the prothrombin-Y510N variant in appropriate coagulation assays and found that the specificity constant for activation of the mutant zymogen by factor Xa is impaired approximately fivefold. Thrombin generation assay using patient's plasma and prothrombin-deficient plasma supplemented with either wild-type or prothrombin-Y510N revealed that both peak height and time to peak for the prothrombin mutant are decreased; however, the endogenous thrombin generation potential is increased. Further analysis indicated that the thrombin mutant exhibits resistance to antithrombin and is inhibited by the serpin with approximately 12-fold slower rate constant. Protein C activation by thrombin-Y510N was also decreased by approximately 10-fold; however, thrombomodulin overcame the catalytic defect. The Na⁺-concentration-dependence of the amidolytic activities revealed that the dissociation constant for the interaction of Na⁺ with the mutant has been elevated approximately 20-fold. These results suggest that Y510 (Y184a in chymotrypsin numbering) belongs to network of residues involved in binding Na⁺. A normal protein C activation by thrombin-Y510N suggests that thrombomodulin modulates the conformation of the Na⁺-binding loop of thrombin. The clotting defect of thrombin-Y510N appears to be compensated by its markedly lower reactivity with antithrombin, explaining patient's normal hemostatic phenotype.

Structure of human factor VIIa-soluble tissue factor with calcium, magnesium and rubidium

Coagulation factor VIIa (FVIIa) consists of a γ-carboxyglutamic acid (GLA) domain, two epidermal growth factor-like (EGF) domains and a protease domain. FVIIa binds three Mg2+ ions and four Ca2+ ions in the GLA domain, one Ca2+ ion in the EGF1 domain and one Ca2+ ion in the protease domain. Further, FVIIa contains an Na+ site in the protease domain. Since Na+ and water share the same number of electrons, Na+ sites in proteins are difficult to distinguish from waters in X-ray structures. Here, to verify the Na+ site in FVIIa, the structure of the FVIIa-soluble tissue factor (TF) complex was solved at 1.8 Å resolution containing Mg2+, Ca2+ and Rb+ ions. In this structure, Rb+ replaced two Ca2+ sites in the GLA domain and occupied three non-metal sites in the protease domain. However, Rb+ was not detected at the expected Na+ site. In kinetic experiments, Na+ increased the amidolytic activity of FVIIa towards the synthetic substrate S-2288 (H-D-Ile-Pro-Arg-p-nitroanilide) by ∼20-fold; however, in the presence of Ca2+, Na+ had a negligible effect. Ca2+ increased the hydrolytic activity of FVIIa towards S-2288 by ∼60-fold in the absence of Na+ and by ∼82-fold in the presence of Na+. In molecular-dynamics simulations, Na+ stabilized the two Na+-binding loops (the 184-loop and 220-loop) and the TF-binding region spanning residues 163-180. Ca2+ stabilized the Ca2+-binding loop (the 70-loop) and Na+-binding loops but not the TF-binding region. Na+ and Ca2+ together stabilized both the Na+-binding and Ca2+-binding loops and the TF-binding region. Previously, Rb+ has been used to define the Na+ site in thrombin; however, it was unsuccessful in detecting the Na+ site in FVIIa. A conceivable explanation for this observation is provided.

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.

Sodium-induced population shift drives activation of thrombin

The equilibrium between active E and inactive E* forms of thrombin is assumed to be governed by the allosteric binding of a Na⁺ ion. Here we use molecular dynamics simulations and Markov state models to sample transitions between active and inactive states. With these calculations we are able to compare thermodynamic and kinetic properties depending on the presence of Na⁺. For the first time, we directly observe sodium-induced conformational changes in long-timescale computer simulations. Thereby, we are able to explain the resulting change in activity. We observe a stabilization of the active form in presence of Na⁺ and a shift towards the inactive form in Na⁺-free simulations. We identify key structural features to quantify and monitor this conformational shift. These include the accessibility of the S1 pocket and the reorientation of W215, of R221a and of the Na⁺ loop. The structural characteristics exhibit dynamics at various timescales: Conformational changes in the Na⁺ binding loop constitute the slowest observed movement. Depending on its orientation, it induces conformational shifts in the nearby substrate binding site. Only after this shift, residue W215 is able to move freely, allowing thrombin to adopt a binding-competent conformation.

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.

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.

Probing light chain mutation effects on thrombin via molecular dynamics simulations and machine learning

Thrombin is a key component for chemotherapeutic and antithrombotic therapy development. As the physiologic and pathologic roles of the light chain still remain vague, here, we continue previous efforts to understand the impacts of the disease-associated single deletion of LYS9 in the light chain. By combining supervised and unsupervised machine learning methodologies and more traditional structural analyses on data from 10 μs molecular dynamics simulations, we show that the conformational ensemble of the ΔK9 mutant is significantly perturbed. Our analyses consistently indicate that LYS9 deletion destabilizes both the catalytic cleft and regulatory functional regions and result in some conformational changes that occur in tens to hundreds of nanosecond scaled motions. We also reveal that the two forms of thrombin each prefer a distinct binding mode of a Na+ ion. We expand our understanding of previous experimental observations and shed light on the mechanisms of the LYS9 deletion associated bleeding disorder by providing consistent but more quantitative and detailed structural analyses than early studies in literature. With a novel application of supervised learning, i.e. the decision tree learning on the hydrogen bonding features in the wild-type and ΔK9 mutant forms of thrombin, we predict that seven pairs of critical hydrogen bonding interactions are significant for establishing distinct behaviors of wild-type thrombin and its ΔK9 mutant form. Our calculations indicate the LYS9 in the light chain has both localized and long-range allosteric effects on thrombin, supporting the opinion that light chain has an important role as an allosteric effector.

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.

Predicting Displaceable Water Sites Using Mixed-Solvent Molecular Dynamics

Water molecules are an important factor in protein-ligand binding. Upon binding of a ligand with a protein's surface, waters can either be displaced by the ligand or may be conserved and possibly bridge interactions between the protein and ligand. Depending on the specific interactions made by the ligand, displacing waters can yield a gain in binding affinity. The extent to which binding affinity may increase is difficult to predict, as the favorable displacement of a water molecule is dependent on the site-specific interactions made by the water and the potential ligand. Several methods have been developed to predict the location of water sites on a protein's surface, but the majority of methods are not able to take into account both protein dynamics and the interactions made by specific functional groups. Mixed-solvent molecular dynamics (MixMD) is a cosolvent simulation technique that explicitly accounts for the interaction of both water and small molecule probes with a protein's surface, allowing for their direct competition. This method has previously been shown to identify both active and allosteric sites on a protein's surface. Using a test set of eight systems, we have developed a method using MixMD to identify conserved and displaceable water sites. Conserved sites can be determined by an occupancy-based metric to identify sites which are consistently occupied by water even in the presence of probe molecules. Conversely, displaceable water sites can be found by considering the sites which preferentially bind probe molecules. Furthermore, the inclusion of six probe types allows the MixMD method to predict which functional groups are capable of displacing which water sites. The MixMD method consistently identifies sites which are likely to be nondisplaceable and predicts the favorable displacement of water sites that are known to be displaced upon ligand binding.

Visualizing correlated motion with HDBSCAN clustering

Correlated motion analysis provides a method for understanding communication between and dynamic similarities of biopolymer residues and domains. The typical equal-time correlation matrices-frequently visualized with pseudo-colorings or heat maps-quickly convey large regions of highly correlated motion but hide more subtle similarities of motion. Here we propose a complementary method for visualizing correlations within proteins (or general biopolymers) that quickly conveys intuition about which residues have a similar dynamic behavior. For grouping residues, we use the recently developed non-parametric clustering algorithm HDBSCAN. Although the method we propose here can be used to group residues using correlation as a similarity matrix-the most straightforward and intuitive method-it can also be used to more generally determine groups of residues which have similar dynamic properties. We term these latter groups "Dynamic Domains", as they are based not on spatial closeness but rather closeness in the column space of a correlation matrix. We provide examples of this method across three human proteins of varying size and function-the Nf-Kappa-Beta essential modulator, the clotting promoter Thrombin and the mismatch repair protein (dimer) complex MutS-alpha. Although the examples presented here are from all-atom molecular dynamics simulations, this visualization technique can also be used on correlations matrices built from any ensembles of conformations from experiment or computation.

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.

Phenotypic and genetic analysis of dysprothrombinemia due to a novel homozygous mutation

OBJECTIVE

We study the phenotype and genotype of a novel gene mutation of factor II (FII) that leads to dysprothrombinemia, and do the meta-analysis to illuminate its molecular pathogenesis. It will further contribute to our comprehension of the pathogenesis of this type of disease.

METHODS

The prothrombin time (PT), activated partial thromboplastin time (APTT) and the activities of other factors were determined by the one-stage clotting method. The prothrombin antigen was measured with enzyme-linked immunosorbent assay (ELISA). Function of the mutant protein was evaluated by thrombin generation tests. Potential mutations in exons, exon-intron boundaries and 5', 3' untranslated sequences of prothrombin gene were screened by polymerase chain reaction and direct sequencing. Suspected mutations were confirmed by reverse sequencing. The structure change of this protein was analyzed by model and bioinformatics analyses.

RESULTS

Phenotypic analysis revealed that the proband had an obviously prolonged PT, APTT, reduced prothrombin activity but normal antigen levels. The other tests were normal. Sequencing analysis detected a homozygous g.26329T>G in the catalytic domain resulting in p.Tyr510Asp. His parents and uncle were heterozygous for this mutation. The thrombin generation test showed that the mutant protein had obstacles in thrombin generation. Bioinformatics and model analyses illuminated that the mutation will be probably damaging and perturbing the structure of Na+-binding site, which will affect the activation of prothrombin.

CONCLUSION

This was the first report of such a mutation in the position which was associated with dysprothrombinemia.

Uncovering Large-Scale Conformational Change in Molecular Dynamics without Prior Knowledge

As the length of molecular dynamics (MD) trajectories grows with increasing computational power, so does the importance of clustering methods for partitioning trajectories into conformational bins. Of the methods available, the vast majority require users to either have some a priori knowledge about the system to be clustered or to tune clustering parameters through trial and error. Here we present non-parametric uses of two modern clustering techniques suitable for first-pass investigation of an MD trajectory. Being non-parametric, these methods require neither prior knowledge nor parameter tuning. The first method, HDBSCAN, is fast-relative to other popular clustering methods-and is able to group unstructured or intrinsically disordered systems (such as intrinsically disordered proteins, or IDPs) into bins that represent global conformational shifts. HDBSCAN is also useful for determining the overall stability of a system-as it tends to group stable systems into one or two bins-and identifying transition events between metastable states. The second method, iMWK-Means, with explicit rescaling followed by K-Means, while slower than HDBSCAN, performs well with stable, structured systems such as folded proteins and is able to identify higher resolution details such as changes in relative position of secondary structural elements. Used in conjunction, these clustering methods allow a user to discern quickly and without prior knowledge the stability of a simulated system and identify both local and global conformational changes.

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.

Loop Electrostatics Asymmetry Modulates the Preexisting Conformational Equilibrium in Thrombin

Thrombin exists as an ensemble of active (E) and inactive (E*) conformations that differ in their accessibility to the active site. Here we show that redistribution of the E*-E equilibrium can be achieved by perturbing the electrostatic properties of the enzyme. Removal of the negative charge of the catalytic Asp102 or Asp189 in the primary specificity site destabilizes the E form and causes a shift in the 215-217 segment that compromises substrate entrance. Solution studies and existing structures of D102N document stabilization of the E* form. A new high-resolution structure of D189A also reveals the mutant in the collapsed E* form. These findings establish a new paradigm for the control of the E*-E equilibrium in the trypsin fold.

Dynamics Govern Specificity of a Protein-Protein Interface: Substrate Recognition by Thrombin

Biomolecular recognition is crucial in cellular signal transduction. Signaling is mediated through molecular interactions at protein-protein interfaces. Still, specificity and promiscuity of protein-protein interfaces cannot be explained using simplistic static binding models. Our study rationalizes specificity of the prototypic protein-protein interface between thrombin and its peptide substrates relying solely on binding site dynamics derived from molecular dynamics simulations. We find conformational selection and thus dynamic contributions to be a key player in biomolecular recognition. Arising entropic contributions complement chemical intuition primarily reflecting enthalpic interaction patterns. The paradigm “dynamics govern specificity” might provide direct guidance for the identification of specific anchor points in biomolecular recognition processes and structure-based drug design.

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.

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.

Why Ser and not Thr brokers catalysis in the trypsin fold

Although Thr is equally represented as Ser in the human genome and as a nucleophile is as good as Ser, it is never found in the active site of the large family of trypsin-like proteases that utilize the Asp/His/Ser triad. The molecular basis of the preference of Ser over Thr in the trypsin fold was investigated with X-ray structures of the thrombin mutant S195T free and bound to an irreversible active site inhibitor. In the free form, the methyl group of T195 is oriented toward the incoming substrate in a conformation seemingly incompatible with productive binding. In the bound form, the side chain of T195 is reoriented for efficient substrate acylation without causing steric clash within the active site. Rapid kinetics prove that this change is due to selection of an active conformation from a preexisting ensemble of reactive and unreactive rotamers whose relative distribution determines the level of activity of the protease. Consistent with these observations, the S195T substitution is associated with a weak yet finite activity that allows identification of an unanticipated important role for S195 as the end point of allosteric transduction in the trypsin fold. The S195T mutation abrogates the Na(+)-dependent enhancement of catalytic activity in thrombin, activated protein C, and factor Xa and significantly weakens the physiologically important allosteric effects of thrombomodulin on thrombin and of cofactor Va on factor Xa. The evolutionary selection of Ser over Thr in trypsin-like proteases was therefore driven by the need for high catalytic activity and efficient allosteric regulation.

Automated identification of elemental ions in macromolecular crystal structures

Many macromolecular model-building and refinement programs can automatically place solvent atoms in electron density at moderate-to-high resolution. This process frequently builds water molecules in place of elemental ions, the identification of which must be performed manually. The solvent-picking algorithms in phenix.refine have been extended to build common ions based on an analysis of the chemical environment as well as physical properties such as occupancy, B factor and anomalous scattering. The method is most effective for heavier elements such as calcium and zinc, for which a majority of sites can be placed with few false positives in a diverse test set of structures. At atomic resolution, it is observed that it can also be possible to identify tightly bound sodium and magnesium ions. A number of challenges that contribute to the difficulty of completely automating the process of structure completion are discussed.

Functional and structural characterization of Vibrio cholerae extracellular serine protease B, VesB

The chymotrypsin subfamily A of serine proteases consists primarily of eukaryotic proteases, including only a few proteases of bacterial origin. VesB, a newly identified serine protease that is secreted by the type II secretion system in Vibrio cholerae, belongs to this subfamily. VesB is likely produced as a zymogen because sequence alignment with trypsinogen identified a putative cleavage site for activation and a catalytic triad, His-Asp-Ser. Using synthetic peptides, VesB efficiently cleaved a trypsin substrate, but not chymotrypsin and elastase substrates. The reversible serine protease inhibitor, benzamidine, inhibited VesB and served as an immobilized ligand for VesB affinity purification, further indicating its relationship with trypsin-like enzymes. Consistent with this family of serine proteases, N-terminal sequencing implied that the propeptide is removed in the secreted form of VesB. Separate mutagenesis of the activation site and catalytic serine rendered VesB inactive, confirming the importance of these features for activity, but not for secretion. Similar to trypsin but, in contrast to thrombin and other coagulation factors, Na(+) did not stimulate the activity of VesB, despite containing the Tyr(250) signature. The crystal structure of catalytically inactive pro-VesB revealed that the protease domain is structurally similar to trypsinogen. The C-terminal domain of VesB was found to adopt an immunoglobulin (Ig)-fold that is structurally homologous to Ig-folds of other extracellular Vibrio proteins. Possible roles of the Ig-fold domain in stability, substrate specificity, cell surface association, and type II secretion of VesB, the first bacterial multidomain trypsin-like protease with known structure, are discussed.

Acidosis, magnesium and acetylsalicylic acid: effects on thrombin

Thrombin, an enzyme from the hydrolase family, is the main component of the blood coagulation system. In ischemic stroke it acts as a serine protease that converts soluble fibrinogen into insoluble strands of fibrin forming blood clots in the brain. It has been found to phosphoresce at room temperature in the millisecond and microsecond ranges. The phosphorescence of thrombin was studied under physiological conditions, in acidosis (decrease of pH from 8.0 to 5.0) and on the addition of salts (magnesium sulfate and sodium chloride) and of acetylsalicylic acid, and its connection with thrombin function is discussed. Acidosis significantly increased the internal dynamics of thrombin. We propose that lactate-acidosis plays a protective role in stroke, preventing the formation of clots. The addition of NaCl and MgSO(4) in different concentrations increased the internal dynamics of thrombin. Also, the addition of MgSO(4) decreased thrombin-induced platelet aggregation. However, magnesium sulfate and acetylsalicylic acid in the therapeutic concentrations used for treatment of ischemic stroke had no effect on thrombin internal dynamics. The data obtained will help to elucidate the conformational stability of thrombin under conditions modulating lactate-acidosis and in the presence of magnesium sulfate.

Diversity of allosteric regulation in proteases

Allostery is a fundamental regulatory mechanism that is based on a functional modulation of a site by a distant site. Allosteric regulation can be triggered by binding of diverse allosteric effectors, ranging from small molecules to macromolecules, and is therefore offering promising opportunities for functional modulation in a wide range of applications including the development of chemical probes or drug discovery. Here, we provide an overview of key classes of allosteric protease effectors, their corresponding molecular mechanisms, and their practical implications.

Evolution of thrombin and other hemostatic proteases by survey of protochordate, hemichordate, and echinoderm genomes

Protochordate genomes enable a prevalence of hemostasis evolution. Broad searches were performed for homologs of human serine proteases of hemostasis on the genomes of Branchiostoma floridae, Saccoglossus kowalevskii, and Strongylocentrotus purpuratus. Sequences were analyzed by multiple bioinformatic tools. The survey revealed numerous homologous components. Amphioxus was rich in some serine proteases not accompanied by gamma-carboxyglutamic or kringle domains similar more to thrombin than to other coagulation factors. The serine proteases found in amphioxus exhibited the attributes similar to those of thrombin by phylogeny relationships, sequence conservation, gene synteny, spatial structure, and ligand docking. A few plasminogen- and plasminogen activators-like proteases with kringles were also present. Those serine proteases demonstrated the greatest proximity rather to plasminogen or plasminogen activators than to thrombin. Searching for homologs of serine protease hemostatic factors in acorn worm and sea urchin revealed several components similar to those found in amphioxus. Hypothetically, the common ancestor of chordates had three separate serine proteases that evolved independently into immunoglobulin-like and kringle proteases in lancelets, and prothrombin, plasminogen activators, and plasminogen in vertebrates. Ancestral proteases evolved in vertebrates into hemostasis factors after merging the proper N-terminal domains and duplications.

Structural biology of factor VIIa/tissue factor initiated coagulation

Factor VII (FVII) consists of an N-terminal gamma-carboxyglutamic acid domain followed by two epidermal growth factor-like (EGF1 and EGF2) domains and the C-terminal protease domain. Activation of FVII results in a two-chain FVIIa molecule consisting of a light chain (Gla-EGF1-EGF2 domains) and a heavy chain (protease domain) held together by a single disulfide bond. During coagulation, the complex of tissue factor (TF, a transmembrane glycoprotein) and FVIIa activates factor IX (FIX) and factor X (FX). FVIIa is structurally "zymogen-like" and when bound to TF, it is more "active enzyme-like." FIX and FX share structural homology with FVII. Three structural biology aspects of FVIIa/TF are presented in this review. One, regions in soluble TF (sTF) that interact with FVIIa as well as mapping of Ca2+, Mg2+, Na+ and Zn2+ sites in FVIIa and their functions; two, modeled interactive regions of Gla and EGF1 domains of FXa and FIXa with FVIIa/sTF; and three, incompletely formed oxyanion hole in FVIIa/sTF and its induction by substrate/inhibitor. Finally, an overview of the recognition elements in TF pathway inhibitor is provided.

Antithrombin effect of polyphenol-rich extracts from black chokeberry and grape seeds

Thrombin is a serine protease that cleaves the peptide bonds in proteins located on the carboxyl side of arginine. Thrombin plays a central role in thromboembolic diseases, which are the major cause of mortality. The aim of the study was to estimate the effects of plant extracts on proteolytic properties of thrombin. Thrombin was incubated with polyphenol-rich extracts from berries of Aronia melanocarpa or seeds of Vitis vinifera (0.5, 5, 50 µg/mL) and with polyphenols ((+)-catechin, (-)-epicatechin, gallic acid, chlorogenic acid, procyanidin B1, cyanidin, cyanidin 3-glucoside, quercetin). The in vitro experiments showed that both extracts in all used concentrations inhibited proteolytic activity of thrombin observed as inhibition of thrombin-induced fibrinogen polymerization, stabilized fibrin formation, and platelet aggregation. Moreover, thrombin amidolytic activity was inhibited by polyphenols belonging to the flavonoid class. Results presented in this study indicate that polyphenol-rich extracts from berries of A. melanocarpa and seeds of V. vinifera may become promising dietary supplements in the prevention of thrombotic states.

Structure and behavior of human α-thrombin upon ligand recognition: thermodynamic and molecular dynamics studies

Thrombin is a serine proteinase that plays a fundamental role in coagulation. In this study, we address the effects of ligand site recognition by alpha-thrombin on conformation and energetics in solution. Active site occupation induces large changes in secondary structure content in thrombin as shown by circular dichroism. Thrombin-D-Phe-Pro-Arg-chloromethyl ketone (PPACK) exhibits enhanced equilibrium and kinetic stability compared to free thrombin, whose difference is rooted in the unfolding step. Small-angle X-ray scattering (SAXS) measurements in solution reveal an overall similarity in the molecular envelope of thrombin and thrombin-PPACK, which differs from the crystal structure of thrombin. Molecular dynamics simulations performed with thrombin lead to different conformations than the one observed in the crystal structure. These data shed light on the diversity of thrombin conformers not previously observed in crystal structures with distinguished catalytic and conformational behaviors, which might have direct implications on novel strategies to design direct thrombin inhibitors.

Thrombin plasticity

Thrombin is the final protease generated in the blood coagulation cascade. It has multiple substrates and cofactors, and serves both pro- and anti-coagulant functions. How thrombin activity is directed throughout the evolution of a clot and the role of conformational change in determining thrombin specificity are issues that lie at the heart of the haemostatic balance. Over the last 20 years there have been a great number of studies supporting the idea that thrombin is an allosteric enzyme that can exist in two conformations differing in activity and specificity. However, recent work has shown that thrombin in its unliganded state is inherently flexible in regions that are important for activity. The effect of flexibility on activity is discussed in this review in context of the zymogen-to-protease conformational transition. Understanding thrombin function in terms of 'plasticity' provides a new conceptual framework for understanding regulation of enzyme activity in general. This article is part of a Special Issue entitled: Proteolysis 50 years after the discovery of lysosome.

Thrombin-aptamer recognition: a revealed ambiguity

Aptamers are structured oligonucleotides that recognize molecular targets and can function as direct protein inhibitors. The best-known example is the thrombin-binding aptamer, TBA, a single-stranded 15-mer DNA that inhibits the activity of thrombin, the key enzyme of coagulation cascade. TBA folds as a G-quadruplex structure, as proved by its NMR structure. The X-ray structure of the complex between TBA and human α-thrombin was solved at 2.9-Å resolution, but did not provide details of the aptamer conformation and the interactions with the protein molecule. TBA is rapidly processed by nucleases. To improve the properties of TBA, a number of modified analogs have been produced. In particular, a modified TBA containing a 5′-5′ polarity inversion site, mTBA, has higher stability and higher affinity toward thrombin with respect to TBA, although it has a lower inhibitory activity. We present the crystal structure of the thrombin–mTBA complex at 2.15-Å resolution; the resulting model eventually provides a clear picture of thrombin–aptamers interaction, and also highlights the structural bases of the different properties of TBA and mTBA. Our findings open the way for a rational design of modified aptamers with improved potency as anticoagulant drugs.

Looking at the proteases from a simple perspective

Proteases have received enormous interest from the research and medical communities because of their significant roles in several human diseases. Some examples include the involvement of thrombin in thrombosis, HIV-1 protease in Acquired Immune Deficiency Syndrome, cruzain in Trypanosoma cruzi infection, and membrane-type 1 matrix metalloproteinase in tumor invasion and metastasis. Many efforts has been undertaken to design effective inhibitors featuring potent inhibitory activity, specificity, and metabolic stability to those proteases involved in such pathologies. Protease inhibitors usually target the active site, but some of them act by other inhibitory mechanisms. The understanding of the structure-function relationships of proteases and inhibitors has an impact on new inhibitor drugs designing. In this paper, the structures of four proteases (thrombin, HIV-protease, cruzain, and a matrix metalloproteinase) are briefly reviewed, and used as examples of the importance of proteases for the development of new treatment strategies, leading to a longer and healthier life.

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.

Redesigning allosteric activation in an enzyme

Enzyme activation by monovalent cations is widely documented in plants and the animal world. In type II enzymes, activation entails two steps: binding of the monovalent cation to its allosteric site and transduction of this event into enhanced catalytic activity. The effect has exquisite specificity for either Na(+) or K(+), the most abundant cations present in physiological environments. Enzymes requiring K(+) such as kinases and molecular chaperones are not activated as well or at all by the larger cation Cs(+) or the smaller cations Na(+) and Li(+). Enzymes requiring Na(+) such as β-galactosidase and clotting proteases are not activated as well by Li(+), or the larger cations K(+), Rb(+), and Cs(+). Efforts to switch specificity between Na(+) and K(+) in this large class of enzymes and completely redesign the mechanism of allosteric transduction leading to enhanced catalytic activity have so far been unsuccessful. Here we show how mutagenesis of two loops defining the Na(+) binding site of thrombin, a Na(+)-activated clotting protease, generates a construct that is most active in the presence of K(+) toward synthetic and physiological substrates. The effect is the result of a higher binding affinity and more efficient allosteric transduction of binding into enhanced catalytic activity for K(+) compared to Na(+), which represents a complete reversal of the properties of wild type. In addition, the construct features altered specificity toward physiological substrates resulting in a significant anticoagulant profile. The findings are relevant to all Na(+)-activated proteases involved in blood coagulation and the complement system.

Complex assemblies of factors IX and X regulate the initiation, maintenance, and shutdown of blood coagulation

Blood hemostasis is accomplished by a complex network of (anti-)coagulatory and fibrinolytic processes. These physiological processes are implemented by the assembly of multiprotein complexes involving both humoral and cellular components. Coagulation factor X, and particularly, factor IX, exemplify the dramatic enhancement that is obtained by the synergistic interaction of cell surface, inorganic and protein cofactors, protease, and substrate. With a focus on structure-function relationship, we review the current knowledge of activity modulation principles in the coagulation proteases factors IX and X and indicate future challenges for hemostasis research. This chapter is organized by describing the principles of hierarchical activation of blood coagulation proteases, including endogenous and exogenous protease activators, cofactor binding, substrate specificities, and protein inhibitors. We conclude by outlining pharmaceutical opportunities for unmet needs in hemophilia and thrombosis.