Computational simulations assessment of mutations impact on streptokinase (SK) from a group G streptococci with enhanced activity - insights into the functional roles of structural dynamics flexibility of SK and stabilization of SK-µplasmin catalytic complex.

Streptokinase (SK), a plasminogen activator (PA) that converts inactive plasminogen (Pg) to plasmin (Pm), is a protein secreted by groups A, C, and G streptococci (GAS, GCS, and GGS, respectively), with high sequence divergence and functional heterogeneity. While roles of some residual changes in altered SK functionality are shown, the underlying structural mechanisms are less known. Herein, using computational approaches, we analyzed the conformational basis for the increased activity of SK from a GGS (SKG132) isolate with four natural residual substitutions (Ile33Phe, Arg45Gln, Asn228Lys, Phe287Ile) compared to the standard GCS (SKC). Using the crystal structure of SK.Pm catalytic complex as main template SKC.µPm catalytic complex was modeled through homology modeling process and validated by several online validation servers. Subsequently, SKG132.µPm structure was constructed by altering the corresponding residual substitutions. Results of three independent MD simulations showed increased RMSF values for SKG132.µPm, indicating the enhanced structural flexibility compared to SKC.µPm, specially in 170 and 250 loops and three regions: R1 (149-161), R2 (182-215) and R3 (224-229). In parallel, the average number of Hydrogen bonds in 170 loop, R2 and R3 (especially for Asn228Lys) of SKG132 compared to that of the SKC was decreased. Accordingly, residue interaction networks (RINs) analyses indicated that Asn228Lys might induce more level of structural flexibility by generation of free Lys256, while Phe287Ile and Ile33Phe enhanced the stabilization of the SKG132.µPm catalytic complex. These results denoted the potential role of the optimal dynamic state and stabilized catalytic complex for increased PA potencies of SK as a thrombolytic drug.

PDB2Graph: A toolbox for identifying critical amino acids map in proteins based on graph theory.

The integrative and cooperative nature of protein structure involves the assessment of topological and global features of constituent parts. Network concept takes complete advantage of both of these properties in the analysis concomitantly. High compatibility to structural concepts or physicochemical properties in addition to exploiting a remarkable simplification in the system has made network an ideal tool to explore biological systems. There are numerous examples in which different protein structural and functional characteristics have been clarified by the network approach. Here, we present an interactive and user-friendly Matlab-based toolbox, PDB2Graph, devoted to protein structure network construction, visualization, and analysis. Moreover, PDB2Graph is an appropriate tool for identifying critical nodes involved in protein structural robustness and function based on centrality indices. It maps critical amino acids in protein networks and can greatly aid structural biologists in selecting proper amino acid candidates for manipulating protein structures in a more reasonable and rational manner. To introduce the capability and efficiency of PDB2Graph in detail, the structural modification of Calmodulin through allosteric binding of Ca(2+) is considered. In addition, a mutational analysis for three well-identified model proteins including Phage T4 lysozyme, Barnase and Ribonuclease HI, was performed to inspect the influence of mutating important central residues on protein activity.