The effects of ligand binding on the backbone dynamics of the kringle 1 domain of human plasminogen

In Vitro Study of the Fibrinolytic Activity via Single Chain Urokinase-Type Plasminogen Activator and Molecular Docking of FGFC1

Fungi fibrinolytic compound 1 (FGFC1) is a rare marine-derived compound that can enhance fibrinolysis both in vitro and in vivo. The fibrinolytic activity characterization of FGFC1 mediated by plasminogen (Glu-/Lys-) and a single-chain urokinase-type plasminogen activator (pro-uPA) was further evaluated. The binding sites and mode of binding between FGFC1 and plasminogen were investigated by means of a combination of in vitro experiments and molecular docking. A 2.2-fold enhancement of fibrinolytic activity was achieved at 0.096 mM FGFC1, whereas the inhibition of fibrinolytic activity occurred when the FGFC1 concentration was above 0.24 mM. The inhibition of fibrinolytic activity of FGFC1 by 6-aminohexanoic acid (EACA) and tranexamic acid (TXA) together with the docking results revealed that the lysine-binding sites (LBSs) play a crucial role in the process of FGFC1 binding to plasminogen. The action mechanism of FGFC1 binding to plasminogen was inferred, and FGFC1 was able to induce plasminogen to exhibit an open conformation by binding through the LBSs. The molecular docking results showed that docking of ligands (EACA, FGFC1) with receptors (KR1–KR5) mainly occurred through hydrophilic and hydrophobic interactions. In addition, the binding affinity values of EACA to KR1–KR5 were −5.2, −4.3, −3.7, −4.5, and −4.3 kcal/moL, respectively, and those of FGFC1 to KR1–KR5 were −7.4, −9.0, −6.3, −8.3, and −6.7 kcal/moL, respectively. The findings demonstrate that both EACA and FGFC1 bound to KR1–KR5 with moderately high affinity. This study could provide a theoretical basis for the clinical pharmacology of FGFC1 and establish a foundation for practical applications of FGFC1.

EVOLUTION OF A CATALYTICALLY EFFECTIVE MODEL ENZYME: THE IMPORTANCE OF TUNED CONFORMATIONAL FLUCTUATIONS

Possible causal connections between the dynamics of a thermally fluctuating model enzyme molecule and catalysis are explored. The model is motivated by observations from experiment and simulation that amino acid residues residing in different enzymatic domains may show markedly different degrees of conformational freedom. Consequently, we are interested in the catalytic efficacy of an enzyme as a function of long-range many-atom cooperative effects resulting from strong, moderate, and weak interactions between enzymatic residues. Here we show and quantify through molecular dynamics simulations how the number and distribution of these interactions affects an enzyme's conformational fluctuation dynamics and its effectiveness as a catalyst. For any given distribution of "stiff" and "loose" enzymatic domains, catalytic fitness is defined as the number of chemical events — specifically the number of times a catalytic residue and substrate surmount a chemical reaction barrier — during molecular dynamics simulation. Through mutation, recombination, and a selection procedure following the ideas of Darwinian evolution, a genetic algorithm drives a population of enzyme molecules to greater catalytic fitness by modifying the mix of stiff and loose interactions. Approximately 30,000 different enzyme molecules are generated by the genetic algorithm — each with a unique number and distribution of strong, moderate, and weak inter-residue interactions. While the catalytically least fit enzyme exhibits 16 chemical events, the fittest boasts 253. That point mutations far from the active-site chemistry in the fittest enzyme have a strong effect on the number of chemical events suggests that catalysis depends, in part, on long-range many-atom globally correlated dynamical fluctuations.

NMR backbone dynamics of VEK-30 bound to the human plasminogen kringle 2 domain

To gain insights into the mechanisms for the tight and highly specific interaction of the kringle 2 domain of human plasminogen (K2(Pg)) with a 30-residue internal peptide (VEK-30) from a group A streptococcal M-like protein, the dynamic properties of free and bound K2(Pg) and VEK-30 were investigated using backbone amide (15)N-NMR relaxation measurements. Dynamic parameters, namely the generalized order parameter, S(2), the local correlation time, tau(e), and the conformational exchange contribution, R(ex), were obtained for this complex by Lipari-Szabo model-free analysis. The results show that VEK-30 displays distinctly different dynamic behavior as a consequence of binding to K2(Pg), manifest by decreased backbone flexibility, particularly at the binding region of the peptide. In contrast, the backbone dynamics parameters of K2(Pg) displayed similar patterns in the free and bound forms, but, nonetheless, showed interesting differences. Based on our previous structure-function studies of this interaction, we also made comparisons of the VEK-30/K2(Pg) dynamics results from different kringle modules complexed with small lysine analogs. The differences in dynamics observed for kringles with different ligands provide what we believe to be new insights into the interactions responsible for protein-ligand recognition and a better understanding of the differences in binding affinity and binding specificity of kringle domains with various ligands.

