Amino Acid Preferences Research Articles

The Wheel of Fortune: Helical Wheel Alanine Scanning of a Spider Venom Antimicrobial Peptide Reveals Residues Involved in Antimicrobial and Cytotoxic Activity.

A preference for several amino acids is observed to occur at particular positions of cationic α-helical antimicrobial peptides (AMPs), which ensures the formation of amphipathic regions once they assume their correct secondary structure in membranes or membrane-mimicking environments and makes them active against pathogens. This study determined the effect of alanine mutations on the secondary structure and bioactivity of lyp1987 (GRLQAFLAKMKEIAAQTL-NH2), a cationic α-helical AMP obtained from the venom of Lycosa poonaensis which exhibits broad range activity against Gram-positive and Gram-negative bacteria with micromolar minimum inhibitory concentrations (MIC). CD spectroscopy revealed no significant difference in the secondary structure, with all alanine-substituted analogs exhibiting predominantly α-helical structure in buffered 2,2,2-trifluoroethanol solution. Alanine substitution at Glu12 and Thr17 increased the activity of lyp1987 against Gram-positive and -negative bacteria, while alanine substitution at Lys9 increased its selectivity against Gram-positive bacteria. Further investigation can be done to determine positions and substitutions that will give less cytotoxic analogs.

Bacterial functional responses to environmental variability: a case study in three distinct mountain lakes within a single watershed

The utilization of carbon substrates (CSs) is crucial for the functioning of heterotrophic bacterial communities and aquatic ecosystems. This study used the Biolog EcoplateTM technique to explore spatial and temporal community-level physiological profiles (CLPPs) of bacteria in three eutrophic-dystrophic mountain lakes in Bulgaria. Bacterial metabolic activity varied among the lakes, being highest in the lake dominated by phytoplankton and lowest in the lake dominated by macrophyte primary producers. CLPP highlighted specific bacterial metabolic preferences to certain carbon guilds (CGs), showing a higher affinity for carbohydrates and polymers and lower preferences for amino acids and amines. The overall metabolic dissimilarity among the lakes was approximately 32%, primarily driven by the CGs of amino acids and amines. The most preferred CSs included the carbohydrates D-Mannitol and N-Acetyl-D-glucosamine, and the carboxylic acid D-Galacturonic acid. Bacterial functional diversity was high (H′: 3.18 − 3.33) due to the high evenness of CS utilization. Multidimensional analyses revealed significant influences of lake-specific characteristics, sampling month, and the interaction between these factors on CLPP patterns. Nutrients and water temperature were identified as the most influential environmental factors affecting bacterial communities. This study enhances our understanding of bacterial functional responses to changing environments, demonstrating the effectiveness of the Biolog EcoplateTM technique in providing insights into these changes.

Genetic features of avian influenza (A/H5N8) clade 2.3.4.4b isolated from quail in Egypt

Several genotypes of the highly pathogenic avian influenza (HPAI) virus H5N8 subtype within clade 2.3.4.4b continue to circulate in different species of domestic birds across Egypt. It is believed that quail contribute to virus replication and adaptation to other gallinaceous poultry species and humans. This study provides genetic characterization of the full genome of HPAI H5N8 isolated from quail in Egypt. The virus was isolated from a commercial quail farm associated with respiratory signs. To characterize the genetic features of the detected virus, gene sequencing via Sanger technology and phylogenetic analysis were performed. The results revealed high nucleotide identity with the HPAI H5N8 virus from Egypt, which has multiple basic amino acid motifs PLREKRRKR/GLF at the hemagglutinin (HA) cleavage site. Phylogenetic analysis of the eight gene segments revealed that the quail isolate is grouped with HPAI H5N8 viruses of clade 2.3.4.4b and closely related to the most recent circulating H5N8 viruses in Egypt. Whole-genome characterization revealed amino acid preferences for avian receptors with few mutations, indicating their affinity for human-like receptors and increased virulence in mammals, such as S123P, S133A, T156A and A263T in the HA gene. In addition, the sequencing results revealed a lack of markers associated with influenza antiviral resistance in the neuraminidase and matrix-2 coding proteins. The results of the present study support the spread of HPAIV H5N8 to species other than chickens in Egypt. Therefore, continuous surveillance of AIV in different bird species in Egypt followed by full genomic characterization is needed for better virus control and prevention.

The Change Rate of the Fbxl21 Gene and the Amino Acid Composition of Its Protein Correlate with the Species-Specific Lifespan in Placental Mammals.

