Crystallographic and kinetic evidence of allostery in a trypsin-like protease

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.

Structure-mechanics statistical learning uncovers mechanical relay in proteins

A protein's adaptive response to its substrates is one of the key questions driving molecular physics and physical chemistry. This work employs the recently developed structure-mechanics statistical learning method to establish a mechanical perspective. Specifically, by mapping all-atom molecular dynamics simulations onto the spring parameters of a backbone-side-chain elastic network model, the chemical moiety specific force constants (or mechanical rigidity) are used to assemble the rigidity graph, which is the matrix of inter-residue coupling strength. Using the S1A protease and the PDZ3 signaling domain as examples, chains of spatially contiguous residues are found to exhibit prominent changes in their mechanical rigidity upon substrate binding or dissociation. Such a mechanical-relay picture thus provides a mechanistic underpinning for conformational changes, long-range communication, and inter-domain allostery in both proteins, where the responsive mechanical hotspots are mostly residues having important biological functions or significant mutation sensitivity.

Protein residues exhibit specific routes of mechanical relay as the adaptive responses to substrate binding or dissociation. On such physically contiguous connections, residues experience prominent changes in their coupling strengths.

How Thrombomodulin Enables W215A/E217A Thrombin to Cleave Protein C but Not Fibrinogen

The W215A/E217A mutant thrombin is called "anticoagulant thrombin" because its activity toward its procoagulant substrate, fibrinogen, is reduced more than 500-fold whereas in the presence of thrombomodulin (TM) its activity toward its anticoagulant substrate, protein C, is reduced less than 10-fold. To understand how these mutations so dramatically alter one activity over the other, we compared the backbone dynamics of wild type thrombin to those of the W215A/E217A mutant thrombin by hydrogen-deuterium exchange coupled to mass spectrometry (HDX-MS). Our results show that the mutations cause the 170s, 180s, and 220s C-terminal β-barrel loops near the sites of mutation to exchange more, suggesting that the structure of this region is disrupted. Far from the mutation sites, residues at the N-terminus of the heavy chain, which need to be buried in the Ile pocket for correct structuring of the catalytic triad, also exchange much more than in wild type thrombin. TM binding causes reduced H/D exchange in these regions and also alters the dynamics of the β-strand that links the TM binding site to the catalytic Asp 102 in both wild type thrombin and in the W215A/E217A mutant thrombin. In contrast, whereas TM binding reduces the dynamics the 170, 180 and 220 s C-terminal β-barrel loops in WT thrombin, this region remains disordered in the W215A/E217A mutant thrombin. Thus, TM partially restores the catalytic activity of W215A/E217A mutant thrombin by allosterically altering its dynamics in a manner similar to that of wild type 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.

Conformational Plasticity-Rigidity Axis of the Coagulation Factor VII Zymogen Elucidated by Atomistic Simulations of the N-Terminally Truncated Factor VIIa Protease Domain

The vast majority of coagulation factor VII (FVII), a trypsin-like protease, circulates as the inactive zymogen. Activated FVII (FVIIa) is formed upon proteolytic activation of FVII, where it remains in a zymogen-like state and it is fully activated only when bound to tissue factor (TF). The catalytic domains of trypsin-like proteases adopt strikingly similar structures in their fully active forms. However, the dynamics and structures of the available corresponding zymogens reveal remarkable conformational plasticity of the protease domain prior to activation in many cases. Exactly how ligands and cofactors modulate the conformational dynamics and function of these proteases is not entirely understood. Here, we employ atomistic simulations of FVIIa (and variants hereof, including a TF-independent variant and N-terminally truncated variants) to provide fundamental insights with atomistic resolution into the plasticity-rigidity interplay of the protease domain conformations that appears to govern the functional response to proteolytic and allosteric activation. We argue that these findings are relevant to the FVII zymogen, whose structure has remained elusive despite substantial efforts. Our results shed light on the nature of FVII and demonstrate how conformational dynamics has played a crucial role in the evolutionary adaptation of regulatory mechanisms that were not present in the ancestral trypsin. Exploiting this knowledge could lead to engineering of protease variants for use as next-generation hemostatic therapeutics.

