Features and structure of a cold active N-acetylneuraminate lyase

N-acetylneuraminate lyases (NALs) are enzymes that catalyze the reversible cleavage and synthesis of sialic acids. They are therefore commonly used for the production of these high-value sugars. This study presents the recombinant production, together with biochemical and structural data, of the NAL from the psychrophilic bacterium Aliivibrio salmonicida LFI1238 (AsNAL). Our characterization shows that AsNAL possesses high activity and stability at alkaline pH. We confirm that these properties allow for the use in a one-pot reaction at alkaline pH for the synthesis of N-acetylneuraminic acid (Neu5Ac, the most common sialic acid) from the inexpensive precursor N-acetylglucosamine. We also show that the enzyme has a cold active nature with an optimum temperature for Neu5Ac synthesis at 20°C. The equilibrium constant for the reaction was calculated at different temperatures, and the formation of Neu5Ac acid is favored at low temperatures, making the cold active enzyme a well-suited candidate for use in such exothermic reactions. The specific activity is high compared to the homologue from Escherichia coli at three tested temperatures, and the enzyme shows a higher catalytic efficiency and turnover number for cleavage at 37°C. Mutational studies reveal that amino acid residue Asn 168 is important for the high k_cat. The crystal structure of AsNAL was solved to 1.65 Å resolution and reveals a compact, tetrameric protein similar to other NAL structures. The data presented provides a framework to guide further optimization of its application in sialic acid production and opens the possibility for further design of the enzyme.

Characterization of the sialic acid synthase from Aliivibrio salmonicida suggests a novel pathway for bacterial synthesis of 7-O-acetylated sialic acids

Resolving the enzymatic pathways leading to sialic acids (Sias) in bacteria are vitally important for understanding their roles in pathogenesis and for subsequent development of tools to combat infections. A detailed characterization of the involved enzymes is also essential due to the highly applicable properties of Sias, i.e., as used in a wide range of medical applications and human nutrition. Bacterial strains that produce Sias display them mainly on their cell surface to mimic animal cells thereby evading the host's immune system. Despite several studies, little is known about the virulence mechanisms of the fish pathogen Aliivibrio salmonicida. The genome of A. salmonicida LFI1238 contains a gene cluster homologous to the Escherichia coli neuraminic acid (Neu) gene cluster involved in biosynthesis of Sias found in the E. coli capsule. This cluster is probably responsible for the biosynthesis of Neu found in A. salmonicida. In this work, we have produced and characterized the sialic acid (Sia) synthase NeuB1, the key enzyme in the pathway. The Sia synthase is an enzyme producing N-acetylneuraminic acid by the condensation of N-acetylmannosamine and phosphoenolpyruvate. Genome content, kinetic data obtained, together with structural considerations, have led us to the prediction that the substrate for NeuB1 from A. salmonicida, E. coli and Streptococcus agalactiae among others, is 4-O-acetyl-N-acetylmannosamine. This means that the product of its enzymatic reaction is 7-O-acetyl-N-acetylneuraminic acid. We propose a pathway for production of this Sia in A. salmonicida, and present evidence for the presence of diacetylated Neu in the bacterium.

Thermodynamics and structure of a salmon cold active goose-type lysozyme

Atlantic salmon goose-type lysozyme (SalG) was previously shown to display features of cold-adaptation as well as renaturation following heat treatment. In this study differential scanning calorimetry (DSC) was carried out to investigate unfolding and potential refolding, while X-ray crystallography was used to study structural factors contributing to the temperature-related characteristics. The recombinant SalG has a melting temperature (T(m)) of 36.8 degrees C under thermal denaturation conditions and regains activity after returning to permissive (low) temperature. Furthermore, refolding is dramatically reduced in solutions with high SalG concentrations, coupled with significant protein precipitation. The structural features of SalG closely resemble those of other g-type lysozymes. However, the N-terminal region of SalG is less anchored to the rest of the molecule due to the absence of disulphide bonds, thus, contributing significantly to the low T(m) of SalG. The absence of disulphide bonds and the distribution of salt bridges may at the same time ease refolding leading to renaturation.

Structure of uracil-DNA N-glycosylase (UNG) from Vibrio cholerae: mapping temperature adaptation through structural and mutational analysis

The crystal structure of Vibrio cholerae uracil-DNA N-glycosylase (vcUNG) has been determined to 1.5 A resolution. Based on this structure, a homology model of Aliivibrio salmonicida uracil-DNA N-glycosylase (asUNG) was built. A previous study demonstrated that asUNG possesses typical cold-adapted features compared with vcUNG, such as a higher catalytic efficiency owing to increased substrate affinity. Specific amino-acid substitutions in asUNG were suggested to be responsible for the increased substrate affinity and the elevated catalytic efficiency by increasing the positive surface charge in the DNA-binding region. The temperature adaptation of these enzymes has been investigated using structural and mutational analyses, in which mutations of vcUNG demonstrated an increased substrate affinity that more resembled that of asUNG. Visualization of surface potentials revealed a more positive potential for asUNG compared with vcUNG; a modelled double mutant of vcUNG had a potential around the substrate-binding region that was more like that of asUNG, thus rationalizing the results obtained from the kinetic studies.

Structure of a highly stable mutant of human fibroblast growth factor 1

Fibroblast growth factors (FGFs) are involved in diverse cellular processes such as cell migration, angiogenesis, osteogenesis, wound healing and embryonic and foetal development. Human acidic fibroblast growth factor (FGF-1) is the only member of the FGF family that binds with high affinity to all four FGF receptors and thus is considered to be the human mitogen with the broadest specificity. However, pharmacological applications of FGF-1 are limited owing to its low stability. It has previously been reported that the introduction of single mutations can significantly improve the stability of FGF-1 and its resistance to proteolytic degradation. Here, the structure of the Q40P/S47I/H93G triple mutant of FGF-1, which exhibits much higher stability, a prolonged half-life and enhanced mitogenic activity, is presented. Compared with the wild-type structure, three localized conformational changes in the stable triple mutant were observed, which is in agreement with the perfect energetic additivity of the single mutations described in a previous study. The huge change in FGF-1 stability (the denaturation temperature increased by 21.5 K, equivalent to DeltaDeltaG(den) = 24.3 kJ mol(-1)) seems to result from the formation of a short 3(10)-helix (position 40), an improvement in the propensity of amino acids to form beta-sheets (position 47) and the rearrangement of a local hydrogen-bond network (positions 47 and 93).

Structural evidence for lack of inhibition of fish goose-type lysozymes by a bacterial inhibitor of lysozyme

It is known that bacteria contain inhibitors of lysozyme activity. The recently discovered Escherichia coli inhibitor of vertebrate lysozyme (Ivy) and its potential interactions with several goose-type (g-type) lysozymes from fish were studied using functional enzyme assays, comparative homology modelling, protein-protein docking, and molecular dynamics simulations. Enzyme assays carried out on salmon g-type lysozyme revealed a lack of inhibition by Ivy. Detailed analysis of the complexes formed between Ivy and both hen egg white lysozyme (HEWL) and goose egg white lysozyme (GEWL) suggests that electrostatic interactions make a dominant contribution to inhibition. Comparison of three dimensional models of aquatic g-type lysozymes revealed important insertions in the beta domain, and specific sequence substitutions yielding altered electrostatic surface properties and surface curvature at the protein-protein interface. Thus, based on structural homology models, we propose that Ivy is not effective against any of the known fish g-type lysozymes. Docking studies suggest a weaker binding mode between Ivy and GEWL compared to that with HEWL, and our models explain the mechanistic necessity for conservation of a set of residues in g-type lysozymes as a prerequisite for inhibition by Ivy.