This article proposes a methodology for establishing a relationship between the change rate of a given gene (relative to a given taxon) together with the amino acid composition of the proteins encoded by this gene and the traits of the species containing this gene. The methodology is illustrated based on the mammalian genes responsible for regulating the circadian rhythms that underlie a number of human disorders, particularly those associated with aging. The methods used are statistical and bioinformatic ones. A systematic search for orthologues, pseudogenes, and gene losses was performed using our previously developed methods. It is demonstrated that the least conserved Fbxl21 gene in the Euarchontoglires superorder exhibits a statistically significant connection of genomic characteristics (the median of dN/dS for a gene relative to all the other orthologous genes of a taxon, as well as the preference or avoidance of certain amino acids in its protein) with species-specific lifespan and body weight. In contrast, no such connection is observed for Fbxl21 in the Laurasiatheria superorder. This study goes beyond the protein-coding genes, since the accumulation of amino acid substitutions in the course of evolution leads to pseudogenization and even gene loss, although the relationship between the genomic characteristics and the species traits is still preserved. The proposed methodology is illustrated using the examples of circadian rhythm genes and proteins in placental mammals, e.g., longevity is connected with the rate of Fbxl21 gene change, pseudogenization or gene loss, and specific amino acid substitutions (e.g., asparagine at the 19th position of the CRY-binding domain) in the protein encoded by this gene.

Evaluating the impact of Torulaspora delbrueckii and amino acid concentration on the nitrogen metabolism of Oenococcus oeni

This study investigates the impact of Saccharomyces cerevisiae and Torulaspora delbrueckii inoculation strategies and amino acid supplementation on nitrogen metabolism and malolactic fermentation (MLF) performance in Oenococcus oeni. Oenococcus oeni has low demand of nitrogen and typically prefers peptides over amino acids. In this work we studied the nitrogen metabolism (proteins, peptides, amino acids and biogenic amines (BA)) of three O. oeni strains. Higher initial amino acid concentrations generally accelerated MLF, particularly with the PSU-1 strain, which exhibited stuck fermentation in some conditions with low initial amino acid concentration related with coinoculation of yeasts. Nitrogen metabolism in O. oeni showed a preference for peptide-bound amino acids, and gene expression analyses highlighted a general upregulation in response to increased amino acid concentrations. BA production, particularly cadaverine and putrescine, was strain-dependent and associated with the presence of the odc gene. Besides, the production of 2-phenylethylamine was related with alcoholic fermentation (AF) inoculation strategy. Overall, the study demonstrates the complexity of nitrogen metabolism in O. oeni, emphasizing that it is influenced by both the yeast species used in AF and the initial amino acid composition, highlighting the importance of peptide composition, offering a potential tool for optimizing MLF in winemaking.

Effects of two fillers and process conditions on the water treatment efficiency of a continuous packed bed biofilm reactor.

This study evaluated the treatment efficiency of two selected fillers and their combination for improving the water quality of aquaculture wastewater using a packed bed biofilm reactor (PBBR) under various process conditions. The fillers used were nanosheet (NS), activated carbon (AC), and a combination of both. The results indicated that the use of combined fillers and the hydraulic retention time (HRT) of 4 h significantly enhanced water quality in the PBBR. The removal rates of chemical oxygen demand, NO2-─N, total suspended solids（TSS）, and chlorophyll a were 63.55%, 74.25%, 62.75%, and 92.85%, respectively. The microbiota analysis revealed that the presence of NS increased the abundance of microbial phyla associated with nitrogen removal, such as Nitrospirae and Proteobacteria. The difference between the M1 and M2 communities was minimal. Additionally, the microbiota in different PBBR samples displayed similar preferences for carbon sources, and carbohydrates and amino acids were the most commonly utilized carbon sources by microbiota. These results indicated that the combination of NS and AC fillers in a PBBR effectively enhanced the treatment efficiency of aquaculture wastewater when operated at an HRT of 4 h. The findings provide valuable insights into optimizing the design of aquaculture wastewater treatment systems.

The effectiveness of selection in a species affects the direction of amino acid frequency evolution.

