Antimicrobial Peptide Databases as the Guiding Resource in New Antimicrobial Agent Identification via Computational Methods
1.1. Antimicrobial Peptides (AMPs)—Natural Defenders with Therapeutic Applications Antimicrobial peptides constitute a crucial component of innate immunity, and are present in the majority of living organisms. They possess the ability to induce immune responses and exhibit activity against a wide range of microorganisms, including bacteria, fungi, and viruses. Moreover, AMPs can target cancer cells, making them potential therapeutic agents for cancer treatment . AMPs are peptide compounds composed of amino acids, typically containing 6 to 100 residues, although peptides with longer sequences also exist. They occur in four main structural forms: α-helices, β-sheets, extended conformations, and loops. Under natural conditions, α-helical and β-sheet structures predominate. Linear peptides lose their structure in solution, whereas cyclic peptides form stable β-sheets maintained by disulfide bonds. Antimicrobial peptides are characterized by a cationic charge and a high content of hydrophobic residues, which contribute to their innate activity against various pathogens and potential applications in therapeutic contexts . AMPs are widely distributed in nature and are found in all living organisms. An example of an AMP in humans is the LL-37 peptide, also known as cathelicidin . Other examples of human AMPs include cathelicidins, such as FALL-39 and defensins, which are divided into α-defensins (e.g., HNP-1, HNP-2, HNP-3, and HNP-4, primarily found in neutrophils) and β-defensins (e.g., hBD-1, hBD-2, hBD-3, and hBD-4, expressed in epithelial tissues) . 1.2. Mechanisms of Action and Immunomodulatory Properties of AMPs Antimicrobial peptides interact with the cell membranes of microorganisms, leading to their permeabilization and lysis through mechanisms such as the “barrel-stave” model, the toroidal pore model, and the carpet model . They also act at the intracellular level, disrupting key processes such as protein synthesis, nucleic acid synthesis, and cell wall component biosynthesis, ultimately causing structural destabilization and cell death . Furthermore, AMPs exhibit immunomodulatory activity, supporting the recruitment and activation of immune cells and stimulating the adaptive immune response, thereby enhancing their effectiveness in defending against infections . Their application as antibiotics is more efficient and safer, as they are less toxic and do not cause side effects . 1.3. Classification and Sources of Antimicrobial Peptides AMPs are classified based on their source, activity, structural characteristics, and amino acid content. Regarding their sources, AMPs are divided into peptides derived from mammals, amphibians, insects, and microorganisms . In mammals, the primary AMP families are cathelicidins and defensins, which protect the organism against infections and exhibit varying expression levels depending on the stage of life, such as in breast milk, which supports infant health . Amphibian-derived peptides, particularly from frogs, play a key role in pathogen defense, with magainin being one of the most well-known examples . In insects, AMPs such as cecropins are synthesized in fat bodies and hemocytes, demonstrating both anti-inflammatory and anticancer activities . Meanwhile, microorganisms, including bacteria and fungi, are sources of AMPs like nisin and gramicidin, which have industrial applications, although their production can be expensive. AMPs possess broad protective properties and potential applications in various fields, including medicine . 1.4. Antimicrobial Peptide Databases Antimicrobial peptide databases are a crucial tool in supporting research on their properties, structure, and function. They provide quick access to data such as amino acid sequences, biological sources, physicochemical properties, antibacterial activity, and three-dimensional peptide structures. These data are sourced from the scientific literature, databases like UniProt and NCBI, and experimental results. The database creation process involves collecting, filtering, and verifying information, using tools such as CD-HIT to ensure high data quality . AMPs are classified based on function (e.g., antibacterial, antiviral), structure, or source of origin, facilitating more precise analysis. Some popular databases, such as dbAMP, DBAASP, and DRAMP, offer advanced analytical tools, including machine learning (ML) algorithms. ML models, trained on data from these databases, enable the prediction of peptide activity, the identification of new sequences, and the design of peptides with desired properties, which contributes to combating bacterial resistance to antibiotics. The effectiveness of these models is assessed using metrics like accuracy and specificity, enabling the further optimization and refinement of the research processes. The types of databases are presented in . 1.5. Antimicrobial Peptide Databases—Availability The selection of AMP databases for analysis was based on the article “A Review on Antimicrobial Peptides Databases and the Computational Tools”. The article presents AMPs as a promising alternative to traditional antibiotics. It includes a review of antimicrobial peptide databases. The databases discussed in the article are presented in . 1.6. Antimicrobial Peptide Databases—Reference Database The dbAMP is a comprehensive database containing information on antimicrobial peptides, including their sequences, biological activity, post-translational modifications (PTM), structure, and physicochemical properties . In 2024, the database was updated to version 3.0, enriching it with new analytical tools, including a text analysis system based on natural language processing, which enables the automatic identification of publications related to AMPs . This update also introduced extensive data on 53 functional activities of peptides, and information on peptides designed based on proteomic and transcriptomic data . The dbAMP database was selected as the reference database for the analysis. It stands out from the others primarily due to having the largest number of registered antimicrobial peptides, making it one of the most comprehensive available databases in this field. Regular updates, such as the addition of 9078 AMPs compared to the year 2022, indicate the dynamic growth of this platform. An additional advantage is the precise taxonomy of peptides, which allows for detailed analysis of their diversity, as well as a significant percentage of peptides derived from vertebrates ( ). Thanks to these features, the database is ideally suited to be a reference database. 1.7. Antimicrobial Peptide Databases—Review The selection of the databases listed below for the analysis described in the article is based on their diversity in terms of the number of registered peptides and their types, which allows for a more comprehensive picture of the compared sequences ( ). Each of these databases specializes in different groups of peptides, such as synthetic peptides, antiparasitic peptides, or those exhibiting activity against biofilms, making the comparison of their contents more thorough. As a result, the analysis of these databases facilitates a better understanding of peptide diversity and enables the formulation of broader conclusions. APD is one of the first databases on natural antimicrobial peptides, containing information on the sequences and activities of peptides derived from various organisms. Peptides are classified according to their biological properties, such as antibacterial, antiviral, or anticancer activities ; BaAMPs focuses on peptides tested against biofilms, providing verified experimental data ; CAMP is a database that collects data on AMP sequences, their origin, and biological activity, including synthetic AMPs ; CyBase offers data on cyclic proteins, supporting research on their structures and functions ; dadp collects data on antimicrobial peptides, focusing on precursor sequences and their bioactive fragments ; DBAASP collects data on antimicrobial peptides, providing information on their structures, conditions of action, and molecular targets. It also includes predictive tools supporting peptide design ; DRAMP is a database containing peptides with defined sequences, categorized into general and patent sets, with data on toxicity and hemolytic activity ; InverPep focuses on peptides derived from invertebrates, offering data on sequences, structures, and physicochemical properties ; ParaPep specializes in antiprotozoal peptides, offering information on their structures and mechanisms of action ; SATPdb contains data on therapeutic peptides, enabling sequence and structure similarity searches . 1.8. Diamond–High-Throughput Protein Alignment Among the various amino acid sequence comparison tools used in bioinformatics analyses, this study utilized diamond software as the main tool. It is an advanced open-source tool designed for the rapid alignment of DNA and protein sequences within large databases. Diamond delivers a performance significantly superior to that of the traditional methods, such as BLASTX, by employing a simplified protein alphabet and advanced optimization technologies that accelerate the analysis process with minimal loss of precision . The tool enables the efficient processing of large datasets, though it requires substantial computational resources. Diamond is particularly valuable in metagenomic projects and the identification of potential antimicrobial peptides (AMPs), though it necessitates users possessing the necessary expertise for the proper configuration and interpretation of results . This study aims to evaluate the quality of information within selected AMP databases, considering their accessibility, content, and research potential. Below, we present the results obtained via analysis and a comparison of per-database and cross-database peptide sequences. A diamond , high-throughput protein alignment program was used to compare the degree of sequence similarity among peptides across the individual databases. The redundancy of the data was also evaluated. As a proof of concept, the proposed methodology was used for an in silico evaluation of the selected species proteomes in order to identify putative antimicrobial peptide candidates. An example candidate was further evaluated via a combination of structural analysis based on the computed homology-based structural model, the in silico digestion of the source protein, and the antimicrobial potential. 