Structural adaptation of endonuclease I from the cold-adapted and halophilic bacterium Vibrio salmonicida

The crystal structure of the periplasmic/extracellular endonuclease I from Vibrio salmonicida has been solved to 1.5 A resolution and, in comparison to the corresponding endonucleases from V. cholerae and V. vulnificus, serves as a model system for the investigation of the structural determinants involved in the temperature and NaCl adaptation of this enzyme class. The overall fold of the three enzymes is essentially similar, but the V. salmonicida endonuclease displays a significantly more positive surface potential than the other two enzymes owing to the presence of ten more Lys residues. However, if the optimum salt concentrations for the V. salmonicida and V. cholerae enzymes are taken into consideration in the electrostatic surface-potential calculation, the potentials of the two enzymes become surprisingly similar. The higher number of basic residues in the V. salmonicida protein is therefore likely to be a result, at least in part, of adaptation to the more saline habitat of V. salmonicida (seawater) than V. cholerae (brackish water). The hydrophobic core of all three enzymes is almost identical, but the V. salmonicida endonuclease has a slightly lower number of internal hydrogen bonds. This, together with repulsive forces between the basic residues on the protein surface of V. salmonicida endonuclease I and differences in the distribution of salt bridges, probably results in higher flexibility of regions of the V. salmonicida protein. This is likely to influence both the catalytic activity and the stability of the protein.

Electrostatic interactions play an essential role in DNA repair and cold-adaptation of uracil DNA glycosylase

Life has adapted to most environments on earth, including low and high temperature niches. The increased catalytic efficiency and thermoliability observed for enzymes from organisms living in constantly cold regions when compared to their mesophilic and thermophilic cousins are poorly understood at the molecular level. Uracil DNA glycosylase (UNG) from cod (cUNG) catalyzes removal of uracil from DNA with an increased k(cat) and reduced K(m) relative to its warm-active human (hUNG) counterpart. Specific issues related to DNA repair and substrate binding/recognition (K(m)) are here investigated by continuum electrostatics calculations, MD simulations and free energy calculations. Continuum electrostatic calculations reveal that cUNG has surface potentials that are more complementary to the DNA potential at and around the catalytic site when compared to hUNG, indicating improved substrate binding. Comparative MD simulations combined with free energy calculations using the molecular mechanics-Poisson Boltzmann surface area (MM-PBSA) method show that large opposing energies are involved when forming the enzyme-substrate complexes. Furthermore, the binding free energies obtained reveal that the Michaelis-Menten complex is more stable for cUNG, primarily due to enhanced electrostatic properties, suggesting that energetic fine-tuning of electrostatics can be utilized for enzymatic temperature adaptation. Energy decomposition pinpoints the residual determinants responsible for this adaptation.

Sequence comparison and environmental adaptation of a bacterial endonuclease

The periplasmic/extracellular bacterial enzyme endonuclease I was chosen as a model system to identify features that might be responsible for temperature- and salt adaptation. A statistical study of amino acid sequence properties belonging to endonuclease I enzymes from three mesophilic habitats (non-marine, brackish water and marine), and three marine temperature groups (psychrophile, intermediate and mesophile) has been conducted. Ten new endonuclease I genes have been sequenced in order to increase the sample size. A bioinformatical method of property dependent statistical analysis of alignments has been applied. To our knowledge this is the first time these methods have been used in order to investigate environmental adaptation of enzymes. Adaptation to low temperature seems to involve increased surface isoelectric point and hydrophobicity in contrast to salt adaptation in which the isoelectric point and hydrophobicity at the surface decreases. Redistribution of charge and hydrophobicity might be the most important signature for cold adaptation and salt adaptation of this enzyme class. The results indicate that general trends of adaptation are possible to elucidate from the amino acid sequences. Also in this paper a new scale of stratified B-factors, derived from the Protein Data Bank, is presented.

Comparative studies of endonuclease I from cold-adapted Vibrio salmonicida and mesophilic Vibrio cholerae

Endonuclease I is a periplasmic or extracellular enzyme present in many different Proteobacteria. The endA gene encoding endonuclease I from the psychrophilic and mildly halophilic bacterium Vibrio salmonicida and from the mesophilic brackish water bacterium Vibrio cholerae have been cloned, over-expressed in Escherichia coli, and purified. A comparison of the enzymatic properties shows large differences in NaCl requirements, optimum pH, temperature stability and catalytic efficiency of the two proteins. The V. salmonicida EndA shows typical cold-adapted features such as lower unfolding temperature, lower temperature optimum for activity, and higher specific activity than V. cholerae EndA. The thermodynamic activation parameters confirm the psychrophilic nature of V. salmonicida EndA with a much lower activation enthalpy. The optimal conditions for enzymatic activity coincide well with the corresponding optimal requirements for growth of the organisms, and the enzymes function predominantly as DNases at physiological concentrations of NaCl. The periplasmic or extracellular localization of the enzymes, which renders them constantly exposed to the outer environment of the cell, may explain this fine-tuning of biochemical properties.

The first structure of a cold-active catalase fromVibrio salmonicidaat 1.96 Å reveals structural aspects of cold adaptation

The cold-adapted catalase from the fish-pathogenic bacterium Vibrio salmonicida (VSC) has recently been characterized and shown to be two times more catalytically efficient compared with catalase from the mesophilic human pathogen Proteus mirabilis [PMC; Lorentzen et al. (2006), Extremophiles, 10, 427-440]. VSC is also less temperature-stable, with a half-life of 5 min at 333 K compared with 50 min for PMC. This was the background for solving the crystal structure of the cold-adapted VSC to 1.96 A and performing an extensive structural comparison of VSC and PMC. The comparison revealed that the entrance (the major channel) leading to the catalytically essential haem group, is locally more flexible and slightly wider in VSC. This might explain the enhanced catalytic efficiency of the nearly diffusion-controlled degradation of hydrogen peroxide into water and molecular oxygen in VSC. The reduced thermal stability of the cold-adapted VSC may be explained by a reduced number of ion-pair networks. The four C-terminal alpha-helices are displaced in the structures, probably owing to missing ionic interactions in VSC compared with PMC, and this is postulated as an initiation site for unfolding the cold-adapted enzyme. VSC is the first crystal structure reported of a cold-adapted monofunctional haem-containing catalase.

Effect of fish skin mucus on the soluble proteome of Vibrio salmonicida analysed by 2-D gel electrophoresis and tandem mass spectrometry