Role of the activation peptide in the mechanism of protein C activation

Protein C is a natural anticoagulant activated by thrombin in a reaction accelerated by the cofactor thrombomodulin. The zymogen to protease conversion of protein C involves removal of a short activation peptide that, relative to the analogous sequence present in other vitamin K-dependent proteins, contains a disproportionately high number of acidic residues. Through a combination of bioinformatic, mutagenesis and kinetic approaches we demonstrate that the peculiar clustering of acidic residues increases the intrinsic disorder propensity of the activation peptide and adversely affects the rate of activation. Charge neutralization of the acidic residues in the activation peptide through Ala mutagenesis results in a mutant activated by thrombin significantly faster than wild type. Importantly, the mutant is also activated effectively by other coagulation factors, suggesting that the acidic cluster serves a protective role against unwanted proteolysis by endogenous proteases. We have also identified an important H-bond between residues T176 and Y226 that is critical to transduce the inhibitory effect of Ca2+ and the stimulatory effect of thrombomodulin on the rate of zymogen activation. These findings offer new insights on the role of the activation peptide in the function of protein C.

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.

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.

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.

Membrane interactions of intrinsically disordered proteins: The example of alpha-synuclein

Peripheral membrane proteins associate reversibly with biological membranes that, compared to protein binding partners, are structurally labile and devoid of specific binding pockets. Membranes in different subcellular compartments vary primarily in their chemical composition and physical properties, and recognition of these features is therefore critical for allowing such proteins to engage their proper membrane targets. Intrinsically disordered proteins (IDPs) are well-suited to accomplish this task using highly specific and low- to moderate-affinity interactions governed by recognition principles that are both similar to and different from those that mediate the membrane interactions of rigid proteins. IDPs have also evolved multiple mechanisms to regulate membrane (and other) interactions and achieve their impressive functional diversity. Moreover, IDP-membrane interactions may have a kinetic advantage in fast processes requiring rapid control of such interactions, such as synaptic transmission or signaling. Herein we review the biophysics, regulation and functional implications of IDP-membrane interactions and include a brief overview of some of the methods that can be used to study such interactions. At each step, we use the example of alpha-synuclein, a protein involved in the pathogenesis of Parkinson's disease and one of the best characterized membrane-binding IDP, to illustrate some of the principles discussed.

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.

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.

Discovery of a novel conformational equilibrium in urokinase-type plasminogen activator

Although trypsin-like serine proteases have flexible surface-exposed loops and are known to adopt higher and lower activity conformations, structural determinants for the different conformations have remained largely obscure. The trypsin-like serine protease, urokinase-type plasminogen activator (uPA), is central in tissue remodeling processes and also strongly implicated in tumor metastasis. We solved five X-ray crystal structures of murine uPA (muPA) in the absence and presence of allosteric molecules and/or substrate-like molecules. The structure of unbound muPA revealed an unsuspected non-chymotrypsin-like protease conformation in which two β-strands in the core of the protease domain undergoes a major antiparallel-to-parallel conformational transition. We next isolated two anti-muPA nanobodies; an active-site binding nanobody and an allosteric nanobody. Crystal structures of the muPA:nanobody complexes and hydrogen-deuterium exchange mass spectrometry revealed molecular insights about molecular factors controlling the antiparallel-to-parallel equilibrium in muPA. Together with muPA activity assays, the data provide valuable insights into regulatory mechanisms and conformational flexibility of uPA and trypsin-like serine proteases in general.

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.

NMR reveals a dynamic allosteric pathway in thrombin

Although serine proteases are found ubiquitously in both eukaryotes and prokaryotes, and they comprise the largest of all of the peptidase families, their dynamic motions remain obscure. The backbone dynamics of the coagulation serine protease, apo-thrombin (S195M-thrombin), were compared to the substrate-bound form (PPACK-thrombin). R₁, R₂, ¹⁵N-{¹H}NOEs, and relaxation dispersion NMR experiments were measured to capture motions across the ps to ms timescale. The ps-ns motions were not significantly altered upon substrate binding. The relaxation dispersion data revealed that apo-thrombin is highly dynamic, with μs-ms motions throughout the molecule. The region around the N-terminus of the heavy chain, the Na⁺-binding loop, and the 170 s loop, all of which are implicated in allosteric coupling between effector binding sites and the active site, were dynamic primarily in the apo-form. Most of the loops surrounding the active site become more ordered upon PPACK-binding, but residues in the N-terminal part of the heavy chain, the γ-loop, and anion-binding exosite 1, the main allosteric binding site, retain μs-ms motions. These residues form a dynamic allosteric pathway connecting the active site to the main allosteric site that remains in the substrate-bound form.

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.

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.

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.

