Comparative analysis of HiSeq3000 and BGISEQ-500 sequencing platform over whole genome sequencing metagenomics data

Recent advances in sequencing technologies and platforms have enabled to generate metagenomics sequences using different sequencing platforms. In this study, we analyzed and compared shotgun metagenomic sequences generated by HiSeq3000 and BGISEQ-500 platforms from 12 sediment samples collected across the Norwegian coast. Metagenomics DNA sequences were normalized to an equal number of bases for both platforms and further evaluated by using different taxonomic classifiers, reference databases, and assemblers. Normalized BGISEQ-500 sequences retained more reads and base counts after preprocessing, while a slightly higher fraction of HiSeq3000 sequences were taxonomically classified. Kaiju classified a higher percentage of reads relative to Kraken2 for both platforms, and comparison of reference database for taxonomic classification showed that MAR database outperformed RefSeq. Assembly using MEGAHIT produced longer assemblies and higher total contigs count in majority of HiSeq3000 samples than using metaSPAdes, but the assembly statistics notably improved with unprocessed or normalized reads. Our results indicate that both platforms perform comparably in terms of the percentage of taxonomically classified reads and assembled contig statistics for metagenomics samples. This study provides valuable insights for researchers in selecting an appropriate sequencing platform and bioinformatics pipeline for their metagenomics studies.

Phylogenetic Revision of the Genus Aliivibrio: Intra- and Inter-Species Variance Among Clusters Suggest a Wider Diversity of Species

Genus Aliivibrio is known to harbor species exhibiting bioluminescence as well as pathogenic behavior affecting the fish farming industry. Current phylogenetic understanding of Aliivibrio has largely remained dormant after reclassification disentangled it from the Vibrio genus in 2007. There is growing evidence of wider diversity, but until now the lack of genomes and selective use of type strains have limited the ability to compare and classify strains firmly. In this study, a total of 143 bacterial strains, including 51 novel sequenced strains, were used to strengthen phylogenetic relationships in Aliivibrio by exploring intra-species and inter-species relations. Multilocus sequence analysis (MLSA), applying the six housekeeping genes 16S ribosomal RNA (rRNA), gapA, gyrB, pyrH, recA, and rpoA, inferred 12 clades and a singular branch in Aliivibrio. Along with four new phylogenetic clades, the MLSA resolved prior inconsistencies circumscribing Aliivibrio wodanis and formed a unique clade we propose as the novel species Aliivibrio sp. “friggae.” Furthermore, phylogenetic assessment of individual marker genes showed gyrB, pyrH, and recA superior to the 16S rRNA gene, resolving accurately for most species clades in Aliivibrio. In this study, we provide a robust phylogenetic groundwork for Aliivibrio as a reference point to classification of species.

The MAR databases: development and implementation of databases specific for marine metagenomics