Vibrio salmonicida is the causative agent of cold-water vibriosis in farmed marine fish species. Adherence of pathogenic bacteria to mucosal surfaces is considered to be the first steps in the infective processes, and proteins involved are regarded as virulence factors. The global protein expression profile of V. salmonicida, grown with and without the presence of fish skin mucus in the synthetic media, was compared. Increased levels of proteins involved in motility, oxidative stress responses, and general stress responses were demonstrated as an effect of growth in the presence of mucus compared to non-mucus containing media. Enhanced levels of the flagellar proteins FlaC, FlaD and FlaE indicate increased motility capacity, while enhanced levels of the heat shock protein DnaK and the chaperonin GroEL indicate a general stress response. In addition, we observed that peroxidases, TPx.Grx and AhpC, involved in the oxidative stress responses, were induced by mucus proteins. The addition of mucus to the culture medium did not significantly alter the growth rate of V. salmonicida. An analysis of mucus proteins suggests that the mucus layer harbours a protein species that potentially possesses catalytic activity against DNA, and a protein with iron chelating activity. This study represents the first V. salmonicida proteomic analysis, and provides specific insight into the proteins necessary for the bacteria to challenge the skin mucus barrier of the fish.

The structure of Vibrio cholerae extracellular endonuclease I reveals the presence of a buried chloride ion

The crystal structure of a periplasmic/extracellular endonuclease from Vibrio cholerae has been solved at low and at neutral pH. Crystals grown at pH 4.6 and 6.9 diffracted to 1.6 A (on BM01A at the ESRF) and 1.95 A (on a rotating-anode generator), respectively. The structures of the endonuclease were compared with the structure of a homologous enzyme in V. vulnificus. The structures of the V. cholerae enzyme at different pH values are essentially identical to each other and to the V. vulnificus enzyme. However, interesting features were observed in the solvent structures. Both V. cholerae structures reveal the presence of a chloride ion completely buried within the core of the protein, with the nearest solvent molecule approximately 7 A away. Magnesium, which is essential for catalysis, is present in the structure at neutral pH, but is absent at low pH, and may partly explain the inactivity of the enzyme at lower pH.

Free energy calculations show that acidic P1 variants undergo large pKa shifts upon binding to trypsin

Serine proteinases and their protein inhibitors belong to one of the most comprehensively studied models of protein-protein interactions. It is well established that the narrow trypsin specificity is caused by the presence of a negatively charged aspartate at the specificity pocket. X-ray crystallography as well as association measurements revealed, surprisingly, that BPTI with glutamatic acid as the primary binding (P1) residue was able to bind to trypsin. Previous free energy calculations showed that there was a substantially unfavorable binding free energy associated with accommodation of ionized P1 Glu at the S1-site of trypsin. In this study, the binding of P1 Glu to trypsin has been systematically investigated in terms of the protonation states of P1 Glu and Asp189, the orientation of Gln192, as well as the possible presence of counterions using the linear interaction energy (LIE) approach and the free energy perturbation (FEP) method. Twenty-four conceivable binding arrangements were evaluated and quantitative agreement with experiments is obtained when the P1 Glu binds in its protonated from. The results suggest that P1 Glu is one of the variants of BPTI that inhibit trypsin strongest at low pH, contrary to the specificity profile of trypsin, suggesting a new regulation mechanism of trypsin-like enzymes.

Predicting proteinase specificities from free energy calculations

The role of the primary binding residue (P1) in complexes between three different subtilases (subtilisin Carlsberg, thermitase and proteinase K) and their canonical protein inhibitor eglin c have been studied by free energy calculations. Based on the crystal structures of eglin c in complex with subtilisin Carlsberg and thermitase, and a homology model of the eglin c-proteinase K complex, a total of 57 mutants have been constructed and docked into their host proteins. The binding free energy was then calculated using molecular dynamics (MD) simulations combined with the linear interaction energy (LIE) method for all complexes differing only in the nature of the amino acid at the P1 position. LIE calculations for 19 different complexes for each subtilase were thus carried out excluding proline. The effects of substitutions at the P1 position on the binding free energies are found to be very large, and positively charged residues (Arg, Lys and His) are particularly deleterious for all three enzymes. The charged variants of the acidic side chains are found to bind more favorably as compared to their protonated states in all three subtilases. Furthermore, hydrophobic amino acids are accommodated most favorably at the S1-site in all three enzymes. Comparison of the three series of binding free energies shows only minor differences in the 19 computed relative binding free energies among these subtilases. This is further reflected in the correlation coefficient between the 23 relative binding free energies obtained, including the possible protonation states of ionizable side chains, but excluding the P1 Pro, for subtilisin Carlsberg versus thermitase (0.95), subtilisin versus proteinase K (0.94) and thermitase versus proteinase K (0.96).

The Crystal Structure of Mismatch-specific Uracil-DNA Glycosylase (MUG) from Deinococcus radiodurans Reveals a Novel Catalytic Residue and Broad Substrate Specificity

Deinococcus radiodurans is extremely resistant to the effects of ionizing radiation. The source of the radiation resistance is not known, but an expansion of specific protein families related to stress response and damage control has been observed. DNA repair enzymes are among the expanded protein families in D. radiodurans, and genes encoding five different uracil-DNA glycosylases are identified in the genome. Here we report the three-dimensional structure of the mismatch-specific uracil-DNA glycosylase (MUG) from D. radiodurans (drMUG) to a resolution of 1.75 angstroms. Structural analyses suggest that drMUG possesses a novel catalytic residue, Asp-93. Activity measurements show that drMUG has a modified and broadened substrate specificity compared with Escherichia coli MUG. The importance of Asp-93 for activity was confirmed by structural analysis and abolished activity for the mutant drMUGD93A. Two other microorganisms, Bradyrhizobium japonicum and Rhodopseudomonas palustris, possess genes that encode MUGs with the highest sequence identity to drMUG among all of the bacterial MUGs examined. A phylogenetic analysis indicates that these three MUGs form a new MUG/thymidine-DNA glycosylase subfamily, here called the MUG2 family. We suggest that the novel catalytic residue (Asp-93) has evolved to provide drMUG with broad substrate specificity to increase the DNA repair repertoire of D. radiodurans.

Probing the effect of point mutations at protein-protein interfaces with free energy calculations

We have studied the effect of point mutations of the primary binding residue (P1) at the protein-protein interface in complexes of chymotrypsin and elastase with the third domain of the turkey ovomucoid inhibitor and in trypsin with the bovine pancreatic trypsin inhibitor, using molecular dynamics simulations combined with the linear interaction energy (LIE) approach. A total of 56 mutants have been constructed and docked into their host proteins. The free energy of binding could be reliably calculated for 52 of these mutants that could unambiguously be fitted into the binding sites. We find that the predicted binding free energies are in very good agreement with experimental data with mean unsigned errors between 0.50 and 1.03 kcal/mol. It is also evident that the standard LIE model used to study small drug-like ligand binding to proteins is not suitable for protein-protein interactions. Three different LIE models were therefore tested for each of the series of protein-protein complexes included, and the best models for each system turn out to be very similar. The difference in parameterization between small drug-like compounds and protein point mutations is attributed to the preorganization of the binding surface. Our results clearly demonstrate the potential of free energy calculations for probing the effect of point mutations at protein-protein interfaces and for exploring the principles of specificity of hot spots at the interface.

Structure of the uracil-DNAN-glycosylase (UNG) fromDeinococcus radiodurans

Uracil-DNA glycosylases are DNA-repair enzymes that catalyse the removal of promutagenic uracil from single- and double-stranded DNA, thereby initiating the base-excision repair (BER) pathway. Uracil in DNA can occur by mis-incorporation of dUMP in place of dTMP during DNA synthesis or by deamination of cytosine, resulting in U-A or U-G mispairs. The radiation-resistant bacterium Deinococcus radiodurans has an elevated number of uracil-DNA glycosylases compared with most other organisms. The crystal structure of dr0689 (uracil-DNA N-glycosylase), which has been shown to be the major contributor to the removal of mis-incorporated uracil bases in crude cell extracts of D. radiodurans, is reported.