Discrimination between conformational selection and induced fit protein-ligand binding using Integrated Global Fit analysis

Molecular recognition between proteins and small molecule ligands is at the heart of biological function in cellular systems and the basis of modern rational drug development. Therefore, the mechanisms governing protein-ligand interaction have been objects of research for many decades. The last 15 years has seen a revival of a discussion whether conformational selection (CS) or induced fit (IF) is the most relevant binding mechanism. A decreasing observed rate constant, k obs, with increasing ligand concentration was considered to be a hallmark of CS, but according to contemporary knowledge, a positive saturating behavior of k obs can be explained by both CS and IF mechanisms. The only currently recognized kinetic method to differentiate between both binding mechanisms includes the measurement of two separate series of binding kinetics with variation of either protein or ligand under pseudo-first-order conditions. This study avoids the disadvantage of high protein concentrations and provides evidence that a comprehensive Integrated Global Fit analysis of sets of binding kinetics with just varied ligand concentration in combination with equilibrium data and optional displacement kinetics can effectively differentiate between CS and IF binding mechanisms. The limiting situation, when physical binding dominates over the previous (CS) or subsequent (IF) conformational changes, is carefully analyzed. Finally, the relevance of kinetic methods and the elucidation of more complex binding mechanisms are discussed for advanced rational selection and optimization of drug candidates.

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.

Protein conformational plasticity and complex ligand-binding kinetics explored by atomistic simulations and Markov models

Understanding the structural mechanisms of protein–ligand binding and their dependence on protein sequence and conformation is of fundamental importance for biomedical research. Here we investigate the interplay of conformational change and ligand-binding kinetics for the serine protease Trypsin and its competitive inhibitor Benzamidine with an extensive set of 150 μs molecular dynamics simulation data, analysed using a Markov state model. Seven metastable conformations with different binding pocket structures are found that interconvert at timescales of tens of microseconds. These conformations differ in their substrate-binding affinities and binding/dissociation rates. For each metastable state, corresponding solved structures of Trypsin mutants or similar serine proteases are contained in the protein data bank. Thus, our wild-type simulations explore a space of conformations that can be individually stabilized by adding ligands or making suitable changes in protein sequence. These findings provide direct evidence of conformational plasticity in receptors.

Conformational plasticity influences several aspects of protein function. Here the authors combine extensive MD simulations with Markov state models—using trypsin as model—to reveal new mechanistic details of how conformational plasticity influence ligand-receptors interactions.

Structure-function analyses of human kallikrein-related peptidase 2 establish the 99-loop as master regulator of activity

Human kallikrein-related peptidase 2 (KLK2) is a tryptic serine protease predominantly expressed in prostatic tissue and secreted into prostatic fluid, a major component of seminal fluid. Most likely it activates and complements chymotryptic KLK3 (prostate-specific antigen) in cleaving seminal clotting proteins, resulting in sperm liquefaction. KLK2 belongs to the “classical” KLKs 1–3, which share an extended 99- or kallikrein loop near their non-primed substrate binding site. Here, we report the 1.9 Å crystal structures of two KLK2-small molecule inhibitor complexes. In both structures discontinuous electron density for the 99-loop indicates that this loop is largely disordered. We provide evidence that the 99-loop is responsible for two biochemical peculiarities of KLK2, i.e. reversible inhibition by micromolar Zn²⁺ concentrations and permanent inactivation by autocatalytic cleavage. Indeed, several 99-loop mutants of KLK2 displayed an altered susceptibility to Zn²⁺, which located the Zn²⁺ binding site at the 99-loop/active site interface. In addition, we identified an autolysis site between residues 95e and 95f in the 99-loop, whose elimination prevented the mature enzyme from limited autolysis and irreversible inactivation. An exhaustive comparison of KLK2 with related structures revealed that in the KLK family the 99-, 148-, and 220-loop exist in open and closed conformations, allowing or preventing substrate access, which extends the concept of conformational selection in trypsin-related proteases. Taken together, our novel biochemical and structural data on KLK2 identify its 99-loop as a key player in activity regulation.