We introduce the marine databases; MarRef, MarDB and MarCat (https://mmp.sfb.uit.no/databases/), which are publicly available resources that promote marine research and innovation. These data resources, which have been implemented in the Marine Metagenomics Portal (MMP) (https://mmp.sfb.uit.no/), are collections of richly annotated and manually curated contextual (metadata) and sequence databases representing three tiers of accuracy. While MarRef is a database for completely sequenced marine prokaryotic genomes, which represent a marine prokaryote reference genome database, MarDB includes all incomplete sequenced prokaryotic genomes regardless level of completeness. The last database, MarCat, represents a gene (protein) catalog of uncultivable (and cultivable) marine genes and proteins derived from marine metagenomics samples. The first versions of MarRef and MarDB contain 612 and 3726 records, respectively. Each record is built up of 106 metadata fields including attributes for sampling, sequencing, assembly and annotation in addition to the organism and taxonomic information. Currently, MarCat contains 1227 records with 55 metadata fields. Ontologies and controlled vocabularies are used in the contextual databases to enhance consistency. The user-friendly web interface lets the visitors browse, filter and search in the contextual databases and perform BLAST searches against the corresponding sequence databases. All contextual and sequence databases are freely accessible and downloadable from https://s1.sfb.uit.no/public/mar/.

Effects of Probiotic Supplementation on the Gut Microbiota and Antibiotic Resistome Development in Preterm Infants

Objectives: In 2014 probiotic supplementation (Lactobacillus acidophilus and Bifidobacterium longum subspecies infantis; Infloran^Ⓡ) was introduced as standard of care to prevent necrotizing enterocolitis (NEC) in extremely preterm infants in Norway. We aimed to evaluate the influence of probiotics and antibiotic therapy on the developing gut microbiota and antibiotic resistome in extremely preterm infants, and to compare with very preterm infants and term infants not given probiotics. Study design: A prospective, observational multicenter study in six tertiary-care neonatal units. We enrolled 76 infants; 31 probiotic-supplemented extremely preterm infants <28 weeks gestation, 35 very preterm infants 28-31 weeks gestation not given probiotics and 10 healthy full-term control infants. Taxonomic composition and collection of antibiotic resistance genes (resistome) in fecal samples, collected at 7 and 28 days and 4 months age, were analyzed using shotgun-metagenome sequencing. Results: Median (IQR) birth weight was 835 (680-945) g and 1,290 (1,150-1,445) g in preterm infants exposed and not exposed to probiotics, respectively. Two extremely preterm infants receiving probiotic developed NEC requiring surgery. At 7 days of age we found higher median relative abundance of Bifidobacterium in probiotic supplemented infants (64.7%) compared to non-supplemented preterm infants (0.0%) and term control infants (43.9%). Lactobacillus was only detected in small amounts in all groups, but the relative abundance increased up to 4 months. Extremely preterm infants receiving probiotics had also much higher antibiotic exposure, still overall microbial diversity and resistome was not different than in more mature infants at 4 weeks and 4 months. Conclusion: Probiotic supplementation may induce colonization resistance and alleviate harmful effects of antibiotics on the gut microbiota and antibiotic resistome. Clinical Trial Registration: Clinicaltrials.gov: NCT02197468. https://clinicaltrials.gov/ct2/show/NCT02197468.

Presence of acyl-homoserine lactones in 57 members of the Vibrionaceae family

AIMS

The aim of this study was to use a sensitive method to screen and quantify 57 Vibrionaceae strains for the production of acyl-homoserine lactones (AHLs) and map the resulting AHL profiles onto a host phylogeny.

METHODS AND RESULTS

We used a high-performance liquid chromatography-tandem mass spectrometry (HPLC-MS/MS) protocol to measure AHLs in spent media after bacterial growth. First, the presence/absence of AHLs (qualitative analysis) was measured to choose internal standard for subsequent quantitative AHL measurements. We screened 57 strains from three genera (Aliivibrio, Photobacterium and Vibrio) of the same family (i.e. Vibrionaceae). Our results show that about half of the isolates produced multiple AHLs, typically at 25-5000 nmol l(-1) .

CONCLUSIONS

This work shows that production of AHL quorum sensing signals is found widespread among Vibrionaceae bacteria and that closely related strains typically produce similar AHL profiles.

SIGNIFICANCE AND IMPACT OF THE STUDY

The AHL detection protocol presented in this study can be applied to a broad range of bacterial samples and may contribute to a wider mapping of AHL production in bacteria, for example, in clinically relevant strains.

Identification of type VI secretion systems in Moritella viscosa

The study describes the identification of type VI secretion systems (T6SSs) in Moritella viscosa, the aetiological agent of winter ulcer disease. Despite the availability of commercial vaccines, M. viscosa causes significant financial losses in salmonid farming. The T6SS transports bacterial proteins from the cell into the environment or directly into host cells, and has been implicated with bacterial virulence. The aim of the study was to identify potential T6SSs in M. viscosa and to determine whether it possesses active T6S, providing further insight into the biology of the bacterium. The genome of M. viscosa 06/09/139 was screened for homology with known T6SS encoding genes. Two genetically distinct loci, termed Moritella Type Six Secretion 1 and 2 (mts1 and mts2), were identified as encoding putative T6SSs. Each locus contained known T6S core genes. The mts2 locus contained species specific genes, some of which have not previously been connected with T6S. The mts1 locus showed sequence homology and synteny to T6SSs of the fish pathogen Aliivibrio salmonicida and a non-pathogenic Moritella sp. PE36. The mts2 locus was more similar to a Vibrio parahaemolyticus T6SS. A functional T6SS was confirmed through identification of secreted Mts1-M, a hemolysin coregulated protein (Hcp) which is a part of the secretion system. Both virulent and avirulent M. viscosa isolates expressed two genes encoding Hcp, mts1-M and mts2-M. The results show that M. viscosa has a functional T6S, but the role of the secretion system and possible connections with virulence need further examination.

DeltaProt: a software toolbox for comparative genomics

Background

Statistical bioinformatics is the study of biological data sets obtained by new micro-technologies by means of proper statistical methods. For a better understanding of environmental adaptations of proteins, orthologous sequences from different habitats may be explored and compared. The main goal of the DeltaProt Toolbox is to provide users with important functionality that is needed for comparative screening and studies of extremophile proteins and protein classes. Visualization of the data sets is also the focus of this article, since visualizations can play a key role in making the various relationships transparent. This application paper is intended to inform the reader of the existence, functionality, and applicability of the toolbox.

Results

We present the DeltaProt Toolbox, a software toolbox that may be useful in importing, analyzing and visualizing data from multiple alignments of proteins. The toolbox has been written in MATLAB™ to provide an easy and user-friendly platform, including a graphical user interface, while ensuring good numerical performance. Problems in genome biology may be easily stated thanks to a compact input format. The toolbox also offers the possibility of utilizing structural information from the SABLE or other structure predictors. Different sequence plots can then be viewed and compared in order to find their similarities and differences. Detailed statistics are also calculated during the procedure.

Conclusions

The DeltaProt package is open source and freely available for academic, non-commercial use. The latest version of DeltaProt can be obtained from http://services.cbu.uib.no/software/deltaprot/. The website also contains documentation, and the toolbox comes with real data sets that are intended for training in applying the models to carry out bioinformatical and statistical analyses of protein sequences.

Equipped with the new algorithms proposed here, DeltaProt serves as an auxiliary analysis tool capable of visualizing and comparing orthologus protein sequences. The framework of the algorithms also enables easy incorporation of extra information on structure, if such data is available.

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.

Molecular characterization of cold adaptation based on ortholog protein sequences from Vibrionaceae species

A set of 298 protein families from psychrophilic Vibrio salmonicida was compiled to identify genotypic characteristics that discern it from orthologous sequences from the mesophilic Vibrio/Photobacterium branch of the gamma-Proteobacteria (Vibrionaceae family). In our comparative exploration we employed alignment based bioinformatical and statistical methods. Interesting information was found in the substitution matrices, and the pattern of asymmetries in the amino acid substitution process. Together with the compositional difference, they identified the amino acids Ile, Asn, Ala and Gln as those having the most psycrophilic involvement. Ile and Asn are enhanced whereas Gln and Ala are suppressed. The inflexible Pro residue is also suppressed in loop regions, as expected in a flexible structure. The dataset were also classified and analysed according to the predicted subcellular location, and we made an additional study of 183 intracellular and 65 membrane proteins. Our results revealed that the psychrophilic proteins have similar hydrophobic and charge contributions in the core of the protein as mesophilic proteins, while the solvent-exposed surface area is significantly more hydrophobic. In addition, the psychrophilic intracellular (but not the membrane) proteins are significantly more negatively charged at the surface. Our analysis supports the hypothesis of preference for more flexible amino acids at the molecular surface. Life in cold climate seems to be obtained through many minor structural modifications rather than certain amino acids substitutions.

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 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.

Comparative expression study to increase the solubility of cold adapted Vibrio proteins in Escherichia coli

Functional and structural studies require gene overexpression and purification of soluble proteins. We wanted to express proteins from the psychrophilic bacterium Vibrio salmonicida in Escherichia coli, but encountered solubility problems. To improve the solubility of the proteins, we compared the effects of six N-terminal fusion proteins (Gb1, Z, thioredoxin, GST, MBP and NusA) and an N-terminal His6-tag. The selected test set included five proteins from the fish pathogen V. salmonicida and two related products from the mesophilic human pathogen Vibrio cholerae. We tested the expression in two different expression strains and at three different temperatures (16, 23 and 37 degrees C). His6-tag was the least effective tag, and these vector constructs were also difficult to transform. MBP and NusA performed best, expressing soluble proteins with all fusion partners in at least one of the cell types. In some cases MBP, GST and thioredoxin fusions resulted in products of incorrect size. The effect of temperature is complex: in most cases level of expression increased with temperature, whereas the effect on solubility was opposite. We found no clear connection between the preferred expression temperature of the protein and the temperature of the original host organism's natural habitat.

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.

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.

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).

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.