Increased Flexibility as a Strategy for Cold Adaptation

Uracil DNA glycosylase (UDG) is a DNA repair enzyme in the base excision repair pathway and removes uracil from the DNA strand. Atlantic cod UDG (cUDG), which is a cold-adapted enzyme, has been found to be up to 10 times more catalytically active in the temperature range 15-37 degrees C as compared with the warm-active human counterpart. The increased catalytic activity of cold-adapted enzymes as compared with their mesophilic homologues are partly believed to be caused by an increase in the structural flexibility. However, no direct experimental evidence supports the proposal of increased flexibility of cold-adapted enzymes. We have used molecular dynamics simulations to gain insight into the structural flexibility of UDG. The results from these simulations show that an important loop involved in DNA recognition (the Leu(272) loop) is the most flexible part of the cUDG structure and that the human counterpart has much lower flexibility in the Leu(272) loop. The flexibility in this loop correlates well with the experimental k(cat)/K(m) values. Thus, the data presented here add strong support to the idea that flexibility plays a central role in adaptation to cold environments.

Structure of anionic salmon trypsin in a second crystal form

Anionic salmon trypsin in a second crystal form (ST-IIB) has been refined at 1.83 A, resolution. The crystals are orthorhombic and belong to space group P2(1)2(1)2 with lattice parameters a = 77.09, b = 82.33 and c = 31.16 A. The present structure has been compared to salmon trypsin as it appears in a previously reported crystal form (ST-IIA) with cell dimensions a = 61.95, b = 84.33 and c = 39.11 A [Smalås & Hordvik (1993). Acta Cryst. D49, 318-330]. The presence of a sulfate group involved in several hydrogen bonds to active-site residues, and the location of an additional benzamidine site in the crystal lattice, are the most striking differences between the present and the previous structure. Superposition of main-chain atoms in the two structures give an overall r.m.s. difference of 0.26 A, with the main differences located to areas with different molecular packing. The overall coordinate error is estimated to be between 0.20 and 0.25 A, by the method of Luzzati.

Crystal structures of five bovine chymotrypsin complexes with P1 BPTI variants

The bovine chymotrypsin-bovine pancreatic trypsin inhibitor (BPTI) interaction belongs to extensively studied models of protein-protein recognition. The accommodation of the inhibitor P1 residue in the S1 binding site of the enzyme forms the hot spot of this interaction. Mutations introduced at the P1 position of BPTI result in a more than five orders of magnitude difference of the association constant values with the protease. To elucidate the structural aspects of the discrimination between different P1 residues, crystal structures of five bovine chymotrypsin-P1 BPTI variant complexes have been determined at pH 7.8 to a resolution below 2 A. The set includes polar (Thr), ionizable (Glu, His), medium-sized aliphatic (Met) and large aromatic (Trp) P1 residues and complements our earlier studies of the interaction of different P1 side-chains with the S1 pocket of chymotrypsin. The structures have been compared to the complexes of proteases with similar and dissimilar P1 preferences, including Streptomyces griseus proteases B and E, human neutrophil elastase, crab collagenase, bovine trypsin and human thrombin. The S1 sites of these enzymes share a common general shape of significant rigidity. Large and branched P1 residues adapt in their complexes similar conformations regardless of the polarity and size differences between their S1 pockets. Conversely, long and flexible residues such as P1 Met are present in the disordered form and display a conformational diversity despite similar inhibitory properties with respect to most enzymes studied. Thus, the S1 specificity profiles of the serine proteases appear to result from the precise complementarity of the P1-S1 interface and minor conformational adjustments occurring upon the inhibitor binding.

Trypsin specificity as elucidated by LIE calculations, X-ray structures, and association constant measurements

The variation in inhibitor specificity for five different amine inhibitors bound to CST, BT, and the cold-adapted AST has been studied by use of association constant measurements, structural analysis of high-resolution crystal structures, and the LIE method. Experimental data show that AST binds the 1BZA and 2BEA inhibitors 0.8 and 0.5 kcal/mole more strongly than BT. However, structural interactions and orientations of the inhibitors within the S1 site have been found to be virtually identical in the three enzymes studied. For example, the four water molecules in the inhibitor-free structures of AST and BT are channeled into similar positions in the S1 site, and the nitrogen atom(s) of the inhibitors are found in two cationic binding sites denoted Position1 and Position2. The hydrophobic binding contributions for all five inhibitors, estimated by the LIE calculations, are also in the same order (-2.1 +/- 0.2 kcal/mole) for all three enzymes. Our hypothesis is therefore that the observed variation in inhibitor binding arises from different electrostatic interactions originating from residues outside the S1 site. This is well illustrated by AST, in which Asp 150 and Glu 221B, despite some distance from the S1 binding site, lower the electrostatic potential of the S1 site and thus enhance substrate binding. Because the trends in the experimentally determined binding energies were reproduced by the LIE calculations after adding the contribution from long-range interactions, we find this method very suitable for rational studies of protein-substrate interactions.

Structural consequences of accommodation of four non-cognate amino acid residues in the S1 pocket of bovine trypsin and chymotrypsin

Crystal structures of P1 Gly, Val, Leu and Phe bovine pancreatic trypsin inhibitor (BPTI) variants in complex with two serine proteinases, bovine trypsin and chymotrypsin, have been determined. The association constants for the four mutants with the two enzymes show that the enlargement of the volume of the P1 residue is accompanied by an increase of the binding energy, which is more pronounced for bovine chymotrypsin. Since the conformation of the P1 side-chains in the two S1 pockets is very similar, we suggest that the difference in DeltaG values between the enzymes must arise from the more polar environment of the S1 site of trypsin. This results mainly from the substitutions of Met192 and Ser189 observed in chymotrypsin with Gln192 and Asp189 present in trypsin. The more polar interior of the S1 site of trypsin is reflected by a much higher order of the solvent network in the empty pocket of the enzyme, as is observed in the complexes of the two enzymes with the P1 Gly BPTI variant. The more optimal binding of the large hydrophobic P1 residues by chymotrypsin is also reflected by shrinkage of the S1 pocket upon the accommodation of the cognate residues of this enzyme. Conversely, the S1 pocket of trypsin expands upon binding of such side-chains, possibly to avoid interaction with the polar residues of the walls. Further differentiation between the two enzymes is achieved by small differences in the shape of the S1 sites, resulting in an unequal steric hindrance of some of the side-chains, as observed for the gamma-branched P1 Leu variant of BPTI, which is much more favored by bovine chymotrypsin than trypsin. Analysis of the discrimination of beta-branched residues by trypsin and chymotrypsin is based on the complexes with the P1 Val BPTI variant. Steric repulsion of the P1 Val residue by the walls of the S1 pocket of both enzymes prevents the P1 Val side-chain from adopting the most optimal chi1 value.

The structure of uracil-DNA glycosylase from Atlantic cod (Gadus morhua) reveals cold-adaptation features