The linker connecting the two kringles plays a key role in prothrombin activation

The zymogen prothrombin is proteolytically converted by factor Xa to the active protease thrombin in a reaction that is accelerated >3,000-fold by cofactor Va. This physiologically important effect is paradigmatic of analogous cofactor-dependent reactions in the coagulation and complement cascades, but its structural determinants remain poorly understood. Prothrombin has three linkers connecting the N-terminal Gla domain to kringle-1 (Lnk1), the two kringles (Lnk2), and kringle-2 to the C-terminal protease domain (Lnk3). Recent developments indicate that the linkers, and particularly Lnk2, confer on the zymogen significant flexibility in solution and enable prothrombin to sample alternative conformations. The role of this flexibility in the context of prothrombin activation was tested with several deletions. Removal of Lnk2 in almost its entirety (ProTΔ146-167) drastically reduces the enhancement of thrombin generation by cofactor Va from >3,000-fold to 60-fold because of a significant increase in the rate of activation in the absence of cofactor. Deletion of Lnk2 mimics the action of cofactor Va and offers insights into how prothrombin is activated at the molecular level. The crystal structure of ProTΔ146-167 reveals a contorted architecture where the domains are not vertically stacked, kringle-1 comes within 9 Å of the protease domain, and the Gla-domain primed for membrane binding comes in contact with kringle-2. These findings broaden our molecular understanding of a key reaction of the blood coagulation cascade where cofactor Va enhances activation of prothrombin by factor Xa by compressing Lnk2 and morphing prothrombin into a conformation similar to the structure of ProTΔ146-167.

Multiple conformational selection and induced fit events take place in allosteric propagation

The fact that we observe a single conformational selection event during binding does not necessarily mean that only a single conformational selection event takes place, even though this is the common assumption. Here we suggest that conformational selection takes place not once in a given binding/allosteric event, but at every step along the allosteric pathway. This view generalizes conformational selection and makes it applicable also to other allosteric events, such as post-translational modifications (PTMs) and photon absorption. Similar to binding, at each step along a propagation pathway, conformational selection is coupled with induced fit which optimizes the interactions. Thus, as in binding, the allosteric effects induced by PTMs and light relate not only to population shift; but to conformational selection as well. Conformational selection and population shift take place conjointly.

Histone H4 promotes prothrombin autoactivation

Recent studies have documented the ability of prothrombin to spontaneously convert to the mature protease thrombin when Arg-320 becomes exposed to solvent for proteolytic attack upon mutation of residues in the activation domain. Whether prothrombin autoactivation occurs in the wild-type under conditions relevant to physiology remains unknown. Here, we report that binding of histone H4 to prothrombin under physiological conditions generates thrombin by autoactivation. The effect is abrogated by mutation of the catalytic Ser-525 and requires the presence of the Gla domain. Fluorescence titrations document direct binding of histone H4 to prothrombin with an affinity in the low nm range. Stopped flow data and luminescence resonance energy transfer measurements indicate that the binding mechanism obeys conformational selection. Among the two conformations of prothrombin, collapsed and fully extended, histone H4 binds selectively to the collapsed form and induces a transition toward a new conformation where the distance between Ser-101 in kringle-1 and Ser-210 in kringle-2 increases by 13 Å. These findings confirm the molecular plasticity of prothrombin emerged from recent structural studies and suggest that different conformations of the inter-kringle linker domain determine the functional behavior of prothrombin. The results also broaden our mechanistic understanding of the prothrombotic phenotype observed during cellular damage due to the release of histones in the blood stream. Prothrombin autoactivation induced by histone H4 emerges as a mechanism of pathophysiological relevance through which thrombin is generated independently of activation of the coagulation cascade.

Essential role of conformational selection in ligand binding