Nearly neutral theory predicts that species with higher effective population size (N e ) are better able to purge slightly deleterious mutations. We compare evolution in high-N e vs. low-N e vertebrates to reveal which amino acid frequencies are subject to subtle selective preferences. We take three complementary approaches, two measuring flux and one measuring outcomes. First, we fit non-stationary substitution models of amino acid flux using maximum likelihood, comparing the high-N e clade of rodents and lagomorphs to its low-N e sister clade of primates and colugos. Second, we compare evolutionary outcomes across a wider range of vertebrates, via correlations between amino acid frequencies and N e . Third, we dissect the details of flux in human, chimpanzee, mouse, and rat, as scored by parsimony - this also enables comparison to a historical paper. All three methods agree on which amino acids are preferred under more effective selection. Preferred amino acids tend to be smaller, less costly to synthesize, and to promote intrinsic structural disorder. Parsimony-induced bias in the historical study produces an apparent reduction in structural disorder, perhaps driven by slightly deleterious substitutions. Within highly exchangeable pairs of amino acids, arginine is strongly preferred over lysine, and valine over isoleucine, consistent with more effective selection preferring a marginally larger free energy of folding. These two preferences match differences between thermophiles and mesophilic relatives. These results reveal the biophysical consequences of mutation-selection-drift balance, and demonstrate the utility of nearly neutral theory for understanding protein evolution.

Chloroplast Genomes of Vitis flexuosa and Vitis amurensis: Molecular Structure, Phylogenetic, and Comparative Analyses for Wild Plant Conservation.

The chloroplast genome plays a crucial role in elucidating genetic diversity and phylogenetic relationships. Vitis vinifera L. (grapevine) is an economically important species, prompting exploration of wild genetic resources to enhance stress resilience. We meticulously assembled the chloroplast genomes of two Korean Vitis L. species, V. flexuosa Thunb. and V. amurensis Rupr., contributing valuable data to the Korea Crop Wild Relatives inventory. Through exhaustive specimen collection spanning diverse ecological niches across South Korea, we ensured comprehensive representation of genetic diversity. Our analysis, which included rigorous codon usage bias assessment and repeat analysis, provides valuable insights into amino acid preferences and facilitates the identification of potential molecular markers. The assembled chloroplast genomes were subjected to meticulous annotation, revealing divergence hotspots enriched with nucleotide diversity, thereby presenting promising candidates for DNA barcodes. Additionally, phylogenetic analysis reaffirmed intra-genus relationships and identified related crops, shedding light on evolutionary patterns within the genus. Comparative examination with chloroplast genomes of other crops uncovered conserved sequences and variable regions, offering critical insights into genetic evolution and adaptation. Our study advances the understanding of chloroplast genomes, genetic diversity, and phylogenetic relationships within Vitis species, thereby laying a foundation for enhancing grapevine genetic diversity and resilience to environmental challenges.

Exploration of the binding determinants of protein phosphatase 5 (PP5) reveals a chaperone-independent activation mechanism

The protein phosphatase 5 (PP5) is normally recruited to its substrates by the molecular chaperones, heat shock protein 70 (Hsp70) and heat shock protein 90 (Hsp90). This interaction requires the tetratricopeptide repeat (TPR) domain of PP5, which binds to an EEVD motif at the extreme C-termini of cytosolic Hsp70 and Hsp90 isoforms. In addition to bringing PP5 into proximity with chaperone-bound substrates, this interaction also relieves auto-inhibition in PP5’s catalytic domain, promoting its phosphatase activity. To better understand the molecular determinants of this process, we screened a large, pentapeptide library for binding to PP5. This screen identified the amino acid preferences at each position, which we validated by showing that the optimal sequences bind 4- to 7-fold tighter than the natural EEVD motifs and stimulate PP5’s enzymatic activity. The enhanced affinity for PP5’s TPR domain was confirmed using a protein-adaptive differential scanning fluorimetry (paDSF) assay. Using this increased knowledge of structure-activity relationships, we re-examined affinity proteomics results to look for potential EEVD-like motifs in the C-termini of known PP5-binding partners. This search identified elongator acetyltransferase complex subunit 1 (ELP1/IKBKAP) as a putative partner and, indeed, we found that its C-terminal sequence, LSLLD, binds directly to PP5’s TPR domain in vitro. Consistent with this idea, mutation of ELP1’s terminal aspartate was sufficient to interrupt the interaction with PP5 in vitro and in cells. Together, these findings reveal the sequence preferences of PP5’s TPR domain and expand the scope of PP5’s functions to include chaperone-independent complexes.

Improving aromatic higher alcohol acetates in wines by co-fermentation of Pichia kluyveri and Saccharomyces cerevisiae: growth interaction and amino acid competition.