Uracil-DNA glycosylase (UDG; EC 3.2.2.3) is a DNA-repair protein that catalyses the hydrolysis of promutagenic uracil residues from single- or double-stranded DNA, generating free uracil and abasic DNA. The crystal structure of the catalytic domain of cod uracil-DNA glycosylase (cUDG) has been determined to 1.9 A resolution, with final R factors of 18.61 and 20.57% for the working and test sets of reflections, respectively. This is the first crystal structure of a uracil-DNA glycosylase from a cold-adapted species and a detailed comparison with the human enzyme is performed in order to rationalize the cold-adapted behaviour of the cod enzyme at the structural level. The catalytic domain of cUDG comprises 223 residues, with a sequence identity to the human UDG of 75%. The tertiary structures of the two enzymes are also similar, with an overall displacement in main-chain atomic positions of 0.63 A. The amino-acid substitutions and the differences in intramolecular hydrogen bonds, hydrophobic interactions, ion-pair interactions and electrostatic potentials are compared and discussed in order to gain insight into the factors that cause the increased activity and reduced thermostability of the cod enzyme. In particular, the reduced number of strong ion-pair interactions in the C-terminal half of cUDG is believed to greatly affect the flexibility and/or stability. Increased positive electrostatic surface potential on the DNA-facing side of cUDG seems to be responsible for increasing the affinity for the negatively charged DNA compared with that of hUDG.

Identification, cloning, and expression of uracil-DNA glycosylase from Atlantic cod (Gadus morhua): characterization and homology modeling of the cold-active catalytic domain

Two distinct forms of the highly conserved uracil-DNA glycosylase (UNG) have been isolated from Atlantic cod (Gadus morhua) liver cDNA by rapid amplification of cDNA ends (RACE). From the cDNA sequences, both forms were deduced to encode an open reading frame of 301 amino acids, with an identical 267-amino-acid C-terminal region and different N-terminal regions of 34 amino acids. By comparison with the human UNG sequences, the two forms were identified as possible mitochondrial (cUNG1) and nuclear (cUNG2) forms. Several constructs of recombinant cUNG (rcUNG) were expressed in Escherichia coli in order to optimize the yield. The recombinant enzyme was purified to apparent homogeneity as determined by sodium dodecyl sulfate-polyacrylamide gel electrophoresis. Activity and stability experiments showed that rcUNG was similar to cUNG previously purified from Atlantic cod liver, and was more pH- and temperature labile than a recombinant human UNG (rhUNG). Under optimal assay conditions for both rcUNG and rhUNG, the turnover number (k(cat)) was three times higher for rcUNG compared with rhUNG, with an identical K(M), resulting in a threefold higher catalytic efficiency (k(cat)/K(M)) for rcUNG. These activity and stability experiments reveal cold-adapted features in rcUNG. Homology models of the catalytic domains of Atlantic cod (cUNG) and mouse uracil-DNA glycosylase (mUNG) were built using the human UNG (hUNG) crystal structure as a template. The unique amino acid substitutions observed in cod UNG were mainly located in the N- and C-terminal parts of the sequence. The analysis indicated a more stable N-terminal, a more flexible C-terminal, and a less stabilized core in cUNG as compared with the mammalian UNGs. Substitution of several amino acids in or near the DNA-binding site in cUNG could give rise to a more positively charged surface and a higher electrostatic potential near the active site compared with the mammalian UNGs. The higher potential may increase the electrostatic interactions between the enzyme and DNA, and may explain the increased substrate affinity and, in combination with the higher flexibility, the higher catalytic efficiency observed for rcUNG.

Crystallization and preliminary X-ray diffraction analysis of a cold-adapted uracil-DNA glycosylase from Atlantic cod (Gadus morhua)

Uracil-DNA glycosylase (UDG) is a DNA-repair enzyme involved in the removal of uracil from DNA. The Atlantic cod UDG (cUDG) possesses typical cold-adaptation features, with higher catalytic efficiency and lower thermal stability than the mammalian counterparts. cUDG has been crystallized by the vapour-diffusion method using sodium citrate as the precipitant at pH 7.5. The crystals are monoclinic and belong to space group P2(1), with unit-cell parameters a = 68.58, b = 67.19, c = 68.64 A, beta = 119.85 degrees. There are two molecules in the asymmetric unit, with a corresponding V(M) value of 2.71 A(3) Da(-1) and a solvent content of 54.7%. Synchrotron diffraction data have been collected to 1.9 A resolution using cryogenic conditions (120 K).

Computational analysis of binding of P1 variants to trypsin

The binding of P1 variants of bovine pancreatic trypsin inhibitor (BPTI) to trypsin has been investigated by means of molecular dynamics simulations. The specific interaction formed between the amino acid at the primary binding (P1) position of the binding loop of BPTI and the specificity pocket of trypsin was estimated by use of the linear interaction energy (LIE) method. Calculations for 13 of the naturally occurring amino acids at the P1 position were carried out, and the results obtained were found to correlate well with the experimental binding free energies. The LIE calculations rank the majority of the 13 variants correctly according to the experimental association energies and the mean error between calculated and experimental binding free energies is only 0.38 kcal/mole, excluding the Glu and Asp variants, which are associated with some uncertainties regarding protonation and the possible presence of counter-ions. The three-dimensional structures of the complex with three of the P1 variants (Asn, Tyr, and Ser) included in this study have not at present been solved by any experimental techniques and, therefore, were modeled on the basis of experimental data from P1 variants of similar size. Average structures were calculated from the MD simulations, from which specific interactions explaining the broad variation in association energies were identified. The present study also shows that explicit treatment of the complex water-mediated hydrogen bonding network at the protein-protein interface is of crucial importance for obtaining reliable binding free energies. The successful reproduction of relative binding energies shows that this type of methodology can be very useful as an aid in rational design and redesign of biologically active macromolecules.

Electrostatic effects play a central role in cold adaptation of trypsin

Organisms that live in constantly cold environments have to adapt their metabolism to low temperatures, but mechanisms of enzymatic adaptation to cold environments are not fully understood. Cold active trypsin catalyses reactions more efficiently and binds ligands more strongly in comparison to warm active trypsin. We have addressed this issue by means of comparative free energy calculations studying the binding of positively charged ligands to two trypsin homologues. Stronger inhibition of the cold active trypsin by benzamidine and positively charged P1-variants of BPTI is caused by rather subtle electrostatic effects. The different affinity of benzamidine originates solely from long range interactions, while the increased binding of P1-Lys and -Arg variants of BPTI is attributed to both long and short range effects that are enhanced in the cold active trypsin compared to the warm active counterpart. Electrostatic interactions thus provide an efficient strategy for cold adaptation of trypsin.

Environmental fate of chlorinated bornanes estimated by theoretical descriptors

Theoretical molecular descriptors have been calculated for 36 polychlorinated bornanes, the majority compound class of the insecticide Toxaphene. The results demonstrate that thermodynamic stability by the use of molecular structural energies can be used as a general parameter for persistence. Since these descriptors agree well with polychlorinated bornanes found in the environment, these compounds should be included as important indicator compounds in future trace analytical investigations of polychlorinated bornanes and also within experimental metabolism studies to investigate potential toxic metabolites. Reactivity descriptors such as electronaffinity, hardness, LUMO location and atomic charges may guide to potential chemical reactions like the dechlorination of polychlorinated bornanes in reductive environment. Further it is advised to use these descriptors and other new potential ones in combination with experimental degradation and toxicology studies to explore the relationship between molecular structure and biological effects of chlorobornanes.

