Lysine residue 114 in human antithrombin III is required for heparin pentasaccharide-mediated activation

Estimating glycosaminoglycan-protein interaction affinity: water dominates the specific antithrombin-heparin interaction

Glycosaminoglycan (GAG)-protein interactions modulate many important biological processes. Structure-function studies on GAGs may reveal probes and drugs, but their structural complexity and highly acidic nature confound such work. Productivity will increase if we are able to identify tight-binding oligosaccharides in silico. An extension of the CHARMM force field is presented to enable modeling of polysaccharides containing sulfamate functionality, and is used to develop a reliable alchemical free-energy perturbation protocol that estimates changes in affinity for the prototypical heparin-antithrombin system to within 2.3 kcal/mol using modest simulation times. Inclusion of water is crucial during simulation as solvation energy was equal in magnitude to the sum of all other thermodynamic factors. In summary, we have identified and optimized a reliable method for estimation of GAG-protein binding affinity, and shown that solvation is a crucial component in GAG-protein interactions.

Insights into the induced fit mechanism in antithrombin–heparin interaction using molecular dynamics simulations

Heparin was isolated in the beginning of the 20th century and until today remains as one of the most important drugs able to interfere with the haemostatic process. Due to the side effects produced by heparin therapy, new promising drugs have been developed, as the synthetic pentasaccharide (synthetically derived from the sequence GlcN-GlcA-GlcN-IdoA-GlcN). The anticoagulant activity of this compound is based on potentiation of antithrombin (AT) inhibitory activity upon serine proteinases of clotting cascade, a mechanism based on the conformational modification of AT. In this context, we present here a molecular dynamics (MD) study of the interaction between the synthetic pentasaccharide and AT. The obtained data correctly predicted an induced fit mechanism in AT-pentasaccharide interaction, showing a solvent-exposed P1 residue instead of a hided conformation. Also, the specific contribution of important amino acid residues to the overall process was also characterized, both in (2)S(0) and (1)C(4) conformations of IdoA residue, suggesting that there is no conformational requirement to the interaction of this residue with AT. Altogether, the results show that MD simulations could be used to characterize and quantify the interaction of synthetic compounds with AT, predicting its specific capacity to induce conformational changes in AT structure. Thus, MD simulations of heparin (and heparin-derived)-AT interactions are proposed here as a powerful tool to assist and support drug design of new antithrombotic agents.

Antithrombin III phenylalanines 122 and 121 contribute to its high affinity for heparin and its conformational activation

The dissociation equilibrium constant for heparin binding to antithrombin III (ATIII) is a measure of the cofactor's binding to and activation of the proteinase inhibitor, and its salt dependence indicates that ionic and non-ionic interactions contribute approximately 40 and approximately 60% of the binding free energy, respectively. We now report that phenylalanines 121 and 122 (Phe-121 and Phe-122) together contribute 43% of the total binding free energy and 77% of the energy of non-ionic binding interactions. The large contribution of these hydrophobic residues to the binding energy is mediated not by direct interactions with heparin, but indirectly, through contacts between their phenyl rings and the non-polar stems of positively charged heparin binding residues, whose terminal amino and guanidinium groups are thereby organized to form extensive and specific ionic and non-ionic contacts with the pentasaccharide. Investigation of the kinetics of heparin binding demonstrated that Phe-122 is critical for promoting a normal rate of conformational change and stabilizing AT*H, the high affinity-activated binary complex. Kinetic and structural considerations suggest that Phe-122 and Lys-114 act cooperatively through non-ionic interactions to promote P-helix formation and ATIII binding to the pentasaccharide. In summary, although hydrophobic residues Phe-122 and Phe-121 make minimal contact with the pentasaccharide, they play a critical role in heparin binding and activation of antithrombin by coordinating the P-helix-mediated conformational change and organizing an extensive network of ionic and non-ionic interactions between positively charged heparin binding site residues and the cofactor.

Contribution of basic residues of the 70-80-loop to heparin binding and anticoagulant function of activated protein C