Solution structure and functional characterization of human plasminogen kringle 5

The ligand binding properties of the kringle 5 (K5) domain of human plasminogen have been investigated via intrinsic tryptophan fluorescence. The oleic acid (OA) affinity for K5 was quantified, yielding an association constant K(a) approximately 2.08 x 10(4) mM(-1). Simultaneously, it was determined that OA and trans-4-(aminomethyl)cyclohexanecarboxylic acid (AMCHA) (K(a) approximately 50 mM(-1)) compete for binding to K5. The solution structure of K5 in the presence of 11 mM AMCHA was solved via NMR spectroscopy (protein heavy atom RMSD approximately 0.93 +/- 0.12 A). The AMCHA binding site was localized via (1)H/(15)N chemical shift perturbation mapping assisted by in silico docking. We have found that AMCHA binds at the canonical kringle lysine binding site (LBS), structured by the Pro54-Gly60 segment plus the neighboring Phe36, Thr37, Trp62, Leu71, and Tyr72 residues. The segment 30-42, encompassing LBS residues, appears to be endowed with a higher degree of structural flexibility as suggested by the relatively lower value of S(2), the generalized order parameter, consistent with a higher backbone heavy atom RMSD of approximately 1.22 A (vs 0.84 A overall) between the two monomeric units in the crystal unit cell, of potential significance for ligand binding. OA was found to perturb the same area of the protein, namely, the LBS, as well as Tyr74. Combined with previous studies, the observation of OA binding expands the range of ligands that interact with kringle 5 while it widens the scope of potential biological functions for kringle domains.

Activation of hepatocyte growth factor by urokinase-type plasminogen activator is ionic strength-dependent

The hepatocyte growth factor (HGF) is a multifunctional cytokine that is produced as latent scHGF (single chain HGF). Various proteases reportedly cleave scHGF to generate the active two-chain form (HGF), including u-PA (urokinase-type plasminogen activator), t-PA (tissue-type plasminogen activator), kallikrein, Factor XIa, Factor XIIa, HGF activator and matriptase. Considerable evidence indicates that, in vivo, u-PA activates scHGF in the liver; however, the in vivo results have not been uniformly supported by in vitro experiments. We now report that cleavage of scHGF by high-molecular-mass u-PA (abbreviated u-PA throughout) is sensitive to ionic strength. scHGF cleavage by u-PA was accelerated as the ionic strength was decreased. This result was equivalent irrespective of whether the predominant anion was chloride or acetate. Lmw-u-PA (low-molecular-mass u-PA) was ineffective at cleaving scHGF, regardless of ionic strength. Although scHGF shares homology with plasminogen, EACA (-amino-caproic acid) did not regulate u-PA-mediated scHGF cleavage. Soluble HGF receptor (MET) and soluble u-PAR (u-PA receptor) inhibited the scHGF cleavage. These results support a model in which the ability of u-PA to activate scHGF in vivo may be highly dependent on local conditions within the extracellular space.

A fast-acting, modular-structured staphylokinase fusion with Kringle-1 from human plasminogen as the fibrin-targeting domain offers improved clot lysis efficacy

To develop a fast-acting clot dissolving agent, a clot-targeting domain derived from the Kringle-1 domain in human plasminogen was fused to the C-terminal end of staphylokinase with a linker sequence in between. Production of this fusion protein in Bacillus subtilis and Pichia pastoris was examined. The Kringle domain in the fusion protein produced from B. subtilis was improperly folded because of its complicated disulfide-bond profile, whereas the staphylokinase domain produced from P. pastoris was only partially active because of an N-linked glycosylation. A change of the glycosylation residue, Thr-30, to alanine resulted in a non-glycosylated biologically active fusion. The resulting mutein, designated SAKM3-L-K1, was overproduced in P. pastoris. Each domain in SAKM3-L-K1 was functional, and this fusion showed fibrin binding ability by binding directly to plasmin-digested clots. In vitro fibrin clot lysis in a static environment and plasma clot lysis in a flow-cell system demonstrated that the engineered fusion outperformed the non-fused staphylokinase. The time required for 50% clot lysis was reduced by 20 to 500% under different conditions. Faster clot lysis can potentially reduce the degree of damage to occluded heart tissues.

Comparison of (13)C(alpha)H and (15)NH backbone dynamics in protein GB1