Atomic resolution structures of trypsin provide insight into structural radiation damage

Radiation damage is an inherent problem in protein X-ray crystallography and the process has recently been shown to be highly specific, exhibiting features such as cleavage of disulfide bonds, decarboxylation of acidic residues, increase in atomic B factors and increase in unit-cell volume. Reported here are two trypsin structures at atomic resolution (1.00 and 0.95 A), the data for which were collected at a third-generation synchrotron (ESRF) at two different beamlines. Both trypsin structures exhibit broken disulfide bonds; in particular, the bond from Cys191 to Cys220 is very sensitive to synchrotron radiation. The data set collected at the most intense beamline (ID14-EH4) shows increased structural radiation damage in terms of lower occupancies for cysteine residues, more breakage in the six disulfide bonds and more alternate conformations. It appears that high intensity and not only the total X-ray dose is most harmful to protein crystals.

Cold adapted enzymes

The number of reports on enzymes from cold adapted organisms has increased significantly over the past years, and reveals that adaptive strategies for functioning at low temperature varies among enzymes. However, the high catalytic efficiency at low temperature seems, for the majority of cold active enzymes, to be accompanied by a reduced thermal stability. Increased molecular flexibility to compensate for the low working temperature, is therefore still the most dominating theory for cold adaptation, although there also seem to be other adaptive strategies. The number of experimentally determined 3D structures of enzymes possessing cold adaptation features is still limited, and restricts a structural rationalization for cold activity. The present summary of structural characteristics, based on comparative studies on crystal structures (7), homology models (7), and amino acid sequences (24), reveals that there are no common structural feature that can account for the low stability, increased catalytic efficiency, and proposed molecular flexibility. Analysis of structural features that are thought to be important for stability (e.g. intra-molecular hydrogen bonds and ion-pairs, proline-, methionine-, glycine-, or arginine content, surface hydrophilicity, helix stability, core packing), indicates that each cold adapted enzyme or enzyme system use different small selections of structural adjustments for gaining increased molecular flexibility that in turn give rise to increased catalytic efficiency and reduced stability. Nevertheless, there seem to be a clear correlation between cold adaptation and reduced number of interactions between structural domains or subunits. Cold active enzymes also seem, to a large extent, to increase their catalytic activity by optimizing the electrostatics at and around the active site.

Substitutions at the P(1) position in BPTI strongly affect the association energy with serine proteinases

The role of the S(1) subsite in trypsin, chymotrypsin and plasmin has been examined by measuring the association with seven different mutants of bovine pancreatic trypsin inhibitor (BPTI); the mutants contain Gly, Ala, Ser, Val, Leu, Arg, and Trp at the P(1) position of the reactive site. The effects of substitutions at the P(1) position on the association constants are very large, comprising seven orders of magnitude for trypsin and plasmin, and over five orders for chymotrypsin. All mutants showed a decrease of the association constant to the three proteinases in the same order: Ala>Gly>Ser>Arg>Val>Leu>Trp. Calorimetric and circular dichroism methods showed that none of the P1 substitutions, except the P1-Val mutant, lead to destabilisation of the binding loop conformation. The X-ray structure of the complex formed between bovine beta-trypsin and P(1)-Leu BPTI showed that the P(1)-Leu sterically conflicts with the side-chain of P(3)-Ile, which thereby is forced to rotate approximately 90 degrees. Ile18 (P(3)) in its new orientation, in turn interacts with the Tyr39 side-chain of trypsin. Introduction of a large side-chain at the P1' position apparently leads to a cascade of small alterations of the trypsin-BPTI interface that seem to destabilise the complex by it adopting a less optimized packing and by tilting the BPTI molecule up to 15 degrees compared to the native trypsin-BPTI complex.

Electrostatics of mesophilic and psychrophilic trypsin isoenzymes: qualitative evaluation of electrostatic differences at the substrate binding site

A qualitative evaluation of electrostatic features of the substrate binding region of seven isoenzymes of trypsin has been performed by using the continuum electrostatic model for the solution of the Poisson-Boltzmann equation. The sources of the electrostatic differences among the trypsins have been sought by comparative calculations on selective charges: all charges, conserved charges, partial charges, unique cold trypsin charges, and a number of charge mutations. As expected, most of the negative potential at the S(1) region of all trypsins is generated from Asp(189), but the potential varies significantly among the seven trypsin isoenzymes. The three cold active enzymes included in this study possess a notably lower potential at and around the S(1)-pocket compared with the warm active counterparts; this finding may be the main contribution to the increased binding affinity. The source of the differences are nonconserved charged residues outside the specificity pocket, producing electric fields at the S(1)-pocket that are different in both sign and magnitude. The surface charges of the mesophilic trypsins generally induce the S(1) pocket positively, whereas surface charges of the cold trypsins produce a negative electric field of this region. Calculations on mutants, where charged amino acids were substituted between the trypsins, showed that mutations in Loop2 (residues 221B and 224) and residue 175, in particular, were responsible for the low potential of the cold enzymes.

Evaluation of protein-protein association energies by free energy perturbation calculations

The association energy upon binding of different amino acids in the specificity pocket of trypsin was evaluated by free energy perturbation calculations on complexes between bovine trypsin (BT) and bovine pancreatic trypsin inhibitor (BPTI). Three simulations of mutations of the primary binding residue (P(1)) were performed (P(1)-Ala to Gly, P(1)-Met to Gly and P(1)-Met to Ala) and the resulting differences in association energy (DeltaDeltaG(a)) are 2. 28, 5.08 and 2.93 kcal/mol for P(1)-Ala to Gly, P(1)-Met to Gly and to Ala with experimental values of 1.71, 4.62 and 2.91 kcal/mol, respectively. The calculated binding free energy differences are hence in excellent agreement with the experimental binding free energies. The binding free energies, however, were shown to be highly dependent on water molecules at the protein-protein interface and could only be quantitatively estimated if the correct number of such water molecules was included. Furthermore, the cavities that were formed when a large amino acid side-chain is perturbed to a smaller one seem to create instabilities in the systems and had to be refilled with water molecules in order to obtain reliable results. In addition, if the protein atoms that were perturbed away were not replaced by water molecules, the simulations dramatically overestimated the initial state of the free energy perturbations.

Structural comparison of psychrophilic and mesophilic trypsins. Elucidating the molecular basis of cold-adaptation

Structural rationalizations for differences in catalytic efficiency and stability between mesophilic and cold-adapted trypsins have been suggested from a detailed comparison of eight trypsin structures. Two trypsins, from Antarctic fish and Atlantic cod, have been constructed by homology modeling techniques and compared with six existing X-ray structures of both cold-adapted and mesophilic trypsins. The structural analysis focuses on the cold trypsin residue determinants found in a more extensive comparison of 27 trypsin sequences, and reveals a number of structural features unique to the cold-adapted trypsins. The increased substrate affinity of the psychrophilic trypsins is probably achieved by a lower electrostatic potential of the S1 binding pocket particularly arising from Glu221B, and from the lack of five hydrogen bonds adjacent to the catalytic triad. The reduced stability of the cold trypsins is expected to arise from reduced packing in two distinct core regions, fewer interdomain hydrogen bonds and from a destabilized C-terminal alpha-helix. The helices of the cold trypsins lack four hydrogen bonds and two salt-bridges, and they have poorer van der Waals packing interactions to the body of the molecule, compared to the mesophilic counterparts.