The role of basic residues of the 70-80-loop, Arg(74), Arg(75), and Lys(78) (chymotrypsin numbering) in the catalytic function of activated protein C (APC) was investigated by expressing mutants of protein C in which these residues were replaced with Ala in three separate constructs. Following purification to homogeneity and activation by thrombin, the catalytic properties of the mutants were characterized with respect to their ability to cleave the chromogenic substrate Spectrozyme PCa, react with protein C inhibitor (PCI), and inactivate factor Va. Relative to wild-type APC, the mutants cleaved Spectrozyme PCa with identical or improved catalytic efficiencies. Similarly, PCI inhibited mutants with identical or improved second-order rate constants (k(2)) in the absence of heparin. However, the heparin-catalyzed inhibition of mutants by PCI was impaired approximately 10-fold. Analysis of k(2) values by a ternary complex model revealed that the affinities of mutants for heparin were impaired to a similar extent. Moreover, analysis of the NaCl gradient elution profiles of APC derivatives from Heparin-Sepharose supported this conclusion. An oligosaccharide containing 14 residues efficiently catalyzed the PCI inhibition of APC by a template mechanism. Further studies revealed that the ability of Arg(74) and Arg(75) mutants to inactivate factor Va was markedly impaired. We conclude that basic residues of the 70-80-loop are critical for the catalytic function of APC.

Lysine 114 of antithrombin is of crucial importance for the affinity and kinetics of heparin pentasaccharide binding

Lys(114) of the plasma coagulation proteinase inhibitor, antithrombin, has been implicated in binding of the glycosaminoglycan activator, heparin, by previous mutagenesis studies and by the crystal structure of antithrombin in complex with the active pentasaccharide unit of heparin. In the present work, substitution of Lys(114) by Ala or Met was shown to decrease the affinity of antithrombin for heparin and the pentasaccharide by approximately 10(5)-fold at I 0.15, corresponding to a reduction in binding energy of approximately 50%. The decrease in affinity was due to the loss of two to three ionic interactions, consistent with Lys(114) and at least one other basic residue of the inhibitor binding cooperatively to heparin, as well as to substantial nonionic interactions. The mutation minimally affected the initial, weak binding of the two-step mechanism of pentasaccharide binding to antithrombin but appreciably (>40-fold) decreased the forward rate constant of the conformational change in the second step and greatly (>1000-fold) increased the reverse rate constant of this step. Lys(114) is thus of greater importance for the affinity of heparin binding than any of the other antithrombin residues investigated so far, viz. Arg(47), Lys(125), and Arg(129). It contributes more than Arg(47) and Arg(129) to increasing the rate of induction of the activating conformational change, a role presumably exerted by interactions with the nonreducing end trisaccharide unit of the heparin pentasaccharide. However, its major effect, also larger than that of these two residues, is in maintaining antithrombin in the activated state by interactions that most likely involve the reducing end disaccharide unit.

The region of antithrombin interacting with full-length heparin chains outside the high-affinity pentasaccharide sequence extends to Lys136 but not to Lys139

The interaction of a well-defined pentasaccharide sequence of heparin with a specific binding site on antithrombin activates the inhibitor through a conformational change. This change increases the rate of antithrombin inhibition of factor Xa, whereas acceleration of thrombin inhibition requires binding of both inhibitor and proteinase to the same heparin chain. An extended heparin binding site of antithrombin outside the specific pentasaccharide site has been proposed to account for the higher affinity of the inhibitor for full-length heparin chains by interacting with saccharides adjacent to the pentasaccharide sequence. To resolve conflicting evidence regarding the roles of Lys136 and Lys139 in this extended site, we have mutated the two residues to Ala or Gln. Mutation of Lys136 decreased the antithrombin affinity for full-length heparin by at least 5-fold but minimally altered the affinity for the pentasaccharide. As a result, the full-length heparin and pentasaccharide affinities were comparable. The reduced affinity for full-length heparin was associated with the loss of one ionic interaction and was caused by both a lower overall association rate constant and a higher overall dissociation rate constant. In contrast, mutation of Lys139 affected neither full-length heparin nor pentasaccharide affinity. The rate constants for inhibition of thrombin and factor Xa by the complexes between antithrombin and full-length heparin or pentasaccharide were unaffected by both mutations, indicating that neither Lys136 nor Lys139 is involved in heparin activation of the inhibitor. Together, these results show that Lys136 forms part of the extended heparin binding site of antithrombin that participates in the binding of full-length heparin chains, whereas Lys139 is located outside this site.

Analysis of a structural determinant in thrombin-protease nexin 1 complexes that mediates clearance by the low density lipoprotein receptor-related protein