Antimicrobial peptides constitute a crucial component of innate immunity, and are present in the majority of living organisms. They possess the ability to induce immune responses and exhibit activity against a wide range of microorganisms, including bacteria, fungi, and viruses. Moreover, AMPs can target cancer cells, making them potential therapeutic agents for cancer treatment . AMPs are peptide compounds composed of amino acids, typically containing 6 to 100 residues, although peptides with longer sequences also exist. They occur in four main structural forms: α-helices, β-sheets, extended conformations, and loops. Under natural conditions, α-helical and β-sheet structures predominate. Linear peptides lose their structure in solution, whereas cyclic peptides form stable β-sheets maintained by disulfide bonds. Antimicrobial peptides are characterized by a cationic charge and a high content of hydrophobic residues, which contribute to their innate activity against various pathogens and potential applications in therapeutic contexts . AMPs are widely distributed in nature and are found in all living organisms. An example of an AMP in humans is the LL-37 peptide, also known as cathelicidin . Other examples of human AMPs include cathelicidins, such as FALL-39 and defensins, which are divided into α-defensins (e.g., HNP-1, HNP-2, HNP-3, and HNP-4, primarily found in neutrophils) and β-defensins (e.g., hBD-1, hBD-2, hBD-3, and hBD-4, expressed in epithelial tissues) . 

Antimicrobial peptides interact with the cell membranes of microorganisms, leading to their permeabilization and lysis through mechanisms such as the “barrel-stave” model, the toroidal pore model, and the carpet model . They also act at the intracellular level, disrupting key processes such as protein synthesis, nucleic acid synthesis, and cell wall component biosynthesis, ultimately causing structural destabilization and cell death . Furthermore, AMPs exhibit immunomodulatory activity, supporting the recruitment and activation of immune cells and stimulating the adaptive immune response, thereby enhancing their effectiveness in defending against infections . Their application as antibiotics is more efficient and safer, as they are less toxic and do not cause side effects . 