Functional characterization of isoschizomeric His-Cys box homing endonucleases from Naegleria

Several species within the amoeboflagellate genus Naegleria harbor an optional ORF containing group I introns in their nuclear small subunit ribosomal DNA. The different ORFs encode homing endonucleases with 65 to 95% identity at the amino-acid level. I-NjaI, I-NanI and I-NitI, from introns in Naegleria jamiesoni, N. andersoni and N. italica, respectively, were analyzed in more detail and found to be isoschizomeric endonucleases that recognize and cleave an approximal 19-bp partially symmetrical sequence, creating a pentanucleotide 3' overhang upon cleavage. The optimal conditions for cleavage activity with respect to temperature, pH, salt and divalent metal ions were investigated. The optimal cleavage temperature for all three endonucleases was found to be 37 degrees C and the activity was dependent on the concentration of NaCl with an optimum at 200 mM. Divalent metal ions, primarily Mg2+, are essential for Naegleria endonuclease activity. Whereas both Mn2+ and Ca2+ could substitute for Mg2+, but with a slower cleavage rate, Zn2+ was unable to support cleavage. Interestingly, the pH dependence of DNA cleavage was found to vary significantly between the I-NitI and I-NjaI/I-NanI endonucleases with optimal pH values at 6.5 and 9, respectively. Site-directed mutagenesis of conserved I-NjaI residues strongly supports the hypothesis that Naegleria homing endonucleases share a similar zinc-binding structure and active site with the His-Cys box homing endonuclease I-PpoI.