We recently identified a synthetic peptide, Pro47-Ile58, derived from the mature protease nexin 1 (PN1) sequence, that inhibited the low density lipoprotein receptor-related protein (LRP)-mediated internalization of thrombin-PN1 (Th-PN1) complexes. Presently, we have analyzed this sequence in Th-PN1 complex catabolism using two independent approaches: 1) An antibody was generated against Pro47-Ile58, which inhibited complex degradation by 70% but had no effect on the binding of the complexes to cell surface heparins. This places the structural determinant in PN1 mediating complex internalization by the LRP outside of the heparin-binding site. 2) Site-directed genetic variants of PN1 with a single Ala substitution at His48, or two Ala substitutions, one at His48 and another at Asp49, were expressed in Sf9 insect cells. The catabolic rate of complexes formed between Th and the singly substituted and doubly substituted variants was lowered to 50 and 15%, respectively, when compared with the catabolic rate of native Th-PN1 complexes. This is the first analysis of a structural determinant in a serine protease inhibitor (SERPIN) required for LRP-mediated internalization and in part may explain the cryptic nature of this site in the unreacted serine protease inhibitor.

Deconvolution of the fluorescence emission spectrum of human antithrombin and identification of the tryptophan residues that are responsive to heparin binding

Heparin causes an allosterically transmitted conformational change in the reactive center loop of antithrombin and a 40% enhancement of tryptophan fluorescence. We have expressed four human antithrombins containing single Trp --> Phe mutations and determined that the fluorescence of antithrombin is a linear combination of the four tryptophans. The contributions to the spectrum of native antithrombin at 340 nm were 8% for Trp-49, 10% for Trp-189, 19% for Trp-225, and 63% for Trp-307. Trp-225 and Trp-307 accounted for the majority of the heparin-induced fluorescence enhancement, contributing 37 and 36%, respectively. Trp-49 and Trp-225 underwent spectral shifts of 15 nm to blue and 5 nm to red, respectively, in the antithrombin-heparin complex. The blue shift for Trp-49 is consistent with partial burial by contact with heparin, whereas the red shift for Trp-225 and large enhancement probably result from increased solvent access upon heparin-induced displacement of the contact residue Ser-380. The enhancement for Trp-307 may result from the heparin-induced movement of helix H seen in the crystal structure. The time-resolved fluorescence properties of individual tryptophans of wild-type antithrombin were also determined using the four variants and showed that Trp-225 and Trp-307 experienced the largest change in lifetime upon heparin binding, providing support for the steady-state fluorescence deconvolution.

Molecular Modeling of Ligand and Mutation Sites of the Type A Domains of Human von Willebrand Factor and Their Relevance to von Willebrand's Disease

von Willebrand factor (vWF) is a large multimeric, multidomain glycoprotein found in platelets, endothelial cells and plasma. The A1, A2, and A3 domains in vWF mediate binding to glycoprotein Ib, ristocetin, botrocetin, collagen, sulphatides, and heparin and provide a protease cleavage site. Mutations causing types 2B, 2M, and 2A von Willebrand's disease (vWD) are located in the A1 and A2 domains. Homology modeling was performed to provide a molecular interpretation of vWF function and mutation sites. This was based on our previous alignment of 75 vWF-A sequences, the doubly wound α/β fold seen in recent vWF-A crystal structures from complement receptor type 3 and lymphocyte function-associated antigen-1, and our new alignment of 28 vWF A1 and A2 sequences from different species. The active site in doubly-wound α/β folds forms a crevice that is located at the switch point between the two halves of the central β-sheet, and usually contains two metal-binding Asp residues in the vWF-A superfamily. Although one of these Asp residues is absent from the A1, A2, and A3 domains, this crevice is shown to correspond to the ristocetin binding site in the A1 domain and the protease cleavage site in the A2 domain. The residues R571-K572-R578-R579-K585 are found to be conserved in 28 A1 sequences and are predicted to constitute the heparin binding site in the A1 domain. Inspection of the type 2M vWD mutation sites that are involved in downregulation of glycoprotein Ib (GpIb) binding to vWF shows that these are spatially clustered at the carboxyl-edge of the β-sheet and above it in the A1 domain and may directly perturb GpIb binding. In contrast, the type 2B vWD mutation sites that are involved in upregulation of GpIb binding to vWF are spatially clustered at the amino edge of this β-sheet and below it and are located on the opposite side of the A1 domain from the type 2M mutation sites. The type 2B mutations are located between the heparin and GpIb binding sites. Because heparin binding inhibits the interaction with GpIb, this provides an explanation of vWF upregulation. The type 2A vWD mutation sites in the A2 domain correspond to buried residues that are otherwise 100% conserved across all 28 species, and are likely to be important for the correct folding of the A2 domain and its physiologically important protease site.

