Evolution of prokaryotic subtilases: genome-wide analysis reveals novel subfamilies with different catalytic residues

Subtilisin-like serine proteases (subtilases) are a very diverse family of serine proteases with low sequence homology, often limited to regions surrounding the three catalytic residues. Starting with different Hidden Markov Models (HMM), based on sequence alignments around the catalytic residues of the S8 family (subtilisins) and S53 family (sedolisins), we iteratively searched all ORFs in the complete genomes of 313 eubacteria and archaea. In 164 genomes we identified a total of 567 ORFs with one or more of the conserved regions with a catalytic residue. The large majority of these contained all three regions around the "classical" catalytic residues of the S8 family (Asp-His-Ser), while 63 proteins were identified as S53 (sedolisin) family members (Glu-Asp-Ser). More than 30 proteins were found to belong to two novel subsets with other evolutionary variations in catalytic residues, and new HMMs were generated to search for them. In one subset the catalytic Asp is replaced by an equivalent Glu (i.e. Glu-His-Ser family). The other subset resembles sedolisins, but the conserved catalytic Asp is not located on the same helix as the nucleophile Glu, but rather on a beta-sheet strand in a topologically similar position, as suggested by homology modeling. The Prokaryotic Subtilase Database (www.cmbi.ru.nl/subtilases) provides access to all information on the identified subtilases, the conserved sequence regions, the proposed family subdivision, and the appropriate HMMs to search for them. Over 100 proteins were predicted to be subtilases for the first time by our improved searching methods, thereby improving genome annotation.     

Cloning of usp45, a gene encoding a secreted protein from Lactococcus lactis subsp. lactis MG1363

We have cloned usp45, a gene encoding an extracellular secretory protein of Lactococcus lactis subsp. lactis strain MG1363. Unidentified secreted 45-kDa protein (Usp45) is secreted by every mesophilic L. lactis strain we tested so far and it is chromosomally encoded. The nucleotide sequence of the usp45 gene revealed an open reading frame of 1383 bp encoding a protein of 461 amino acids (aa), composed of a 27-aa signal peptide and a mature protein initiated at Asp28. The gene contains a consensus promoter sequence and a weak ribosome-binding site; the latter is rather uncommon for Gram-positive bacteria. Expression studies in Escherichia coli showed efficient synthesis and secretion of the protein. Usp45 has an unusual aa composition and distribution, and it is predicted to be structurally homologous with P54 of Enterococcus faecium. Up to now, no biological activity could be postulated for this secreted protein.     

Analysis of growth of Lactobacillus plantarum WCFS1 on a complex medium using a genome-scale metabolic model

A genome-scale metabolic model of the lactic acid bacterium Lactobacillus plantarum WCFS1 was constructed based on genomic content and experimental data. The complete model includes 721 genes, 643 reactions, and 531 metabolites. Different stoichiometric modeling techniques were used for interpretation of complex fermentation data, as L. plantarum is adapted to nutrient-rich environments and only grows in media supplemented with vitamins and amino acids. (i) Based on experimental input and output fluxes, maximal ATP production was estimated and related to growth rate. (ii) Optimization of ATP production further identified amino acid catabolic pathways that were not previously associated with free-energy metabolism. (iii) Genome-scale elementary flux mode analysis identified 28 potential futile cycles. (iv) Flux variability analysis supplemented the elementary mode analysis in identifying parallel pathways, e.g. pathways with identical end products but different co-factor usage. Strongly increased flexibility in the metabolic network was observed when strict coupling between catabolic ATP production and anabolic consumption was relaxed. These results illustrate how a genome-scale metabolic model and associated constraint-based modeling techniques can be used to analyze the physiology of growth on a complex medium rather than a minimal salts medium. However, optimization of biomass formation using the Flux Balance Analysis approach, reported to successfully predict growth rate and by product formation in Escherichia coli and Saccharomyces cerevisiae, predicted too high biomass yields that were incompatible with the observed lactate production. The reason is that this approach assumes optimal efficiency of substrate to biomass conversion, and can therefore not predict the metabolically inefficient lactate formation.     

Accelerating the reconstruction of genome-scale metabolic networks

Background  

The genomic information of a species allows for the genome-scale reconstruction of its metabolic capacity. Such a metabolic reconstruction gives support to metabolic engineering, but also to integrative bioinformatics and visualization. Sequence-based automatic reconstructions require extensive manual curation, which can be very time-consuming. Therefore, we present a method to accelerate the time-consuming process of network reconstruction for a query species. The method exploits the availability of well-curated metabolic networks and uses high-resolution predictions of gene equivalency between species, allowing the transfer of gene-reaction associations from curated networks.     

Epochal evolution of GGII.4 norovirus capsid proteins from 1995 to 2006

Noroviruses are the causative agents of the majority of viral gastroenteritis outbreaks in humans. During the past 15 years, noroviruses of genotype GGII.4 have caused four epidemic seasons of viral gastroenteritis, during which four novel variants (termed epidemic variants) emerged and displaced the resident viruses. In order to understand the mechanisms and biological advantages of these epidemic variants, we studied the genetic changes in the capsid proteins of GGII.4 strains over this period. A representative sample was drawn from 574 GGII.4 outbreak strains collected over 15 years of systematic surveillance in The Netherlands, and capsid genes were sequenced for a total of 26 strains. The three-dimensional structure was predicted by homology modeling, using the Norwalk virus (Hu/NoV/GGI.1/Norwalk/1968/US) capsid as a reference. The highly significant preferential accumulation and fixation of mutations (nucleotide and amino acid) in the protruding part of the capsid protein provided strong evidence for the occurrence of genetic drift and selection. Although subsequent new epidemic variants differed by up to 25 amino acid mutations, consistent changes were observed in only five positions. Phylogenetic analyses showed that each variant descended from its chronologic predecessor, with the exception of the 2006b variant, which is more closely related to the 2002 variant than to the 2004 variant. The consistent association between the observed genetic findings and changes in epidemiology leads to the conclusion that population immunity plays a role in the epochal evolution of GGII.4 norovirus strains.     

Deletion of various carboxy-terminal domains of Lactococcus lactis SK11 proteinase: effects on activity, specificity, and stability of the truncated enzyme

The Lactococcus lactis SK11 cell envelope proteinase is an extracellular, multidomain protein of nearly 2,000 residues consisting of an N-terminal serine protease domain, followed by various other domains of largely unknown function. Using a strategy of deletion mutagenesis, we have analyzed the function of several C-terminal domains of the SK11 proteinase which are absent in cell envelope proteinases of other lactic acid bacteria. The various deletion mutants were functionally expressed in L. lactis and analyzed for enzyme stability, activity, (auto)processing, and specificity toward several substrates. C-terminal deletions of first the cell envelope W (wall) and AN (anchor) domains and then the H (helix) domain leads to fully active, secreted proteinases of unaltered specificity. Gradually increasing the C-terminal deletion into the so-called B domain leads to increasing instability and autoproteolysis and progressively less proteolytic activity. However, the mutant with the largest deletion (838 residues) from the C terminus and lacking the entire B domain still retains proteolytic activity. All truncated enzymes show unaltered proteolytic specificity toward various substrates. This suggests that the main role played by these domains is providing stability or protection from autoproteolysis (B domain), spacing away from the cell (H domain), and anchoring to the cell envelope (W and AN domains). In addition, this study allowed us to more precisely map the main C-terminal autoprocessing site of the SK11 proteinase and the epitope for binding of group IV monoclonal antibodies.