Two competing and mutually exclusive mechanisms of ligand recognition - conformational selection and induced fit - have dominated our interpretation of ligand binding in biological macromolecules for almost six decades. Conformational selection posits the pre-existence of multiple conformations of the macromolecule from which the ligand selects the optimal one. Induced fit, on the other hand, postulates the existence of conformational rearrangements of the original conformation into an optimal one that are induced by binding of the ligand. In the former case, conformational transitions precede the binding event; in the latter, conformational changes follow the binding step. Kineticists have used a facile criterion to distinguish between the two mechanisms based on the dependence of the rate of relaxation to equilibrium, kobs, on the ligand concentration, [L]. A value of kobs decreasing hyperbolically with [L] has been seen as diagnostic of conformational selection, while a value of kobs increasing hyperbolically with [L] has been considered diagnostic of induced fit. However, this simple conclusion is only valid under the rather unrealistic assumption of conformational transitions being much slower than binding and dissociation events. In general, induced fit only produces values of kobs that increase with [L] but conformational selection is more versatile and is associated with values of kobs that increase with, decrease with or are independent of [L]. The richer repertoire of kinetic properties of conformational selection applies to kinetic mechanisms with single or multiple saturable relaxations and explains the behavior of nearly all experimental systems reported in the literature thus far. Conformational selection is always sufficient and often necessary to account for the relaxation kinetics of ligand binding to a biological macromolecule and is therefore an essential component of any binding mechanism. On the other hand, induced fit is never necessary and only sufficient in a few cases. Therefore, the long assumed importance and preponderance of induced fit as a mechanism of ligand binding should be reconsidered.

The serine protease domain of MASP-3: enzymatic properties and crystal structure in complex with ecotin

Mannan-binding lectin (MBL), ficolins and collectin-11 are known to associate with three homologous modular proteases, the MBL-Associated Serine Proteases (MASPs). The crystal structures of the catalytic domains of MASP-1 and MASP-2 have been solved, but the structure of the corresponding domain of MASP-3 remains unknown. A link between mutations in the MASP1/3 gene and the rare autosomal recessive 3MC (Mingarelli, Malpuech, Michels and Carnevale,) syndrome, characterized by various developmental disorders, was discovered recently, revealing an unexpected important role of MASP-3 in early developmental processes. To gain a first insight into the enzymatic and structural properties of MASP-3, a recombinant form of its serine protease (SP) domain was produced and characterized. The amidolytic activity of this domain on fluorescent peptidyl-aminomethylcoumarin substrates was shown to be considerably lower than that of other members of the C1r/C1s/MASP family. The E. coli protease inhibitor ecotin bound to the SP domains of MASP-3 and MASP-2, whereas no significant interaction was detected with MASP-1, C1r and C1s. A tetrameric complex comprising an ecotin dimer and two MASP-3 SP domains was isolated and its crystal structure was solved and refined to 3.2 Å. Analysis of the ecotin/MASP-3 interfaces allows a better understanding of the differential reactivity of the C1r/C1s/MASP protease family members towards ecotin, and comparison of the MASP-3 SP domain structure with those of other trypsin-like proteases yields novel hypotheses accounting for its zymogen-like properties in vitro.

Differential protease activity augments polyphagy in Helicoverpa armigera

Helicoverpa armigera (Lepidoptera: Noctuidae) and other polyphagous agricultural pests are extending their plant host range and emerging as serious agents in restraining crop productivity. Dynamic regulation, coupled with a diversity of digestive and detoxifying enzymes, play a crucial role in the adaptation of polyphagous insects. To investigate the functional intricacy of serine proteases in the development and polyphagy of H. armigera, we profiled the expression of eight trypsin-like and four chymotrypsin-like phylogenetically diverse mRNAs from different life stages of H. armigera reared on nutritionally distinct host plants. These analyses revealed diet- and stage-specific protease expression patterns. The trypsins expressed showed structural variations, which might result in differential substrate specificity and interaction with inhibitors. Protease profiles in the presence of inhibitors and their mass spectrometric analyses revealed insight into their differential activity. These findings emphasize the differential expression of serine proteases and their consequences for digestive physiology in promoting polyphagy in H. armigera.

