Ser(214) is crucial for substrate binding to serine proteases

Codon switching of conserved Ser residues in coagulation and fibrinolytic proteases

BACKGROUND

Unique among all amino acids, Ser is encoded by 2 sets of codons, TCN and AGY (N = any nucleotide, Y = pyrimidine), that cannot interconvert through single nucleotide substitutions. Both codons are documented at the essential residues S195 and S214 within the active site of serine proteases. However, it is not known how the codons interconverted during evolution because replacement of S195 or S214 by other amino acids typically results in loss of activity.

OBJECTIVE

To characterize the prevalence of codon switching among essential and non-essential Ser residues in coagulation and fibrinolytic proteases from different vertebrate lineages.

METHODS

TCN and AGY codon usage was analyzed in >550 sequences.

RESULTS

Evolutionary pressure to preserve the codon of S195 is absolute, with no evidence of interconversion. Pressure to preserve the codon of S214 is also strong, but an AGY↔TCN interconversion is observed in factor VII-inactive and protein C from ray-finned fish. In both cases, the interconversion occurred in genes that were rapidly evolving. In contrast, codon switching at nonessential Ser residues in the kringle domains of coagulation and fibrinolytic proteases is quite common and could be identified in half of the kringles analyzed.

CONCLUSION

Codon interconversion of essential Ser residues of coagulation and fibrinolytic proteases only occurred in genes that were rapidly evolving and that-at least in some cases-evolved following genome duplication. Interconversion is common at nonessential Ser residues as found in kringle domains.

Modern non-polarizable force fields diverge in modeling the enzyme-substrate complex of a canonical serine protease

Classical molecular dynamics simulation is a powerful and established method of modern computational chemistry. Being able to obtain accurate information on molecular behavior is crucial to get valuable insights into structure-function relationships that translate into fundamental findings and practical applications. Active sites of enzymes are known to be particularly intricate, therefore, simpler non-polarizable force fields may provide an inaccurate description. In this work, we addressed this hypothesis in a case of a canonical serine triad protease trypsin in its complex with a substrate-mimicking inhibitor. We tested six modern and popular force fields to find that significantly diverging results may be obtained. Amber FB-15 and OPLS-AA/M turned out to model the active site incorrectly. Amber ff19sb and ff15ipq demonstrated mixed performance. The best performing force fields were CHARMM36m and Amber ff99sb-ildn, therefore, they are recommended for use with this and related systems. We speculate that a similar lack of cross-force field convergence may be characteristic of other enzymatic systems. Therefore, we advocate for careful consideration of different force fields in any study within the field of computational enzymology.

Legumain Activity Is Controlled by Extended Active Site Residues and Substrate Conformation

Legumain is a lysosomal cysteine protease with strict specificity for cleaving after asparagine residues. By sequence comparison, legumain belongs to MEROPS clan CD of the cysteine proteases, which indicates its structural and mechanistic relation to caspases. Contrasting caspases, legumain harbors a pH-dependent ligase activity in addition to the protease activity. Although we already have a significant body of knowledge on the catalytic activities of legumain, many mechanistic details are still elusive. In this study, we provide evidence that extended active site residues and substrate conformation are steering legumain activities. Biochemical experiments and bioinformatics analysis showed that the catalytic Cys189 and His148 residues are regulated by sterically close Glu190, Ser215 and Asn42 residues. While Glu190 serves as an activity brake, Ser215 and Asn42 have a favorable effect on legumain protease activity. Mutagenesis studies using caspase-9 as model enzyme additionally showed that a similar Glu190 activity brake is also implemented in the caspases. Furthermore, we show that the substrate's conformational flexibility determines whether it will be hydrolyzed or ligated by legumain. The functional understanding of the extended active site residues and of substrate prerequisites will allow us to engineer proteases with increased enzymatic activity and better ligase substrates, with relevance for biotechnological applications.

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.

Mechanical couplings of protein backbone and side chains exhibit scale-free network properties and specific hotspots for function

•

Statistical learning from protein dynamics unravels rigidities in interaction network.

•

Backbone and side-chain mechanical couplings exhibit scale-free network properties.

•

Graphical depiction of network rigidities captures sequence co-evolution patterns.

•

Functional sites at secondary structure peripheries are mechanical hotspots.

•

Our rigidity scores are compelling metrics for residue biological significance.

A backbone-side-chain elastic network model (bsENM) is devised in this contribution to decipher the network of molecular interactions during protein dynamics. The chemical details in 5 μs all-atom molecular dynamics (MD) simulation are mapped onto the bsENM spring constants by self-consistent iterations. The elastic parameters obtained by this structure-mechanics statistical learning are then used to construct inter-residue rigidity graphs for the chemical components in protein amino acids. A key discovery is that the mechanical coupling strengths of both backbone and side chains exhibit heavy-tailed distributions and scale-free network properties. In both rat trypsin and PDZ3 proteins, the statistically prominent modes of rigidity graphs uncover the sequence-specific coupling patterns and mechanical hotspots. Based on the contributions to graphical modes, our residue rigidity scores in backbone and side chains are found to be very useful metrics for the biological significance. Most functional sites have high residue rigidity scores in side chains while the biologically important glycines are generally next to mechanical hotspots. Furthermore, prominent modes in the rigidity graphs involving side chains oftentimes coincide with the co-evolution patterns due to evolutionary restraints. The bsENM specifically devised to resolve the protein chemical character thus provides useful means for extracting functional information from all-atom MD.