Purification and characterization of a cold-adapted uracil-DNA glycosylase from Atlantic cod (Gadus morhua)

Uracil-DNA glycosylase (UDG; UNG) has been purified 17000-fold from Atlantic cod liver (Gadus morhua). The enzyme has an apparent molecular mass of 25 kDa, as determined by gel filtration, and an isoelectric point above 9.0. Atlantic cUNG is inhibited by the specific UNG inhibitor (Ugi) from the Bacillus subtilis bacteriophage (PBS2), and has a 2-fold higher activity for single-stranded DNA than for double-stranded DNA. cUNG has an optimum activity between pH 7.0-9.0 and 25-50 mM NaCl, and a temperature optimum of 41 degrees C. Cod UNG was compared with the recombinant human UNG (rhUNG), and was found to have slightly higher relative activity at low temperatures compared with their respective optimum temperatures. Cod UNG is also more pH- and temperature labile than rhUNG. At pH 10.0, the recombinant human UNG had 66% residual activity compared with only 0.4% for the Atlantic cUNG. At 50 degrees C, cUNG had a half-life of 0.5 min compared with 8 min for the rhUNG. These activity and stability experiments reveal cold-adapted features in cUNG.

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.

Arctic charr (Salvelinus alpinus) vitellogenin: development and validation of an enzyme-linked immunosorbent assay

Vitellogenin (Vtg) was isolated from plasma of estradiol-17 beta-treated Arctic charr males by double precipitation with MgCl2-EDTA and distilled water, followed by ion-exchange chromatography. The monomeric form of Vtg, as revealed by sodium dodecyl sulfate-polyacrylamide gel electrophoresis, was 158 kDa. The purified Vtg was used to raise a polyclonal antibody for Vtg (AbVtg), and the specificity of the AbVtg was assessed by Western blot analysis. No cross-reactivity was observed with plasma from control males. Using this AbVtg, a competitive enzyme-linked immunosorbent assay was developed. The detection limit of the assay was 2 ng ml-1, and the intra- and inter-assay variations determined from plasma samples were 8.6 and 13.3%, respectively. The assay was validated by quantification of Vtg in plasma samples obtained during a reproductive cycle of Arctic charr. Vtg of females increased gradually from 3 mg ml-1 in early March to a peak value of 22 mg ml-1 in late August, followed by a rapid drop to 2 mg ml-1 at the time of spawning in mid-October. The temporal changes in plasma Vtg of females correlate well with the reproductive cycle. Vtg was undetectable in males, except on some sampling dates during July-September when minute amounts (3-13 micrograms ml-1) were detected in some individuals.

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.

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.