This study presents a site-resolved experimental view of backbone C(alpha)H and NH internal motions in the 56-residue immunoglobulin-binding domain of streptococcal protein G, GB1. Using (13)C(alpha)H and (15)NH NMR relaxation data [T(1), T(2), and NOE] acquired at three resonance frequencies ((1)H frequencies of 500, 600, and 800 MHz), spectral density functions were calculated as F(omega) = 2omegaJ(omega) to provide a model-independent way to visualize and analyze internal motional correlation time distributions for backbone groups in GB1. Line broadening in F(omega) curves indicates the presence of nanosecond time scale internal motions (0.8 to 5 nsec) for all C(alpha)H and NH groups. Deconvolution of F(omega) curves effectively separates overall tumbling and internal motional correlation time distributions to yield more accurate order parameters than determined by using standard model free approaches. Compared to NH groups, C(alpha)H internal motions are more broadly distributed on the nanosecond time scale, and larger C(alpha)H order parameters are related to correlated bond rotations for C(alpha)H fluctuations. Motional parameters for NH groups are more structurally correlated, with NH order parameters, for example, being larger for residues in more structured regions of beta-sheet and helix and generally smaller for residues in the loop and turns. This is most likely related to the observation that NH order parameters are correlated to hydrogen bonding. This study contributes to the general understanding of protein dynamics and exemplifies an alternative and easier way to analyze NMR relaxation data.

Heat capacities and a snapshot of the energy landscape in protein GB1 from the pre-denaturation temperature dependence of backbone NH nanosecond fluctuations

Protein stability is usually characterized calorimetrically by a melting temperature and related thermodynamic parameters. Despite its importance, the microscopic origin of the melting transition and the relationship between thermodynamic stability and dynamics remains a mystery. Here, NMR relaxation parameters were acquired for backbone 15NH groups of the 56 residue immunoglobulin-binding domain of streptococcal protein G over a pre-denaturation temperature range of 5-50 degrees C. Relaxation data were analyzed using three methods: the standard three-Lorentzian model free approach; the F(omega)=2omegaJ(omega) spectral density approach that yields motional correlation time distributions, and a new approach that determines frequency-dependent order parameters. Regardless of the method of analysis, the temperature dependence of internal motional correlation times and order parameters is essentially the same. Nanosecond time-scale internal motions are found for all NHs in the protein, and their temperature dependence yields activation energies ranging up to about 33kJ/mol residue. NH motional barrier heights are structurally correlated, with the largest energy barriers being found for residues in the most "rigid" segments of the fold: beta-strands 1 and 4 and the alpha-helix. Trends in this landscape also parallel the free energy of folding-unfolding derived from hydrogen-deuterium (H-D) exchange measurements, indicating that the energetics for internal motions occurring on the nanosecond time-scale mirror those occurring on the much slower time-scale of H-D exchange. Residual heat capacities, derived from the temperature dependence of order parameters, range from near zero to near 100J/mol K residue and correlate with this energy landscape. These results provide a unique picture of this protein's energy landscape and a relationship between thermodynamic stability and dynamics that suggests thermosensitive regions in the fold that could initiate the melting process.

Trifluoroethanol increases the stability of Delta(5)-3-ketosteroid isomerase. 15N NMR relaxation studies

In the equilibrium unfolding process of Delta(5)-3-ketosteroid isomerase from Pseudomonas testosteroni by urea, it was observed that the enzyme stability increases by 2.5 kcal/mol in the presence of 5% trifluoroethanol (TFE). To elucidate the increased enzyme stability by TFE, the backbone dynamics of Delta(5)-3-ketosteroid isomerase were studied in the presence and absence of 5% TFE by (15)N NMR relaxation measurements, and the motional parameters (S(2), tau(e), and R(ex)) were extracted from the relaxation data using the model-free formalism. The presence of 5% TFE causes little change or a slight increase in the order parameters (S(2)) for a number of residues, which are located mainly in the dimer interface region. However, the majority of the residues exhibit reduced order parameters in the presence of 5% TFE, indicating that high frequency (pico- to nanosecond) motions are generally enhanced by TFE. The results suggest that the entropy can be an important factor for the enzyme stability, and the increase in entropy by TFE is partially responsible for the increased stability of Delta(5)-3-ketosteroid isomerase.

Backbone dynamics of the regulatory domain of calcium vector protein, studied by (15)N relaxation at four fields, reveals unique mobility characteristics of the intermotif linker