A distance geometry-based description and validation of protein main-chain conformation

Understanding the protein main-chain conformational space forms the basis for the modelling of protein structures and for the validation of models derived from structural biology techniques. Presented here is a novel idea for a three-dimensional distance geometry-based metric to account for the fine details of protein backbone conformations. The metrics are computed for dipeptide units, defined as blocks of Cαi-1-O i-1-Cαi -O i -Cαi+1 atoms, by obtaining the eigenvalues of their Euclidean distance matrices. These were computed for ∼1.3 million dipeptide units collected from nonredundant good-quality structures in the Protein Data Bank and subjected to principal component analysis. The resulting new Euclidean orthogonal three-dimensional space (DipSpace) allows a probabilistic description of protein backbone geometry. The three axes of the DipSpace describe the local extension of the dipeptide unit structure, its twist and its bend. By using a higher-dimensional metric, the method is efficient for the identification of Cα atoms in an unlikely or unusual geometrical environment, and its use for both local and overall validation of protein models is demonstrated. It is also shown, for the example of trypsin proteases, that the detection of unusual conformations that are conserved among the structures of this protein family may indicate geometrically strained residues of potentially functional importance.

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.

Effects of single N-glycosylation site knockout on folding and defibrinogenating activities of acutobin recombinants from HEK293T

Acutobin, the α-fibrinogenase from Deinagkistrodon acutus venom, contains four N-glycosylation sites with disialylated complex-typed glycans. Here, we explore the functional roles of each of the N-glycan by site-directed mutagenesis. The wild-type (ATB-wt) and single glycan-knockout mutants of recombinant acutobin were prepared from HEK293T, demonstrating that mutations at Asn(77), Asn(81) and Asn(100) impaired the folding while the S79A mutant and various Asn(229)-deglycosylated mutants were correctly folded. Based on homology modeling of acutobin and multiple sequence alignment with various venom thrombin-like enzymes, the importance of a hydrophilic environment at each glycosylation site to the enzyme folding could be rationalized. Remarkably, all the mutants showed similar catalytic activities for the chromogenic substrate and similar thermal stabilities as ATB-wt, suggesting that the glycan knockout did not affect the gross conformation and stability of the active sites. Although SDS-PAGE analyses revealed that ATB-wt and the D229-mutant degraded all human fibrinogen subunits faster but less specifically in vitro as compared with other mutants that cleaved only the α-subunit, ATB-wt and D229-mutant were not able to release fibrinogen-peptide A and thus coagulated human plasma slower than the other mutants did. In the mice model, the defibrinogenating effect of ATB-wt was stronger and lasting-longer than those of all the mutants. Taken together, all the glycans contribute to the pharmacokinetics of acutobin and ATB-wt in vivo, and the microenvironment around the Asn(229)-glycan appears to regulate the fibrinogen-chain specificity of acutobin while the N-glycans at positions 77, 81 and 100 are crucial for its folding.

trans-Protease activity and structural insights into the active form of the alphavirus capsid protease