Higher alcohol acetates (HAAs) are potent aroma-active esters that impart desirable fruity and floral aromas. However, the conversion of higher alcohol precursors into HAAs is extremely low in winemaking. To investigate the underlying yeast-yeast interaction on targeted improvement of aromatic HAAs, we evaluated fermentation activity, cell viability, amino acid consumption and HAA production when Pichia kluyveri and Saccharomyces cerevisiae were inoculated concurrently or sequentially. Pichia kluyveri PK-21 possessed the ability to survive and increased HAA level up to 5.2-fold in mixed fermentation. Such an increment may benefit from the efficient conversion of higher alcohol precursors into HAAs (>27-fold higher than S. cerevisiae). During mixed fermentation, the two yeasts exhibited crucial interactions regarding cell growth and amino acid competition. Saccharomyces cerevisiae dominated over the co-inoculated P. kluyveri by efficient uptake of amino acids and biomass production. However, this dominance decreased in sequential fermentation, where P. kluyveri growth increased due to the consumption of preferred amino acids prior to S. cerevisiae. Pearson correlation analysis indicated that phenylalanine and aspartic acid may act as positive amino acids in boosting P. kluyveri growth and HAA production. Laboratory-scale winemaking validated the fermentation performance of P. kluyveri in sequential inoculum, resulting in a balanced aroma profile with enhanced floral and tropical fruity characteristics in the final wines. This study proposes a microbial, non-genetically engineered approach for targeted increase of HAA production in winemaking and the findings provide new insights into yeast-yeast interactions. © 2024 Society of Chemical Industry.

Conformational Preferences of the Non-Canonical Amino Acids (2S,4S)-5-Fluoroleucine, (2S,4R)-5-Fluoroleucine, and 5,5'-Difluoroleucine in a Protein.

Proteins produced with leucine analogues, where CH2F groups substitute specific methyl groups, can readily be probed by 19F NMR spectroscopy. As CF and CH groups are similar in hydrophobicity and size, fluorinated leucines are expected to cause minimal structural perturbation, but the impact of fluorine on the rotational freedom of CH2F groups is unclear. We present high-resolution crystal structures of Escherichia coli peptidyl-prolyl cis-trans isomerase B (PpiB) prepared with uniform high-level substitution of leucine by (2S,4S)-5-fluoroleucine, (2S,4R)-5-fluoroleucine, or 5,5'-difluoroleucine. Apart from the fluorinated leucine residues, the structures show complete structural conservation of the protein backbone and the amino acid side chains except for a single isoleucine side chain located next to a fluorine atom in the hydrophobic core of the protein. The carbon skeletons of the fluorinated leucine side chains are also mostly conserved. The CH2F groups show a strong preference for staggered rotamers and often appear locked into single rotamers. Substitution of leucine CH3 groups for CH2F groups is thus readily tolerated in the three-dimensional (3D) structure of a protein, and the rotation of CH2F groups can be halted at cryogenic temperatures.

Changing selection on amino acid substitutions in Gag protein between major HIV-1 subtypes.

Amino acid preferences at a protein site depend on the role of this site in protein function and structure as well as on external constraints. All these factors can change in the course of evolution, making amino acid propensities of a site time-dependent. When viral subtypes divergently evolve in different host subpopulations, such changes may depend on genetic, medical, and sociocultural differences between these subpopulations. Here, using our previously developed phylogenetic approach, we describe sixty-nine amino acid sites of the Gag protein of human immunodeficiency virus type 1 (HIV-1) where amino acids have different impact on viral fitness in six major subtypes of the type M. These changes in preferences trigger adaptive evolution; indeed, 32 (46 per cent) of these sites experienced strong positive selection at least in one of the subtypes. At some of the sites, changes in amino acid preferences may be associated with differences in immune escape between subtypes. The prevalence of an amino acid in a protein site within a subtype is only a poor predictor for whether this amino acid is preferred in this subtype according to the phylogenetic analysis. Therefore, attempts to identify the factors of viral evolution from comparative genomics data should integrate across multiple sources of information.

P1' specificity of the S219V/R203G mutant tobacco etch virus protease.