Comparative molecular dynamics of mesophilic and psychrophilic protein homologues studied by 1.2 ns simulations

It is well established that the dynamic motion of proteins plays an important functional role, and that the adaptation of a protein molecule to its environment requires optimization of internal non-covalent interactions and protein-solvent interactions. Serine proteinases in general, and trypsin in particular has been used as a model system in exploring possible structural features for cold adaptation. In this study, a 500 p.s. and a 1200 p.s. molecular dynamics (MD) simulation at 300 K of both anionic salmon trypsin and cationic bovine trypsin are analyzed in terms of molecular flexibility, internal non-covalent interactions and protein-solvent interactions. The present MD simulations do not indicate any increased flexibility of the cold adapted enzyme on an overall basis. However, the apparent higher flexibility and deformability of the active site of anionic salmon trypsin may lower the activation energy for ligand binding and for catalysis, and might be a reason for the increased binding affinity and catalytic efficiency compared to cationic bovine trypsin.

Residue determinants and sequence analysis of cold-adapted trypsins

The digestive enzyme trypsin is among the most extensively studied proteins, and its structure has been reported from a large number of organisms. This article focuses on the trypsins from vertebrates adapted to life at low temperatures. Cold-adapted organisms seem to have compensated for the reduced reaction rates at low temperatures by evolving more active and less temperature-stable enzymes. We have analyzed 27 trypsin sequences from a variety of organisms to find unique attributes for the cold-adapted trypsins, comparing trypsins from salmon, Antarctic fish, cod, and pufferfish to other vertebrate trypsins. Both the "cold" and the "warm" active trypsins have about 50 amino acids that are unique and conserved within each class. The main unique features of the cold-adapted trypsins attributable to low-temperature adaptation seem to be (1) reduced hydrophobicity and packing density of the core, mainly because of a lower (Ile + Leu)/(Ile + Leu + Val) ratio, (2) reduced stability of the C-terminal, (3) lack of one warm trypsin conserved proline residue and one proline tyrosine stacking, (4) difference in charge and flexibility of loops extending the binding pocket, and (5) different conformation of the "autolysis" loop that is likely to be involved in substrate binding.

The crystal structures of the complexes between bovine beta-trypsin and ten P1 variants of BPTI

The high-resolution X-ray structures have been determined for ten complexes formed between bovine beta-trypsin and P1 variants (Gly, Asp, Glu, Gln, Thr, Met, Lys, His, Phe, Trp) of bovine pancreatic trypsin inhibitor (BPTI). All the complexes were crystallised from the same conditions. The structures of the P1 variants Asp, Glu, Gln and Thr, are reported here for the first time in complex with any serine proteinase. The resolution of the structures ranged from 1.75 to 2.05 A and the R-factors were about 19-20 %. The association constants of the mutants ranged from 1.5x10(4) to 1.7x10(13) M-1. All the structures could be fitted into well-defined electron density, and all had very similar global conformations. All the P1 mutant side-chains could be accomodated at the primary binding site, but relative to the P1 Lys, there were small local changes within the P1-S1 interaction site. These comprised: (1) changes in the number and dynamics of water molecules inside the pocket; (2) multiple conformations and non-optimal dihedral angles for some of the P1 side-chains, Ser190 and Gln192; and (3) changes in temperature factors of the pocket walls as well as the introduced P1 side-chain. Binding of the cognate P1 Lys is characterised by almost optimal dihedral angles, hydrogen bonding distances and angles, in addition to considerably lower temperature factors. Thus, the trypsin S1 pocket seems to be designed particularly for lysine binding.

High-resolution structures of three new trypsin-squash-inhibitor complexes: a detailed comparison with other trypsins and their complexes

An anionic trypsin from Atlantic salmon and bovine trypsin have been complexed with the squash-seed inhibitors, CMTI-I (Cucurbita maxima trypsin inhibitor I, P1 Arg) and CPTI-II (Cucurbita pepo trypsin inhibitor II, P1 Lys). The crystal structures of three such complexes have been determined to 1.5-1.8 A resolution and refined to crystallographic R factors ranging from 17.6 to 19.3%. The two anionic salmon-trypsin complexes (ST-CPTI and ST-CMTI) and the bovine-trypsin complex (BT-CPTI) have been compared to other trypsin-inhibitor complexes by means of general structure and primary and secondary binding features. In all three new structures, the primary binding residue of the inhibitor binds to trypsin in the classical manner, but with small differences in the primary and secondary binding patterns. Lysine in CPTI-II binds deeper in the specificity pocket of bovine trypsin than lysine in other known lysine-bovine-trypsin complexes, and anionic salmon trypsin lacks some of the secondary binding interactions found in the complexes formed between squash inhibitors and bovine trypsin. The ST-CMTI complex was formed from the reactive-site-cleaved form of the inhibitor. However, well defined electron density was observed for the P1-P1' peptide bond, together with a hydrogen-bonding pattern virtually identical to those of all serine-protease-protein-inhibitor complexes, indicating a resynthesis of the scissile bond.

Structure of a non-psychrophilic trypsin from a cold-adapted fish species

The crystal structure of cationic trypsin (CST) from the Atlantic salmon (Salmo salar) has been refined at 1.70 A resolution. The crystals are orthorhombic, belong to space group P212121, with lattice parameters a = 65.91, b = 83.11 and c = 154.79 A, and comprise four molecules per asymmetric unit. The structure was solved by molecular replacement with AMoRe and refined with X-PLOR to an R value of 17.4% and Rfree of 21.5% for reflections |F| > 3sigmaF between 8.0 and 1.7 A resolution. The four non-crystallographic symmetry (NCS) related molecules in the asymmetric unit display r.m.s. deviations in the range 0.31-0.74 A for main-chain atoms, with the largest differences confined to two loops. One of these is the calcium-binding loop where the electron-density indicates a calcium ion for only one of the four molecules. In order to find structural rationalizations for the observed difference in thermostability and catalytic efficiency of CST, anionic salmon trypsin (AST) and bovine trypsin (BT), the three structures have been extensively compared. The largest deviations for the superimposed structures occur in the surface loops and particularly in the so-called 'autolysis loop'. Both the salmon enzymes possess a high methionine content, lower overall hydrophobicity and enhanced surface hydrophilicity, compared with BT. These properties have so far been correlated to cold-adaptation features, while in this work it is shown that the non-psychrophilic cationic salmon trypsin shares these features with the psychrophilic anionic salmon trypsin.

The crystal structure of anionic salmon trypsin in complex with bovine pancreatic trypsin inhibitor

The complex formed between anionic salmon trypsin (ST) and bovine pancreatic trypsin inhibitor (BPTI) has been crystallised, and the X-ray structure has been solved using the molecular replacement method. The crystals are hexagonal and belong to space group P6(1)22 with lattice parameters of a = b = 83.12 A and c = 222.15 A. Data have been collected to 2.1 A and the structure has been refined to a crystallographic R-factor of 20.6%. Catalysis by salmon trypsin is distinguished by a Km value 20-fold lower than that for mammalian trypsins, and a k(cat) twice as high. The present ST-BPTI complex serves as a model for the Michaelis-Menten complex, and has been compared with corresponding bovine and rat trypsin (RT) complexes. The binding of BPTI to salmon trypsin is characterised by stronger primary interactions in the active site, and a somewhat looser secondary binding.