A specific antidote for dabigatran: functional and structural characterization

We present an antidote for dabigatran that effectively reverses its anticoagulative effect in human plasma in vitro and in rats in vivo. The antidote shares structural features with thrombin in the mode of binding but has no activity in coagulation tests.

Autoactivation of thrombin precursors

Trypsin-like proteases are synthesized as inactive zymogens and convert to the mature form upon activation by specific enzymes, often assisted by cofactors. Central to this paradigm is that the zymogen does not convert spontaneously to the mature enzyme, which in turn does not feed back to activate its zymogen form. In the blood, the zymogens prothrombin and prethrombin-2 require the prothrombinase complex to be converted to the mature protease thrombin, which is unable to activate prothrombin or prethrombin-2. Here, we show that replacement of key residues within the activation domain causes these zymogens to spontaneously convert to thrombin. The conversion is started by the zymogen itself, which is capable of binding ligands at the active site, and is abrogated by inactivation of the catalytic residue Ser-195. The product of autoactivation is functionally and structurally equivalent to wild-type thrombin. Zymogen autoactivation is explained by conformational selection, a basic property of the trypsin fold uncovered by structural and rapid kinetics studies. Both the zymogen and protease undergo a pre-existing equilibrium between active and inactive forms. The equilibrium regulates catalytic activity in the protease and has the potential to unleash activity in the zymogen to produce autoactivation. A new strategy emerges for the facile production of enzymes through zymogen autoactivation that is broadly applicable to trypsin-like proteases of biotechnological and clinical interest.

The dynamic structure of thrombin in solution

The backbone dynamics of human α-thrombin inhibited at the active site serine were analyzed using R(1), R(2), and heteronuclear NOE experiments, variable temperature TROSY 2D [(1)H-(15)N] correlation spectra, and R(ex) measurements. The N-terminus of the heavy chain, which is formed upon zymogen activation and inserts into the protein core, is highly ordered, as is much of the double beta-barrel core. Some of the surface loops, by contrast, remain very dynamic with order parameters as low as 0.5 indicating significant motions on the ps-ns timescale. Regions of the protein that were thought to be dynamic in the zymogen and to become rigid upon activation, in particular the γ-loop, the 180s loop, and the Na(+) binding site have order parameters below 0.8. Significant R(ex) was observed in most of the γ-loop, in regions proximal to the light chain, and in the β-sheet core. Accelerated molecular dynamics simulations yielded a molecular ensemble consistent with measured residual dipolar couplings that revealed dynamic motions up to milliseconds. Several regions, including the light chain and two proximal loops, did not appear highly dynamic on the ps-ns timescale, but had significant motions on slower timescales.

Investigating the role of T7 and T12 residues on the biological properties of thrombin-binding aptamer: enhancement of anticoagulant activity by a single nucleobase modification

An acyclic pyrimidine analogue, containing a five-member cycle fused on the pyrimidine ring, was synthesized and introduced at position 7 or 12 of the 15-mer oligodeoxynucleotide GGTTGGTGTGGTTGG, known as thrombin-binding aptamer (TBA). Characterization by 1H NMR and CD spectroscopies of the resulting aptamers, TBA-T7b and TBA-T12b, showed their ability to fold into the typical antiparallel chairlike G-quadruplex structure formed by TBA. The apparent CD melting temperatures indicated that the introduction of the acyclic residue, mainly at position 7, improves the thermal stability of resulting G-quadruplexes with respect to TBA. The anticoagulant activity of the new molecules was then valued in PT assay, and it resulted that TBA-T7b is more potent than TBA in prolonging clotting time. On the other hand, in purified fibrinogen assay the thrombin inhibitory activity of both modified sequences was lower than that of TBA using human enzyme, whereas the potency trend was again reversed using bovine enzyme. Obtained structure-activity relationships were investigated by structural and computational studies. Taken together, these results reveal the active role of TBA residues T7 and T12 and the relevance of some amino acids located in the anion binding exosite I of the protein in aptamer-thrombin interaction.

Cloning of a trypsin-like serine protease and expression patterns during Plasmodium falciparum invasion in the mosquito, Anopheles dirus (Peyton and Harrison)