The alphavirus capsid protein (CP) is a serine protease that possesses cis-proteolytic activity essential for its release from the nascent structural polyprotein. The released CP further participates in viral genome encapsidation and nucleocapsid core formation, followed by its attachment to glycoproteins and virus budding. Thus, protease activity of the alphavirus capsid is a potential antialphaviral target to arrest capsid release, maturation, and structural polyprotein processing. However, the discovery of capsid protease inhibitors has been hampered due to the lack of a suitable screening assay and of the crystal structure in its active form. Here, we report the development of a trans-proteolytic activity assay for Aura virus capsid protease (AVCP) based on fluorescence resonance energy transfer (FRET) for screening protease inhibitors. Kinetic parameters using fluorogenic peptide substrates were estimated, and the K(m) value was found to be 2.63 ± 0.62 μM while the k(cat)/K(m) value was 4.97 × 10(4) M(-1) min(-1). Also, the crystal structure of the trans-active form of AVCP has been determined to 1.81-Å resolution. Structural comparisons of the active form with the crystal structures of available substrate-bound mutant and inactive blocked forms of the capsid protease identify conformational changes in the active site, the oxyanion hole, and the substrate specificity pocket residues, which could be critical for rational drug design. IMPORTANCE The alphavirus capsid protease is an attractive antiviral therapeutic target. In this study, we have described the formerly unappreciated trans-proteolytic activity of the enzyme and for the first time have developed a FRET-based protease assay for screening capsid protease inhibitors. Our structural studies unveil the structural features of the trans-active protease, which has been previously proposed to exist in the natively unfolded form (M. Morillas, H. Eberl, F. H. Allain, R. Glockshuber, and E. Kuennemann, J. Mol. Biol. 376:721-735, 2008, doi:http://dx.doi.org/10.1016/j.jmb.2007.11.055). The different enzymatic forms have been structurally compared to reveal conformational variations in the active and substrate binding sites. The flexible active-site residue Ser218, the disordered C-terminal residues after His261, and the presence of a water molecule in the oxyanion hole of AVCPΔ2 (AVCP with a deletion of the last two residues at the C terminus) reveal the effect of the C-terminal Trp267 deletion on enzyme structure. New structural data reported in this study along with the fluorogenic assay will be useful in substrate specificity characterization, high-throughput protease inhibitor screening, and structure-based development of antiviral drugs.

Structure of a proteolytically resistant analogue of (NLys)5SFTI-1 in complex with trypsin: evidence for the direct participation of the Ser214 carbonyl group in serine protease-mediated proteolysis

Peptide-peptoid hybrids are found to be potent inhibitors of serine proteases. These engineered peptidomimetics benefit from both types of units of the biopolymeric structure: the natural inhibitor part serves as a good binding template, while the P1-positioned peptoid component provides complete resistance towards proteolysis. In this report, the mechanism of proteolytic resistance of a P1 peptoid-containing analogue is postulated based on the crystal structure of the (NLys)(5)-modified sunflower trypsin inhibitor SFTI-1 in complex with bovine trypsin solved at 1.29 Å resolution. The structural differences between the (NLys)(5)SFTI-1-trypsin complex and the native SFTI-1-trypsin complex are surprisingly small and reveal the key role of the carbonyl group of the Ser214 residue of the enzyme, which is crucial for binding of the inhibitor and plays a crucial role in proteolysis mediated by serine proteases. The incorporated NLys5 peptoid residue prevents Ser214 from forming a hydrogen bond to the P1 residue, and in turn Gln192 does not form a hydrogen bond to the carbonyl group of the P2 residue. It also increases the distance between the Ser214 carbonyl group and the Ser195 residue, thus preventing proteolysis. The hybrid inhibitor structure reported here provides insight into protein-protein interaction, which can be efficiently and selectively probed with the use of peptoids incorporated within endogenous peptide ligands.

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.

The role of factor XIa (FXIa) catalytic domain exosite residues in substrate catalysis and inhibition by the Kunitz protease inhibitor domain of protease nexin 2

To select residues in coagulation factor XIa (FXIa) potentially important for substrate and inhibitor interactions, we examined the crystal structure of the complex between the catalytic domain of FXIa and the Kunitz protease inhibitor (KPI) domain of a physiologically relevant FXIa inhibitor, protease nexin 2 (PN2). Six FXIa catalytic domain residues (Glu(98), Tyr(143), Ile(151), Arg(3704), Lys(192), and Tyr(5901)) were subjected to mutational analysis to investigate the molecular interactions between FXIa and the small synthetic substrate (S-2366), the macromolecular substrate (factor IX (FIX)) and inhibitor PN2KPI. Analysis of all six Ala mutants demonstrated normal K(m) values for S-2366 hydrolysis, indicating normal substrate binding compared with plasma FXIa; however, all except E98A and K192A had impaired values of k(cat) for S-2366 hydrolysis. All six Ala mutants displayed deficient k(cat) values for FIX hydrolysis, and all were inhibited by PN2KPI with normal values of K(i) except for K192A, and Y5901A, which displayed increased values of K(i). The integrity of the S1 binding site residue, Asp(189), utilizing p-aminobenzamidine, was intact for all FXIa mutants. Thus, whereas all six residues are essential for catalysis of the macromolecular substrate (FIX), only four (Tyr(143), Ile(151), Arg(3704), and Tyr(5901)) are important for S-2366 hydrolysis; Glu(98) and Lys(192) are essential for FIX but not S-2366 hydrolysis; and Lys(192) and Tyr(5901) are required for both inhibitor and macromolecular substrate interactions.

Molecular dissection of Streptomyces trypsin on substrate recognition

We recently identified residue 71 of two homologous serine proteases from Streptomyces omiyaensis (SOT) and Streptomyces griseus (SGT) as a crucial residue for differences in their topological specificities, i.e. recognition of a distinct three-dimensional structure. To study the role of this key residue in substrate recognition, we used surface plasmon resonance analysis to evaluate the affinities of inactive mutants, in which residues 71 of SOT and SGT were substituted respectively with Leu and Tyr, toward different types of collagens. We identified another amino acid residue involved in the interaction with collagens from analyses of inactive chimeras between SOT and SGT using an in vivo DNA shuffling system. Results showed that residue 72 contributes to collagen binding. By substituting Leu71 and Gln72 with Tyr and Arg, respectively, SGT mutant showed a change in topological specificity and high hydrolytic activity toward type IV collagen comparable to SOT. We demonstrated that the neighboring residues 71 and 72 in the N-terminal β-barrel domain of the enzyme synergistically play an important role in substrate recognition.

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.

Engineering thrombin for selective specificity toward protein C and PAR1

Thrombin elicits functional responses critical to blood homeostasis by interacting with diverse physiological substrates. Ala-scanning mutagenesis of 97 residues covering 53% of the solvent accessible surface area of the enzyme identifies Trp(215) as the single most important determinant of thrombin specificity. Saturation mutagenesis of Trp(215) produces constructs featuring k(cat)/K(m) values for the hydrolysis of fibrinogen, protease-activated receptor PAR1, and protein C that span five orders of magnitude. Importantly, the effect of Trp(215) replacement is context dependent. Mutant W215E is 10-fold more specific for protein C than fibrinogen and PAR1, which represents a striking shift in specificity relative to wild-type that is 100-fold more specific for fibrinogen and PAR1 than protein C. However, when the W215E mutation is combined with deletion of nine residues in the autolysis loop, which by itself shifts the specificity of the enzyme from fibrinogen and PAR1 to protein C, the resulting construct features significant activity only toward PAR1. These findings demonstrate that thrombin can be re-engineered for selective specificity toward protein C and PAR1. Mutations of Trp(215) provide important reagents for dissecting the multiple functional roles of thrombin in the blood and for clinical applications.

Active site prediction using evolutionary and structural information

Motivation: The identification of catalytic residues is a key step in understanding the function of enzymes. While a variety of computational methods have been developed for this task, accuracies have remained fairly low. The best existing method exploits information from sequence and structure to achieve a precision (the fraction of predicted catalytic residues that are catalytic) of 18.5% at a corresponding recall (the fraction of catalytic residues identified) of 57% on a standard benchmark. Here we present a new method, Discern, which provides a significant improvement over the state-of-the-art through the use of statistical techniques to derive a model with a small set of features that are jointly predictive of enzyme active sites.

Results: In cross-validation experiments on two benchmark datasets from the Catalytic Site Atlas and CATRES resources containing a total of 437 manually curated enzymes spanning 487 SCOP families, Discern increases catalytic site recall between 12% and 20% over methods that combine information from both sequence and structure, and by ≥50% over methods that make use of sequence conservation signal only. Controlled experiments show that Discern's improvement in catalytic residue prediction is derived from the combination of three ingredients: the use of the INTREPID phylogenomic method to extract conservation information; the use of 3D structure data, including features computed for residues that are proximal in the structure; and a statistical regularization procedure to prevent overfitting.

Contact:kimmen@berkeley.edu

Supplementary information:Supplementary data are available at Bioinformatics online.

Prostate-specific antigen is a "chymotrypsin-like" serine protease with unique P1 substrate specificity

Prostate-specific antigen (PSA), a serine protease belonging to the human kallikrein family, is best known as a prostate cancer biomarker. Emerging evidence suggests that PSA may also play a salient role in prostate cancer development and progression. With large amounts of enzymatically active PSA continuously and selectively produced by all stages of prostate cancer, PSA is an attractive target. PSA inhibitors, therefore, may represent a promising class of therapeutics and/or imaging agents. PSA displays chymotrypsin-like specificity, cleaving after hydrophobic residues, in addition to possessing a unique ability to cleave after glutamine in the P1 position. In this study, we investigated the structural motifs of the PSA S1 pocket that give it a distinct architecture and specificity when compared to the S1 pocket of chymotrypsin. Using the previously described PSA substrate Ser-Ser-Lys-Leu-Gln (SSKLQ) as a template, peptide aldehyde based inhibitors containing novel P1 aldehydes were made and tested against both proteases. Glutamine derivative aldehydes were highly specific for PSA while inhibitors with hydrophobic P1 aldehydes were potent inhibitors of both proteases with K(i) values <500 nM. The crystal structure of PSA was used to generate a model that allowed GOLD docking studies to be performed to further understand the critical interactions required for inhibitor binding to the S1 pockets of PSA and chymotrypsin. In conclusion, these results provide experimental and structural evidence that the S1 specificity pocket of PSA is distinctly different from that of chymotrypsin and that the development of highly specific PSA inhibitors is feasible.

A one-pot, simple methodology for cassette randomisation and recombination for focused directed evolution

Protein engineering is currently performed either by rational design, focusing in most cases on only a few positions modified by site-directed mutagenesis, or by directed molecular evolution, in which the entire protein-encoding gene is subjected to random mutagenesis followed by screening or selection of desired phenotypes. A novel alternative is focused directed evolution, in which only fragments of a protein are randomised while the overall scaffold of a protein remains unchanged. For this purpose, we developed a PCR technique using long, spiked oligonucleotides, which allow randomising of one or several cassettes in any given position of a gene. This method allows over 95% incorporation of mutations independently of their position within the gene, yielding sufficient product to generate large libraries, and the possibility of simultaneously randomising more than one locus at a time, thus originating recombination. The high efficiency of this method was verified by creating focused mutant libraries of Pseudomonas fluorescens esterase I (PFEI), screening for altered substrate selectivity and validating against libraries created by error-prone PCR. This led to the identification of two mutants within the OSCARR library with a 10-fold higher catalytic efficiency towards p-nitrophenyl dodecanoate. These PFEI variants were also modelled in order to explain the observed effects.

Role of the A chain in thrombin function

The A chain of thrombin is covalently linked to the catalytic B chain but is separate from any known epitope for substrate recognition. In this study we present the results of the Ala replacement of 12 charged residues controlling the stability of the A chain and its interaction with the B chain. Residues Arg4 and Glu8 play a significant role in substrate recognition, even though they are located > 20 A away from residues of the catalytic triad, the primary specificity pocket and the Na+ site. The R4A mutation causes significant perturbation of Na+ binding, fibrinogen clotting and PAR1 cleavage, but modest reduction of protein C activation in the presence of thrombomodulin. These findings challenge our current paradigm of thrombin structure-function relations focused exclusively on the properties of the catalytic B chain, and explain why certain naturally occurring mutations of the A chain cause serious bleeding.

Structure of the glucanase inhibitor protein (GIP) family from phytophthora species suggests coevolution with plant endo-beta-1,3-glucanases

During invasion of their plant hosts, species of the oomycete genus Phytophthora secrete glucanase inhibitor proteins (GIPs) into the plant apoplast, which bind and inhibit the activity of plant extracellular endo-beta-1,3-glucanases (EGases). GIPs show structural homology to the chymotrypsin class of serine proteases (SP) but lack proteolytic activity due to the absence of an intact catalytic triad and, thus, belong to a broader class of proteins called serine protease homologs (SPH). To study the evolutionary relationship between GIPs and functional SP, database searches were used to identify 48 GIP homologs in the P. sojae, P. ramorum, and P. infestans genomes, composing GIPs, SPH, and potentially functional SP. Analyses of P. infestans-inoculated tomato leaves showed that P. infestans GIPs and tomato EGases are present in the apoplast and form stable complexes in planta. Studies of the temporal expression of a four-membered GIP family from P. infestans (PiGIP1 to PiGIP4) further revealed that the genes show distinctly different patterns during an infection timecourse. Codon evolution analyses of GIP homologs identified several positively selected peptide sites and structural modeling revealed them to be in close proximity to rapidly evolving EGase residues, suggesting that the interaction between GIPs and EGases has the hallmarks of a coevolving molecular arms race.

Thrombin

Thrombin is a Na+-activated, allosteric serine protease that plays opposing functional roles in blood coagulation. Binding of Na+ is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme, but is dispensable for cleavage of the anticoagulant protein C. The anticoagulant function of thrombin is under the allosteric control of the cofactor thrombomodulin. Much has been learned on the mechanism of Na+ binding and recognition of natural substrates by thrombin. Recent structural advances have shed light on the remarkable molecular plasticity of this enzyme and the molecular underpinnings of thrombin allostery mediated by binding to exosite I and the Na+ site. This review summarizes our current understanding of the molecular basis of thrombin function and allosteric regulation. The basic information emerging from recent structural, mutagenesis and kinetic investigation of this important enzyme is that thrombin exists in three forms, E*, E and E:Na+, that interconvert under the influence of ligand binding to distinct domains. The transition between the Na+ -free slow from E and the Na+ -bound fast form E:Na+ involves the structure of the enzyme as a whole, and so does the interconversion between the two Na+ -free forms E* and E. E* is most likely an inactive form of thrombin, unable to interact with Na + and substrate. The complexity of thrombin function and regulation has gained this enzyme pre-eminence as the prototypic allosteric serine protease. Thrombin is now looked upon as a model system for the quantitative analysis of biologically important enzymes.

Dielectric studies of water clusters in cyclodextrins: Relevance to the transition between slow and fast forms of thrombin

Cyclodextrins are useful models in the study of hydrogen bonded water clusters. In alpha-cyclodextrin hexahydrate (alpha-CD.6H2O), water molecules are ordered and occupy well-defined positions whereas in the larger beta-cyclodextrin dodecahydrate (beta-CD.12H2O), there is considerable disorder with water molecules freely arranged over several possible sites. Here it is shown that beta-CD exhibits substantial structural flexibility and proton mobility compared with alpha-CD which is relatively very rigid and exhibits negligible short-range protonic conduction. These properties are directly controlled by the effective dielectric constant of the molecule, which is determined by the rotational freedom of water molecules in the hydrogen bond network. This model may be relevant to proteins where water clusters of this kind are found on the protein surface and occasionally in the protein interior. The case of thrombin, an allosteric enzyme incorporating a network of 20 internal hydrogen bonded water molecules, is discussed.

Thrombin allostery

Thrombin is a Na(+)-activated, allosteric serine protease that plays multiple functional roles in blood pathophysiology. Binding of Na(+) is the major driving force behind the procoagulant, prothrombotic and signaling functions of the enzyme. This review summarizes our current understanding of the molecular basis of thrombin allostery with special emphasis on the kinetic aspects of Na(+) activation. The molecular mechanism of thrombin allostery is a remarkable example of long-range communication that offers a paradigm for many other biological systems.

Rapid kinetics of Na+ binding to thrombin

The kinetic mechanism of Na(+) binding to thrombin was resolved by stopped-flow measurements of intrinsic fluorescence. Na(+) binds to thrombin in a two-step mechanism with a rapid phase occurring within the dead time of the spectrometer (<0.5 ms) followed by a single-exponential slow phase whose k(obs) decreases hyperbolically with increasing [Na(+)]. The rapid phase is due to Na(+) binding to the enzyme E to generate the E:Na(+) form. The slow phase is due to the interconversion between E(*) and E, where E(*) is a form that cannot bind Na(+). Temperature studies in the range from 5 to 35 degrees C show significant enthalpy, entropy, and heat capacity changes associated with both Na(+) binding and the E to E(*) transition. As a result, under conditions of physiologic temperature and salt concentrations, the E(*) form is negligibly populated (<1%) and thrombin is almost equally partitioned between the E (40%) and E:Na(+) (60%) forms. Single-site Phe mutations of all nine Trp residues of thrombin enabled assignment of the fluorescence changes induced by Na(+) binding mainly to Trp-141 and Trp-215, and to a lesser extent to Trp-148, Trp-207, and Trp-237. However, the fast phase of fluorescence increase is influenced to different extents by all Trp residues. The distribution of these residues over the entire thrombin surface demonstrates that Na(+) binding induces long-range effects on the structure of the enzyme as a whole, contrary to the conclusions drawn from recent structural studies. These findings elucidate the mechanism of Na(+) binding to thrombin and are relevant to other clotting factors and enzymes allosterically activated by monovalent cations.

Role of Na+and K+in Enzyme Function

Metal complexation is a key mediator or modifier of enzyme structure and function. In addition to divalent and polyvalent metals, group IA metals Na+and K+play important and specific roles that assist function of biological macromolecules. We examine the diversity of monovalent cation (M+)-activated enzymes by first comparing coordination in small molecules followed by a discussion of theoretical and practical aspects. Select examples of enzymes that utilize M+as a cofactor (type I) or allosteric effector (type II) illustrate the structural basis of activation by Na+and K+, along with unexpected connections with ion transporters. Kinetic expressions are derived for the analysis of type I and type II activation. In conclusion, we address evolutionary implications of Na+binding in the trypsin-like proteases of vertebrate blood coagulation. From this analysis, M+complexation has the potential to be an efficient regulator of enzyme catalysis and stability and offers novel strategies for protein engineering to improve enzyme function.

Crystal structure of thrombin in a self-inhibited conformation

The activating effect of Na(+) on thrombin is allosteric and depends on the conformational transition from a low activity Na(+)-free (slow) form to a high activity Na(+)-bound (fast) form. The structures of these active forms have been solved. Recent structures of thrombin obtained in the absence of Na(+) have also documented inactive conformations that presumably exist in equilibrium with the active slow form. The validity of these inactive slow form structures, however, is called into question by the presence of packing interactions involving the Na(+) site and the active site regions. Here, we report a 1.87A resolution structure of thrombin in the absence of inhibitors and salts with a single molecule in the asymmetric unit and devoid of significant packing interactions in regions involved in the allosteric slow --> fast transition. The structure shows an unprecedented self-inhibited conformation where Trp-215 and Arg-221a relocate >10A to occlude the active site and the primary specificity pocket, and the guanidinium group of Arg-187 penetrates the protein core to fill the empty Na(+)-binding site. The extreme mobility of Trp-215 was investigated further with the W215P mutation. Remarkably, the mutation significantly compromises cleavage of the anticoagulant protein C but has no effect on the hydrolysis of fibrinogen and PAR1. These findings demonstrate that thrombin may assume an inactive conformation in the absence of Na(+) and that its procoagulant and anticoagulant activities are closely linked to the mobility of residue 215.

Discrimination of Esterase and Peptidase Activities of Acylaminoacyl Peptidase from Hyperthermophilic Aeropyrum pernix K1 by a Single Mutation

It has been shown that highly conserved residues that form crucial structural elements of the catalytic apparatus may be used to account for the evolutionary history of enzymes. Using saturation mutagenesis, we investigated the role of a conserved residue (Arg(526)) at the active site of acylaminoacyl peptidase from hyperthermophilic Aeropyrum pernix K1 in substrate discrimination and catalytic mechanism. This enzyme has both peptidase and esterase activities. The esterase activity of the wild-type enzyme with p-nitrophenyl caprylate as substrate is approximately 7 times higher than the peptidase activity with Ac-Leu-p-nitroanilide as substrate. However, with the same substrates, this difference was increased to approximately 150-fold for mutant R526V. A more dramatic effect occurred with mutant R526E, which essentially completely abolished the peptidase activity but decreased the esterase activity only by a factor of 2, leading to a 785-fold difference in the enzyme activities. These results provide rare examples that illustrate how enzymes can be evolved to discriminate their substrates by a single mutation. The possible structural and energetic effects of the mutations on k(cat) and K(m) of the enzyme were discussed based on molecular dynamics simulation studies.

High resolution crystal structures of free thrombin in the presence of K+ reveal the molecular basis of monovalent cation selectivity and an inactive slow form

Structural biology has recently advanced our understanding of the molecular mechanisms of activation and selectivity in monovalent cation activated enzymes. Here we report a 1.9 Angstrom resolution crystal structure of free thrombin, a Na(+) selective enzyme, in the presence of KCl. There are two molecules in the asymmetric unit, one with the cation site bound to K(+) and the other with this site free. The K(+)-bound form shows key differences compared with the Na(+)-bound structure that explain the different kinetics of activation. The cation-free form, on the other hand, assumes a conformation where the monovalent cation binding site is completely disordered, the S1 pocket is inaccessible to substrate and binding to exosite I is compromised by an unprecedented >20 Angstrom shift in the position of the autolysis loop. This form, named S(*), corresponds to the inactive Na(+)-free slow form identified by early kinetic studies. A simple model of thrombin allostery that incorporates the contribution of S(*) is proposed.

Functional requirements for the optimal catalytic configuration of the AChE active center

Functional analysis of the HuAChE active center architecture revealed that accommodation of structurally diverse substrates and other ligands is achieved through interactions with specific subsites such as the acyl pocket, cation binding site, hydrophobic site or the oxyanion hole. Recent studies have begun to unravel the role of this active center architecture in maintaining the optimal catalytic facility of the enzyme through inducing proper alignment of the catalytic triad. The exact positioning of the catalytic glutamate (Glu334) seems to be determined by a hydrogen bond network including several polar residues and water molecules. Disruption of this network by replacement of Ser229 by alanine is thought to remove the Glu334 carboxylate from the vicinity of His447 abolishing catalytic activity. The proper orientation of the catalytic histidine side chain is maintained by these polar interactions as well as through "aromatic trapping" by residues lining the HuAChE active center gorge. Thus, replacement of aromatic residues in the vicinity of His447, as in the F295A/F338A or in the Y72N/Y124Q/W286A/F295L/F297V/Y337A (hexamutant which mimicks the aromatic lining of HuBChE) enzymes, resulted in a dramatic decrease in catalytic activity, which was proposed to originate from catalytically nonproductive mobility of His447. Yet, HuBChE is catalytically efficient indicating that "aromatic trapping" is not the only way to conformationally stabilize the His447 side chain. A possible restriction of this mobility in a series of F295X/F338A HuAChEs was examined in silico followed by site-directed mutagenesis. Both simulations and reactivities of the actual F295X/F338A enzymes, carrying various aliphatic residues at position 295, indicate that of the bulky amino acids, like leucine or isoleucine, only methionine was capable of maintaining the catalytically viable conformation of His447. The F295M/F338A HuAChE was only two-fold less reactive than the F338A enzyme toward acetylthiocholine, and exhibited wild type-like reactivity toward covalent modifiers of the catalytic Ser203. The findings are consistent with the notion that different combinations of steric interference and specific polar interactions serve to maintain the position of His447 and thereby the high efficiency of the catalytic machinery. The two seemingly conflicting demands on the architecture of the active center-flexible accommodation of substrate and optimal juxtaposition of residues of the catalytic triad, demonstrate the truly amazing molecular design of the AChE active center.

Mannan-binding-lectin-associated serine proteases, characteristics and disease associations

Mannan-binding lectin (MBL)-associated serine proteases (MASPs) circulate in plasma as zymogens in complexes with MBL and with L- and H-ficolin. Upon binding of MBL or ficolin to pathogen-associated molecular patterns, the MASPs are activated. MASP-2 can now cleave C4 and C2 to generate the C3 convertase, C4bC2b. The functions of the other two MASPs, MASP-1 and MASP-3 have not been elucidated. MASP-1 can cleave C2, and with low efficiency also C3, and may serve a function through direct C3 activation. No natural substrate for MASP-3 has been identified. MBL deficiency, occurring at a frequency of about 10%, is the most common congenital immunodeficiency and is associated with susceptibility to infections and autoimmune disorders. Inherited MASP-2 deficiency has been described as the result of a mutation causing the exchange of aspartic acid with a glycine at position 105, a position in the first domain, CUB1, involved in calcium binding. This mutation abolishes the binding to MBL and ficolins, and deprives MASP-2 of functional activity. The index case suffered from recurrent severe infections and autoimmune reactions. The gene frequency of the mutation among Caucasians is 3.6%. It is not found in Chinese, who present a different mutation also associated with MASP-2 deficiency.

Thrombin: a paradigm for enzymes allosterically activated by monovalent cations

Enzymes activated by monovalent cations are abundantly represented in plants and in the animal world. The mechanism, of activation involves formation of a ternary intermediate with the enzyme-substrate complex, or binding of the cation to an allosteric site in the protein. Thrombin is a Na+-activated enzyme with procoagulant, anticoagulant and signaling roles. The binding of Na+ influences allosterically thrombin function and offers a paradigm for regulatory control of protease activity and specificity. Here we review the molecular basis of thrombin allostery as recently emerged from mutagenesis and structural studies. The role of Na+ in blood coagulation and the evolution of serine proteases are also discussed.

Hirudin Binding Reveals Key Determinants of Thrombin Allostery

Thrombin exists in two allosteric forms, slow (S) and fast (F), that recognize natural substrates and inhibitors with significantly different affinities. Because under physiologic conditions the two forms are almost equally populated, investigation of thrombin function must address the contribution from the S and F forms and the molecular origin of their differential recognition of ligands. Using a panel of 79 Ala mutants, we have mapped for the first time the epitopes of thrombin recognizing a macromolecular ligand, hirudin, in the S and F forms. Hirudin binding is a relevant model for the interaction of thrombin with fibrinogen and PAR1 and is likewise influenced by the allosteric S-->F transition. The epitopes are nearly identical and encompass two hot spots, one in exosite I and the other in the Na+ site at the opposite end of the protein. The higher affinity of the F form is due to the preferential interaction of hirudin with Lys-36, Leu-65, Thr-74, and Arg-75 in exosite I; Gly-193 in the oxyanion hole; and Asp-221 and Asp-222 in the Na+ site. Remarkably, no correlation is found between the energetic and structural involvements of thrombin residues in hirudin recognition, which invites caution in the analysis of protein-protein interactions in general.

Energetic and structural consequences of perturbing Gly-193 in the oxyanion hole of serine proteases

The oxyanion hole of serine proteases is formed by the backbone N atoms of the catalytic Ser-195 and Gly-193 and engages the backbone O atom of the P1 residue of substrate in an important H-bonding interaction. The energetic contribution of this interaction in the ground and transition states is presently unknown. Measurements of the individual rate constants defining the catalytic mechanism of substrate hydrolysis for wild-type thrombin and trypsin and their G193A and G193P mutants reveal that Gly-193 is required for optimal substrate binding and acylation. Crystal structures of the G193A and G193P mutants of thrombin bound to the active site inhibitor H-d-Phe-Pro-Arg-CH2Cl document the extent of perturbation induced by the replacement of Gly-193. The Ala mutant weakens the H-bonding interaction of the N atom of residue 193, whereas the Pro substitution abrogates it altogether with additional small shifts of the protein backbone. From the kinetic and structural data, we estimate that the H-bonding interaction in the oxyanion hole contributes a stabilization of the ground and transition states of > 1.5 kcal/mol but < 3.0 kcal/mol. These results shed light on a basic aspect of the enzyme-substrate interaction in the entire family of trypsin-like serine proteases.

Exosite-driven substrate specificity and function in coagulation

Macromolecular substrate recognition and serine proteinase specificity lie at the heart of the tightly regulated hemostatic response. Mechanisms established for the less specific serine proteinases of digestion have played a dominant role in guiding investigations of the basis for the narrow specificities exhibited by the coagulation enzymes. These concepts have also dominated the development of specific inhibitors of coagulation for therapeutic purposes. Studies of the enzymology and physical biochemistry of prothrombinase challenge these prevailing ideas by establishing a principal role for exosites within the enzyme in determining substrate recognition and directing the action of the enzyme on its biological substrate. Mechanisms by which narrow protein substrate specificity is achieved by prothrombinase also apply to several other reactions of coagulation. These strategies are increasingly evident in the action of other families of enzymes that act with high specificity on protein substrates. Exosite-driven enzymic function probably represents a widely employed biological strategy for the achievement of high macromolecular substrate specificity.

Thrombomodulin changes the molecular surface of interaction and the rate of complex formation between thrombin and protein C

The interaction of thrombin with protein C triggers a key down-regulatory process of the coagulation cascade. Using a panel of 77 Ala mutants, we have mapped the epitope of thrombin recognizing protein C in the absence or presence of the cofactor thrombomodulin. Residues around the Na(+) site (Thr-172, Lys-224, Tyr-225, and Gly-226), the aryl binding site (Tyr-60a), the primary specificity pocket (Asp-189), and the oxyanion hole (Gly-193) hold most of the favorable contributions to protein C recognition by thrombin, whereas a patch of residues in the 30-loop (Arg-35 and Pro-37) and 60-loop (Phe-60h) regions produces unfavorable contributions to binding. The shape of the epitope changes drastically in the presence of thrombomodulin. The unfavorable contributions to binding disappear and the number of residues promoting the thrombin-protein C interaction is reduced to Tyr-60a and Asp-189. Kinetic studies of protein C activation as a function of temperature reveal that thrombomodulin increases >1,000-fold the rate of diffusion of protein C into the thrombin active site and lowers the activation barrier for this process by 4 kcal/mol. We propose that the mechanism of thrombomodulin action is to kinetically facilitate the productive encounter of thrombin and protein C and to allosterically change the conformation of the activation peptide of protein C for optimal presentation to the thrombin active site.

Investigation of the role of a second conserved serine in carboxylesterases via site-directed mutagenesis

Carboxylesterases are enzymes that catalyze the hydrolysis of ester and amide moieties. These enzymes have an active site that is composed of a nucleophile (Ser), a base (His), and an acid (Glu) that is commonly known as a catalytic triad. It has previously been observed that the majority of carboxylesterases and lipases contain a second conserved serine in their active site [Proteins, 34 (1999) 184]. To investigate whether this second serine is also involved in the catalytic mechanism, it was mutated to an alanine, a glycine or a cysteine. Site-directed mutagenesis of this conserved serine resulted in a loss of specific activity, in both the S247G and S247A mutants (5- to 15-fold), which was due to a decrease in the rate of catalysis (kcat). Due to the instability of the S247C mutant no reliable data could be attained. A carbamate inhibitor, carbaryl, was then employed to investigate whether this decrease in the kcat was due to the rate of formation of the acyl-enzyme intermediate (k2) or the rate of deacylation (k3). The S247A mutant was found only to alter k2 (2.5-fold decrease), with no effect on k3. Together with information inferred from a human carboxylesterase crystal structure, it was concluded that this serine provides an important structural support for the spatial orientation of the glutamic acid, stabilizing the catalytic triad so that it can perform the hydrolysis.

Molecular dissection of Na+ binding to thrombin

Na(+) binding near the primary specificity pocket of thrombin promotes the procoagulant, prothrombotic, and signaling functions of the enzyme. The effect is mediated allosterically by a communication between the Na(+) site and regions involved in substrate recognition. Using a panel of 78 Ala mutants of thrombin, we have mapped the allosteric core of residues that are energetically linked to Na(+) binding. These residues are Asp-189, Glu-217, Asp-222, and Tyr-225, all in close proximity to the bound Na(+). Among these residues, Asp-189 shares with Asp-221 the important function of transducing Na(+) binding into enhanced catalytic activity. None of the residues of exosite I, exosite II, or the 60-loop plays a significant role in Na(+) binding and allosteric transduction. X-ray crystal structures of the Na(+)-free (slow) and Na(+)-bound (fast) forms of thrombin, free or bound to the active site inhibitor H-d-Phe-Pro-Arg-chloromethyl-ketone, document the conformational changes induced by Na(+) binding. The slow --> fast transition results in formation of the Arg-187:Asp-222 ion pair, optimal orientation of Asp-189 and Ser-195 for substrate binding, and a significant shift of the side chain of Glu-192 linked to a rearrangement of the network of water molecules that connect the bound Na(+) to Ser-195 in the active site. The changes in the water network and the allosteric core explain the thermodynamic signatures linked to Na(+) binding and the mechanism of thrombin activation by Na(+). The role of the water network uncovered in this study establishes a new paradigm for the allosteric regulation of thrombin and other Na(+)-activated enzymes involved in blood coagulation and the immune response.

Crystal structure of the anticoagulant slow form of thrombin

Using the thrombin mutant R77aA devoid of the site of autoproteolytic degradation at exosite I, we have solved for the first time the structure of thrombin free of any inhibitors and effector molecules and stabilized in the Na(+)-free slow form. The slow form shows subtle differences compared with the currently available structures of the Na(+)-bound fast form that carry inhibitors at the active site or exosite I. The most notable differences are the displacement of Asp-189 in the S1 specificity pocket, a downward shift of the 190-193 strand, a rearrangement of the side chain of Glu-192, and a significant shift in the position of the catalytic Ser-195 that is no longer within H-bonding distance from His-57. The structure of the slow form explains the reduced specificity toward synthetic and natural substrates and suggests a molecular basis for its anticoagulant properties.