Proteases that recognize linear amino acid sequences with high specificity became indispensable tools of recombinant protein technology for the removal of various fusion tags. Due to its stringent sequence specificity, the catalytic domain of the nuclear inclusion cysteine protease of tobacco etch virus (TEV PR) is also a widely applied reagent for enzymatic removal of fusion tags. For this reason, efforts have been made to improve its stability and modify its specificity. For example, P1' autoproteolytic cleavage-resistant mutant (S219V) TEV PR was found not only to be nearly impervious to self-inactivation, but also exhibited greater stability and catalytic efficiency than the wild-type enzyme. An R203G substitution has been reported to further relax the P1' specificity of the enzyme, however, these results were obtained from crude intracellular assays. Until now, there has been no rigorous comparison of the P1' specificity of the S219V and S219V/R203G mutants in vitro, under carefully controlled conditions. Here, we compare the P1' amino acid preferences of these single and double TEV PR mutants. The in vitro analysis was performed by using recombinant protein substrates representing 20 P1' variants of the consensus TENLYFQ*SGT cleavage site, and synthetic oligopeptide substrates were also applied to study a limited set of the most preferred variants. In addition, the enzyme-substrate interactions were analyzed in silico. The results indicate highly similar P1' preferences for both enzymes, many side-chains can be accommodated by the S1' binding sites, but the kinetic assays revealed lower catalytic efficiency for the S219V/R203G than for the S219V mutant.

Amino acid preference and fermentation performance of Pichia kluyveri strains in a synthetic wort

Pichia kluyveri is a maltose-negative yeast applied for non- and low-alcoholic beer (NABLAB) production due to its fruity aroma attributes, however, related literature is still limited. Thus, this study aims to investigate basic fermentation parameters, such as growth, sugar consumption, and production of primary fermentation metabolites, but also the amino acid preference of three different P. kluyveri strains in an agitated model brewers’ wort substrate. Furthermore, the effect of agitation on these parameters was studied for one of the strains. The three strains assessed, PKCH_01, PKCH_02 and PKCH_03, shared a glucophilic nature, a cell growth from 4E+07 to 6E+07 E+07 CFU/mL, similar amino acid preference, and equivalent ethanol yield factor (Yse) from 0.45 to 0.5 g ethanol/g glucose. In addition, agitation increased the cell growth yield factor by 78% and favored the uptake of amino acids, minimizing their release to the medium. Interestingly, a strong preference for methionine (above 80% consumption) by P. kluyveri in a beer-related medium was reported for the first time, whereas glutamine and glutamate were within the low-preference group (below 15% consumption). This study provides a better understanding about fermentation and physiological attributes of P. kluyveri, thereby revealing its potential for brewing applications.

How Do Molecular Tweezers Bind to Proteins? Lessons from X-ray Crystallography.

To understand the biological relevance and mode of action of artificial protein ligands, crystal structures with their protein targets are essential. Here, we describe and investigate all known crystal structures that contain a so-called "molecular tweezer" or one of its derivatives with an attached natural ligand on the respective target protein. The aromatic ring system of these compounds is able to include lysine and arginine side chains, supported by one or two phosphate groups that are attached to the half-moon-shaped molecule. Due to their marked preference for basic amino acids and the fully reversible binding mode, molecular tweezers are able to counteract pathologic protein aggregation and are currently being developed as disease-modifying therapies against neurodegenerative diseases such as Alzheimer's and Parkinson's disease. We analyzed the corresponding crystal structures with 14-3-3 proteins in complex with mono- and diphosphate tweezers. Furthermore, we solved crystal structures of two different tweezer variants in complex with the enzyme Δ1-Pyrroline-5-carboxyl-dehydrogenase (P5CDH) and found that the tweezers are bound to a lysine and methionine side chain, respectively. The different binding modes and their implications for affinity and specificity are discussed, as well as the general problems in crystallizing protein complexes with artificial ligands.

Residue K28 of Zika Virus NS5 Protein Is Implicated in Virus Replication and Antagonism of STAT2.

The identification of four potential nonstructural 5 (NS5) residues-K28, K45, V335, and S749-that share the same amino acid preference in STAT2-interacting flaviviruses [Dengue virus (DENV) and Zika virus (ZIKV)], but not in STAT2-non-interacting flaviviruses [West Nile virus (WNV) and/or Yellow fever virus (YFV)] from an alignment of multiple flavivirus NS5 sequences, implied a possible association with the efficiency of ZIKV to antagonize the human signal transducer and activator of transcription factor 2 (STAT2). Through site-directed mutagenesis and reverse genetics, mutational impacts of these residues on ZIKV growth in vitro and STAT2 antagonism were assessed using virus growth kinetics assays and STAT2 immunoblotting. The results showed that mutations at the residue K28 significantly reduced the efficiency of ZIKV to antagonize STAT2. Further investigation involving residue K28 demonstrated its additional effects on the phenotypes of ZIKV-NS5 nuclear bodies. These findings demonstrate that K28, identified from sequence alignment, is an important determinant of replication and STAT2 antagonism by ZIKV.