Molecular Modeling of Ligand and Mutation Sites of the Type A Domains of Human von Willebrand Factor and Their Relevance to von Willebrand's Disease

von Willebrand factor (vWF) is a large multimeric, multidomain glycoprotein found in platelets, endothelial cells and plasma. The A1, A2, and A3 domains in vWF mediate binding to glycoprotein Ib, ristocetin, botrocetin, collagen, sulphatides, and heparin and provide a protease cleavage site. Mutations causing types 2B, 2M, and 2A von Willebrand's disease (vWD) are located in the A1 and A2 domains. Homology modeling was performed to provide a molecular interpretation of vWF function and mutation sites. This was based on our previous alignment of 75 vWF-A sequences, the doubly wound α/β fold seen in recent vWF-A crystal structures from complement receptor type 3 and lymphocyte function-associated antigen-1, and our new alignment of 28 vWF A1 and A2 sequences from different species. The active site in doubly-wound α/β folds forms a crevice that is located at the switch point between the two halves of the central β-sheet, and usually contains two metal-binding Asp residues in the vWF-A superfamily. Although one of these Asp residues is absent from the A1, A2, and A3 domains, this crevice is shown to correspond to the ristocetin binding site in the A1 domain and the protease cleavage site in the A2 domain. The residues R571-K572-R578-R579-K585 are found to be conserved in 28 A1 sequences and are predicted to constitute the heparin binding site in the A1 domain. Inspection of the type 2M vWD mutation sites that are involved in downregulation of glycoprotein Ib (GpIb) binding to vWF shows that these are spatially clustered at the carboxyl-edge of the β-sheet and above it in the A1 domain and may directly perturb GpIb binding. In contrast, the type 2B vWD mutation sites that are involved in upregulation of GpIb binding to vWF are spatially clustered at the amino edge of this β-sheet and below it and are located on the opposite side of the A1 domain from the type 2M mutation sites. The type 2B mutations are located between the heparin and GpIb binding sites. Because heparin binding inhibits the interaction with GpIb, this provides an explanation of vWF upregulation. The type 2A vWD mutation sites in the A2 domain correspond to buried residues that are otherwise 100% conserved across all 28 species, and are likely to be important for the correct folding of the A2 domain and its physiologically important protease site.

Conserved basic residues in the C-type lectin and short complement repeat domains of the G3 region of proteoglycans

Aggrecan is the major proteoglycan of the extracellular matrix in cartilage. It contains two N-terminal globular regions, G1 and G2, and one C-terminal globular region, G3. G3 is implicated in the intracellular processing of aggrecan and contains a C-type lectin carbohydrate recognition domain (CRD), frequent occurrences of a C-terminal short complement repeat (SCR) domain, and occasionally an N-terminal epidermal growth factor domain. The CRD and SCR domains in 13 G3 sequences were each subjected to structural analysis. Alignment of 131 sequences from all seven groups in the CRD superfamily defined a consensus length of 136 residues, in which 32% of residues were conserved. Although the G3 CRD sequences agreed with this consensus, they also contained five fully conserved basic residues that are atypical of the CRD superfamily. Homology modelling showed that four of these residues are located on a surface region on the CRD that is separate from the Ca2+-binding residues involved in carbohydrate interactions. One conserved basic residue is identical in position with that of a conserved basic residue that mediates hyaluronate binding in the structurally related proteoglycan tandem repeat (PTR) domain in G1 and in link protein. The alignment of 13 G3 SCR sequences with 101 sequences in the SCR superfamily showed good agreement with conserved residues in the SCR superfamily. There are also five conserved basic residues in the G3 SCR that are atypical of the SCR superfamily, and homology modelling showed that all five were located on one surface of the SCR. It is concluded that both the CRD and SCR domains in G3 possess basic residues that are atypical of their superfamilies and might be related to function, and that the G3 CRD domain shows an evolutionary relationship to the PTR domain in G1.