AMPs are classified based on their source, activity, structural characteristics, and amino acid content. Regarding their sources, AMPs are divided into peptides derived from mammals, amphibians, insects, and microorganisms . In mammals, the primary AMP families are cathelicidins and defensins, which protect the organism against infections and exhibit varying expression levels depending on the stage of life, such as in breast milk, which supports infant health . Amphibian-derived peptides, particularly from frogs, play a key role in pathogen defense, with magainin being one of the most well-known examples . In insects, AMPs such as cecropins are synthesized in fat bodies and hemocytes, demonstrating both anti-inflammatory and anticancer activities . Meanwhile, microorganisms, including bacteria and fungi, are sources of AMPs like nisin and gramicidin, which have industrial applications, although their production can be expensive. AMPs possess broad protective properties and potential applications in various fields, including medicine . 

Antimicrobial peptide databases are a crucial tool in supporting research on their properties, structure, and function. They provide quick access to data such as amino acid sequences, biological sources, physicochemical properties, antibacterial activity, and three-dimensional peptide structures. These data are sourced from the scientific literature, databases like UniProt and NCBI, and experimental results. The database creation process involves collecting, filtering, and verifying information, using tools such as CD-HIT to ensure high data quality . AMPs are classified based on function (e.g., antibacterial, antiviral), structure, or source of origin, facilitating more precise analysis. Some popular databases, such as dbAMP, DBAASP, and DRAMP, offer advanced analytical tools, including machine learning (ML) algorithms. ML models, trained on data from these databases, enable the prediction of peptide activity, the identification of new sequences, and the design of peptides with desired properties, which contributes to combating bacterial resistance to antibiotics. The effectiveness of these models is assessed using metrics like accuracy and specificity, enabling the further optimization and refinement of the research processes. The types of databases are presented in . 

The selection of AMP databases for analysis was based on the article “A Review on Antimicrobial Peptides Databases and the Computational Tools”. The article presents AMPs as a promising alternative to traditional antibiotics. It includes a review of antimicrobial peptide databases. The databases discussed in the article are presented in . 

The dbAMP is a comprehensive database containing information on antimicrobial peptides, including their sequences, biological activity, post-translational modifications (PTM), structure, and physicochemical properties . In 2024, the database was updated to version 3.0, enriching it with new analytical tools, including a text analysis system based on natural language processing, which enables the automatic identification of publications related to AMPs . This update also introduced extensive data on 53 functional activities of peptides, and information on peptides designed based on proteomic and transcriptomic data . The dbAMP database was selected as the reference database for the analysis. It stands out from the others primarily due to having the largest number of registered antimicrobial peptides, making it one of the most comprehensive available databases in this field. Regular updates, such as the addition of 9078 AMPs compared to the year 2022, indicate the dynamic growth of this platform. An additional advantage is the precise taxonomy of peptides, which allows for detailed analysis of their diversity, as well as a significant percentage of peptides derived from vertebrates ( ). Thanks to these features, the database is ideally suited to be a reference database. 

The selection of the databases listed below for the analysis described in the article is based on their diversity in terms of the number of registered peptides and their types, which allows for a more comprehensive picture of the compared sequences ( ). Each of these databases specializes in different groups of peptides, such as synthetic peptides, antiparasitic peptides, or those exhibiting activity against biofilms, making the comparison of their contents more thorough. As a result, the analysis of these databases facilitates a better understanding of peptide diversity and enables the formulation of broader conclusions. APD is one of the first databases on natural antimicrobial peptides, containing information on the sequences and activities of peptides derived from various organisms. Peptides are classified according to their biological properties, such as antibacterial, antiviral, or anticancer activities ; BaAMPs focuses on peptides tested against biofilms, providing verified experimental data ; CAMP is a database that collects data on AMP sequences, their origin, and biological activity, including synthetic AMPs ; CyBase offers data on cyclic proteins, supporting research on their structures and functions ; dadp collects data on antimicrobial peptides, focusing on precursor sequences and their bioactive fragments ; DBAASP collects data on antimicrobial peptides, providing information on their structures, conditions of action, and molecular targets. It also includes predictive tools supporting peptide design ; DRAMP is a database containing peptides with defined sequences, categorized into general and patent sets, with data on toxicity and hemolytic activity ; InverPep focuses on peptides derived from invertebrates, offering data on sequences, structures, and physicochemical properties ; ParaPep specializes in antiprotozoal peptides, offering information on their structures and mechanisms of action ; SATPdb contains data on therapeutic peptides, enabling sequence and structure similarity searches . 

Among the various amino acid sequence comparison tools used in bioinformatics analyses, this study utilized diamond software as the main tool. It is an advanced open-source tool designed for the rapid alignment of DNA and protein sequences within large databases. Diamond delivers a performance significantly superior to that of the traditional methods, such as BLASTX, by employing a simplified protein alphabet and advanced optimization technologies that accelerate the analysis process with minimal loss of precision . The tool enables the efficient processing of large datasets, though it requires substantial computational resources. Diamond is particularly valuable in metagenomic projects and the identification of potential antimicrobial peptides (AMPs), though it necessitates users possessing the necessary expertise for the proper configuration and interpretation of results . This study aims to evaluate the quality of information within selected AMP databases, considering their accessibility, content, and research potential. Below, we present the results obtained via analysis and a comparison of per-database and cross-database peptide sequences. A diamond , high-throughput protein alignment program was used to compare the degree of sequence similarity among peptides across the individual databases. The redundancy of the data was also evaluated. As a proof of concept, the proposed methodology was used for an in silico evaluation of the selected species proteomes in order to identify putative antimicrobial peptide candidates. An example candidate was further evaluated via a combination of structural analysis based on the computed homology-based structural model, the in silico digestion of the source protein, and the antimicrobial potential. 