Understanding specific gene regulation during responses to malaria infection is key to dissecting the mosquito defense mechanisms and host/parasite interactions. A full-length serine protease cDNA was isolated from the fat body of female Anopheles dirus, a major malaria vector in Thailand. The predicted amino acid sequence of SERF4 identifies it as a member of the serine protease family containing a single trypsin-like protease domain. Digestive trypsins function in the female mosquito midgut and are inducible in two phases in this tissue upon blood intake. However, the gene was highly upregulated in the midgut at day 3 postinfection by Plasmodium falciparum. In situ hybridization confirmed that SERF4 transcripts were located in the midgut epithelial cells rather than hemocytes or other tissues associated with the midgut. SERF4 was also strongly downregulated in the whole insects at day 16 after infection in comparison with the blood-fed control. Changes in the expression of the SERF4 gene in response to infection with this human malaria parasite suggest a role in parasite-specific innate immunity.

Conformational selection or induced fit? A critical appraisal of the kinetic mechanism

For almost five decades, two competing mechanisms of ligand recognition, conformational selection and induced fit, have dominated our interpretation of ligand binding in biological macromolecules. When binding-dissociation events are fast compared to conformational transitions, the rate of approach to equilibrium, k(obs), becomes diagnostic of conformational selection or induced fit based on whether it decreases or increases, respectively, with the ligand concentration, [L]. However, this simple conclusion based on the rapid equilibrium approximation is not valid in general. Here we show that conformational selection is associated with a rich repertoire of kinetic properties, with k(obs) decreasing or increasing with [L] depending on the relative magnitude of the rate of ligand dissociation, k(off), and the rate of conformational isomerization, k(r). We prove that, even for the simplest two-step mechanism of ligand binding, a decrease in k(obs) with [L] is unequivocal evidence of conformational selection, but an increase in k(obs) with [L] is not unequivocal evidence of induced fit. Ligand binding to glucokinase, thrombin, and its precursor prethrombin-2 are used as relevant examples. We conclude that conformational selection as a mechanism for a ligand binding to its target may be far more common than currently believed.

Conformational selection in trypsin-like proteases

For over four decades, two competing mechanisms of ligand recognition--conformational selection and induced-fit--have dominated our interpretation of protein allostery. Defining the mechanism broadens our understanding of the system and impacts our ability to design effective drugs and new therapeutics. Recent kinetics studies demonstrate that trypsin-like proteases exist in equilibrium between two forms: one fully accessible to substrate (E) and the other with the active site occluded (E*). Analysis of the structural database confirms existence of the E* and E forms and vouches for the allosteric nature of the trypsin fold. Allostery in terms of conformational selection establishes an important paradigm in the protease field and enables protein engineers to expand the repertoire of proteases as therapeutics.

Crystal structures of prethrombin-2 reveal alternative conformations under identical solution conditions and the mechanism of zymogen activation

Prethrombin-2 is the immediate zymogen precursor of the clotting enzyme thrombin, which is generated upon cleavage at R15 and separation of the A chain and catalytic B chain. The X-ray structure of prethrombin-2 determined in the free form at 1.9 Å resolution shows the 215-217 segment collapsed into the active site and occluding 49% of the volume available for substrate binding. Remarkably, some of the crystals harvested from the same crystallization well, under identical solution conditions, diffract to 2.2 Å resolution in the same space group but produce a structure in which the 215-217 segment moves >5 Å and occludes 24% of the volume available for substrate binding. The two alternative conformations of prethrombin-2 have the side chain of W215 relocating >9 Å within the active site and are relevant to the allosteric E*-E equilibrium of the mature enzyme. Another unanticipated feature of prethrombin-2 bears on the mechanism of prothrombin activation. R15 is found buried within the protein in ionic interactions with E14e, D14l, and E18, thereby making its exposure to solvent necessary for proteolytic attack and conversion to thrombin. On the basis of this structural observation, we constructed the E14eA/D14lA/E18A triple mutant to reduce the level of electrostatic coupling with R15 and promote zymogen activation. The mutation causes prethrombin-2 to spontaneously convert to thrombin, without the need for the snake venom ecarin or the physiological prothrombinase complex.