Purification and characterization of pancreatic elastase from North Atlantic salmon (Salmo salar)

An elastase I-like enzyme was purified to homogeneity from the pyloric caeca of North Atlantic salmon (Salmo salar) and compared with porcine elastase I. The molecular weight and isoelectric point were estimated to be 27 kDa and over 9.3, respectively. The pH optimum was between 8.0 and 9.5, and the enzyme was unstable at pH values below 4. Kinetic properties examined using Suc-(Ala)3-p-nitroanilide showed that the catalytic efficiency of salmon elastase was about 2.5 times higher than that of porcine elastase. Furthermore, the salmon enzyme was less stable at lower pH values and temperatures than the porcine enzyme. The preference for amino acids at the primary binding site was found to be different from that of the porcine elastase. The salmon elastase binding pocket seems to prefer more branched aliphatic residues than the porcine elastase.

Structure of native pancreatic elastase from North Atlantic salmon at 1.61 Å resolution

The crystal structure of native salmon pancreatic elastase (SPE) has been solved by molecular-replacement methods, and refined by conventional conjugate-gradient methods and simulated-annealing techniques. The final R value is 17.2% for 21 389 reflections between 8.0 and 1.61 A, and the corresponding free R value is 23.9%. The overall tertiary structure of SPE is remarkably similar to that of porcine pancreatic elastase I (PPE), to which it shows about 67% sequence identity. The primary structure of SPE is determined from the electron-density maps, and only about 15 side chains are somewhat uncertain. Interesting differences between SPE and PPE, are one sequence deletion assigned to position 186, the residue 192 at the entrance of the specificity pocket is substituted from a Gln in PPE to Asn in SPE, and one of the calcium ligands is different. Furthermore, electron density is missing in SPE for the last three residues of the C-terminal helix. A comparison of the present amino-acid sequence of SPE with other sequences available indicates that SPE belongs to the class 1 pancreatic elastases.

Molecular cloning and characterization of anionic and cationic variants of trypsin from Atlantic salmon

Pancreatic cDNA libraries from Atlantic salmon (Salmo salar) were constructed and screened with salmon trypsin-specific probes. Five clones containing near full-length transcripts were selected for further characterization. Comparison of deduced amino acid sequences revealed that all variants possessed the canonical serine protease catalytic triad, consisting of histidine, aspartic acid and serine residues, a substrate-binding pocket with aspartic acid at the bottom, and 12 cysteine residues comprising six disulphide bridges. Translation in vitro of one of the trypsin clones produced a protein with the expected molecular mass of 24.5 kDa. Three of the Atlantic salmon trypsins (SalTRP-I, SalTRP-IA and SalTRP-IB) possessed very similar sequences and may represent allelic variants encoded by the same gene focus; however, existence as tetraploid loci or isoloci where disomic inheritance is incomplete may also exist in Atlantic salmon and cannot be excluded. Two other trypsin clones (SalTRP-II and SalTRP-III) are probably encoded by separate gene loci. Analysis of genomic DNA by Southern blotting and hybridization to a trypsin probe showed a complex pattern, indicative of a large number of gene loci for trypsin in Atlantic salmon. The charged amino acid distribution showed that four of the Atlantic salmon trypsin clones encoded anionic forms of the enzyme, while the fifth clone represented a cationic variant. Multiple alignments of the Atlantic salmon trypsin sequences with trypsin, chymotrypsin and elastase from different species placed all Atlantic salmon sequences approximately equidistant from trypsins of other species. Interestingly, the distance between the anionic and cationic variants from Atlantic salmon was similar to the distance between salmon and mammalian trypsins, revealing an early separation of these two types of trypsin, possibly prior to the derivation of fish during evolution. A structural model based on X-ray diffraction studies of the salmon trypsin protein was very similar to that of the mammalian enzyme. All residues which differ in charge between anionic and cationic trypsins were located at exposed regions of the proteins.

Crystallization and preliminary X-ray crystallographic studies of native elastase from North Atlantic salmon (Salmo salar)

Crystals of elastase from North Atlantic salmon have been grown from 2-methyl-2,4-pentanediol by the hanging-drop vapour-diffusion method at room temperature. They grow to dimensions of 0.7 x 0.4 x 0.3 mm in three weeks. The crystals belong to the tetragonal space group P4(1)2(1)2 or P4(3)2(1)2 with cell dimensions a = b = 68.0 A and c = 84.0 A. There are eight molecules in the unit cell. The crystals diffract to at least 1.6 A resolution and are suitable for a high-resolution crystal structure determination.

Refined crystal structure of lysozyme from the rainbow trout (Oncorhynchus mykiss)

Lysozymes (E.C. 3.2.1.17) are well characterized ubiquitous enzymes that have an antibacterial effect. The lysozymes from rainbow trout (RBTL) (Oncorhynchus mykiss) could be particularly interesting in aquaculture since they show higher activity than egg-white lysozyme and lysozymes from other fish species against a variety of pathogenic bacteria. Two lysozymes, I and II, differing only in a single amino acid, were purified from the kidney of rainbow trout and shown to belong to the c-type class of lysozymes. The type II form was shown to be much more potent against a variety of bacteria than the type I enzyme. We have grown crystals from a mixture containing about 80% type I and 20% type II lysozyme from rainbow trout, and solved the X-ray crystal structure. The crystals are trigonal with a = 76.68, c = 54.46 A and space group P3(1)21. The phase problem was solved by the molecular-replacement method, and the structure was refined to an R-factor of 17.4% using data to 1.8 A resolution. The crystal structure shows that the three-dimensional structure of rainbow trout lysozyme is very similar to the previously solved structures of other c-type lysozymes. The single polypeptide of 129 amino acids is folded into two domains separated by a deep cleft which contains the active site. Secondary-structure elements, four alpha-helices and a three-stranded beta-sheet, are located in the same sequential positions as in the hen, turkey and human enzymes. The beta-sheet is found to be common for structures of both c- and g-type lysozymes. We suggest that differences in antibiotic activity of the two forms of RBTL are probably due to small differences in the hydophobicity of a small surface region.

Comparative molecular dynamics simulation studies of salmon and bovine trypsins in aqueous solution

The flexibility and conformational behaviour of salmon and bovine trypsins were modelled with a 300 ps molecular dynamics simulation in aqueous solution. Trajectories from both trypsins were analysed to eventually detect differences in mobility that could explain observed variations in stability and activity. The simulations were performed at 300 K with all the acidic groups deprotonated and the basic groups protonated. The radius of gyration, the overall r.m.s. deviation from the starting structure as a function of time, together with the r.m.s. deviation from the starting structures as a function of residue number, demonstrated that the simulations were stable and representative of the X-ray structures of both enzymes. Isotropic Debye-Waller factors were calculated from the fluctuations for main-chain atoms and were in good agreement with experimental values. The overall dynamic properties of the two enzymes were similar. Based on the present 300 ps molecular dynamics simulation, it cannot be concluded that either of the two enzymes is more 'flexible' than the other. However, there are clearly differences in mobility on a more detailed level and for particular regions.