During the analysis, in all cases, a greater number of comparisons were obtained than the sum of sequences in the compared databases. This results from the methodology used for sequence comparison. When a reference database, containing, for example, 5000 sequences, is compared with a target database with 3000 sequences, each comparison involves matching each sequence from the reference database with every sequence in the target database. For example, sequence A1 from the reference database is compared sequentially with all sequences in the target database: A1 with B1 (100% match), A1 with B2 (60% match), and so on. Then, this process is repeated for the next sequence from the reference database, for instance, A2 with B1 (20% match), A2 with B2 (25% match), etc., ( ). As a result, for each of the 5000 sequences in the reference database, up to 3000 comparisons with the target database are generated. Therefore, the final number of results in the tsv file will be equal to the product of the number of sequences in the reference database and in the target database (in this case, up to 15 million possible comparisons, although in practice, only the significant results are often recorded). Thus, the number of obtained comparisons significantly exceeds the sum of the number of sequences in both databases, as each match between sequences is treated as a separate result ( ). The discrepancy between the theoretical and obtained sequence values arises from the use of a reduced alphabet in diamond , which enhances processing speed but leads to the exclusion of sequences containing non-standard symbols, such as atypical amino acids or errors in the source data. Based on the conducted research, it is possible to define and standardize indices that describe the repeatability of peptide sequences in databases. The introduced metrics, namely the Database Absolute-Identity Repeatability Index diamond (DAIRI d ) and the Inter-Database Absolute-Identity Repeatability Index diamond (IDAIRI d ), allow for assessing sequence redundancy both within a single database and across multiple databases. The compared value (redundancy) refers to the number of repeating peptide sequences whose similarity reaches 100%, reflecting the percentage of identical peptides within a given database or between databases. These indices were specifically developed to quantify this percentage, based on data obtained through comparisons using the diamond tool. The DAIRI d and IDAIRI d metrics provide a standardized approach in order to evaluate how frequently specific peptide sequences reoccur within a given database or across multiple databases. The results of DAIRI d and IDAIRI d relative to the reference database (dbAMP) are presented in . Based on the obtained data, it is possible to observe the extent to which the analyzed databases are similar to themselves (self-comparison) as well as to the reference database. The analysis of the obtained indices from self-comparison and inter-database comparisons with the reference database was extended by an additional assessment of the similarity distribution in relation to peptide sequence length. The graphical representation of the analysis was performed using the RStudio environment, which facilitated a more precise presentation and interpretation of the data. The results of the self-comparison (DAIRI d ) of the analyzed databases reveal the proportion of peptide sequences exhibiting 100% similarity ( ). The highest percentage of identical sequences is found in the BaAMPs (82.66%) and CancerPPD (70.23%) databases, while the lowest is observed in CyBase (8.99%). The DAIRI d value for the reference database is 17.21%, indicating the relatively low similarity of the database to itself, as this value is less than half of the obtained results and nearly twice as low as the average DAIRI d , which is 33.85%. The database Absolute-Identity Repeatability Index diamond results are divided into three categories: Databases with similar sequence counts relative to the reference database dbAMP ( ). The degree of absolute similarity between the dbAMP and CAMP databases, as determined by IDAIRI d at 27.80%, suggests that more than a quarter of all peptides in both databases share identical sequences. This is represented on the graph by a signal bar at 100% similarity, with minimal scatter for other signals, indicating low peptide diversity in the reference database. A similar result is observed for the DBAASP database, which contains the highest number of identical peptides when compared to the other databases, despite a greater range of values in the dataset. The signal bar at 100% similarity is particularly intense for peptides of up to 25 amino acids in length, showing that this database contains a large number of peptides identical to those in dbAMP, although the length distribution is more varied. In comparison, the DRAMP database shows only 17.12% of identical sequences relative to dbAMP, the lowest value in this comparison, despite a similarly intense signal bar at 100% similarity. The scatter of similarity is small for peptides shorter than 100 amino acids, but increases for longer peptides, suggesting a decrease in the number of similar peptides as the sequence length increases. In the SATPdb database, identical peptides are found only for sequences of up to 50 amino acids in length, with similarity sharply decreasing for longer peptides. This suggests the greater diversity of longer sequences in SATPdb compared to dbAMP as peptide length increases. Databases with low sequence similarity relative to the reference database dbAMP ( )–lower IDAIRI d values arise from the normalization of this index. Lower values are characteristic of databases where the ratio of the number of peptides in the compared database to the reference database is low. In the APD database, which is the least diverse in this comparison relative to dbAMP, 100% similarity is evenly distributed across peptides of up to 50 amino acids in length, with similarity increasing for longer sequences. The CancerPPD database exhibits significant diversity, as similarity is observed only for peptides of up to 50 amino acids long and for a single group of peptides of around 70 amino acids in length. This suggests relatively high sequence diversity in this database compared to dbAMP. On the other hand, the CyBase database, which contains cyclic peptides, shows the lowest percentage of identical sequences relative to dbAMP (0.28%). Despite this, the scatter of data around 100 amino acids indicates low sequence diversity, and the absence of 100% similarity signals for certain peptide lengths suggests the uniqueness of this database. Databases with the least similar sequence counts compared to the reference database ( )–IDAIRI d decrease to values below 1%. The BaAMPs database contains the fewest identical sequences when compared to dbAMP. Despite its low IDAIRI d value, the scatter of signals in this database closely matches the signal values of the reference database, indicating low peptide diversity within this database. Similarly, the dadp database, although containing the highest number of identical sequences relative to dbAMP among the databases with the least similarity, shows minimal signal scatter, much like BaAMPs. The InverPep database, despite its low IDAIRI d value, exhibits a scatter of data similar to other databases in this group, suggesting that databases with fewer sequences may rely on more extensive counterparts. Lastly, the results for the ParaPep database align with previous analyses, indicating a similar lack of diversity when compared to the reference database. In order to demonstrate the applicability of the whole database as a query pattern for the identification of novel, putative AMPs, we have set up an simulation where a selected reference database (dbAMP) is used to assess the potential of a whole proteome. For that purpose, we have selected the proteome of Naja naja ( N. naja , indian cobra). It is one of very few complete venomous snake proteomes available in UniProt resources. Furthermore, it is also one of the limited number of species where the complete genome is also available, opening up possibilities for further experimentation and scientific work . The sequence similarity search of the dbAMP peptide collection against the N. naja proteome provided an extensive interaction network consisting of 12,815 records, with the sequence similarity ranging from 100.0 to 19.0%. The complexity of the obtained data file made it impossible to analyze the interaction in the raw format. To interpret the obtained data, they were passed into the Cytoscape software for visualization ( ) . For the proof of concept analysis, we have selected a relatively simple interaction network represented by graph in f. In this particular case, sequence similarity matching correlated two AMP peptides dbAMP_07416, 159aa ( https://awi.cuhk.edu.cn/dbAMP/information.php?db=dbAMP_07416 (accessed on 17 January 2025)) and dbAMP_10197, 135aa ( https://awi.cuhk.edu.cn/dbAMP/information.php?db=dbAMP_10197 (accessed on 17 January 2025)) with total of nineteen proteins from the N. naja proteome. The UniProt identifiers, as well as the respective percentage identity toward the AMPs and the beginning and end of the peptide matching sequence, are presented in . In order to further assess the potential of the identified targets A0A8C64F2, A0A8C6XI40, A0A8C6XPH8, and A0A8C6X4F2, we have analyzed structure corresponding to the related peptides. For all four targets, no experimental models were available in the PDB database, nor in AlphaFold respositories. We have built respective 3D models using the SWISS-MODEL webserver as an automated homology-modeling server. Models were built in the automated mode. Models with highest sequence coverage are analyzed in the manuscript. In each case, a high degree of structural similarity was observed . As an example, presents the combined results for the structural investigation of query peptides, as well as for the four cases indicated in in light gray. In the selected cases from the N. naja proteome, we wished to demonstrate the presence of the homologous domain in the larger structures. As indicated in , we have tested selected cases for putative cleavage sites that would allow for the release of the desired domain from the protein. In each tested case, such sites were found either directly flanking the AMP domain or within 10aa from the flanking amino acid. Additionally, for both the query and the putative subject sequences, we have generated electrostatic surface representations. 

The results of the self-comparison (DAIRI d ) reveal that the highest level of redundancy is observed in the BaAMPs (82.66%) and CancerPPD (70.23%) databases, suggesting a large number of repeating sequences. Conversely, the lowest DAIRI d values were recorded for the CyBase (8.99%) and APD (11.70%) databases, indicating a greater uniqueness of the sequences collected within them. The reference database dbAMP exhibits a moderate level of redundancy (17.21%), placing it below the average for all analyzed databases (33.85%). Inter-database analysis (IDAIRI d ) reveals that the databases most similar to the reference dbAMP are CAMP (27.80%) and DBAASP (27.92%), which may result from shared data sources or similar sequence selection criteria. Relatively high IDAIRI d values were also obtained for the SATPdb (20.05%) and DRAMP (17.12%) databases. The lowest similarity with dbAMP was recorded for CyBase (0.28%), which contained the relatively fewest sequences similar to the reference database, suggesting that it is the most unique compared to the reference database. In contrast, the DBAASP database exhibits the greatest similarity with dbAMP, making it the most identical of all the compared databases, which may indicate a high sequence homology between these databases. Analysis of the similarity distribution in relation to sequence length reveals that greater similarity is observed for shorter peptides, while longer sequences exhibit a greater spread of values. In many databases, such as SATPdb, DRAMP, and CancerPPD, identical sequences dominate among peptides shorter than 50 amino acids, with their numbers decreasing as the length increases. Conversely, some databases, such as CyBase, are characterized by greater sequence uniqueness across various length intervals, indicating their specific nature. The obtained results suggest that the developed indices may serve as valuable tools for assessing the quality and uniqueness of peptide databases. In the future, these indices may be used to optimize bioinformatics analysis and improve the classification and selection of sequences within peptide databases. In our proof-of-concept analysis, we have demonstrated that the approach similar to high-throughput processing is able to guide the potentially successful identification of new, putative AMP sequences. In our case, we were able to demonstrate that relatively low sequence similarity hits still represent viable targets. All the identified cases retained the key structural aspects of the query peptides (fold, secondary structure domain organization), while varied sequences influenced parameters, potentially improving AMP properties. 

4.1. Diamond— A Useful Tool in Database Analysis Diamond (v. 2.2.10.164) is an advanced, open-source software designed for matching DNA and protein sequences with large databases such as NCBI-nr and KEGG. It is characterized by exceptional performance—up to 20,000 times faster than BLASTX—allowing for the processing of vast datasets on standard servers, thus eliminating the need for supercomputers in metagenomic and evolutionary analyses . Although this tool requires significant memory and computational resources, for smaller datasets, tools such as BLASTX may be sufficient, making diamond particularly optimal for large-scale projects . One of the key factors accelerating diamond ’s performance is the use of a reduced protein alphabet. The traditional set of 20 amino acids has been replaced with an 11-letter set ([KREDQN], [C], [G], [H], [ILV], [M], [F], [Y], [W], [P], [STA]), enabling faster sequence matching with minimal loss of sensitivity . While this approach may affect precision in cases requiring highly specific matches, the differences are negligible in most applications. Additionally, the use of spaced seed technology allows for the analysis of selected positions in longer sequence fragments, enhancing sensitivity without significantly increasing analysis time. However, for more demanding tasks, such as detecting rare mutations, alternative tools may provide higher accuracy . Diamond also introduces double indexing, optimizing data locality and reducing the number of memory access operations. These innovations enable the processing of large datasets much faster than tools such as BLAST or MMSeqs2. For instance, in testing, diamond processed 246 million reads in just 2.3 hours, whereas BLASTX would require approximately 800,000 processor hours . Despite this, such high speeds demand substantial memory resources, which may pose a challenge on less powerful computational machines. While performance is one of diamond ’s key strengths, its application in biological analyses is equally important. The tool enables the fast and accurate comparison of protein sequences across various databases, which is crucial for identifying potential antimicrobial peptide (AMP) candidates. Diamond ’s ability to process large datasets in a short time is particularly valuable in metagenomic projects that require the analysis of hundreds of thousands or millions of sequences. Thus, diamond can be used to evaluate the proteomes of various species to identify potential AMP candidates, with results being integrated with structural analysis and antibacterial potential assessments . It is important to note that users must possess the necessary expertise to properly configure the tool and interpret the results, especially in more complex analyses requiring parameter adjustments. 4.2. Data Preprocessing and Sequence Compatibility for Diamond Tool Analysis After downloading sequences from the selected database, the analysis proceeded using the diamond tool. The first step was to preprocess the data to remove sequences that were not compatible with the diamond framework. This program utilizes a reduced alphabet , which significantly optimizes processing speed but leads to the rejection of sequences containing non-standard symbols, such as unusual amino acids or characters resulting from errors in the source data. This step is crucial, as it ensures that all analyzed sequences are comprehensible to the tool and can be further processed. Scripts used for data generation are included as a . 4.3. Efficient Sequence Comparison with Diamond: Converting Fasta to dmnd and Interpreting Results in tsv In the next step, the data stored in the fasta format were converted into the dmnd format. The dmnd format is a binary structure optimized for fast homology searching in large databases. Compared to traditional text formats such as fasta , dmnd is much more efficient in terms of data storage and processing. Its compactness reduces file sizes, which shortens read and write times. Furthermore, this format includes built-in indexes that enable rapid searching and sequence comparison, making it extremely useful in bioinformatics analyses involving large datasets. After converting the files into the dmnd format, the actual comparative analysis was performed. The diamond tool compared sequences from the created database with other datasets, generating results in tsv (Tab Separated Values) format. This format, due to its tabular structure, facilitates the later interpretation of the results. The tsv file contains key information about sequence matches, such as query sequence identifiers (qseqid) and database sequence identifiers (sseqid), the percentage of identical residues in the match (pident), match length (length), the E value indicating the statistical significance of the match (evalue), and the bit score describing the quality of the match (bitscore) ( ). 4.4. Calculation of DAIRI d and IDAIRI d Indices The Database Absolute-Identity Repeatability Index diamond (DAIRI d ) represents the percentage of comparisons performed using diamond in which 100% sequence identity was observed between databases. It is calculated according to the following formula (Equation (1)). (1) D A I R I d = n d i a m o n d c o m p a r i s o n s w i t h 100 % s i m i l a r i t y n t o t a l d i a m o n d c o m p a r i s o n s · 100 % where n diamond _ comparisons _ with _100%_ similarity represents the number of comparisons yielding 100% sequence identity and n total _ diamond _ comparisons denotes the total number of diamond -based comparisons performed. The resulting value indicates the proportion of comparisons within a given database that correspond to completely identical sequences. The Inter-Database Absolute-Identity Repeatability Index diamond (IDAIRI d ) extends the concept of DAIRI d by incorporating differences in the number of sequences between the compared and reference databases. It is defined as follows (Equation (2)). (2) I D A I R I d = n d i a m o n d c o m p a r i s o n s w i t h 100 % s i m i l a r i t y n t o t a l d i a m o n d c o m p a r i s o n s · n s e q u e n c e s i n t h e c o m p a r e d d a t a b a s e n s e q u e n c e s i n t h e r e f e r e n c e d a t a b a s e · 100 % The first component of this equation represents the proportion of 100% identical comparisons, while the second factor, the ratio of sequence counts in the compared database to those in the reference database, introduces a normalization element that accounts for differences in database sizes. Since databases vary in the number of sequences they contain, direct comparisons may lead to biased interpretations when one database is significantly larger or smaller than another. Even if a smaller database contains all the sequences present in a larger database, their absolute count will be lower, necessitating a normalization adjustment. The normalization factor ensures that the results are not artificially inflated or diminished due to database size disparities. When both databases contain an equal number of sequences, the normalization factor equals 1, making IDAIRI d equivalent to DAIRI d . If the compared database contains fewer sequences than the reference database, the normalization factor is less than 1, reducing the IDAIRI d value and preventing overestimation due to the larger reference dataset. Conversely, when the compared database is larger, the normalization factor exceeds 1, adjusting IDAIRI d to reflect the broader sequence repertoire in the compared dataset. 4.5. Graphical Analysis in Rstudio The data obtained from sequence comparisons between databases, processed using the tool and stored in tsv files, were subjected to graphical analysis using the Rstudio (v. 2024.12.1) environment. Three types of plots were prepared. The first type consisted of self-comparison scatter plots, which depicted the distribution of peptide similarity ( Y -axis) in relation to peptide sequence length ( X -axis). The second type of plots were scatter plots of similarity, identical to the previous type but differing in the data they represented, as they included information on comparisons between the reference database and the other databases. The final type of plots were histograms, illustrating the number of peptides of a given length with a specific degree of similarity. 4.6. Naja naja Proteome as an AMP Potential Source Analysis The proteome of Naja naja ( N. naja , Indian cobra) was downloaded directly from the UniProt database as a plain text fasta file . The reference proteome UP000694559 contained 29, 718 entries. After transformation into a diamomd database file, it was used during the search, with the complete dbAMP database as a query. The resulting tsv file was then inspected and adjusted for import into Cytoscape (v. 3.10.3), a network analysis and visualization software . In short, the file was imported into a spreadsheet software and modified. Suitable column labels based on the diamond nomenclature were introduced. Additionally, helper spreadsheets were prepared based on the data obtained from UniProt. In particular, information allowing for the identification of toxins in the proteome were extracted and formatted according to the Cytoscape requirements. The content of the prepared spreadsheets was imported into Cytoscape in the following sequence: interaction network (calculated similarity between N. naja proteins and dbAMP peptides), source and target identifiers allowing for visual differentiation between network elements, identifiers related to the protein type allowing for toxins, and non-toxin diffentiation. 4.7. Naja naja Proteome Hits—Structure Analysis and Visualization Structural models of the N. naja targets selected for proof of concept analysis were built using the SWISS-MODEL, a fully automated protein structure homology-modeling server in the automated mode . The predicted structure of the AMP peptides was obtained from the dbAMP database webpage and used without modification. Electrostatic potential representation was generated using the PyMOL APBS Electrostatics plug-in. The obtained models were visualized in the PyMOL software (opensource v. 3.0.0) . 

Diamond (v. 2.2.10.164) is an advanced, open-source software designed for matching DNA and protein sequences with large databases such as NCBI-nr and KEGG. It is characterized by exceptional performance—up to 20,000 times faster than BLASTX—allowing for the processing of vast datasets on standard servers, thus eliminating the need for supercomputers in metagenomic and evolutionary analyses . Although this tool requires significant memory and computational resources, for smaller datasets, tools such as BLASTX may be sufficient, making diamond particularly optimal for large-scale projects . One of the key factors accelerating diamond ’s performance is the use of a reduced protein alphabet. The traditional set of 20 amino acids has been replaced with an 11-letter set ([KREDQN], [C], [G], [H], [ILV], [M], [F], [Y], [W], [P], [STA]), enabling faster sequence matching with minimal loss of sensitivity . While this approach may affect precision in cases requiring highly specific matches, the differences are negligible in most applications. Additionally, the use of spaced seed technology allows for the analysis of selected positions in longer sequence fragments, enhancing sensitivity without significantly increasing analysis time. However, for more demanding tasks, such as detecting rare mutations, alternative tools may provide higher accuracy . Diamond also introduces double indexing, optimizing data locality and reducing the number of memory access operations. These innovations enable the processing of large datasets much faster than tools such as BLAST or MMSeqs2. For instance, in testing, diamond processed 246 million reads in just 2.3 hours, whereas BLASTX would require approximately 800,000 processor hours . Despite this, such high speeds demand substantial memory resources, which may pose a challenge on less powerful computational machines. While performance is one of diamond ’s key strengths, its application in biological analyses is equally important. The tool enables the fast and accurate comparison of protein sequences across various databases, which is crucial for identifying potential antimicrobial peptide (AMP) candidates. Diamond ’s ability to process large datasets in a short time is particularly valuable in metagenomic projects that require the analysis of hundreds of thousands or millions of sequences. Thus, diamond can be used to evaluate the proteomes of various species to identify potential AMP candidates, with results being integrated with structural analysis and antibacterial potential assessments . It is important to note that users must possess the necessary expertise to properly configure the tool and interpret the results, especially in more complex analyses requiring parameter adjustments. 

After downloading sequences from the selected database, the analysis proceeded using the diamond tool. The first step was to preprocess the data to remove sequences that were not compatible with the diamond framework. This program utilizes a reduced alphabet , which significantly optimizes processing speed but leads to the rejection of sequences containing non-standard symbols, such as unusual amino acids or characters resulting from errors in the source data. This step is crucial, as it ensures that all analyzed sequences are comprehensible to the tool and can be further processed. Scripts used for data generation are included as a . 

In the next step, the data stored in the fasta format were converted into the dmnd format. The dmnd format is a binary structure optimized for fast homology searching in large databases. Compared to traditional text formats such as fasta , dmnd is much more efficient in terms of data storage and processing. Its compactness reduces file sizes, which shortens read and write times. Furthermore, this format includes built-in indexes that enable rapid searching and sequence comparison, making it extremely useful in bioinformatics analyses involving large datasets. After converting the files into the dmnd format, the actual comparative analysis was performed. The diamond tool compared sequences from the created database with other datasets, generating results in tsv (Tab Separated Values) format. This format, due to its tabular structure, facilitates the later interpretation of the results. The tsv file contains key information about sequence matches, such as query sequence identifiers (qseqid) and database sequence identifiers (sseqid), the percentage of identical residues in the match (pident), match length (length), the E value indicating the statistical significance of the match (evalue), and the bit score describing the quality of the match (bitscore) ( ). 

d and IDAIRI d Indices The Database Absolute-Identity Repeatability Index diamond (DAIRI d ) represents the percentage of comparisons performed using diamond in which 100% sequence identity was observed between databases. It is calculated according to the following formula (Equation (1)). (1) D A I R I d = n d i a m o n d c o m p a r i s o n s w i t h 100 % s i m i l a r i t y n t o t a l d i a m o n d c o m p a r i s o n s · 100 % where n diamond _ comparisons _ with _100%_ similarity represents the number of comparisons yielding 100% sequence identity and n total _ diamond _ comparisons denotes the total number of diamond -based comparisons performed. The resulting value indicates the proportion of comparisons within a given database that correspond to completely identical sequences. The Inter-Database Absolute-Identity Repeatability Index diamond (IDAIRI d ) extends the concept of DAIRI d by incorporating differences in the number of sequences between the compared and reference databases. It is defined as follows (Equation (2)). (2) I D A I R I d = n d i a m o n d c o m p a r i s o n s w i t h 100 % s i m i l a r i t y n t o t a l d i a m o n d c o m p a r i s o n s · n s e q u e n c e s i n t h e c o m p a r e d d a t a b a s e n s e q u e n c e s i n t h e r e f e r e n c e d a t a b a s e · 100 % The first component of this equation represents the proportion of 100% identical comparisons, while the second factor, the ratio of sequence counts in the compared database to those in the reference database, introduces a normalization element that accounts for differences in database sizes. Since databases vary in the number of sequences they contain, direct comparisons may lead to biased interpretations when one database is significantly larger or smaller than another. Even if a smaller database contains all the sequences present in a larger database, their absolute count will be lower, necessitating a normalization adjustment. The normalization factor ensures that the results are not artificially inflated or diminished due to database size disparities. When both databases contain an equal number of sequences, the normalization factor equals 1, making IDAIRI d equivalent to DAIRI d . If the compared database contains fewer sequences than the reference database, the normalization factor is less than 1, reducing the IDAIRI d value and preventing overestimation due to the larger reference dataset. Conversely, when the compared database is larger, the normalization factor exceeds 1, adjusting IDAIRI d to reflect the broader sequence repertoire in the compared dataset. 

The data obtained from sequence comparisons between databases, processed using the tool and stored in tsv files, were subjected to graphical analysis using the Rstudio (v. 2024.12.1) environment. Three types of plots were prepared. The first type consisted of self-comparison scatter plots, which depicted the distribution of peptide similarity ( Y -axis) in relation to peptide sequence length ( X -axis). The second type of plots were scatter plots of similarity, identical to the previous type but differing in the data they represented, as they included information on comparisons between the reference database and the other databases. The final type of plots were histograms, illustrating the number of peptides of a given length with a specific degree of similarity. 

The proteome of Naja naja ( N. naja , Indian cobra) was downloaded directly from the UniProt database as a plain text fasta file . The reference proteome UP000694559 contained 29, 718 entries. After transformation into a diamomd database file, it was used during the search, with the complete dbAMP database as a query. The resulting tsv file was then inspected and adjusted for import into Cytoscape (v. 3.10.3), a network analysis and visualization software . In short, the file was imported into a spreadsheet software and modified. Suitable column labels based on the diamond nomenclature were introduced. Additionally, helper spreadsheets were prepared based on the data obtained from UniProt. In particular, information allowing for the identification of toxins in the proteome were extracted and formatted according to the Cytoscape requirements. The content of the prepared spreadsheets was imported into Cytoscape in the following sequence: interaction network (calculated similarity between N. naja proteins and dbAMP peptides), source and target identifiers allowing for visual differentiation between network elements, identifiers related to the protein type allowing for toxins, and non-toxin diffentiation. 

Structural models of the N. naja targets selected for proof of concept analysis were built using the SWISS-MODEL, a fully automated protein structure homology-modeling server in the automated mode . The predicted structure of the AMP peptides was obtained from the dbAMP database webpage and used without modification. Electrostatic potential representation was generated using the PyMOL APBS Electrostatics plug-in. The obtained models were visualized in the PyMOL software (opensource v. 3.0.0) . 