I-NjaI, a nuclear intron-encoded homing endonuclease from Naegleria, generates a pentanucleotide 3' cleavage-overhang within a 19 base-pair partially symmetric DNA recognition site

Different species of the amoebo-flagellate Naegleria harbor optional group I introns in the nuclear ribosomal DNA that contain open reading frames. Intron proteins from Naegleria jamiesoni, Naegleria andersoni, and Naegleria italica (named I-NjaI, I-NanI and I-NitI, respectively) were expressed in Escherichia coli and found to be isoschizomeric homing endonucleases that specifically recognize and cleave intron-lacking homologous alleles of ribosomal DNA. The I-NjaI endonuclease was affinity purified, characterized in more detail, and found to generate five-nucleotide 3' staggered ends at the intron insertion site which differs from the ends generated by all other known homing endonucleases. The recognition site was delimited and found to cover an approximately 19 base-pair partially symmetric sequence spanning both the cleavage site and the intron insertion site. The palindromic feature was supported by mutational analysis of the target DNA. All single-site substitutions within the recognition sequence were cleaved by the purified I-NjaI endonuclease, but at different efficiencies. The center of symmetry and cleavage was found to be completely degenerate in specificity, which resembles that of the subclass IIW bacterial restriction enzymes.

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.

Temperature and pH sensitivity of trypsins from Atlantic salmon (Salmo salar) in comparison with bovine and porcine trypsin

