A structural comparison of 21 inhibitor complexes of the aspartic proteinase from Endothia parasitica.

The aspartic proteinases are an important family of enzymes associated with several pathological conditions such as hypertension (renin), gastric ulcers (pepsin), neoplastic disease (cathepsins D and E), and AIDS (HIV proteinase). Studies of inhibitor binding are therefore of great importance for design of novel inhibitors for potential therapeutic applications. Numerous X-ray analyses have shown that transition-state isostere inhibitors of aspartic proteinases bind in similar extended conformations in the active-site cleft of the target enzyme. Upon comparison of 21 endothiapepsin inhibitor complexes, the hydrogen bond lengths were found to be shortest where the isostere (P1-P'1) interacts with the enzyme's catalytic aspartate pair. Hydrogen bonds with good geometry also occur at P'2, and more so at P3, where a conserved water molecule is involved in the interactions. Weaker interactions also occur at P2, where the side-chain conformations of the inhibitors appear to be more variable than at the more tightly held positions. At P2 and, to a lesser extent, P3, the side-chain conformations depend intriguingly on interactions with spatially adjacent side chains, namely P'1 and P1, respectively. The tight binding at P1-P'1, P3, and P'2 is also reflected in the larger number of van der Waals contacts and the large decreases in solvent-accessible area at these positions, as well as their low temperature factors. Our analysis substantiates earlier proposals for the locations of protons in the transition-state complex. Aspartate 32 is probably ionized in the complexes, its charge being stabilized by 1, or sometimes 2, hydrogen bonds from the transition-state analogues at P1. The detailed comparison also indicates that the P1 and P2 residues of substrate in the ES complex may be strained by the extensive binding interactions at P3, P'1, and P'2 in a manner that would facilitate hydrolysis of the scissile peptide bond.     

Direct observation by X-ray analysis of the tetrahedral "intermediate" of aspartic proteinases.

We report the X-ray analysis at 2.0 A resolution for crystals of the aspartic proteinase endothiapepsin (EC 3.4.23.6) complexed with a potent difluorostatone-containing tripeptide renin inhibitor (CP-81,282). The scissile bond surrogate, an electrophilic ketone, is hydrated in the complex. The pro-(R) (statine-like) hydroxyl of the tetrahedral carbonyl hydrate is hydrogen-bonded to both active-site aspartates 32 and 215 in the position occupied by a water in the native enzyme. The second hydroxyl oxygen of the hydrate is hydrogen-bonded only to the outer oxygen of Asp 32. These experimental data provide a basis for a model of the tetrahedral intermediate in aspartic proteinase-mediated cleavage of the amide bond. This indicates a mechanism in which Asp 32 is the proton donor and Asp 215 carboxylate polarizes a bound water for nucleophilic attack. The mechanism involves a carboxylate (Asp 32) that is stabilized by extensive hydrogen bonding, rather than an oxyanion derivative of the peptide as in serine proteinase catalysis.     

Crystal structures of Aspergillus oryzae aspartic proteinase and its complex with an inhibitor pepstatin at 1.9A resolution

The X-ray structures of Aspergillus oryzae aspartic proteinase (AOAP) and its complex with inhibitor pepstatin have been determined at 1.9A resolution. AOAP was crystallized in an orthorhombic system with the space group P2(1)2(1)2(1) and cell dimensions of a=49.4A, b=79.4A, and c=93.6A. By the soaking of pepstatin, crystals are transformed into a monoclinic system with the space group C2 and cell dimensions of a=106.8A, b=38.6A, c=78.7A, and beta=120.3 degrees. The structures of AOAP and AOAP/pepstatin complex were refined to an R-factor of 0.177 (R(free)=0.213) and of 0.185 (0.221), respectively. AOAP has a crescent-shaped structure with two lobes (N-lobe and C-lobe) and the deep active site cleft is constructed between them. At the center of the active site cleft, two Asp residues (Asp33 and Asp214) form the active dyad with a hydrogen bonding solvent molecule between them. Pepstatin binds to the active site cleft via hydrogen bonds and hydrophobic interactions with the enzyme. The structures of AOAP and AOAP/pepstatin complex including interactions between the enzyme and pepstatin are very similar to those of other structure-solved aspartic proteinases and their complexes with pepstatin. Generally, aspartic proteinases cleave a peptide bond between hydrophobic amino acid residues, but AOAP can also recognize the Lys/Arg residue as well as hydrophobic amino acid residues, leading to the activation of trypsinogen and chymotrypsinogen. The X-ray structure of AOAP/pepstatin complex and preliminary modeling show two possible sites of recognition for the positively charged groups of Lys/Arg residues around the active site of AOAP.     

A neutron Laue diffraction study of endothiapepsin: implications for the aspartic proteinase mechanism.

Current proposals for the catalytic mechanism of aspartic proteinases are largely based on X-ray structures of bound oligopeptide inhibitors possessing nonhydrolyzable analogues of the scissile peptide bond. However, the positions of protons on the catalytic aspartates and the ligand in these complexes have not been determined with certainty. Thus, our objective was to locate crucial protons at the active site of an inhibitor complex since this will have major implications for a detailed understanding of the mechanism of action. We have demonstrated that high-resolution neutron diffraction data can be collected from crystals of the fungal aspartic proteinase endothiapepsin bound to a transition state analogue (H261). The neutron structure of the complex has been refined at a resolution of 2.1 A to an R-factor of 23.5% and an R(free) of 27.4%. This work represents the largest protein structure studied to date by neutron crystallography at high resolution. The neutron data demonstrate that 49% of the main chain nitrogens have exchanged their hydrogen atoms with D2O in the mother liquor. The majority of residues resisting exchange are buried within core beta-sheet regions of the molecule. The neutron maps confirm that the protein has a number of buried ionized carboxylate groups which are likely to give the molecule a net negative charge even at very low pH, thereby accounting for its low pI. The functional groups at the catalytic center have clearly undergone H-D exchange despite being buried by the inhibitor occupying the active site cleft. Most importantly, the data provide convincing evidence that Asp 215 is protonated and that Asp 32 is the negatively charged residue in the transition state complex. This has an important bearing on mechanistic proposals for this class of proteinase.     

Veni,vidi, vici - atomic resolution unravelling the mysteries of protein function

Atomic resolution macromolecular crystallography has become a powerful and versatile tool in structural biology; the number of atomic resolution structures is steadily increasing. Novel techniques are being developed and the use of complementary methods that span the field from sample preparation to validation and analysis of the resulting models has emerged. These allow the fuller exploitation of the information stored in crystal structures and reveal a depth of structural detail that was unattainable in the recent past.     

Atomic resolution analysis of the catalytic site of an aspartic proteinase and an unexpected mode of binding by short peptides

The X-ray structures of native endothiapepsin and a complex with a hydroxyethylene transition state analog inhibitor (H261) have been determined at atomic resolution. Unrestrained refinement of the carboxyl groups of the enzyme by using the atomic resolution data indicates that both catalytic aspartates in the native enzyme share a single negative charge equally; that is, in the crystal, one half of the active sites have Asp 32 ionized and the other half have Asp 215 ionized. The electron density map of the native enzyme refined at 0.9 A resolution demonstrates that there is a short peptide (probably Ser-Thr) bound noncovalently in the active site cleft. The N-terminal nitrogen of the dipeptide interacts with the aspartate diad of the enzyme by hydrogen bonds involving the carboxyl of Asp 215 and the catalytic water molecule. This is consistent with classical findings that the aspartic proteinases can be inhibited weakly by short peptides and that these enzymes can catalyze transpeptidation reactions. The dipeptide may originate from autolysis of the N-terminal Ser-Thr sequence of the enzyme during crystallization.     

Isolation and sequence analysis of the genomic DNA fragment encoding an aspartic proteinase inhibitor homologue from potato (Solanum tuberosum L.)

A genomic DNA clone encoding an aspartic proteinase inhibitor of potato was isolated from a lambda EMBL3 phage library using the aspartic proteinase inhibitor cDNA as a hybridization probe. The gene has all characteristic sequences normally found in eucaryotic genes. Typical CAAT and TATA box sequences were found in the 5-upstream region. In this part are also two putative regulatory AGGA box sequences located. In the genomic sequence there are no intron sequences interrupting the coding region. An open reading frame of the gene encodes a precursor protein of 217 amino acids which shows high percent identity with the aspartic proteinase inhibitor cDNA.     

X-ray, neutron and NMR studies of the catalytic mechanism of aspartic proteinases

Current proposals for the catalytic mechanism of aspartic proteinases are largely based on X-ray structures of bound oligopeptide inhibitors possessing non-hydrolysable analogues of the scissile peptide bond. Until recent years, the positions of protons on the catalytic aspartates and the ligand in these complexes had not been determined with certainty due to the inadequate resolution of these analyses. There has been much interest in locating the catalytic protons at the active site of aspartic proteinases since this has major implications for detailed understanding of the mechanism of action and the design of improved transition state mimics for therapeutic applications. In this review we discuss the results of studies which have shed light on the locations of protons at the catalytic centre. The first direct determination of the proton positions stemmed from neutron diffraction data collected from crystals of the fungal aspartic proteinase endothiapepsin bound to a transition state analogue (H261). The neutron structure of the complex at a resolution of 2.1 A provided evidence that Asp 215 is protonated and that Asp 32 is the negatively charged residue in the transition state complex. Atomic resolution X-ray studies of inhibitor complexes have corroborated this finding. A similar study of the native enzyme established that it, unexpectedly, has a dipeptide bound at the catalytic site which is consistent with classical reports of inhibition by short peptides and the ability of pepsins to catalyse transpeptidation reactions. Studies by NMR have confirmed the findings of low-barrier and single-well hydrogen bonds in the complexes with transition state analogues.     

On the reliability of peptide nonplanarity seen in ultra‐high resolution crystal structures

Ultra‐high resolution protein crystal structures have been considered as relatively reliable sources for defining details of protein geometry, such as the extent to which the peptide unit deviates from planarity. Chellapa and Rose (Proteins 2015; 83:1687) recently called this into question, reporting that for a dozen representative protein structures determined at ～1 Å resolution, the diffraction data could be equally well fit with models restrained to have highly planar peptides, i.e. having a standard deviation of the ω torsion angles of only ～1° instead of the typically observed value of ～6°. Here, we document both conceptual and practical shortcomings of that study and show that the more tightly restrained models are demonstrably incorrect and do not fit the diffraction data equally well. We emphasize the importance of inspecting electron density maps when investigating the agreement between a model and its experimental data. Overall, this report reinforces that modern standard refinement protocols have been well‐conceived and that ultra‐high resolution protein crystal structures, when evaluated carefully and used with an awareness of their levels of coordinate uncertainty, are powerful sources of information for providing reliable information about the details of protein geometry.     

Sub-atomic resolution crystal structure of cholesterol oxidase: what atomic resolution crystallography reveals about enzyme mechanism and the role of the FAD cofactor in redox activity

The crystal structure of cholesterol oxidase, a 56kDa flavoenzyme was anisotropically refined to 0.95A resolution. The final crystallographic R-factor and R(free) value is 11.0% and 13.2%, respectively. The quality of the electron density maps has enabled modeling of alternate conformations for 83 residues in the enzyme, many of which are located in the active site. The additional observed structural features were not apparent in the previous high-resolution structure (1.5A resolution) and have enabled the identification of a narrow tunnel leading directly to the isoalloxazine portion of the FAD prosthetic group. The hydrophobic nature of this narrow tunnel suggests it is the pathway for molecular oxygen to access the isoalloxazine group for the oxidative half reaction. Resolving the alternate conformations in the active site residues provides a model for the dynamics of substrate binding and a potential oxidation triggered gating mechanism involving access to the hydrophobic tunnel. This structure reveals that the NE2 atom of the active site histidine residue, H447, critical to the redox activity of this flavin oxidase, acts as a hydrogen bond donor rather than as hydrogen acceptor. The atomic resolution structure of cholesterol oxidase has revealed the presence of hydrogen atoms, dynamic aspects of the protein and how side-chain conformations are correlated with novel structural features such as the oxygen tunnel. This new structural information has provided us with the opportunity to re-analyze the roles played by specific residues in the mechanism of the enzyme.     

X-ray structures of Sap1 and Sap5: structural comparison of the secreted aspartic proteinases from Candida albicans

Proteolytic activity is an important virulence factor for Candida albicans (C. albicans). It is attributed to the family of the secreted aspartic proteinases (Saps) from C. albicans with a minimum of 10 members. Saps show controlled expression and regulation for the individual stages of the infection process. Distinct isoenzymes can be responsible for adherence and tissue damage of local infections, while others cause systemic diseases. Earlier, only the structures of Sap2 and Sap3 were known. In our research, we have now succeeded in solving the X-ray crystal structures of the apoenzyme of Sap1 and Sap5 in complex with pepstatin A at 2.05 and 2.5 A resolution, respectively. With the structure of Sap1, we have completed the set of structures of isoenzyme subgroup Sap1-3. Of subgroup Sap4-6, the structure of the enzyme Sap5 is the first structure that has been described up to now. This facilitates comparison of structural details as well as inhibitor binding modes among the different subgroup members. Structural analysis reveals a highly conserved overall secondary structure of Sap1-3 and Sap5. However, Sap5 clearly differs from Sap1-3 by its electrostatic overall charge as well as through structural conformation of its entrance to the active site cleft. Design of inhibitors specific for Sap5 should concentrate on the S4 and S3 pockets, which significantly differ from Sap1-3 in size and electrostatic charge. Both Sap1 and Sap5 seem to play a major part in superficial Candida infections. Determination of the isoenzymes' structures can contribute to the development of new Sap-specific inhibitors for the treatment of superficial infections with a structure-based drug design program.     

Atomic resolution structures of Escherichia coli and Bacillus anthracis malate synthase A: Comparison with isoform G and implications for structure‐based drug discovery

Enzymes of the glyoxylate shunt are important for the virulence of pathogenic organisms such as Mycobacterium tuberculosis and Candida albicans. Two isoforms have been identified for malate synthase, the second enzyme in the pathway. Isoform A, found in fungi and plants, comprises ～530 residues, whereas isoform G, found only in bacteria, is larger by ～200 residues. Crystal structures of malate synthase isoform G from Escherichia coli and Mycobacterium tuberculosis were previously determined at moderate resolution. Here we describe crystal structures of E. coli malate synthase A (MSA) in the apo form (1.04 Å resolution) and in complex with acetyl‐coenzyme A and a competitive inhibitor, possibly pyruvate or oxalate (1.40 Å resolution). In addition, a crystal structure for Bacillus anthracis MSA at 1.70 Å resolution is reported. The increase in size between isoforms A and G can be attributed primarily to an inserted α/β domain that may have regulatory function. Upon binding of inhibitor or substrate, several active site loops in MSA undergo large conformational changes. However, in the substrate bound form, the active sites of isoforms A and G from E. coli are nearly identical. Considering that inhibitors bind with very similar affinities to both isoforms, MSA is as an excellent platform for high‐resolution structural studies and drug discovery efforts.     

Atomic resolution structure of prokaryotic phospholipase A2: analysis of internal motion and implication for a catalytic mechanism

We have found a secreted phospholipase A(2) (PLA(2), EC 3.1.1.4) from Streptomyces violaceoruber A-2688, which is the first PLA(2) identified in prokaryote, and determined its tertiary structure by NMR and X-ray analyses. In this study, we collected the X-ray diffraction data of the bacterial PLA(2) at room temperature (297 K) using conventional MoK(alpha) radiation and refined the structure at a 1.05 A resolution. The atomic resolution analysis led us to introduce disordered conformations and hydrogen atoms into a full anisotropic model. The molecular motion, which is expressed as the sum of rigid-body motion and internal motion of protein, is roughly estimated as the thermal motion when the X-ray diffraction data are collected at room temperature. In this study, we applied a TLS (rigid-body motion in terms of translation, libration, and screw motions) model to analyze the rigid-body motion of the bacterial PLA(2) and calculated the internal motion by subtracting the estimate of the rigid-body motion from the observed anisotropic temperature factor. We also subjected the TLS model to estimate the internal motion of the bovine pancreatic PLA(2) using the anisotropic temperature factor deposited in the Protein Data Bank. Both results indicate that the localization of regions exhibiting larger internal motion in the bacterial PLA(2) is almost the same as that in the bovine pancreatic PLA(2), suggesting that although the tertiary structure of the bacterial PLA(2) is strikingly different from that of the bovine pancreatic PLA(2), the internal motion, which is associated with the calcium(II) ion-binding, phospholipid-binding, and allosteric interfacial activation, is commonly observed in both PLA(2)s.     

Ultrahigh resolution drug design. II. Atomic resolution structures of human aldose reductase holoenzyme complexed with Fidarestat and Minalrestat: implications for the binding of cyclic imide inhibitors

The X-ray structures of human aldose reductase holoenzyme in complex with the inhibitors Fidarestat (SNK-860) and Minalrestat (WAY-509) were determined at atomic resolutions of 0.92 A and 1.1 A, respectively. The hydantoin and succinimide moieties of the inhibitors interacted with the conserved anion-binding site located between the nicotinamide ring of the coenzyme and active site residues Tyr48, His110, and Trp111. Minalrestat's hydrophobic isoquinoline ring was bound in an adjacent pocket lined by residues Trp20, Phe122, and Trp219, with the bromo-fluorobenzyl group inside the specificity pocket. The interactions between Minalrestat's bromo-fluorobenzyl group and the enzyme include the stacking against the side-chain of Trp111 as well as hydrogen bonding distances with residues Leu300 and Thr113. The carbamoyl group in Fidarestat formed a hydrogen bond with the main-chain nitrogen atom of Leu300. The atomic resolution refinement allowed the positioning of hydrogen atoms and accurate determination of bond lengths of the inhibitors, coenzyme NADP+ and active-site residue His110. The 1'-position nitrogen atom in the hydantoin and succinimide moieties of Fidarestat and Minalrestat, respectively, form a hydrogen bond with the Nepsilon2 atom of His 110. For Fidarestat, the electron density indicated two possible positions for the H-atom in this bond. Furthermore, both native and anomalous difference maps indicated the replacement of a water molecule linked to His110 by a Cl-ion. These observations suggest a mechanism in which Fidarestat is bound protonated and becomes negatively charged by donating the proton to His110, which may have important implications on drug design.     

Variable metallation of human superoxide dismutase: atomic resolution crystal structures of Cu-Zn, Zn-Zn and as-isolated wild-type enzymes

Human Cu-Zn superoxide dismutase (SOD1) protects cells from the effects of oxidative stress. Mutations in SOD1 are linked to the familial form of amyotrophic lateral sclerosis. Several hypotheses for their toxicity involve the mis-metallation of the enzyme. We present atomic-resolution crystal structures and biophysical data for human SOD1 in three metallation states: Zn-Zn, Cu-Zn and as-isolated. These data represent the first atomic-resolution structures for human SOD1, the first structure of a reduced SOD1, and the first structure of a fully Zn-substituted SOD1 enzyme. Recombinantly expressed as-isolated SOD1 contains a mixture of Zn and Cu at the Cu-binding site. The Zn-Zn structure appears to be at least as stable as the correctly (Cu-Zn) metallated enzyme. These data raise the possibility that in a cellular environment with low availability of free copper, Zn-Zn may be the preferred metallation state of SOD1 prior to its interaction with the copper chaperone.     

Hepatocyte uptake of alpha 1-proteinase inhibitor-trypsin complexes in vitro: evidence for a shared uptake mechanism for proteinase complexes of alpha 1-proteinase inhibitor and antithrombin III

In vivo clearance studies have indicated that the clearance of proteinase complexes of the homologous serine proteinase inhibitors alpha 1-proteinase inhibitor and antithrombin III occurs via a specific and saturable pathway located on hepatocytes. In vitro hepatocyte-uptake studies with antithrombin III-proteinase complexes confirmed the hepatocyte uptake and degradation of these complexes, and demonstrated the formation of a disulfide interchange product between the ligand and a cellular protein. We now report the results of in vitro hepatocyte uptake studies with alpha 1-proteinase inhibitor-trypsin complexes. Trypsin complexes of alpha 1-proteinase inhibitor were prepared and purified to homogeneity. Uptake of these complexes by hepatocytes was time and concentration-dependent. Competition experiments with alpha 1-proteinase inhibitor, alpha 1-proteinase inhibitor-trypsin, and antithrombin III-thrombin indicated that the proteinase complexes of these two inhibitors are recognized by the same uptake mechanism, whereas the native inhibitor is not. Uptake studies were performed at 37 degrees C with 125I-alpha 1-proteinase inhibitor-trypsin and analyzed by sodium dodecyl sulfate-gel electrophoresis in conjunction with autoradiography. These studies demonstrated time-dependent uptake and degradation of the ligand to low molecular weight peptides. In addition, there was a time-dependent accumulation of a high molecular weight complex of ligand and a cellular protein. This complex disappeared when gels were performed under reducing conditions. The sole cysteine residue in alpha 1-proteinase inhibitor was reduced and alkylated with iodoacetamide. Trypsin complexes of the modified inhibitor were prepared and purified to homogeneity. Uptake and degradation studies demonstrated no differences in the results obtained with this modified complex as compared to unmodified alpha 1-proteinase inhibitor-trypsin complex. In addition, the high molecular weight disulfide interchange product was still present on sodium dodecyl sulfate-polyacrylamide gel electrophoresis of solubilized cells. Clearance and clearance competition studies with alpha 1-proteinase inhibitor-trypsin, alkylated alpha 1-proteinase inhibitor-trypsin, antithrombin III-thrombin, and anti-thrombin III-factor IXa further demonstrated the shared hepatocyte uptake mechanism for all these complexes.     

Towards atomic resolution with crystals grown in gel: the case of thaumatin seen at room temperature

One reason for introducing a gel in the crystallization medium of proteins is its ability to reduce convection in solution. This can lead to better nucleation and growth conditions, and to crystals having enhanced diffraction properties. We report here the X-ray characterization at room temperature of high-quality crystals of the intensely sweet thaumatin prepared in a sodium tartrate solution gelified with 0.15% (m/v) agarose. Using a synchrotron radiation, these crystals diffracted to a previously unachieved resolution. A diffraction dataset was collected from four crystals at a resolution of 1.2 A with a R(sym) of 3.6% and a completeness of 99%. Refinement was carried out to a final crystallographic R-factor of 12.0%. The quality of the electron density map allowed for the observation of fine structural details in the protein and its solvation shell. Crystallization in gel might be used more generally to improve the quality of macromolecular crystals. Advantages provided by the gelified medium in the frame of structural studies are emphasized.     

High resolution crystal structures of triosephosphate isomerase complexed with its suicide inhibitors: The conformational flexibility of the catalytic glutamate in its closed, liganded active site

The key residue of the active site of triosephosphate isomerase (TIM) is the catalytic glutamate, which is proposed to be important (i) as a catalytic base, for initiating the reaction, as well as (ii) for the subsequent proton shuttling steps. The structural properties of this glutamate in the liganded complex have been investigated by studying the high resolution crystal structures of typanosomal TIM, complexed with three suicide inhibitors: (S)-glycidol phosphate ((S)-GOP, at 0.99 Å resolution), (R)-glycidol phosphate, ((R)-GOP, at 1.08 Å resolution), and bromohydroxyacetone phosphate (BHAP, at 1.97 Å resolution). The structures show that in the (S)-GOP active site this catalytic glutamate is in the well characterized, competent conformation. However, an unusual side chain conformation is observed in the (R)-GOP and BHAP complexes. In addition, Glu97, salt bridged to the catalytic lysine in the competent active site, adopts an unusual side chain conformation in these two latter complexes. The higher chemical reactivity of (S)-GOP compared with (R)-GOP, as known from solution studies, can be understood: the structures indicate that in the case of (S)-GOP, Glu167 can attack the terminal carbon of the epoxide in a stereoelectronically favored, nearly linear O–C–O arrangement, but this is not possible for the (R)-GOP isomer. These structures confirm the previously proposed conformational flexibility of the catalytic glutamate in its closed, liganded state. The importance of this conformational flexibility for the proton shuttling steps in the TIM catalytic cycle, which is apparently achieved by a sliding motion of the side chain carboxylate group above the enediolate plane, is also discussed.     

Clitocypin, a new cysteine proteinase inhibitor, is monomeric: impact on the mechanism of folding

The molecular mass of clitocypin, a new type of cysteine proteinase inhibitor from the mushroom Clitocybe nebularis, has been determined by analytical ultracentrifugation and gel exclusion chromatography. The result is in agreement with the formula mass of 16.8 kDa, demonstrating that the inhibitor is a monomer in aqueous solution. This enables the kinetics of unfolding and refolding to be interpreted in terms of folding in a kinetically two state, highly cooperative transition from the thermally unfolded state.     

Investigating the effects of double mutation C30A/C75A on onconase structure: Studies at atomic resolution

The structure of onconase C30A/C75A double mutant has been determined at 1.12Å resolution. The structure has high structural homology to other onconase structures. The changes being results of mutation are relatively small, distributed asymmetrically around the two mutated positions, and they are observed not only in the mutation region but expanded to entire molecule. Different conformation of Lys31 side chain that influences the hydrogen bonding network around catalytic triad is probably responsible for lower catalytic efficiency of double mutant. The decrease in thermal stability observed for the onconase variant might be explained by a less dense packing as manifested by the increase of the molecular volume and the solvent accessible surface area. © 2013 Wiley Periodicals, Inc. Biopolymers 101: 454–460, 2014.