Preprint
Review

Small Natural Cyclic Peptides from DBAASP Database

Altmetrics

Downloads

155

Views

45

Comments

0

A peer-reviewed article of this preprint also exists.

Submitted:

22 May 2024

Posted:

23 May 2024

You are already at the latest version

Alerts
Abstract
Antimicrobial peptides (AMPs) are promising tools for combating microbial resistance. However, their therapeutic potential is hindered by two intrinsic drawbacks – low target affinity and poor in vivo stability. Macrocyclization, a process that improves the pharmacological properties and bioactivity of peptides, can address these limitations. As a result, macrocyclic peptides represent attractive drug candidates. Moreover, many drugs are macrocycles that originated from natural product scaffolds, suggesting that nature offers solutions to the challenges faced by AMPs. In this review, we explore natural cyclic peptides from the DBAASP database. DBAASP is a comprehensive repository of data on antimicrobial/cytotoxic activities and structures of peptides. We analyze the data on small (≤25 AA) ribosomal and non-ribosomal cyclic peptides from DBAASP according to their amino acid composition, bonds used for cyclization, targets they act on, and mechanisms of action. This analysis will enhance our understanding of the small cyclic peptides that nature has provided to defend living organisms.
Keywords: 
Subject: Medicine and Pharmacology  -   Medicine and Pharmacology

1. Introduction

Numerous endeavours are underway to explore natural peptides as prototypes for novel drugs. For instance, antimicrobial peptides (AMPs) are considered hopeful tools to solve the problem of resistance [1,2]. Along with antimicrobial properties, AMPs possess anticancer activities [3], cell-penetrating abilities [4], immune system modulating capabilities [5], and other biologically significant properties [6], majorly attributed to their common feature of membrane targeting. Although, in some cases, their activities rely upon the interaction with membrane or cytoplasmic proteins. Peptide-based drugs predominantly work by inhibiting interactions between large surfaces, such as protein-protein, protein-nucleic acids, protein-lipid and lipid-lipid interactions. Simultaneously, there are small molecular weight peptide-based drugs that inhibit interactions with specific enzymes or receptors. It is worth noting, that the current design process for novel peptide-based drugs remains challenging. The main hurdle lies in creating druggable peptide designs [7].
The druggability of a compound refers to its ability to meet specific criteria related to absorption, distribution, metabolism, excretion, and toxicity. Absorption and distribution determine bioavailability, which translates to the capability of a compound to pass the membrane barrier and resist chemical destruction, e.g., be metabolically and chemically stable. Drugs need to be formulated to resist rapid chemical destruction in living organisms and to be easily absorbed and distributed. For traditional small molecular weight drugs, druggability is determined by the Rule of 5 (Ro5) [8], which defines limits for molecular size (MW<500 KDa) and polarity (<5 hydrogen bond donors, <10 hydrogen bond acceptors, cLogP lipophilicity <5) to achieve cell permeability and oral bioavailability. These limitations in drug discovery have constrained the range of targets that can be effectively inhibited by small molecules [9], which means that targets are mostly limited to the active sites of enzymes and receptors. Consequently, the small size of traditional drugs makes it challenging for them to effectively inhibit large surface interactions. This limitation has led to a broad category of target sites referred to as ‘undruggable’ binding sites [10,11,12]. Notably, half of the beyond-Rule-of-5 (bRo5) drugs that target flat or groove-like interfaces are macrocycles [11,12], many of which originate from natural product scaffolds.
The unique physicochemical properties of peptide-based drugs, including their larger size and more flexible backbone, enable them to act as potent inhibitors of ‘undruggable’ binding sites. However, it is worth noting that peptides may face two inherent limitations: membrane impermeability and poor in vivo stability. Macrocyclization improves the pharmacological properties and bioactivity of peptides [13]. Consequently, macrocyclic peptides emerge as promising candidates for drug development, deserving significant attention.
Macrocyclic peptides constitute a highly diverse family of peptides found in various organisms, including bacteria, fungi, protists, etc. Recent studies on natural and synthetic model systems have demonstrated the increasing utilization of cyclic peptide natural products in designing bRo5 molecules, which exhibit exceptional cell permeability and oral absorption [14]. Macrocyclization is employed to enhance metabolic stability. To achieve cell permeability, macrocycles must possess a ‘chameleonic’ ability, allowing them to change their conformation to expose polar groups in an aqueous solution while concealing them during traversal of lipid membranes [15].
Small macrocyclic peptides represent more valuable drug candidates. The Database of Antimicrobial Activity and Structure of Peptides (DBAASP) is a comprehensive repository for small cyclic peptides [16], freely accessible at www.dbaasp.org. Currently, the database contains information on 597 non-ribosomal and 1413 ribosomal cyclic peptides (data was obtained on 13/05/24). The length of non-ribosomal cyclic peptides does not exceed 25 amino acids (AA), making them small peptides. The number of ribosomal cyclic peptides of length less than 26 AA equals 484. In this review, we will explore small cyclic peptides of DBAASP that have been offered by nature to defend living organisms. Volume of DBAASP allows to uncover common peculiarities of the cyclic structures formed in natural peptides.

2. Definition of intrachain bonds and rings of natural cyclic peptides

Various kinds of intrachain bonds, such as amide (AMD), ester (EST), thioether (TIE), amine (AMN), etc., are used to cyclize natural peptides (see Table 1). As a result, many kinds of rings are created.
In general, a certain kind of ring, for example, lactam (LAC), can be formed by different types of bonds and three kinds of chain linking: sidechain-to-sidechain (SSB), sidechain-to-mainchain (SMB), and mainchain-to-mainchain (MMB). Consequently, in this work, each ring formed by a specific intra-chain bond has been characterized by the following three parameters: the type of bond (AMD, for example), the chains participating in the bonding (SMB, for example) and the type of ring formed (LAC, for example). In some cases, special names are employed for widespread types of rings to emphasize the differences between rings formed by similar bonds. For instance, a peptide is a polymer of amino acids with a main chain composed of a string of amide bonds. Consequently, the majority of rings formed by intrachain bonds constitute macrolactams. Amide bonds formed between the N-terminal amine and C-terminal carboxyl groups create a macrolactam ring and can be defined by the following three parameters (AMD, MMB, LAC). An analogous macrolactam ring formed by bonding between sidechain amine (for example, of Lys) and C-terminal carboxyl, is presented in the three-parameter definition as (AMD, SMB, LAC). To emphasize that the first macrolactam corresponds to the case of head-to-tail linking, the ring has been renamed as NCB. Therefore, a ring of head-to-tail cyclized peptides in the three-parameter definition is presented as (AMD, MMB, NCB) instead of the above-mentioned (AMD, MMB, LAC). Another example involves lactams formed in cyclic dipeptides, where the ring name diketopiperazine (DKP) is used instead of lactam.
Thus, rings closed by different kinds of bonds between pairs of residues spaced by one or more amino acids can be defined as lactams (macrolactams, macrocyclic amides). In this work, lactams have been classified based on bonds that close them (Table 1). Above, examples of rings closed by amide bonds, such as NCB and DKP, have been presented. As an additional example, the rings closed by ester bonds between carboxyl groups of side chains (for example, of Asp or Glu) and C-termini of peptides can be considered. This type of ring can be simultaneously classified as both, a lactam and a lactone (cyclic ester). However, for more specificity, it has been defined as LCN (lactone). All bonds used for ring closing in small natural peptides, along with the names of corresponding rings as used in DBAASP [16] and in this work, are presented in Table 1.

3. Small Cyclic Peptides of DBAASP

DBAASP is a manually-curated database that offers detailed information on the physicochemical and structural features, as well as antimicrobial and cytotoxic activities of peptides assessed in vitro. The current version of the database contains data on more than 22000 peptides, the antimicrobial/cytotoxic potency of which has been assessed in vitro [16]. DBAASP provides information on peptides belonging to both Ro5 and bRo5 groups. In terms of size, peptides can be classified into two types: those with a length of ≤5 AA, referred to as ultra-short peptides (USPs), and those with a length of 6-25 AA, referred to as short peptides (SPs).

3.1. Ultra Short Cyclic Peptides