Four differently charged trypsins were purified from pyloric caeca of Atlantic salmon (Salmo salar). The isoelectric points of three anionic isoforms were 4.70, 4.60, and 4.55 (anionic trypsin I, II and III, respectively). And for the first time a cationic isoform (isoelectric point above 9.3) has been isolated from a marine species. The apparent molecular weights of all four isoforms were about 25 kDa as determined by SDS-PAGE. The salmon enzymes were inhibited by serine proteinase inhibitors in general and also by specific trypsin inhibitors. Anionic trypsin I and the cationic isoform were further examined. Anionic trypsin I showed the typical cold-adaptation features, low pH and temperature stability (also lower Gibb's free energy of GdnHCl-induced unfolding) and high catalytic efficiency as compared to the mammalian trypsins. The cationic isoform did not show these features, but resembled the mammalian trypsins.

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.

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.

Cold adaption of enzymes: structural comparison between salmon and bovine trypsins

The crystal structure of an anionic form of salmon trypsin has been determined at 1.82 A resolution. We report the first structure of a trypsin from a phoikilothermic organism in a detailed comparison to mammalian trypsins in order to look for structural rationalizations for the cold-adaption features of salmon trypsin. This form of salmon trypsin (ST II) comprises 222 residues, and is homologous to bovine trypsin (BT) in about 65% of the primary structure. The tertiary structures are similar, with an overall displacement in main chain atomic positions between salmon trypsin and various crystal structures of bovine trypsin of about 0.8 A. Intramolecular hydrogen bonds and hydrophobic interactions are compared and discussed in order to estimate possible differences in molecular flexibility which might explain the higher catalytic efficiency and lower thermostability of salmon trypsin compared to bovine trypsin. No overall differences in intramolecular interactions are detected between the two structures, but there are differences in certain regions of the structures which may explain some of the observed differences in physical properties. The distribution of charged residues is different in the two trypsins, and the impact this might have on substrate affinity has been discussed.

Intragenic rearrangements in the mitochondrial NADH dehydrogenase subunit 6 gene of vertebrates

We have sequenced the mitochondrial-encoded NADH dehydrogenase subunit 6 gene from 19 species of birds. Comparison of the derived amino-acid sequences in 22 avian species, six mammals, and two fishes, reveals an intragenic rearrangement in mammals. The C-terminal half of the mammalian protein includes an internal insertion of 10-15 amino acids and a C-terminal deletion of 8-9 amino acids. Based on comparative sequence alignments and hydropathy profile analysis, five hydrophobic segments (designated I to V) corresponding to transmembrane regions are proposed. In this structural model of NADH dehydrogenase subunit 6, the mammalian insertion is found in a variable loop region between transmembrane segments IV and V.

Fluorescence studies of beta-adrenergic ligand binding to alpha 1-acid glycoprotein with 1-anilino-8-naphthalene sulfonate, isoprenaline, adrenaline and propranolol

The present study shows that ANS (1-anilino-8-naphthalene sulfonate), propranolol, isoprenaline, adrenaline and dopamine have common binding sites on AAG (alpha 1-acid glycoprotein). A fluorescence technique was employed to characterize the interaction between the ligands and AAG at 20-22 degrees. The binding of ANS to AAG caused increased fluorescence intensity at emission and excitation wavelengths of 400 and 470 nm. In this situation, propranolol displaced ANS in a concentration-dependent mode with an apparent dissociation constant of 6.2 +/- 0.01 microM, whereas isoprenaline did not reduce the ANS-AAG fluorescence. However, in the presence of AAG, catecholamines caused a marked increase of fluorescence at excitation and emission wavelengths of 250 and 325 nm, respectively. These wavelengths were employed to characterize the binding of isoprenaline, adrenaline and propranolol to AAG. Two subsets of binding sites were demonstrated. The Kd values were 0.87 +/- 0.03 and 25.1 +/- 10.7 microM for ANS, 0.76 +/- 0.09 and 133 +/- 30.4 microM for propranolol, 140 +/- 14 and 2.18 +/- 0.58 mM for isoprenaline, 137 +/- 24 and 14.8 +/- 0.1 mM for adrenaline, respectively. AAG had identical high affinity binding capacity for these ligands (n approximately 1). However, the second class of binding sites showed ligand-dependent binding capacity: n = 1 for ANS, n approximately 10 for propranolol, n approximately 15 for adrenaline, n approximately 20 for isoprenaline, respectively. ANS, propranolol, dopamine and adrenaline caused concentration-dependent inhibition of isoprenaline binding to AAG with apparent dissociation constants of 5.1 +/- 1.8 microM, 6.4 +/- 1.1 microM, 0.57 +/- 0.13 mM and 1.5 +/- 0.46 mM, respectively.

Inhibition of 7-hydroxymethotrexate formation by amsacrine

The inhibition of methotrexate (MTX) biotransformation to 7-hydroxymethotrexate (7-OH-MTX) by 4'-(9-acridinylamino)-methanesulfon-m-anisidide (mAMSA) was studied in bile-drained rats in vivo and in incubates of isolated rat hepatocytes and rat-liver homogenate in vitro. In vivo, i.v. administration of 10 mg/kg mAMSA prior to [3H]-MTX infusion (50 mg/kg) led to a significant alteration in 7-OH-MTX kinetics. 7-OH-MTX peak concentrations and AUC in bile and serum were reduced by 75% and the recovery of MTX as 7-OH-MTX in bile and urine decreased by 70%, whereas MTX pharmacokinetics remained unaltered. In suspensions of isolated hepatocytes, 10 microM mAMSA led to a 54% decrease in 7-OH-MTX formation. However, the hepatocellular influx and efflux of MTX was not perturbed by mAMSA. Preincubation of rat-liver homogenates with 1.25-10 microM mAMSA reduced the formation of 7-OH-MTX by up to 73%. mAMSA appeared to inhibit MTX hydroxylation competitively, exhibiting a Ki of 3 microM. Due to its inhibition of the MTX-oxidizing system, mAMSA may be beneficial in combination chemotherapy with MTX by reducing 7-OH-MTX-associated toxicity and, possibly, enhancing the cytotoxic effects of MTX.

Effect of 2,3-butanedione on human myeloperoxidase

1. Myeloperoxidase is inhibited by various diketones that are recognized arginine reagents. 2. Although arginine residues in the enzyme were modified in both the light and the dark, enzyme inactivation occurred only in the presence of light. 3. Under conditions where diketones caused inactivation of myeloperoxidase, spectral studies indicated marked damage to the haem residues of the enzyme. 4. It was concluded that diketones serve simply as photosensitizers of visible light-induced inactivation of myeloperoxidase. 5. Studies on other haemoproteins indicated the great ease with which the presence of diketones sensitized haem residues for photodestruction.