Substrate profiling of human transglutaminase 1 using cDNA display and next-generation sequencing.

Human transglutaminase 1 (TG1) modulates skin development, while its involvement in diseases remains poorly understood, necessitating comprehensive exploration of its substrate interactions. To study the substrate profile of TG1, an in vitro selection system based on cDNA display technology was used to screen two peptide libraries with mutations at varying distance from the reactive glutamine. Next-generation sequencing and bioinformatics analysis of the selected DNA pools revealed a detailed TG1 substrate profile, indicating preferred and non-preferred amino acid sequences. The peptide sequence, AEQHKLPSKWPF, was identified showing high reactivity and specificity to TG1. The position weight matrix calculated from the per amino acid enrichment factors was employed to search human proteins using an in-house algorithm, revealing six known TG1 substrate proteins with high scores, alongside a list of candidate substrates currently under investigation. Our findings are expected to assist in future medical diagnoses and development of treatments for skin disorders.

Protein design using structure-based residue preferences

Recent developments in protein design rely on large neural networks with up to 100s of millions of parameters, yet it is unclear which residue dependencies are critical for determining protein function. Here, we show that amino acid preferences at individual residues—without accounting for mutation interactions—explain much and sometimes virtually all of the combinatorial mutation effects across 8 datasets (R2 ~ 78-98%). Hence, few observations (~100 times the number of mutated residues) enable accurate prediction of held-out variant effects (Pearson r > 0.80). We hypothesized that the local structural contexts around a residue could be sufficient to predict mutation preferences, and develop an unsupervised approach termed CoVES (Combinatorial Variant Effects from Structure). Our results suggest that CoVES outperforms not just model-free methods but also similarly to complex models for creating functional and diverse protein variants. CoVES offers an effective alternative to complicated models for identifying functional protein mutations.

Learning the shape of protein microenvironments with a holographic convolutional neural network

Proteins play a central role in biology from immune recognition to brain activity. While major advances in machine learning have improved our ability to predict protein structure from sequence, determining protein function from its sequence or structure remains a major challenge. Here, we introduce holographic convolutional neural network (H-CNN) for proteins, which is a physically motivated machine learning approach to model amino acid preferences in protein structures. H-CNN reflects physical interactions in a protein structure and recapitulates the functional information stored in evolutionary data. H-CNN accurately predicts the impact of mutations on protein stability and binding of protein complexes. Our interpretable computational model for protein structure-function maps could guide design of novel proteins with desired function.

Extended Cleavage Specificity of two Hematopoietic Serine Proteases from a Ray-Finned Fish, the Spotted Gar (Lepisosteus oculatus).

The extended cleavage specificities of two hematopoietic serine proteases originating from the ray-finned fish, the spotted gar (Lepisosteus oculatus), have been characterized using substrate phage display. The preference for particular amino acids at and surrounding the cleavage site was further validated using a panel of recombinant substrates. For one of the enzymes, the gar granzyme G, a strict preference for the aromatic amino acid Tyr was observed at the cleavable P1 position. Using a set of recombinant substrates showed that the gar granzyme G had a high selectivity for Tyr but a lower activity for cleaving after Phe but not after Trp. Instead, the second enzyme, gar DDN1, showed a high preference for Leu in the P1 position of substrates. This latter enzyme also showed a high preference for Pro in the P2 position and Arg in both P4 and P5 positions. The selectivity for the two Arg residues in positions P4 and P5 suggests a highly specific substrate selectivity of this enzyme. The screening of the gar proteome with the consensus sequences obtained by substrate phage display for these two proteases resulted in a very diverse set of potential targets. Due to this diversity, a clear candidate for a specific immune function of these two enzymes cannot yet be identified. Antisera developed against the recombinant gar enzymes were used to study their tissue distribution. Tissue sections from juvenile fish showed the expression of both proteases in cells in Peyer's patch-like structures in the intestinal region, indicating they may be expressed in T or NK cells. However, due to the lack of antibodies to specific surface markers in the gar, it has not been possible to specify the exact cellular origin. A marked difference in abundance was observed for the two proteases where gar DDN1 was expressed at higher levels than gar granzyme G. However, both appear to be expressed in the same or similar cells, having a lymphocyte-like appearance.