The number of natural USPs (peptides with a length of ≤5 AA) in DBAASP is 275, comprising 7 ribosomal and 268 non-ribosomal peptides. Among these, there are 232 ultra-short cyclic peptides (USCPs), including 230 non-ribosomal and 2 ribosomal peptides. The numbers of USCPs of different lengths and intrachain bonds have been presented in Table 2.
Amide bond. According to Table 2, rings of non-ribosomal USCPs are mainly (159 peptide) closed by amide bonds. Among these, 147 are closed by mainchain-to-mainchain amide bonds, including 65 NCB and 82 DKP. In 12 peptides containing a LAC ring, the rings are closed by amide bonds between the side chain of Ornithine or Trp and C-terminal carboxyl group.
Amine and disulfide bonds. A lactam ring can also be formed by amine bonds between the amine group of main chain and the β-carbon of side chain of Cys, as observed in the case of peptides containing a β-lactam ring (see Table 2).
Amine bonds also participate in the formation of five-membered saturated Pyrrolidine ring in the exotic peptide, Gliotoxin (DBAASPN_21167). Gliotoxin is a unique non-ribosomal dipeptide in DBAASP, containing three different rings, including a disulfide bond. In this peptide, the disulfide bond links two α-carbon atoms (creating DSL ring) and so it differs from other disulfide-bonded peptides in DBAASP, which are created through disulfide bonding between Cys residues (CST ring). It is worth noting that Gliotoxin is also cyclized by a head-to-tail amide bond, and therefore, constituting of a diketopiperazine ring.
CST rings, created through disulfide bonding between Cys residues, appear only in three non-ribosomal peptides: PM181110 (DBAASPN_21712), Malformin E (DBAASPN_21150), and Malformin A1 (DBAASPN_21151).
Ester bond. In 65 USCP, comprising only non-ribosomal peptides, LCN rings are closed by ester bonds (EST) through the binding between the side chain of Thr or N-terminal hydroxycarboxylic acids, (such as 5-hydroxy-2,4-dimethyl-3-oxooctanoic acid, 2-hydroxy-3-methylbutyric acid, 3-hydroxy-2-methyloctynoic acid, 3-hydroxy-2,4-dimethyldecanoic acid, 3-hydroxy-2,4,6-trimethyldecanoic acid, 2-hydroxy-4-methylpentanoic acid, 2,6,10,14-methyl-3-hydroxy-elcosadienoic acid, etc.) with the carboxyl group of the C-terminal amino acid.
Thioether (Sulfide) bond. Thioether bonds (TIE) between the side chains of neighboring Cys and Val residues close five or six-membered rings of varying saturation. Saturated five-membered rings are named thiazolinidine (THZD), and occurred in β-lactam antibiotics of the penicillin or methicillin family (DBAASPN_20466, DBAASPN_20468, DBAASPN_20907). In Cephalosporin C (DBAASPN_20467), the cyclized structure is presented by unsaturated six-membered Thiazin (TZ) ring.
Thus, non-ribosomal USCPs are primarily cyclized by amide and ester bonds, while other types (DSB, TIE, and AMN) are used rarely. Among USCPs, only 2 are ribosomal and both use AMD bonds for cyclization (refer to Table 2).
Natural USCPs rarely use disulfide bonds (DSB) for cyclization. Cys is a scarce amino acid in non-ribosomal linear USPs, as depicted in Figure 1A. Non-ribosomal peptides are rich in unusual amino acids (Figure 1). In contrast to ultra-short linear peptides, USCPs are rich in Trp, Phe, and Pro (Figure 1B,C). Pro is especially abundant in the cyclic dipeptides and tripeptides. This fact can be explained by the necessity of Pro to cyclize ultra-short peptides, as it has the propensity to promote sharp turns in the main chain of the polypeptide. For example, it has been shown that cyclic tripeptides contain at least one cis peptide bond, while tetrapeptides are the smallest cyclic peptides that can adopt an all-trans conformation, although they generally adopt alternating cis/trans conformation [17]. Consequently, with the lengthening of peptide chains, the sharp need for Pro in cyclization diminishes.
The results of the analysis indicate that natural USCPs are mainly synthesized non-ribosomally, predominantly by fungi (131 peptides) and bacteria (102 peptides), only 7 peptides originate from archaea, 1 from Animalia, and 2 from plantae. They act on a variety of cells, organisms or viral particles (cancer cells, fungi, gram-positive and gram-negative bacteria, viruses, parasites). Notably, their target objects majorly are proteins (membrane or cytoplasmic, see Figure 2). This observation is explainable. USCPs have the ability to satisfy the Rule of 5 and so have druggable targets, enzymes or receptors. For example, beta-lactam antibiotics targeted the enzyme responsible for the cell wall biosynthesis. They have a common core skeleton composed of an enclosed dipeptide formed by the condensation of L-cysteine and D-valine [18]. While for the head-to-tail cyclized dipeptides, core skeleton of which consist of 2,5-diketopiperazine ring, is suggested that they inhibit regulators of bacterial cell population density (‘quorum sensing’) [19,20]. For some member of cyclic dipeptides (for example, cyclic (L-Pro – L-Leu; L-Ser – L-Trp; L-Tyr – L-Pro), quorum sensing activities has been studied [19]. For other 2,5-diketopiperazines, constituted by hydroxyprolines, an ability to stimulate cytokines or chemokines in macrophages have been supposed [5]. We have to note, that, as a rule, USCPs are not membrane-active peptides (see Figure 2).

3.2. Short cyclic peptides

In contrast to USCPs, a membrane (lipid bilayer) is a main target for short cyclic peptides (SCPs) (see Figure 2). Therefore, they are mainly membrane-active peptides (MAP). It is worth noting that, in addition to the membrane, there are other ‘undruggable’ binding sites targeted by SCPs (Figure 2). The MWs of short peptides (SPs) exceed 600 Da and they are categorized as beyond Rule of 5 compounds. They can interact with ‘undruggable’ binding sites and would be considered as a basis for the development of beyond-Rule-of-5 (bRo5) drugs [11]. SPs can be considered druggable if they are capable to change the surface accessibility of polar groups during the transition from water to the membrane [15]. An additional requirement is resistance to chemical destruction, which can be achieved through cyclization. Consequently, it is interesting to explore the types of natural SCPs that occur.
DBAASP stores data on 1565 ribosomal and 464 non-ribosomal SPs, which are peptides with length in the interval of 6-25 AA. Among these, 482 ribosomal and 367 non-ribosomal peptides are cyclized by intrachain covalent bonds.

3.2.1. Peptides could interact with several ‘undruggable’ binding sites

‘Undruggable’ binding sites can include surfaces of the proteins, DNA, RNA, lipids, etc. Often, undruggable binding sites can consist of surfaces of two or more molecules. For example, in the case of membrane proteins, the binding site can constitute surfaces of the protein and lipids simultaneously. Darobactin (DBAASPR_17389) exclusively targets the bacterial insertase BamA6 – the central unit of the BAM [21], while macrocyclic beta-hairpin-shaped peptide JB-95 interacts with both outer membrane and BamA [22].
In general, if a molecule is a membrane protein, instances where there is a distinct binding site on the protein surface, fully separated from the lipid environment, are rare. Moreover, if the targets of a peptide are inside of the cell (such as cytoplasmic proteins, DNA, RNA, ribosome, mitochondrial membrane, etc.), it has to translocate through the lipid bilayer in an energetically dependent or independent manner and so has to interact with the cytoplasmic membrane. Consequently, it is clear, why the majority of AMPs are considered membrane-active peptides (MAP).

3.2.2. Non-ribosomal short cyclic peptides

Various kinds of intra-chain bonds, such as amide, ester, thioether, amine, etc., are used to cyclize non-ribosomal SCPs (see Table 3). As a result, many kinds of rings are created. The number and size of rings in peptides also vary. The majority of non-ribosomal cyclic peptides consist of non-proteinogenic amino acids and have modifications on termini. Only 24 non-ribosomal SCPs of DBAASP consist exclusively of proteinogenic amino acids and do not have modification on their C- or N-termini.
As mentioned above, DBAASP stores data on 465 non-ribosomal SPs, including 368 cyclic peptides.

3.2.2.1. Intra-chain bonds and created rings in non-ribosomal peptides

Amide bonds. Among non-ribosomal SCPs, more than half (204 peptides) have been cyclized by amide bonds (Table 3).
Rings created by amide bonds can be divided into two groups. The number of peptides with rings created by linking terminal amines of the main chain to the terminal carboxyl of the main chain, that is NCB, equals 138. In another peptide group (66 peptides), rings are closed by amide bonds between side chain amine groups and C-terminal carboxyl groups of the main chain, forming LAC rings (Table 3). For instance, amide bonds between amine group of the side chains of 2,4-diaminobutyrate (DAB), Lys, or 2,3-diaminopropionic acid (DAP), or N-terminal 3-aminotetradecanoic acid and main chain carboxyl groups are cyclized into LAC rings.
Ester bond. Apart from amide bonds, other side chain-to-main chain bonds are also used to cyclize non-ribosomal SCPs. For example, ester bonds (ESTs) between the side chains of Ser, Homo Ser, dimethyl-serine, Trh, Allo-Thr, N-methyl-Trh, Tyr, N-terminal hydroxy-fatty acids and C-terminal, carboxyl groups form rings in 142 SCPs defined as lactone (LCN).
Thioether (Sulfide), amine and carbon bonds. Non-ribosomal SCPs closing their rings by amine, thioether and carbon bonds are also presented in the DBAASP. However, in limited numbers. Their scarcity can be explained by many various factors, including a lack of sufficient data in DBAASP at the current moment. Therefore, it would be difficult to make some generalizations about their exoticism. We will only have to describe them.
Amino acid Cys can participate in the creation of intrachain thioether (sulfide) bonds, resulting in the formation of various kinds of rings. Thioether bonds between the side chain of Cys with the main chain carbonyl group of neighbouring amino acids form an unsaturated Thiazoline (THZN) ring in the SCPs of the Bacitracin family (DBAASPN_17106-DBAASPN_17109). While, in the SCPs of Pagoamide (DBAASPN_18305), Haligramide (DBAASPN_21016, DBAASPN_21017), and Bisebromoamide (DBAASPN_18509) families, bonding between similar groups of atoms cyclize chain into more aromatic (unsaturated) Thiazole (THZ) rings. The same kinds of bonds in Lugdunin (DBAASPN_13304) create a saturated Thiazolidine ring.
LAC ring can also be formed through amine bonds between the side chain carbons of 2,3-diaminoacrylic acids and main chain amine groups, as observed in peptides of the Callyaerin (DBAASPN_9851-9856) family. Amine bonds between the side chain carbons of Leu and main chain amine groups in Ilamycin C (DBAASPN_14863), Ilamycin E (DBAASPN_14866), and Rufomycins (DBAASPN_ 22199-22200) produce six-membered delta lactam (δLAC) rings.
In Vancomycin (DBAASPN_21634), ether bonds between the side chains of beta-oxy tyrosines and Phenylglycine form a macrolactam (LAC) ring.
In a few DBAASP peptides (DBAASPN_20988, DBAASPN_20991-20999, DBAASPN_21001, DBAASPN_21634), side chains are stapled by carbon bond, using the two Aryl groups creating a Biaryl ring (BAR).
Thus, like USCPs, non-ribosomal SCPs are mainly cyclized by amide bonds (about 55% of non-ribosomal SCPs). ESTs are more frequent (about 40% of non-ribosomal SCPs) among other bonds used for cyclization (see Table 3).

3.2.3. Ribosomal short cyclic peptides

Opposed to non-ribosomal SCPs, ribosomal SCPs mainly use Cys residue for cyclization. To create rings, Cys participates in disulfide and sulfide (thioether) bond formation. In over 80% of ribosomal SCPs, rings are formed through the creation of cystines (CST), aminovinyl cysteines (AVC), and lanthionines (LAN). Such rings are not formed in non-ribosomal SCPs (see Table 3).
Peptides with rings, created by amide bonds represent only about 20% of ribosomal SCPs of DBAASP, whereas over half of non-ribosomal SCPs are cyclized through amide bonds. Also, as observed in Table 3, lactone rings (LCNs) are not represented in ribosomal SCPs, while they constitute about 40% of non-ribosomal SCPS.
One additional difference, concerning five-membered rings, is of note. In ribosomal SCPs, five-membered rings (such as OXZ and OXZN) are created by ether bonds, while in non-ribosomal SCPs, five-membered rings are mainly closed by thioether bonds, forming rings such as THZ, THZN, and TZD. Only THZ is a shared five-membered ring, appearing in both, non-ribosomal and ribosomal SCPs.
3.2.3.1. Intrachain bonds and created rings in ribosomal peptides
Disulfide bond. In the majority of ribosomal SCPs, rings are formed through the disulfide bonds. Figure 4 represents the distribution of lengths for ribosomal disulfide-bonded peptides in DBAASP, ranging up to 50 AA. The total number of such peptides is 1100, which constitutes the majority of all cyclic ribosomal peptides (1269) of that length.
The lengths of approximately 360 ribosomal disulfide-bonded peptides fall within the interval of 6-25 AA in length. That is, they are short. In this section, we will focus on disulfide-bonded short ribosomal cyclic peptides (RSCPs).
The amino acid composition of the RSCPs from DBAASP, compared to the compositions of UniProt [23] sequences (Figure 5), reveals that RSCPs are more hydrophobic (with a higher abundance of phenylalanine and isoleucine) and more basic (due to the higher abundance of lysine) compared to the ‘average protein’. RSCPs also contain more cysteines and glycines than the ‘average protein’. The presence of glycine is likely related to the requirements for flexibility to close the ring and, at the same time, glycine can promote ‘chameleonic’ ability of peptide, while the presence of cysteine is crucial for cyclization, which is important for the activity and chemical stability.
The distribution of cysteines along the peptide chain is also an interesting aspect to investigate, in order to assess the size of cycles that can be closed. The analysis of the distribution of i-spaced amino acid pairs (DiSAAP), specifically focusing on cysteine pairs, in RSCPs of DBAASP reveals a statistically significant abundance of pairs with 5 and 8 AA spaces, while pairs with 0 and 1 spaces are scarce compared to random distribution (Figure 6). These findings can be explained by the role of Cys residues. Cys can be the source of intrachain bonds and through post-translational modification, it can stabilize the certain peptide structures. Taking into account the geometry of the side chain, it is problematic to form disulfide bond by pairs of Cys residues spaced by 0 or 1 residues because of steric hindrances [24]. Therefore, the appearance of neighbouring Cys pairs or those separated by one residue is limited. On the other hand, the abundance of Cys pairs with 5 and 8 AA spaces, suggests the presence of cyclized loops consisting of 7 and 10 AA, characterizing many disulfide-bonded RSCPs in DBAASP. For instance, the abundance of Cys pairs spaced by 5 AA points to the prevalence of ‘Rana box’-containing peptides (RBPs) among RSCPs of DBAASP. RBPs, such as Brevinins (DBAASPR_1009), Nigrocins (DBAASPR_1542), etc., contain small C-terminal loops closed into rings by disulfide bridges and N terminal linear tails containing more than 7 AA [25]. Consequently, the overall topology of chains of such peptides resemble a lasso. In the DBAASP, peptides with bigger disulfide-bonded loops that differ in amino acid composition from RBPs [26] or are situated at the N-terminus of peptides [27] can be found.
Moreover, when the linear tail is absent, the peptide chain topology resembles a hairpin. Small size, rigid structure, stability to proteases, and wide spectrum of biological functions make hairpin AMPs an attractive molecular basis for drug design [28]. Many disulfide-bonded RSCPs from DBAASP, such as Protegrins (DBAASPR_672), Tachyplesins (DBAASPR_2229), Arenicins (DBAASPR_2025), Θ-defensins (DBAASPR_823), etc. have hairpin topology of chains.
Figure 13.
Thus, among disulfide-bonded ribosomal peptides of DBAASP, about 290 contain a single disulfide bond, allowing us to envision their structure. These peptides can adopt either a ‘lasso-like’ (LL) structure (Figure 7A), with a C- or N- terminal cyclized loop, or a ‘hairpin-like’ (HL) structure, formed when the Cysteines are located near the peptide termini (Figure 7B).
To define the type of structure, the parameter s and p have been introduced. For the given length of peptide L, positions n1 and n2 (n1 < n2) of first (C1) and second (C2) cysteines, respectively, the value of s and p are calculated as:
s = ( L n 2 + n 1 + 1 ) / L
p = ( L n 2 ) / ( n 1 1 )
Relying on the values of s and p, the correspondence of the peptides to LL or HL structure has been assessed. If s≥0.5 and 0.3≥ p≥3 or n1=1, then the peptide structure is defined as LL, and otherwise (when s<0.5), the peptide is defined as HL. In the case of LL peptides, we distinguish N-LL (when L–n2>n1) and C-LL (when L–n2<n1). Among ribosomal SCPs that are cyclized by one disulfide-bond, the majority have been defined as C-LL, and only one as N-LL. The number of HL peptides is about 60 (see Table 4). The size of loops of LL peptides are about 7-8 AA (including Cys-s), while HL loops are a little larger, about 10-11 AA. Such macrocyclic peptides have some conformational freedom and, consequently, flexibility.
The structures of about 70 peptides among disulfide-bonded ribosomal SCPs are stabilized by two or more disulfide bonds. Various combinations of pairing can be considered for a given number n=2k (k=1, 2, …, m) of cysteines (C) that participate in the formation of k disulfide bonds. The number (N) of variants of pairing to form k disulfide bonds is assessed by the formula:
N = ( 2 k ) ! / ( 2 k   k ! ) ,
while the number of structures created as a result, can be classified into three groups (Figure 8). If cysteines of the peptide with k disulfide bonds are enumerated according to the increasing positions along the chain as C1, C2, C3, …, Cn, then pairing schemes that define three different groups of structures can be presented. For example, the pairing of cysteines according to the scheme C1–Cn, C2–Cn-1, …, Cn/2–Cn/2+1 (Figure 8A) gives the formation of ‘ladder-like’ multicyclic structures (LD). Another scheme C1–C2, C3–C4, …, Cn-1–Cn (Figure 8B) is the formula to form structures consisting of ‘strings of rings’ (SR) closed by disulfide bonds. All other variants of pairing cause the crossing of rings formed by disulfide bonds (Figure 8C), and thus, the formation of structures with ‘crossed-rings’ (CR).
Consequently, disulfide bonds can create three kinds of multicyclic structures: LD, SR, and CR. For k=2, each kind of structure is presented with a single variant of pairing. With the increasing of k, the number of variants of pairing causing CR is increased, while variants that correspond to LD and SR are always equal to one. Therefore, if we consider the process of pairing of cysteines as random, we can suppose that, in disulfide-bonded peptides, CR structures prevail. Indeed, in the majority of peptides in DBAASP with 3 and more disulfide bonds, pairing of cysteines results in crossing of cycles, which means CR structures are formed, while other types of structures appear rarely (see Table 4). Although, opposed to this, in the short peptides with 2 disulfide bonds, the LD structure, having a hairpin topology, markedly prevails over the other two types of structure, and thus, this fact cannot be considered as a result of a random process. We can also note non-randomness concerning the fact that C-LL type of structure prevails over N-LL for peptides with one disulfide bond. Apparently, simple structures, like LL and HP, are self-sufficient to perform defense functions. The abundance of Gly and Pro relative to the ‘average protein’ is apparently linked with the necessity to promote such simple structures, which can be the building blocks to form longer and more structurally complex AMPs.
Thioether bond (Sulfide). DBAASP stores more than 85 ribosomal peptides with thioether bonds, ranging in length from 6 to 38 AA. Peptides are rich in Cys, Ser, and Thr. Thioether bonds are formed between the side chains of Cys and the side chains of Thr or Ser (their derivatives), as in lantibiotics [29], or between the side chains of Cys with the carbonyl group (main chain) of the neighbouring amino acid, as in peptides containinig thiazole ring. Thioether bonds predominantly form the amino acid lanthionine (see Table 1), and consequently, lantibiotics are more widespread among thioether bond-containing peptides. In lantibiotics, residues participating in side chain-to-side chain bonds are spaced by 2-5 AA along the chain, and thus, they form small rings of 4-7 AA (including Cys and Ser or Thr). Such rings are rarely crossed to form CR structures. For example, Nisin contains 5 rings, and among them, 3 do not cross each other. Only the 4-th and the 5-th are crossed. The majority of lantibiotics have one or more inter-crossed rings in their structures. However, there are peptides that contain a string of uncrossed rings, such as Enterocin W-beta (DBAASPR_1475).
Along with the thioether bonds, lantibiotics use disulfide (e.g., Enterocin W-alpha – DBAASPR_1474) and head-to-tail amide bonds (e.g., Subtilosin-A – DBAASPR_6076) for cyclization. Unlike other lantibiotics, Subtilosin-A have a ladder-like structure. The number of rings varies from 3 to 7.
Another posttranslational modification links the side chains of Cys with the side chains of Thr or Ser by TIE bonds and results in the creation of aminovinyl cysteine (AVC) rings (Table 1). Often, such rings appear in combination with lanthionine rings. However, there are some peptides containing AVC rings without LAN rings, such as Microvionin (DBAASPR_21973), Goadvionin B2 (DBAASPR_21974), Thioholgamides (DBAASPR_21975 and DBAASPR_21976), JBIR-140 (DBAASPR_21979), Thioviridamide (DBAASPR_21978). Structures adopted by these peptides resemble LL-type, because the Cys participating in the formation of AviCys is the C-terminal amino acid.
As mentioned above, there is an opportunity to link the side chain of Cys with the carbonyl group of the neighbouring amino acid through a TIE bond to create aromatic, unsaturated thiazole ring (THZ). The number of THZ-containing peptides in DBAASP is not high. Only 8 peptides, such as Micrococcins, Patellamides, Bottromycines, etc., have been identified (see Table 3). The structure of peptides containing more than one thiazole ring can be presented as a ‘string of rings’.
Among the peptides with TIE bonds in DBAASP, about 70 have a length of less than 25 AA and can be included in the set of ribosomal SCPs (Table 3).
Amide bond. Amide bonds are used more rarely than disulfide bonds to close rings in ribosomal SCPs. 66 SCPs are cyclized into NCB rings, while the number of peptides with LAC rings formed by bonding of side chain-to-main chain equals 30 (Table 3). It is worth noting that, although LAC rings in both, ribosomal and non-ribosomal SCPs are closed by isopeptide bonds, they are distinguished by the groups of side- and main chain atoms participating in bond formation. As mentioned above, LAC rings of the non-ribosomal SCPs and USCPs are closed by bonding between the amine group of side chains of basic residues (such as diaminopropionic acid – DAP, diaminobutyric acid – DAB, Lys, ornithine, 3-aminotetradecanoic acid, 3-aminohexadecanoic acid, etc.) and C-terminal carboxyl groups of the main chain. In contrast, in the ribosomal SCPs, the amine group of the N-terminus forms an isopeptide bond with the carboxyl side chain of a glutamic or aspartic acid residues to create a LAC ring of 7-9 AA. Such peptides form a family of ribosomally synthesized and post-translationally modified peptides called ‘lasso’ peptides (LP) [30]. Biosynthetic gene clusters (BGCs) for lasso peptides are presented in many bacterial genomes [31]. Often, the additional ring is formed in the C-terminal tail by disulfide bonds. The C-terminal tail is trapped within the ring either by bulky amino acids or disulfide bridges, or both [30]. Of note, due to the frequent appearance of Asp and Glu in LP, they are mainly anionic.
Amine, imine, ether, ester, and carbonyl bonds. In addition to disulfide, thioether, and amide bonds, other bonds, such as, amine, imine, ester, and carbonyl are used to close rings in ribosomal SCPs. However, cyclic peptides with such rings are scarce in DBAASP. As mentioned above, one of the reasons of scarcity could be insufficient data in DBAASP at the moment.
As in the case of non-ribosomal SCPs, in ribosomal SCPs, amine bonds are used to create LAC rings. Five such SCPs can be found in the DBAASP. They are Dynobactin A (DBAASPR_20296) and type B lantibiotics of Duramycin (DBAASPR_19104) and Cinnamycin (DBAASPR_21774) families. In contrast to non-ribosomal SCPs, where side chain-to-main chain bonding takes place, in ribosomal SCPs, amine bonds link the side chains of pairs of residues Lys and Ser (by forming Lysinoalanine in lantibiotics, e.g., in Duramicyn B – DBAASPR_19105) or His and Tyr (between the imidazole Nε2 of His and β-carbon of Tyr in Dynobactin A).
LAC rings are created in two SCPs of DBAASP (Bottromycin family, DBAASPR_21130) by bonding of N-terminal amine group nitrogen to the main chain carbonyl group carbon of 3-methyl-valine through imine bonds.
Carbon bonds between the 6-th carbon of the indole ring of Trp with the β-carbon of another residue (Lys, Asn, or Trp) create a LAC ring in ribosomal SCPs, such as Dynobactin A (DBAASPR_20296) and peptides of the Darobactin family (DBAASPR_17389, DBAASPR_22221-22222).
In the Ribosomal SCP Micrococcin P1 (DBAASPR_21120) [32], an unsaturated six-membered pyrimidine ring is formed by two carbon bonds: between the β-carbons of the first and 10-th dehydroalanines (DHA) and between the carbonyl group carbon of the 9-th Cys and the α-carbon of the first DHA.
Ether bonds, which are not represented in non-ribosomal SCPs, participate in the creation of five-membered unsaturated rings in ribosomal SCPs. SCPs of Patellamide family (e.g., DBAASPR_21121) consists of less aromatic oxazoline and meth-oxazoline rings, while Microcyclamide (DBAASPR_21139) and Wewakazoles (DBAASPR_21267) involve more aromatic oxazole and meth-oxazole. Additionally, ether bonds between the C7 indole carbon of one tryphtophan and the β-carbon of another tryphtophan form macrolactam (LAC) ring in post-ribosomally synthesized peptides of the Darobactin family (DBAASPR_17389, DBAASPR_22221- 22222).

4. An impact of cyclization on peptide potency

It is generally accepted that cyclization can increase the potency of antimicrobial peptides. It enhances the stability and structural rigidity of AMPs, resulting in improved resistance against enzymatic degradation and better interaction with microbial membranes, often leading to enhanced antimicrobial efficiency in vivo. Moreover, in certain in vitro studies involving Rana box-containing AMPs [33,34] or disulfide-bonded hairpin peptides [35], the doubt is expressed on the key role of cyclization on maintaining the antimicrobial potency. It is important to note that removing the Rana-box from AMPs or making substitutions in peptide sequence alters physicochemical properties of peptides, making it incorrect to compare their potencies under such conditions. To assess the impact of cyclization on peptide potency, it is advisable to compare cyclic and linear peptides with the same amino acid sequence. Furthermore, it is preferable that in vitro assessments of peptide potency intended for comparison are conducted under identical conditions (target bacterial strains, type of culture media, colony forming units (CFU), measures of activity, etc.). Therefore, we believe that the most reliable data for comparison are those obtained in the same laboratory.
In the DBAASP database, data on a series of peptides designed by various laboratories to evaluate the impact of cyclization on peptide activities have been compiled [36,37,38,39,40]. These series consist of pairs of peptides with identical AA sequences, but differing in the type of intrachain bonds used in cyclic counterparts. The cyclic counterparts are generated using either disulfide bonds [36,37,38,39] and head-to-tail amide bonds [40].
It is worth noting that, for a reliable interpretation of the comparison results, the peculiarities of the broth microdilution method (MIC assay) used to assess peptide potencies in vitro must be considered. The precision of the MIC assay allows for the differentiation of the potencies of two compared compounds with high reliability only if their MICs differ by more than one dilution (two or more).
In the Figure 9, the results of the MIC comparison for the pairs of peptides are presented. There are a total of 38 pairs, numerically labeled and arranged along the x-axis of Figure 9. The first 19 pairs belong to the disulfide bond set (DSB set), while the remaining 19 pairs belong to the head-to-tail amide bond set (HT set). For the given MICs assessed for the linear (MICl) and cyclic (MICc) peptides, differences in the MICs (ΔMIC) have been calculated as
M I C = s g n r r , where
r = log 2 M I C c / M I C l
It is worth emphasizing that the DSB set exhibits more heterogeneity in lengths (ranging from 13 to 24 AA) and chain topology (including both ‘lasso-like’ and ‘hairpin-like’ (including ‘ladder-like’) structures) compared to the HT set. The HT set comprises shorter peptides (6-9 AA) with similar chain topologies (ringed main chains).
As illustrated in Figure 9A, in the majority of pairs tested on gram-negative E. coli strains, the differences between the MIC values do not exceed one dilution interval (with only a few exceptions), suggesting that the cyclization does not significantly affect the activities of peptides against E. coli strains. However, in pairs tested on gram-positive S. aureus strains (Figure 9B), the situation is more complex, making it difficult to draw unambiguous conclusions. The comparison results indicate a decrease in MIC caused by cyclization in both the DSB and HT sets, with the tendency being particularly pronounced for peptides from the HT set. Furthermore, in the majority of cases, the differences do not exceed one dilution interval, falling within the margin of experimental error.
Consequently, we can conclude that the activity assessed on gram-negative E. coli strains is independent of peptide cyclization. However, regarding gram-positive S. aureus strains, cyclization may have some impact in specific cases, particularly for HT-cyclized peptides with lengths ranging from 6 to 9 AA. Nevertheless, overall, this influence is not prominently expressed.
In the case of peptides from the DSB set, the results can be explained by considering the significant role of chain topology in the potency of peptides with certain sequences [35,41]. For instance, it has been shown that the hairpin-like topology of a chain can be stabilized by factors other than disulfide bonding, such as cation-π interactions between aromatic and basic residues [42]. Given that AMPs are rich in aromatic and basic residues, which can support the hairpin topology and thus the potency of the peptide, disulfide bonding may have a negligible effect on potency. Moreover, the absence of intrachain bonds leads to increased flexibility of the AMP, which is often crucial for the antimicobial efficiency. However, the necessity of disulfide bonding to resist proteases remains clear.
Regarding the shorter peptides from the HT set, it could be speculated that the rigid cyclic structure of peptides holds particular value against the envelopes of gram-positive strains, where a different mechanism of action operates compared to gram-negative envelopes. However, further data are required to thoroughly examine this hypothesis.

5. Conclusions

DBAASP is a repository of more than 22000 AMPs, including more than 3900 natural peptides, among which over 2000 peptides are cyclic. Nature offers two ways of synthesis of cyclic peptides: post-ribosomal and non-ribosomal [43,44]. There are 769 non-ribosomally synthesized peptides in DBAASP, including 598 cyclic and small peptides (length ≤25 AA). The post-ribosomal method of creating cyclic peptides involves post-translational modifications. The number of post-ribosomally synthesized peptides in DBAASP equals 3192, including 484 small cyclic peptides.
It should be noted that DBAASP exclusively stores data on peptides for which in vitro testing has been performed. Therefore, several kinds of cyclic peptides that have not been tested in vitro, cannot be included in the database at the current time. This implies that DBAASP does not cover the entire spectrum of cyclic natural AMPs, and it is presumed that many types of cyclic structures not represented in DBAASP are available and/or will be revealed later. Nevertheless, the volume of DBAASP allows us to understand the level of variability of cyclic structures and the types of intrachain bonds responsible for chain cyclization. Consequently, our aim was to uncover common, statistically reliable tendencies, while some particular facts revealed by analysis should be specified with the influx of new data on cyclic peptides.
The data analysis reveals that both system of synthesis contribute to the remarkable structural variability of small natural AMPs. Simultaneously, the number of types of bonds providing the variability of cyclic structure is limited (see Table 1). The bonds used to cyclize natural peptides include amide, ester, ether, thioether, amine, imine, and carbon. The majority of rings created during non-ribosomal synthesis are closed by amide (isopeptide) bonds, and amide bonds are also employed for cyclization in post-ribosomal synthesis. Disulfide bonds are predominantly formed during post-ribosomal synthesis (Table 5). Only three fungal peptides with CST rings, two from the Malformin family (DBAASPN_21150 and DBAASPN_21150), and PM181110 (DBAASPN_21712) are synthesized non-ribosomally. Rings cyclized by ester bonds (LCN) are a hallmark of non-ribosomally synthesized peptides (Table 5), which suggests that an effective machinery for post-translational creation of intrachain ester bonds between proteinogenic residues did not evolve. It is worth noting that there is sufficient data to draw reliable conclusions for cyclic peptides where rings are closed by last three kinds of bonds. However, data on other bonds used by natural cyclic peptides for ring closing are insufficient to reveal common tendencies.
The closing of rings through amide bonds is a shared feature in both non-ribosomal and post-ribosomal synthesis (Table 5). However, the chains and residues participating in bond formation vary depending on whether it occures in eukaryotes or prokaryotes, or how they have been synthesized, non- or post- ribosomally. For instance, peptides with rings closed by head-to-tail amide bonds (NCB rings), which are formed between main chain atom groups, are synthesized by both eukaryotes and prokaryotes. Meanwhile, side chain-to-main chain isopeptide bonds (SMIB) for LAC ring closing are preferentially used by bacteria. Currently, only one exception has been identified: certain species of the fungus genus Aspergillus produce SMIB within the USCPs, such as peptides from the Sclerotiotide family (e.g., DBAASPN_18451). It has to be noted that, in post-ribosomal synthesis, SMIBs mainly involve the side chain carboxyl of acidic amino acids and N-terminal amine groups. In non-ribosomal synthesis, SMIBs are formed between the side chain amine group of basic amino acids and the C-terminal carboxyl group.
In general, rings occurring in post-ribosomally and non-ribosomally synthesized cyclic AMPs differ in type, size, participating chains, and bonds used for closing. Consequently, the protein complexes utilized in the two different methods of synthesizing natural cyclic peptides are distinct. However, there are rings that appear in both systems of synthesis used by nature. For instance, NCB rings closed by amide bonds between main chain group of atoms, or five-membered rings such as thiazole closed by thioether bonds between main and side chains, appear in both systems of synthesis. Consequently, it is reasonable to suppose that proteins responsible for creating similar rings in the two different systems of synthesis are homologous. However, this supposition has been disproven based on the comparison of the proteins involved in creating the rings in the two systems of synthesis. The comparison reveals that the proteins engaged in the creation of similar rings during post-ribosomal and non-ribosomal synthesis are not homologous. Therefore, the supposition on convergent evolution arises [44].
The requirement for cyclic structure formation has an impact on the amino acid composition of cyclic AMPs. Post-ribosomally synthesized AMPs are enriched with Cys, which forms the basis for the creation of ring-forming amino acids such as cystines and lanthionines. Ultra-short cyclic AMPs are rich in Pro and Gly, amino acids that promote chain turning and ring-closing. Some cyclic peptides have a non-traditional (for AMPs) composition of amino acids and utilize acidic residues to close the ring by forming amide bonds with N terminal amine groups, as seen in lasso peptides.
Rings in small cyclic peptides (≤25 AA) closed by different bonds vary in size and the topology of chains. Rings closed by amide, disulfide, and ester bonds constitute either micro- (e.g., DKP, DSL) or macrolactams. Amine, imine, carbon, thioether, and ether bonds are apparently responsible for creating micro rings such as four-membered (βLAC), five-membered (THZ, TZN, TZND, OXZ, OXZN, PRL), and six-membered (PYR) rings. The topology of macrolactams, i.e., polypeptide chains closed by bonds such as amide (creating NCB and LAC rings), disulfide (creating CST rings), and ester (creating LCN rings), can adopt four forms: ‘hairpin-like’, ‘lasso-like’, ‘string of rings’, and ‘crossed rings’ structures. The majority of short cyclic ribosomal peptides closing rings by DSB, AMD, and EST bonds adopt hairpin-like or lasso-like structures. These compact yet flexible structures provide membrane permeability, while intrachain bonds ensure chemical stability. Depending on membrane composition, such membrane-active peptides can perturb lipid bilayer leading to either disruption of membrane or their translocation through it to reach the cytoplasmic target.
We hope that the knowledge of the peculiarities of ring closing in natural peptides revealed in this work will contribute to the design of new cyclic compounds with desired pharmacological features.

Supplementary Materials

The following supporting information can be downloaded at the website of this paper posted on Preprints.org. Table S1 contains the raw data corresponding to Figure 9.

Author Contributions

Conceptualization, M.P.; Methodology, M.P., M.G., E.A., and N.S.; Validation, M.P., M.G., E.A., and N.S.; Formal Analysis, M.P., M.G., E.A., and N.S.; Investigation, M.P., M.G., E.A., and N.S.; Data Curation, M.G., E.A., and N.S.; Writing – Original Draft Preparation, M.P., M.G., E.A., and N.S.; Writing – Review & Editing, M.P., M.G., E.A., and N.S.; Visualization, M.P., M.G., E.A., and N.S.; Supervision, M.P.; Project Administration, M.P.; Funding Acquisition, M.P., M.G., E.A., and N.S.

Funding

This research was funded by Shota Rustaveli National Science Foundation of Georgia (SRNSFG) [FR-21-3368]; and International Science and Technology Center provided through National Institute of Allergy and Infectious Diseases / National Institutes of Health [G-2102]. The article processing charges (APC) was funded by SRNSFG.

Data Availability Statement

The data presented in this study are openly available in the DBAASP database at https://www.dbaasp.org; and are included in the article/supplementary material.

Conflicts of Interest

The authors declare no conflict of interest.

References

  1. Antimicrobial Resistance Collaborators. Global burden of bacterial antimicrobial resistance in 2019: A systematic analysis. Lancet 2022, 399, 629–655. [Google Scholar] [CrossRef]
  2. Wu, Q.; Yang, Z.; Nie, Y.; Shi, Y.; Fan, D. Multi-drug resistance in cancer chemotherapeutics: Mechanisms and lab approaches. Cancer Lett. 2014, 347, 159–166. [Google Scholar] [CrossRef]
  3. Fath, K.M.; Babakhaniyan, K.; Zokaei, M.; Yaghoubian, A.; Akbari, S.; Khorsandi, M.; Soofi, A.; Nabi-Afjadi, M.; Zalpoor, H.; Jalalifar, F.; et al. Anti-cancer peptide-based therapeutic strategies in solid tumors. Cell Mol Biol Lett. 2022, 27, 33. [Google Scholar] [CrossRef]
  4. Gori, A.; Lodigiani, G.; Colombarolli, S.G.; Bergamaschi, G.; Vitali, A. Cell penetrating peptides: Classification, mechanisms, methods of study, and applications. ChemMedChem. 2023, 18, e202300236. [Google Scholar] [CrossRef]
  5. Chen, J.H.; Lan, X.P.; Liu, Y.; Jia, A.Q. The effects of diketopiperazines from Callyspongia sp. on release of cytokines and chemokines in cultured J774A.1 macrophages. Bioorg Med Chem Lett. 2012, 22, 3177–3180. [Google Scholar] [CrossRef]
  6. Mookherjee, N.; Anderson, M.A.; Haagsman, H.P.; Davidson, D.J. Antimicrobial host defence peptides: Functions and clinical potential. Nat Rev Drug Discov. 2020, 19, 311–332. [Google Scholar] [CrossRef]
  7. Di, L. Strategic approaches to optimizing peptide ADME properties. AAPS J. 2015, 17, 134–143. [Google Scholar] [CrossRef]
  8. Lipinski, C.A.; Lombardo, F.; Dominy, B.W.; Feeney, P.J. Experimental and computational approaches to estimate solubility and permeability in drug discovery and development settings. Adv. Drug Deliv. Rev. 2001, 46, 3–26. [Google Scholar] [CrossRef]
  9. Griffith, M.; Griffith, O.L.; Coffman, A.C.; Weible, J.V.; McMichael, J.F.; Spies, N.C.; Koval, J.; Das, I.; Callaway, M.B.; Eldred, J.M.; et al. DGIdb: Mining the druggable genome. Nat. Methods 2013, 10, 1209–1210. [Google Scholar] [CrossRef]
  10. Doak, B.C.; Over, B.; Giordanetto, F.; Kihlberg, J. Oral druggable space beyond the rule of 5: Insights from drugs and clinical candidates. Chem. Biol. 2014, 21, 1115–1142. [Google Scholar] [CrossRef]
  11. Doak, B.C.; Zheng, J.; Dobritzsch, D.; Kihlberg, J. How beyond rule of 5 drugs and clinical candidates bind to their targets. J. Med. Chem. 2016, 59, 2312–2327. [Google Scholar] [CrossRef]
  12. Villar, E.A.; Beglov, D.; Chennamadhavuni, S.; Porco Jr, J.A.; Kozakov, D.; Vajda, S.; Whitty, A. How proteins bind macrocycles. Nat. Chem. Biol. 2014, 10, 723–731. [Google Scholar] [CrossRef]
  13. Dougherty, P.G.; Sahni, A.; Pei, D. Understanding cell penetration of cyclic peptides. Chem Rev. 2019, 119, 10241–10287. [Google Scholar] [CrossRef]
  14. Passioura, T.; Katoh, T.; Goto, Y.; Suga, H. Selection-based discovery of druglike macrocyclic peptides. Annu Rev Biochem. 2014, 83, 727–752. [Google Scholar] [CrossRef]
  15. Whitty, A.; Zhong, M.; Viarengo, L.; Beglov, D.; Hall, D.R.; Vajda, S. Quantifying the chameleonic properties of macrocycles and other high-molecular-weight drugs. Drug Discov. Today 2016, 21, 21,712–717. [Google Scholar] [CrossRef]
  16. Pirtskhalava, M.; Amstrong, A.A.; Grigolava, M.; Chubinidze, M.; Alimbarashvili, E.; Vishnepolsky, B.; Gabrielian, A.; Rosenthal, A.; Hurt, D.E.; Tartakovsky, M. DBAASP v3: Database of antimicrobial/cytotoxic activity and structure of peptides as a resource for development of new therapeutics. Nucleic Acids Res. 2021, 49, D288–D297. [Google Scholar] [CrossRef]
  17. Ramakrishnan, C.; Paul, P.K.C.; Ramanarayan, K. Cyclic peptides––small and big and their conformational aspects. J Biosci. 1985, 8, 239–251. [Google Scholar] [CrossRef]
  18. Hamed, R.B.; Gomez-Castellanos, J.R.; Henry, L.; Ducho, C.; McDonough, M.A.; Schofield, C.J. The enzymes of β-lactam biosynthesis. Nat Prod Rep. 2013, 30, 21–107. [Google Scholar] [CrossRef]
  19. Yu, X.; Li, L.; Sun, S.; Chang, A.; Dai, X.; Li, H.; Wang, Y.; Zhu, H. A cyclic dipeptide from marine fungus penicillium chrysogenum dxy-1 exhibits anti-quorum sensing activity. ACS Omega 2021, 6, 7693–7700. [Google Scholar] [CrossRef]
  20. Whiteley, M.; Diggle, S.P.; Greenberg, E.P. Progress in and promise of bacterial quorum sensing research. Nature 2017, 551, 313–320. [Google Scholar] [CrossRef]
  21. Imai, Y.; Meyer, K.J.; Iinishi, A.; Favre-Godal, Q.; Green, R.; Manuse, S.; Caboni, M.; Mori, M.; Niles, S.; Ghiglieri, M.; et al. A new antibiotic selectively kills Gram-negative pathogens. Nature 2019, 576, 459–464. [Google Scholar] [CrossRef] [PubMed]
  22. Urfer, M.; Bogdanovic, J.; Monte, F.L.; Moehle, K.; Zerbe, K.; Omasits, U.; Ahrens, C.H.; Pessi, G.; Eberl, L.; Robinson, J.A. A peptidomimetic antibiotic targets outer membrane proteins and disrupts selectively the outer membrane in escherichia coli. J Biol Chem. 2016, 291, 1921–1932. [Google Scholar] [CrossRef] [PubMed]
  23. UniProt Consortium. UniProt: A worldwide hub of protein knowledge. Nucleic Acids Res. 2019, 47, D506–D515. [Google Scholar] [CrossRef] [PubMed]
  24. Pirtskhalava, M.K.; Poroikov, V.V.; Tumanyan, V.G. Conformational possibilities of peptide-fragments containing disulfide bridge. Russ J Bioorg Chem. 1980, 6, 1125–1132. [Google Scholar]
  25. Jiang, Y.; Wu, Y.; Wang, T.; Chen, X.; Zhou, M.; Ma, C.; Xi, X.; Zhang, Y.; Chen, T.; Shaw, C.; et al. Brevinin-1GHd: A novel Hylarana guentheri skin secretion-derived Brevinin-1 type peptide with antimicrobial and anticancer therapeutic potential. Biosci Rep. 2020, 40, BSR20200019. [Google Scholar] [CrossRef] [PubMed]
  26. Lee, S.H.; Kim, E.H.; O’neal, J.T.; Dale, G.E.; Holthausen, D.J.; Bowen, J.R.; Quicke, K.M.; Skountzou, I.; Gopal, S.; George, S.; et al. The amphibian peptide Yodha is virucidal for Zika and dengue viruses. Sci Rep. 2021, 11, 602. [Google Scholar] [CrossRef] [PubMed]
  27. Wang, X.; Ren, S.; Guo, C.; Zhang, W.; Zhang, X.; Zhang, B.; Li, S.; Ren, J.; Hu, Y.; Wang, H. Identification and functional analyses of novel antioxidant peptides and antimicrobial peptides from skin secretions of four East Asian frog species. Acta Biochim Biophys Sin (Shanghai) 2017, 49, 550–559. [Google Scholar] [CrossRef] [PubMed]
  28. Randall, J.R.; DuPai, C.D.; Cole, T.J.; Davidson, G.; Groover, K.E.; Slater, S.L.; Mavridou, D.A.I.; Wilke, C.O.; Davies, B.W. Designing and identifying β-hairpin peptide macrocycles with antibiotic potential. Sci. Adv. 2023, 9, eade0008. [Google Scholar] [CrossRef] [PubMed]
  29. Willey, J.M.; van der Donk, W.A. Lantibiotics: Peptides of diverse structure and function. Annu Rev Microbiol. 2007, 61, 477–501. [Google Scholar] [CrossRef]
  30. Cheng, C.; Hua, Z.C. Lasso peptides: Heterologous production and potential medical application. Front Bioeng Biotechnol. 2020, 8, 571165. [Google Scholar] [CrossRef]
  31. Mevaere, J.; Goulard, C.; Schneider, O.; Sekurova, O.N.; Ma, H.; Zirah, S.; Afonso, C.; Rebuffat, S.; Zotchev, S.B.; Li, Y. An orthogonal system for heterologous expression of actinobacterial lasso peptides in Streptomyces hosts. Sci Rep. 2018, 8, 8232. [Google Scholar] [CrossRef] [PubMed]
  32. Duquesne, S.; Destoumieux-Garzón, D.; Peduzzi, J.; Rebuffat, S. Microcins, gene-encoded antibacterial peptides from enterobacteria. Nat Prod Rep. 2007, 24, 708–734. [Google Scholar] [CrossRef] [PubMed]
  33. Bao, K.; Yuan, W.; Ma, C.; Yu, X.; Wang, L.; Hong, M.; Xi, X.; Zhou, M.; Chen, T. Modification targeting the ‘Rana Box’ motif of a novel nigrocin peptide from Hylarana latouchii enhances and broadens its potency against multiple bacteria. Front Microbiol. 2018, 9, 2846. [Google Scholar] [CrossRef] [PubMed]
  34. Zou, W.; Zhang, Y.; Zhou, M.; Chen, X.; Ma, C.; Wang, T.; Jiang, Y.; Chen, T.; Shaw, C.; Wang, L. Exploring the active core of a novel antimicrobial peptide, palustrin-2LTb, from the Kuatun frog, Hylarana latouchii, using a bioinformatics-directed approach. Comput Struct Biotechnol J. 2022, 20, 6192–6205. [Google Scholar] [CrossRef] [PubMed]
  35. Sivanesam, K.; Kier, B.L.; Whedon, S.D.; Chatterjee, C.; Andersen, N.H. Hairpin structure stability plays a role in the activity of two antimicrobial peptides. FEBS Lett. 2016, 590, 4480–4488. [Google Scholar] [CrossRef] [PubMed]
  36. Smirnova, M.P.; Kolodkin, N.I.; Kolobov, A.A.; Afonin, V.G.; Afonina, I.V.; Stefanenko, L.I.; Shpen’, V.M.; Shamova, O.V.; Kolobov, A.A. Indolicidin analogs with broad-spectrum antimicrobial activity and low hemolytic activity. Peptides 2020, 132, 170356. [Google Scholar] [CrossRef] [PubMed]
  37. Wu, M.; Hancock, R.E. Improved derivatives of bactenecin, a cyclic dodecameric antimicrobial cationic peptide. Antimicrob Agents Chemother. 1999, 43, 1274–1276. [Google Scholar] [CrossRef] [PubMed]
  38. Antcheva, N.; Morgera, F.; Creatti, L.; Vaccari, L.; Pag, U.; Pacor, S.; Shai, Y.; Sahl, H.G.; Tossi, A. Artificial beta-defensin based on a minimal defensin template. Biochem J. 2009, 421, 435–447. [Google Scholar] [CrossRef] [PubMed]
  39. Li, S.; Lv, X.; Li, F.; Xiang, J. Characterization of a lymphoid organ specific anti-lipopolysaccharide factor from shrimp reveals structure-activity relationship of the lps-binding domain. Front Immunol. 2019, 10, 872. [Google Scholar] [CrossRef]
  40. Lohan, S.; Mandal, D.; Choi, W.; Konshina, A.G.; Tiwari, R.K.; Efremov, R.G.; Maslennikov, I.; Parang, K. Small amphiphilic peptides: Activity against a broad range of drug-resistant bacteria and structural insight into membranolytic properties. J Med Chem. 2022, 65, 665–687. [Google Scholar] [CrossRef]
  41. Mani, R.; Waring, A.J.; Lehrer, R.I.; Hong, M. Membrane-disruptive abilities of beta-hairpin antimicrobial peptides correlate with conformation and activity: A 31P and 1H NMR study. Biochim Biophys Acta. 2005, 1716, 11–18. [Google Scholar] [CrossRef] [PubMed]
  42. Stanojlovic, V.; Müller, A.; Moazzam, A.; Hinterholzer, A.; Ożga, K.; Berlicki, Ł.; Schubert, M.; Cabrele, C. A conformationally stable acyclic β-hairpin scaffold tolerating the incorporation of poorly β-sheet-prone amino acids. Chembiochem. 2022, 23, e202100604. [Google Scholar] [CrossRef] [PubMed]
  43. Zuber, P. Non-ribosomal peptide synthesis. Curr Opin Cell Biol. 1991, 3, 1046–1050. [Google Scholar] [CrossRef] [PubMed]
  44. McIntosh, J.A.; Donia, M.S.; Schmidt, E.W. Ribosomal peptide natural products: Bridging the ribosomal and nonribosomal worlds. Nat Prod Rep. 2009, 26, 537–559. [Google Scholar] [CrossRef] [PubMed]
Figure 1. (A) The composition of amino acids for non-ribosomal linear peptides of length 2-5.(https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24). (B) The composition of amino acids for non-ribosomal cyclic peptides of length 2-3. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24). (C) The composition of amino acids for non-ribosomal cyclic peptides of length 4-5. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24).
Figure 1. (A) The composition of amino acids for non-ribosomal linear peptides of length 2-5.(https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24). (B) The composition of amino acids for non-ribosomal cyclic peptides of length 2-3. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24). (C) The composition of amino acids for non-ribosomal cyclic peptides of length 4-5. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24).
Preprints 107225 g001aPreprints 107225 g001b
Figure 2. The distribution of peptides acting on specific cellular components (target objects): lipid bilayer, membrane protein, cytoplasmic protein, and DNA/RNA for non-ribosomal USCPs, non- and post-ribosomal SCPs from DBAASP. USCPs – ultra-short cyclic peptides (1-5 AA), SCPs – short cyclic peptides, (6-25 AA); Analysis performed on 13/05/24.
Figure 2. The distribution of peptides acting on specific cellular components (target objects): lipid bilayer, membrane protein, cytoplasmic protein, and DNA/RNA for non-ribosomal USCPs, non- and post-ribosomal SCPs from DBAASP. USCPs – ultra-short cyclic peptides (1-5 AA), SCPs – short cyclic peptides, (6-25 AA); Analysis performed on 13/05/24.
Preprints 107225 g002
Figure 4. Length distribution of ribosomal disulfide-bonded peptides up to 50 AA in DBAASP. (https://dbaasp.org/statistics?page=physicochemical-statistics, analysis performed on 13/05/24).
Figure 4. Length distribution of ribosomal disulfide-bonded peptides up to 50 AA in DBAASP. (https://dbaasp.org/statistics?page=physicochemical-statistics, analysis performed on 13/05/24).
Preprints 107225 g004
Figure 5. Amino acid composition of ribosomal short cyclic peptides (RSCPs) of length 6-25 AA from DBAASP relative to UniProt. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24).
Figure 5. Amino acid composition of ribosomal short cyclic peptides (RSCPs) of length 6-25 AA from DBAASP relative to UniProt. (https://dbaasp.org/statistics?page=compositional-statistics, analysis performed on 13/05/24).
Preprints 107225 g005
Figure 6. Frequencies of the appearance of i-spaced pairs of Cys residues (i = 0, 1, 2,…, 10) in disulfide-bonded ribosomal SCPs.
Figure 6. Frequencies of the appearance of i-spaced pairs of Cys residues (i = 0, 1, 2,…, 10) in disulfide-bonded ribosomal SCPs.
Preprints 107225 g006
Figure 7. Schematic representation of structures of single disulfide bond-containing peptides: (A) ‘lasso-like’ and (B) ‘hairpin-like’. Circles correspond to amino acids.
Figure 7. Schematic representation of structures of single disulfide bond-containing peptides: (A) ‘lasso-like’ and (B) ‘hairpin-like’. Circles correspond to amino acids.
Preprints 107225 g007
Figure 8. Schematic representation of structures of two disulfide bond-containing peptides: (A) ‘ladder- like’; (B) ‘string of rings’; (C) ‘crossed rings’. Circles correspond to amino acids.
Figure 8. Schematic representation of structures of two disulfide bond-containing peptides: (A) ‘ladder- like’; (B) ‘string of rings’; (C) ‘crossed rings’. Circles correspond to amino acids.
Preprints 107225 g008
Figure 9. Differences in the MICs (ΔMIC) for pairs of cyclic and linear AMPs with identical sequences, expressed in units of dilution: (A) Tested on E. coli strains; (B) Tested on S. aureus strains. Green rhombi – DSB set; pink rhombi – HT set.
Figure 9. Differences in the MICs (ΔMIC) for pairs of cyclic and linear AMPs with identical sequences, expressed in units of dilution: (A) Tested on E. coli strains; (B) Tested on S. aureus strains. Green rhombi – DSB set; pink rhombi – HT set.
Preprints 107225 g009
Table 1. Rings and ring-closing intra-chain bonds (RCBs) that formed post-ribosomally or non-ribosomally in small AMPs.
Table 1. Rings and ring-closing intra-chain bonds (RCBs) that formed post-ribosomally or non-ribosomally in small AMPs.
RCBs Participating Chains and Types of Rings In Non-Ribosomal Cyclic AMPs In Ribosomal Cyclic
AMPs
SMB** MMB** SSB**
TIE**
Preprints 107225 i001
THZN***
Preprints 107225 i002
+
THZD***
Preprints 107225 i003
THZD
Preprints 107225 i004
+
THZ***
Preprints 107225 i005
+ +
TZ***Preprints 107225 i006 +
LAN***
Preprints 107225 i007
+
MeLAN***
Preprints 107225 i008
+
AVC***
Preprints 107225 i009
+
MeAVC***
Preprints 107225 i010
+
ETH*
Preprints 107225 i011
OXZN****
Preprints 107225 i012
+
OXZ***
Preprints 107225 i013
+
MeOXZ***
Preprints 107225 i014
+
MeOXZN***
Preprints 107225 i015
+
LAC***
Preprints 107225 i016
+ +
EST*
Preprints 107225 i017
LCN***
Preprints 107225 i018
+
LCN
Preprints 107225 i019
+
Preprints 107225 i020
DSB*
CST***
Preprints 107225 i021
+ +
DSL***
Preprints 107225 i022
+
AMD*
Preprints 107225 i023
NCB***Preprints 107225 i024 + +
LAC
Preprints 107225 i025
+
LAC
Preprints 107225 i026
+
DKP***
Preprints 107225 i027
+
Preprints 107225 i028
IMN*
LAC
Preprints 107225 i029
+
AMN*
Preprints 107225 i030
βLAC***
Preprints 107225 i031
+
δLAC ***
Preprints 107225 i032
+
LAC
Preprints 107225 i033
+
LAC
Preprints 107225 i034
+
PRL***
Preprints 107225 i035
+
Preprints 107225 i036
CAR*
BAR***
Preprints 107225 i037
+
LAC
Preprints 107225 i038
+
LAC
Preprints 107225 i039
+
PYR***
Preprints 107225 i040
PYR
Preprints 107225 i041
+
* RCBs: AMD – amide; AMN – amine; IMN – imine; CAR – carbon; ETH – ether; TIE – thioether (sulfide); DSB – disulfide; EST – ester.** Participating chains: SMB – sidechain-mainchain bond; MMB – mainchain-mainchain bond; SSB – sidechain-sidechain bond.*** Types of rings: THZN – thiazoline (cyclic imino thioethers having an unsaturated five-membered structure, 1 double bond); THZD – thiazolidine (cyclic amino thioethers having a saturated five-membered structure); THZ – thiazole (cyclic imino thioethers having an unsaturated five-membered structure, 2 double bonds); TZ – thiazine (cyclic imino thioethers having an unsaturated six-membered structure, 2 double bonds); OXZN – oxazoline (cyclic imino ethers having an unsaturated five-membered structure, 1 double bond); OXZ – oxazole (cyclic imino ethers having an unsaturated five-membered structure, 2 double bonds); MeOXZN – methyloxazoline; MeOXZ – methyloxazole; LCN – macrocyclic ester, lactone; LAC – macrocyclic amide, macrolactam; βLAC – cyclic amide, β-lactam; δLAC – cyclic amide, δ-lactam; DKP – cyclic amide, head-to-tail cyclized dipeptide, 2,5-diketopiperazine; NCB – macrocyclic amide, head-to-tail cyclized macrolactam; CST – macrolactam closed by cysteine; DSL – disulfide linkage between two Cα; LAN – macrolactam closed by lanthionine; MeLAN – macrolactam closed by methyllanthionine; BAR – macrolactam stapled by biaryl group; PYR – pyridine (a benzene core in which one -CH group is replaced by a nitrogen atom); PRL – cyclic secondary amine, pyrrolidine; AVC – aminovinyl cysteine; MeAVC – aminovinyl methyl cysteine.Rx – ‘R’ denotes an atom or group of atoms, and ‘X’ denotes their number; RCT – side chain of C-terminal amino acids; RNT – side chain of N-terminal amino acids. ‘+’ – appeared, ‘–’ – not appeared.Analysis performed on 13/05/24.
Table 2. The numbers of ultra-short cyclic peptides (USCPs) with defined rings, ring-closing bonds (RCBs) and lengths.
Table 2. The numbers of ultra-short cyclic peptides (USCPs) with defined rings, ring-closing bonds (RCBs) and lengths.
Type of Synthesis Peptide Length Number of Peptides
Total Total Cyclic Certain RCBs, Participating Chains and Types of Rings
AMD*
MMB**
NCB***
AMD
MMB
DKP*
AMD
SMB**
LAC***
AMN*
SMB
βLAC***
AMN
SMB
PRL***
DSB
SSB
CST***
DSB*
SSB**
DSL***
EST*
SMB
LCN***
TIE*
SSB
THZD***
TIE
SSB
TZ***
Non-Ribosomal 2 93 89 0 79 0 4 2 1 1 6 3 1
3 47 32 13 3 12 0 0 0 0 4 0 0
4 55 43 31 0 0 0 0 0 0 10 0 0
5 68 61 21 0 0 0 0 2 0 39 0 0
Post-Ribosomal 2 0 0 0 0 0 0 0 0 0 0 0 0
3 0 0 0 0 0 0 0 0 0 0 0 0
4 1 1 1 0 0 0 0 0 0 0 0 0
5 6 1 1 0 0 0 0 0 0 0 0 0
* RCBs: AMD – amide; AMN – amine; TIE – thioether (sulfide); DSB – disulfide; EST – ester.** Parcipitating chains: MMB – mainchain-mainchain bond; SMB – sidechain-mainchain bond; SSB – sidechain-sidechain bond. *** Types of rings: NCB – macrocyclic amide, head-to-tail cyclized macrolactam; DKP – cyclic amide, head-to-tail cyclized dipeptide; LAC – macrocyclic amide, macrolactam; βLAC – cyclic amide, β-lactam; PRL – cyclic secondary amine, pyrrolidine; CST – macrolactam closed by cystine; DSL – disulfide linkage between two Cα; LCN – macrocyclic ester, lactone; THZD – thiazolidine; TZ – thiazine.Analysis performed on 13/05/24.
Table 3. The numbers of natural short cyclic peptides (SCPs) with defined rings and ring-closing intrachain bonds (RCBs) formed post-ribosomally or non-ribosomally.
Table 3. The numbers of natural short cyclic peptides (SCPs) with defined rings and ring-closing intrachain bonds (RCBs) formed post-ribosomally or non-ribosomally.
RCBs Participating Chains Types of Rings Non-Ribosomal Post-Ribosomal
DSB* SSB** CST*** 0 360
TIE* SSB LAN*** 0 22
MeLAN*** 0 22
AVC*** 0 12
MeAVC*** 0 5
SMB** THZ*** 4 8
THZN*** 4 0
THZD*** 1 0
AMD* MMB* NCB *** 138 66
SMB LAC*** 66 30
AMN* SMB LAC 8 0
δLAC 4 0
SSB LAC 0 5
EST* SMB LCN*** 142 0
ETH* SMB OXZ *** 0 4
MeOXZ*** 0 5
OXZN*** 0 2
MeOXZN*** 0 1
SSB LAC 3 1
IMN* MMB LAC 0 2
CAR* SMB LAC 0 6
SSB LAC 0 4
BAR*** 12 0
PYR*** 0 1
MMB PYR 0 1
Total Number 367 482
* RCBs: DSB – disulfide; TIE – thioether (sulfide); AMD – amide; AMN – amine; EST – ester; ETH – ether; IMN – imine; CAR – carbon.** Parcipitating chains: SSB – sidechain-sidechain bond; SMB – sidechain-mainchain bond; MMB – mainchain-mainchain bond.*** Types of rings: CST – macrolactam closed by cystine; LAN – macrolactam closed by lanthionine; MeLAN – macrolactam closed by methyllanthionine; AVC – aminovinyl cysteine; MeAVC – aminovinyl methyl cysteine; THZ – thiazole; THZN – thiazoline; THZD – thiazolidine; NCB – macrocyclic amide, head-to-tail cyclized macrolactam; LAC – macrocyclic amide, macrolactam; δLAC – cyclic amide, δ-lactam; LCN – macrocyclic ester, lactone; OXZ – oxazole; MeOXZ – methyloxazole; OXZN – oxazoline; MeOXZN – methyloxazoline; BAR – macrolactam stapled by biaryl group; PYR – pyridine.Analysis performed on 13/05/24.
Table 4. The types of structures of disulfide-bonded short cyclic peptides (SCPs).
Table 4. The types of structures of disulfide-bonded short cyclic peptides (SCPs).
Number of DSBs Structural Types
C-LL N-LL HP LD CR
1 227 1 60 NA NA
2 NA NA 29 29 8
3 NA NA 3 3 10
4 NA NA 4 4 16
C-LL and N-LL – C- and N-terminally looped ‘lasso-like’ structures, respectively; HP – ‘hairpin-like’; CR – crossed rings; LD – ‘ladder-like’; NA – not appeared; LD can simultaneously be considered as HP.
Table 5. The numbers of small peptides synthesized post-ribosomally or non-ribosomally in specific source organisms cyclized by amide, disulfide, and ester intrachain bonds (RCBs).
Table 5. The numbers of small peptides synthesized post-ribosomally or non-ribosomally in specific source organisms cyclized by amide, disulfide, and ester intrachain bonds (RCBs).
Type of Synthesis Non-ribosomal Post-ribosomal
RCBs AMD* DSB* EST* AMD DSB EST
Participating Chains MMB** SMB** SSB** SMB MMB SMB SSB SMB
Types of Rings NCB*** LAC*** CST*** LCN*** NCB LAC CST LCN
Sources Animalia 33 0 0 18 18 0 931 0
Plantae 3 0 0 1 139 0 238 0
Protista 0 0 0 0 0 0 2 0
Fungi 74 10 3 61 0 0 24 0
Bacteria 88 66 0 122 19 32 38 0
Archaea 0 0 0 2 0 0 1 0
* RCBs: AMD – amide; DSB – disulfide; EST – ester.** Participating chains: MMB – mainchain-mainchain bond; SMB – sidechain-mainchain bond; SSB – sidechain-sidechain bond.*** Types of rings: NCB – macrocyclic amide, head-to-tail cyclized macrolactam; LAC – macrocyclic amide, macrolactam; CST – macrolactam closed by cystine; LCN – macrocyclic ester, lactone.Analysis performed on 13/05/24.
Disclaimer/Publisher’s Note: The statements, opinions and data contained in all publications are solely those of the individual author(s) and contributor(s) and not of MDPI and/or the editor(s). MDPI and/or the editor(s) disclaim responsibility for any injury to people or property resulting from any ideas, methods, instructions or products referred to in the content.
Copyright: This open access article is published under a Creative Commons CC BY 4.0 license, which permit the free download, distribution, and reuse, provided that the author and preprint are cited in any reuse.
Prerpints.org logo

Preprints.org is a free preprint server supported by MDPI in Basel, Switzerland.

Subscribe

© 2024 MDPI (Basel, Switzerland) unless otherwise stated