CaVP is a calcium-binding protein from amphioxus. It has a modular composition with two domains, but only the two EF-hand motifs localized in the C-terminal domain are functional. We recently determined the solution structure of this regulatory half (C-CaVP) in the Ca(2+)-saturated form and characterized the stepwise ion binding. This paper reports the (15)N nuclear relaxation rates of the Ca(2+)-saturated C-CaVP, measured at four different NMR fields (9.39, 11.74, 14.1, and 18.7 T), which were used to map the spectral density function for the majority of the amide H(N)-N vectors. Fitting the spectral density values at eight frequencies by a model-free approach, we obtained the microdynamic parameters characterizing the global and internal movements of the polypeptide backbone. The two EF-hand motifs, including the ion binding loops, behave like compact structural units with restricted mobility as reflected in the quite uniform order parameter and short internal correlation time (< 20 nsec). Comparative analysis of the two Ca(2+) binding sites shows that site III, having a larger affinity for the metal ion, is generally more rigid, and the amide vector in the second residue of each loop is significantly less restricted. The linker fragment is animated simultaneously by a larger amplitude fast motion and a slow conformational exchange on a microsecond to millisecond time scale. The backbone dynamics of C-CaVP characterized here is discussed in relation with other well-characterized Ca(2+)-binding proteins.

SUPPLEMENTAL MATERIAL

See www.proteinscience.org

Structure and binding determinants of the recombinant kringle-2 domain of human plasminogen to an internal peptide from a group A Streptococcal surface protein

The X-ray crystal structure of a complex of a modified recombinant kringle-2 domain of human plasminogen, K2Pg[C4G/E56D/L72Y] (mK2Pg), containing an upregulated lysine-binding site, bound to a functional 30 residue internal peptide (VEK-30) from an M-type protein of a group A Streptococcus surface protein, has been determined by molecular replacement methods using K4Pg as a model, and refined at 2.7 A resolution to a R-factor of 19.5 %. The X-ray crystal structure shows that VEK-30 exists as a nearly end-to-end alpha-helix in the complex with mK2Pg. The final structure also revealed that Arg17 and His18 of VEK-30 served as cationic loci for Asp54 and Asp56 of the consensus lysine-binding site of mK2Pg, while Glu20 of VEK-30 coordinates with Arg69 of the cationic binding site of mK2Pg. The hydrophobic ligand-binding pocket in mK2Pg, consisting primarily of Trp60 and Trp70, situated between the positive and negative centers of the lysine-binding site, is utilized in a novel manner in stabilizing the interaction with VEK-30 by forming a cation-pi-electron-mediated association with the positive side-chain of Arg17 of this peptide. Additional lysine-binding sites, as well as exosite electrostatic and hydrogen bonding interactions involving Glu9 and Lys14 of VEK-30, were observed in the structural model. The importance of these interactions were tested in solution by investigating the binding constants of synthetic variants of VEK-30 to mK2Pg, and it was found that, Lys14, Arg17, His18, and Glu20 of VEK-30 were the most critical amino acid binding determinants. With regard to the solution studies, circular dichroism analysis of the titration of VEK-30 with mK2Pg demonstrated that the peptidic alpha-helical structure increased substantially when bound to the kringle module, in agreement with the X-ray results. This investigation is the first to delineate structurally the mode of interaction of the lysine-binding site of a kringle with an internal pseudo-lysine residue of a peptide or protein that functionally interacts with a kringle module, and serves as a paradigm for this important class of interactions.

15N NMR relaxation studies of backbone dynamics in free and steroid-bound Delta 5-3-ketosteroid isomerase from Pseudomonas testosteroni

The backbone dynamics of Delta(5)-3-ketosteroid isomerase (KSI) from Pseudomonas testosteroni has been studied in free enzyme and its complex with a steroid ligand, 19-nortestosterone hemisuccinate (19-NTHS), by (15)N relaxation measurements. The relaxation data were analyzed using the model-free formalism to extract the model-free parameters (S(2), tau(e), and R(ex)) and the overall rotational correlation time (tau(m)). The rotational correlation times were 19.23 +/- 0.08 and 17.08 +/- 0.07 ns with the diffusion anisotropies (D( parallel)/D( perpendicular)) of 1.26 +/- 0.03 and 1.25 +/- 0.03 for the free and steroid-bound KSI, respectively. The binding of 19-NTHS to free KSI causes a slight increase in the order parameters (S(2)) for a number of residues, which are located mainly in helix A1 and strand B4. However, the majority of the residues exhibit reduced order parameters upon ligand binding. In particular, strands B3, B5, and B6, which have most of the residues involved in the dimer interaction, have the reduced order parameters in the steroid-bound KSI, indicating the increased high-frequency (pico- to nanosecond) motions in the intersubunit region of this homodimeric enzyme. Our results differ from those of previous studies on the backbone dynamics of monomeric proteins, in which high-frequency internal motions are typically restricted upon ligand binding.