Loading [MathJax]/extensions/MathML/content-mathml.js
info@auctorespublishing.com
+1-(302)-520-2644
+1-(302)-520-2644

Comparative Investigation of some Bioeffective short Peptides of Marine Derivation

  1. Home
  2. Journals
  3. Biotechnology and Bioprocessing
  4. Archives
Submit Guidelines

Research Article | DOI: https://doi.org/10.31579/2766-2314/130

Comparative Investigation of some Bioeffective short Peptides of Marine Derivation

  • Nawal Abd El-Baky 1*
  • Amro A. Amara 1

Protein Research Department, Genetic Engineering and Biotechnology Research Institute (GEBRI), City of Scientific Research and Technological Applications (SRTA-City), Alexandria, Egypt.

*Corresponding Author: Nawal Abd El-Baky, Protein Research Department, Genetic Engineering and Biotechnology Research Institute (GEBRI), City of Scientific Research and Technological Applications (SRTA-City), Alexandria, Egypt.

Citation: Nawal Abd El-Baky, Amro A. Amara, (2024), Comparative Investigation of some Bioeffective Short Peptides of Marine Derivation, J, Biotechnology and Bioprocessing, 5(1); DOI:10.31579/2766-2314/130

Copyright: © 2024, Nawal Abd El-Baky. This is an open access article distributed under the Creative Commons Attribution License, which permits unrestricted use, distribution, and reproduction in any medium, provided the original work is properly cited.

Received: 05 December 2023 | Accepted: 20 December 2023 | Published: 16 January 2024

Keywords: cyclic peptides; hybrids or consensus; intrinsic disorder; marine bioeffective peptides; physicochemical properties

Abstract

Short peptides especially marine ones are unique in their amino acids constituents that build their structure/function/specificity. Their physical and chemical properties make them diverse in their bioeffectiveness. In this study, different bioeffective peptides from marine bio-resources were investigated to map distinctive properties in their structure that could affect their function and specificity. The investigated short peptides were Pardaxins (1-5), Arenicin (1-2), Clavanin-(A-E), Aurelin, Dicynthaurin, Halocidin subunit (A-B), Tachyplesin (II-III), and Polyphemusin (I-II). Their different sequences were analyzed comparatively through alignment, their clustering (phylogeny tree), amino acids-%, and some other physical and chemical parameters. For each of them, cyclic structure based on disulfide bonds if present, and intrinsic disorder, were inspected. First alignment revealed significant diversion among their sequences, while removing unshaded amino acids from the alignment made them closer in their phylogenetic relationships and more evident to design hybrids or consensus of sequences and functions. Polyphemusin I and II had the highest charges (net) with values of +7, yet Pardaxins 1, 3, and 5 had 0 charges. All of them had higher numbers of hydrophilic amino acids than hydrophobic ones except for Aurelin, Tachyplesin II and III, and Polyphemusin I and II. Seven from twenty studied peptides had cyclic structures via disulfide bond between their cysteine residues. One amino acid was not found in all sequences; methionine. All sequences displayed other missed amino acids than methionine. Dicynthaurin is a naturally ordered peptide (via PONDR score), but Halocidin subunit B was totally disordered. This study paves the way for designing hybrids or consensus of these peptides in sequences and functions depending on their distinctive features.

Introduction

Marine organisms are a vast cradle of unusual bioeffective peptides. Peptides from several aquatic species with extraordinary impending nutraceutical and medicinal values are the emphasis of up-to-date research. Their crucial biological roles including microbicidal or antioxidative action, antitumor efficiency, cardioprotective properties, immunomodulatory role, etc. lead to continuous efforts and attempts to design and formulate these peptides (Ngo et al., 2012, Cheung et al., 2015, Ovchinnikova, 2019). Marine bioeffective peptides are distinguished from non-marine-derived ones by their resistance to proteolysis performed by proteases of gastrointestinal tract due to steric factors caused by their cyclic nature (Pangestuti and Kim, 2017, Saint Jean et al., 2017, Ibrahim et al., 2018, Phyo et al., 2018).

Marine bioeffective peptides are diverse and have different characteristics. Pardaxin possess antimicrobial actions and comprises 33 amino acids (Kolusheva et al., 2008).  It is also utilized for repelling sharks (Shai, 1994). Pardaxins 1, 2 and 3 were produced from Pardachirus pavoninus, and Pardaxins 4 and 5 were released by Pardachirus marmoratus (Shai et al., 1988). Pardaxin resembles melittin in its ability to lyse cells of bacteria and mammals (Oren and Shai, 1996). It might prompt apoptosis, leading to antitumor efficacy (Hsu et al., 2011). Arenicin is bactericidal to Gram-negative bacteria and induced apoptosis (Cho and Lee, 2011, Krenev et al., 2020). It also showed antifungal effects (Park et al., 2011). Arenicin 1 and 2 were isolated from Arenicola marina (Ovchinnikova et al., 2004, Umnyakova et al., 2018). 

Aurelin was isolated from Aurelia aurita (jellyfish) and exhibited antimicrobial functions (Ovchinnikova et al., 2006, Shenkarev et al., 2012). Aurelin shared similar features in its structure with defensins and toxins that block channels (Ovchinnikova et al., 2006). It displayed bactericidal effectiveness on gram-negative and -positive bacteria (Ovchinnikova et al., 2006). 

Clavanins (Clavanin-A, -B, -C, -D and –E) were obtained from Styela clava (Tunicate) hemocytes and revealed antimicrobial efficiency (Juliano et al., 2017, Lehrer et al., 2001, Menzel et al., 2002). Their antimicrobial functions, which rely on their α-helical structures, were thoroughly investigated against various microbes and even viruses such as herpes simplex virus 1 (HSV-1) (Lee et al., 1997, Lehrer et al., 2001, Carriel-Gomes et al., 2007). Clavanin-A antimicrobial efficiency contributes to innate immune responses of Styela clava (Juliano et al., 2017). Besides, it caused in vitro viral inactivation of non-enveloped adenovirus and rotavirus (Carriel-Gomes et al., 2007). Clavanin-B was proved to have antiviral effectiveness on HIV (Ireland et al., 2008). 

Dicynthaurin, another marine bioeffective peptide possessing antimicrobial activity, has been isolated from Halocynthia aurantium (Lee et al., 2001a). It demonstrated expansive action on Gram +ve and Gram -ve bacteria nonetheless was not active against Candida albicans (Lee et al., 2001a). Halocidin was also obtained from Halocynthia aurantium and has bactericidal and anti-Candida actions (Jang et al., 2003, Jang et al., 2006, Jang et al., 2002). A peptide derived from Halocidin could effectively treat oral infection with Candida albicans in mice (Shin et al., 2013). 

Polyphemusins from Limulus polyphemus have antimicrobial efficacy (Marggraf et al., 2018, Lee et al., 2022).  They are potent against various microbes and viruses (Lee et al., 2022, Panteleev et al., 2020, Tziveleka et al., 2003). Tachyplesins antimicrobial peptides were released from hemocytes of horseshoe crab (Ohta et al., 1992).Tachyplesin I could penetrate cells and has anticancer properties (Vernen et al., 2019a).

In the present investigation, twenty bioeffective peptides (Pardaxins (1-5), Arenicin (1-2), Clavanin-(A-E), Aurelin, Dicynthaurin, Halocidin subunit (A-B), Tachyplesin (II-III), and Polyphemusin (I-II)) from marine bio-resources were explored to record distinctive properties in their structures affecting their antimicrobial functions. Their different sequences were compared to each other, in addition to amino acids-%, and some other chemical and physical parameters. For each of them, intrinsic disorder and cyclic structures based on disulfide bonds were also inspected. The obtained data can lead to designing of hybrids or consensus of these peptides that may achieve more potent actions or even new functions. 

Material and Methods

Investigated marine bioeffective peptides

The twenty marine bioeffective peptides included in the present work are found in Table 1. 

Pardaxins (1-5), 1, Arenicin 2, Clavanin-A, Clavanin-B, Clavanin-E, Clavanin-C, Clavanin-D, Aurelin, Dicynthaurin, Halocidin subunit A, Halocidin subunit B, Tachyplesin-II, Polyphemusin I, Tachyplesin III, Polyphemusin II have lengths of 33, 21, 21, 23, 23, 22, 23, 22, 40, 30, 18, 15, 17, 18, 16, and 18 amino acids, respectively.

Peptide parent nameSpecific nameSequence
PardaxinsPardaxins lGFFALIPKIISSPLFKTLLSAVGSALSSSGEQE
Pardaxins 2GFFALIPKIISSPIFKTLLSAVGSALSSSGGQE
Pardaxins 3GFFAFIPKIISSPLFKTLLSAVGSALSSSGEQE
Pardaxins 4GFFALIPKIISSPLFKTLLSAVGSALSSSGGQE
Pardaxins 5GFFALIPKIISSPLFKTLLSAVGSALSSSGDQE
ArenicinArenicin 1RWCVYAYVRVRGVLVRYRRCW
Arenicin 2RWCVYAYVRIRGVLVRYRRCW
ClavaninsClavanin-AVFQFLGKIIHHVGNFVHGFSHVF
Clavanin-BVFQFLGRIIHHVGNFVHGFSHVF
Clavanin-CVFHLLGKIIHHVGNFVYGFSHVF
Clavanin-DFKLLGRIIHHVGNFVYGFSHVF
Clavanin-EFKLLGKIIHHVGNFVHGFSHVF
AurelinAurelinAACSDRAHGHICESFKSFCKDSGRNGVKLRANCKKTCGLC
DicynthaurinDicynthaurinILQKAVLDCLKAAGSSLSKAAITAIYNKIT
HalocidinHalocidin subunit AWLNALLHHGLNCAKGVLA
Halocidin subunit BALLHHGLNCAKGVLA
TachyplesinTachyplesin IIRWCFRVCYRGICYRKCR
Tachyplesin IIIWCFRVCYRGICYRKCR
PolyphemusinPolyphemusin IRRWCFRVCYRGFCYRKCR
Polyphemusin IIRRWCFRVCYKGFCYRKCR

Table 1: The studied peptides sequences.

Peptide sequences comparison with Ugene

Alignment No. 1

The up-to-date free version of Ugene 49.1 (Portable) was used to do the peptides sequences comparison (alignment) (Okonechnikov et al., 2012, Rose et al., 2018). This alignment was carried out by the alignment selection with Clustal Omega in the software since the investigated peptides had different lengths. The shading was done for amino acids with analogous physicochemical features. 

Alignment No. 2

Aiming to increase the probability of finding close associations between the twenty peptides sequences to facilitate planning of peptides hybrids, the unshaded amino acids after using the physicochemical shading selection were eliminated. The obtained sequences were aligned once more and even adjusted manually afterwards. This created reduced hypothetical peptides sequences. 

After running each of the fore-mentioned alignments, the phylogeny tree was generated. The involved mode in tree building was Neighbor Joining and distance matrix was Jones-Taylor-Thornton. 

Recording the properties of the twenty peptides’ complete sequences

The peptides general features, which included mass (Da), isoelectric point, net charge, hydrophobicity, and extinction coefficients were attained from the PepDraw web site (https://pepdraw.com).

Drawing cyclic peptides structures based on disulfide bonds 

If studied peptides had cyclic structures via disulfide bond between cysteine residues was detected via PepSMI (https://www.novoprolabs.com/tools/convert-peptide-to-smiles-string).

Concluding the amino acids-% in the investigated peptides

DAMBE 7.2.102 was utilized to estimate the amino acids-% in each of the twenty peptides sequences (complete ones not reduced hypothetical peptides sequences) under investigation (Xia, 2018). 

Data plotting of amino acids-%

The data of amino acids-% were considered in PAST 4.14. software using the clustering analysis option (Hammer et al., 2001). Neighbor joining was chosen with Hoen similarity index to cluster comparable peptide groups. 

The software plotted these data to elucidate the variations within amino acids or within the peptides. It also enabled contrasts between the diverse amino acids as well as viewing the standard deviations of their existence % among peptides, which mirrored the distinction within different peptides constituent from each amino acid. 

Helical Wheel (HW)

Anthreprot 6.9.3. determined the HW of each peptide sequence to record the hydrophobic amino acids and hydrophilic ones on the surface of HW and their distribution through the complete peptides sequences and reduced hypothetical ones (Deléage et al., 1989, Deléage et al., 1988, Deléage et al., 2001, Deléage, 2012). 

Specifying the amino acids that had disorder

All sequences were loaded in the PrDOS server to find any amino acid that has disorder or form with other ones, a disorder region. PrDOS accepted all the sequences lengths involved in this study (Ishida and Kinoshita, 2007). In addition, PONDR (www.pondr.com) was used to find natural disorder compared to crystallization data of natural proteins but in peptides that have sequence length of ≥ 30 amino acids. flDPnn was used to map propensity of disorder for investigated peptides (Hu et al., 2021). 

Results and Discussion

Peptides and specifically peptides from marine derivation showed unlimited potentials (Rawat et al., 2006, Liu et al., 2015, Simonetti et al., 2009, De Clercq, 2000). They found their way in different applications including fungicidal (Jang et al., 2006, Fan et al., 2023, Han et al., 2016), bactericial (Ohta et al., 1992, Panteleev et al., 2017, Saito et al., 1995, Serna et al., 2021, Li et al., 2020), virucidal (even to SARS-CoV-2) (Gogineni et al., 2015, Xie et al., 2016, Riccio et al., 2020, Raihan et al., 2021), anticancer (Vernen et al., 2019a, Ahmed et al., 2021, Vernen et al., 2019b, Zhang et al., 2021, Han et al., 2015, Lin et al., 2013), etc.

This study involved random selection of twenty bioeffective peptides from marine derivation; Arenicin (Arenicin 1 and 2) (Cho and Lee, 2011, Krenev et al., 2020, Park et al., 2011, Ovchinnikova et al., 2004, Umnyakova et al., 2018), Aurelin (Ovchinnikova et al., 2006, Shenkarev et al., 2012), Clavanins (-A, -B, -C, -D and -E) (Silva et al., 2015), Dicynthaurin (Bringezu et al., 2007b, Lee et al., 2001b, Bringezu et al., 2007a, Majerowicz et al., 2007, Wen et al., 2007), Halocidin subunit A and subunit B (Shin et al., 2010, Han et al., 2016), Pardaxins 1-5 (Oren and Shai, 1996), Polyphemusin I and II (Marggraf et al., 2018, Lee et al., 2022), and Tachyplesin II and III (Ohta et al., 1992, Vernen et al., 2019a). Their sequences were compared to each other, in addition to amino acids-%, and other chemical and physical parameters. Their disorder properties and cyclic structures based on disulfide bonds were also checked.

Alignment No. 1

Firstly, we performed the classical step to compare their sequences; alignment. Apparently, the peptides sequences were confusing to the alignment softwares because of their different lengths, and none of them gave the same result (data not shown). Meanwhile, Clustal-Omega gave more appropriate results (Figure 1), where similarity could be observed, even if the amino acids was still scattered. This gave us the idea to use the physicochemical properties mode to shade the similar amino acids and to eliminate the unshaded ones in the following run. 

Figure 1 revealed that peptides from the same source were aligned together. For example, Pardaxins 1 to 5 have been all aligned as a separate cluster. Besides, Tachyphemusin II and III were aligned with Polyphemusin I and II along with Arenicin 1 and 2. In the middle, Clavanin A-E, Aurelin, Halocidin subunit A and B showed similarity. In fact, Clustal Omega exhibited perfect grouping for analogous regions. 

Furthermore, the shadding was based on the physiochemical properties of the amino acids (Figure 1). 

Figure 1: Multiple alignment of the complete sequences of the investigated peptides with Clustal Omega (Ugene software).

Amino acids that appeared without coloring (either did not have unique physiochemical characteristics at their position or did not align with a region of distinctive physciochemical features) were removed and the reduced hypothetical sequences were realigned.

Alignment No. 2

After removing the unshaded amino acids and performing a second alignment run, the similarity became more obvious (Figure 2 and Figure 3). If Arenicin 1 and 2 was moved down manually, all of the hypothetical sequences would align together in different regions as in Figure 3. Such approach might be useful in designing peptides de novo or hybrids.

Figure 2: Multiple alignment of the reduced hypothetical sequences of the investigated peptides.

Figure 3: Multiple alignment of the reduced hypothetical sequences of the investigated peptides after adjustment manually.The phylogeny tree obtained from Ugene of the aligned twenty peptides sequences confirmed that removing unshaded amino acids increased the clustering of peptides of different classes together (Figure 4, and Figure 5). The cadre showed grouping of the peptides (Arenicin 2, Halocidin subunit A and B, Clavanin-B, Tachyplesin III, Clavanin-E, Polyphemusin I, Clavanin-A, Tachyplesin II, Clavanin-C, Polyphemusin II, and Clavanin- D) due to similarity after reducing the original sequence (Figure 5).

Figure 4: Phylogeny tree from the alignment of complete peptides sequences.

Figure 5: Phylogeny tree from the alignment of reduced hypothetical peptides sequences. The cadre showed grouping of the peptides due to similarity after reducing the original sequence.

Recording the properties of the twenty peptides’ complete sequences

Table 2 illustrated the properties of the twenty peptides amino acids sequences. Aurelin had the highest mass of 4300 Da and uppermost hydrophobicity of +45.67 Kcal * mol-1 among all peptides (Table 2). Polyphemusin I and II had the utmost charges (net) with values of +7, yet Pardaxins 1, 3, and 5 had 0 charges (Table 2). Arenicin 1 and 2 exhibited the leading isoelectric point of 11.14 (Table 2).

Peptide Name            Length

Mass

(Da)

isoelectric

point (pI)

Net chargeHydrophobicity Kcal * mol-1

Extinction coefficient

M-1* cm-I

Extinction coefficient2

M-1* cm-I

Arenicin 1212758.4611.14+6+10.511547015595
Arenicin 2212772.4811.14+6 +9.8515595 15470 
Aurelin404300.008.94+5+45.67375 0 
Clavanin-A232665.409.93+1+11.67 0 0 
Clavanin-B232693.4010.90+1+10.60 0 
Clavanin-C232666.429.59+1+10.65 1490 1490 
Clavanin-D222586.4010.38+2+10.591490 1490 
Clavanin-E222532.3910.59+2+14.62 0 0
Dicynthaurin 303103.76 9.91+3+19.72 1490 1490 
Halocidin subunit A 181929.04 9.04+1+12.04 5500 5500 
Halocidin subunit B151515.839.04+1+14.53 00 
Pardaxins 1333393.846.970+15.0300 
Pardaxins 2333321.829.95+1+12.68 0 0 
Pardaxins 3333427.826.970+14.5700
Pardaxins 4333321.829.95+1+12.55 0 0 
Pardaxins 5333379.826.930+15.04 00 
Polyphemusin I182457.1810.54+7+15.2487308480
Polyphemusin II182429.1710.33+7+16.2387308480
Tachyplesin II172267.0910.01+6+14.0287308480
Tachyplesin III162110.999.48+5+12.2187308480

Table 2: Properties of the different peptides’ amino acids sequences.

Drawing cyclic peptides structures based on disulfide bonds

Seven studied peptides (Arenicin 1 and 2, Polyphemusin I and II, Tachyplesin II and III, and Aurelin) had cyclic structures via disulfide bond between cysteine residues (Figure 6). The other thirteen sequences formed only linear structures due to the absence of cysteine residues (Pardaxins l-5 and Clavanin-A-E) or presence of only one cysteine (Dicynthaurin, and 

Halocidin subunit A and B). The most unique cyclic structure was for Aurelin as the peptide contains 6 cysteine residues (Figure 6). These cyclic structures could make them more unaffected by proteolysis performed by proteases of gastrointestinal tract due to steric factors (Pangestuti and Kim, 2017, Saint Jean et al., 2017, Ibrahim et al., 2018, Phyo et al., 2018).

Figure 6: Unique cyclic structures of seven from twenty studied peptides via disulphide bond between cysteines

Concluding the amino acids-% in the investigated peptides

To deepen our understanding to the 20 selected marine bioeffective peptides and their amino acids constituents, amino acid-% distribution was checked. 

Stacked chart of the amino acids-% of investigated peptides in Figure 7 displayed that methionine did not exist at all in all of the 20 peptides. 

Figure 7: Stacked chart of the amino acids-% of investigated peptides.

 

Figure 8: Amino acids-% in investigated peptides with standard deviation

Figure 9: Cluster analysis for amino acids-% in investigated peptides

Demonstrated amino acids-% in investigated peptides with standard deviation of their existence % among peptides, which mirrored the distinction within different peptides constituent from each amino acid. This brief statistical estimation revealed the discrepancy among the amino acids as well as within each peptide class.

The composition of amino acids is strictly associated with bioeffectiveness of these bioeffective peptides. For instance, Arenicin 1 was found to have fungicidal efficacy that is dependent on tryptophan in its N terminus and arginine in its β-turn (Park et al., 2011). 

Cluster analysis for amino acids-% in investigated peptides in Figure 9 indicated that the twenty amino acids-% classified peptides to six major groups, where Pardaxins 1-5 were clustered in one group. The same case was found with Clavanin-A, -B, -C, -D and -E.  Dicynthaurin and Aurelin each could be considered as separate groups (Figure 9). Polyphemusin I and II, Tachyplesin II and III along with Arenicin 1 and 2 are grouped in one cluster. While, Halocidin subunit A and B were clustered together as in Figure 9.

Helical Wheel (HW)

The HW analysis investigated the distribution of the hydrophobic and hydrophilic amino acids on helical folding form of peptides. The results were summarized in Table 3. Table 3 showed that all peptides had higher numbers of hydrophilic amino acids than hydrophobic ones except for Aurelin, Tachyplesin II and III, and Polyphemusin I and II. Hydrophobic amino acids in peptides have critical contribution to their bactericidal effectiveness as they enhance peptides ability to bind and disintegrate membranes of bacteria (Saint Jean et al., 2017). The overall differences between the numbers of the hydrophobic and hydrophilic amino acids ranged from two to four (Table 3). Meanwhile, the oddest difference was in case of Aurelin, which showed differences equal to 10 amino acids (Table 3)

NAME

HYDROPHOBIC AND HYDROPHILIC AMINO ACIDS DISTRIBUTIONNO. OF HYDROPHOBIC AMINO ACIDSNO. OF HYDROPHILIC AMINO ACIDS
PARDAXINS 1GFFALIPKIISSPLFKTLLSAVGSALSSSGEQE1315
PARDAXINS 2GFFALIPKIISSPIFKTLLSAVGSALSSSGGQE1215
PARDAXINS 3GFFAFIPKIISSPLFKTLLSAVGSALSSSGEQE1315
PARDAXINS 4GFFALIPKIISSPLFKTLLSAVGSALSSSGGQE1215
PARDAXINS 5GFFALIPKIISSPLFKTLLSAVGSALSSSGDQE1315
ARENICIN 1RWCVYAYVRVRGVLVRYRRCW69
ARENICIN 2RWCVYAYVRIRGVLVRYRRCW69
CLAVANIN-AVFQFLGKIIHHVGNFVHGFSHVF811
CLAVANIN-BVFQFLGRIIHHVGNFVHGFSHVF811
CLAVANIN-EFKLLGKIIHHVGNFVHGFSHVF811
CLAVANIN-CVFHLLGKIIHHVGNFVYGFSHVF711
CLAVANIN-DFKLLGRIIHHVGNFVYGFSHVF711
AURELINAACSDRAHGHICESFKSFCKDSGRNGVKLRANCKKTCGLC2010
DICYNTHAURINILQKAVLDCLKAAGSSLSKAAITAIYNKIT1215
HALOCIDIN SUBUNIT AWLNALLHHGLNCAKGVLA510
HALOCIDIN SUBUNIT BALLHHGLNCAKGVLA48
TACHYPLESIN IIRWCFRVCYRGICYRKCR64
POLYPHEMUSIN IRRWCFRVCYRGFCYRKCR74
TACHYPLESIN IIIWCFRVCYRGICYRKCR54
POLYPHEMUSIN IIRRWCFRVCYKGFCYRKCR74

Table 3:The hydrophobic amino acids (red) and hydrophilic ones (blue) distribution through the complete peptides sequences.

Table 4 demonstrated the hydrophobic amino acids (red) and hydrophilic ones (blue) distribution through reduced hypothetical peptides sequences. The numbers of hydrophobic amino acids was significantly reduced after sequence reductions. 

Name

Hydrophobic and hydrophilic amino acids distributionNo. of hydrophobic amino acidsNo. of hydrophilic amino acids
Pardaxins 1GFFLIPIISSPLFKTLLSAVGSALSS814
Pardaxins 2GFFLIPIISSPIFKTLLSAVGSALSS814
Pardaxins 3GFFFIPIISSPLFKTLLSAVGSALSS814
Pardaxins 4GFFLIPIISSPLFKTLLSAVGSALSS814
Pardaxins 5GFFLIPIISSPLFKTLLSAVGSALSS814
Arenicin 1RWCYAYRVRGVYRRC54
Arenicin 2RWCYAYRIRGVYRRC54
Clavanin-AVFFLGIIHHVGNVHVF48
Clavanin-BVFFLGIIHHVGNVHVF49
Clavanin-EFLLGIIHHVGNVHVF49
Clavanin-CVFLLGIIHHVGNVHVF410
Clavanin-DFLLGIIHHVGNVHVF49
AurelinAACDHICESFKSFCDSGNGVLRKKTCLC138
DicynthaurinILQKAVLDCLKAAGSSLSAIYNKIT1012
Halocidin subunit AWLNLLHHGLNAA47
Halocidin subunit BLLHHGLNAA35
Tachyplesin IIRWCRVCRGYRKC52
Polyphemusin IRRWCRVCRGYRKC62
Tachyplesin IIIWCRVCRGYRKC42
Polyphemusin IIRRWCRVCKGYRKC52

Table 4:The hydrophobic amino acids (red) and hydrophilic ones (blue) distribution through reduced hypothetical peptides sequences.

Specifying the amino acids that had disorder

The amino acids that had disorder in the twenty investigated complete sequences of marine bioeffective peptides were specified from PrDOS as demonstrated in Figure 10. The red color in the sequences signified amino acids with disorder. According to PrDOS, Halocidin subunit B (15 amino acids length) was totally disordered (Figure 10). None of the peptides was 

totally ordered, nevertheless Polyphemusin II (18 amino acids) had only one disordered amino acid (Figure 10). 

PrDOS could act on all investigated sequences whatever their lengths are (Ishida and Kinoshita, 2007). This algorithm anticipates disorder from sequences of amino acids. The probabilities below 0.5 in the curve showed ordered residues, while those above 0.5 showed disordered ones (Figure 10). 

Figure 10: Disorder probability for 20 peptides’ complete sequences obtained from PrDOS server.

The red colored residues in the sequences represent amino acids with intrinsic disorder.

Figure 11 revealed disorder in seven peptides complete sequences (that comprise 30 or more amino acids) obtained from PONDR. The red color in the sequences exemplifies amino acids with disorder. Dicynthaurin is a naturally ordered peptide (via PONDR score). Aurelin had the highest numbers of disordered residues (16 of 40), Pardaxins 1, 3, and 5 had 11 disordered residues, and Pardaxins 2 and 4 had 12 ones (Figure 11).  The scores below 0.5 showed ordered residues, while those above 0.5 showed disordered ones (Figure 11). 

Figure 11: Intrinsic disorder for seven peptides’ complete sequences (that comprise 30 or more amino acids) obtained from PONDR. The red color in the sequences represents amino acids with disorder.

As can be concluded from Figure 10 and Figure 11, there are differences in the records obtained from both algorithms. These differences originated from the variation in their anticipation method since PONDR depend on matching queries to crystallization data of natural proteins not the sequence itself (Li et al., 1999).

Figure 12 presented flDPnn propensity of disorder in investigated peptides sequences. Records of disorder are vital in designing of de novo or hybrid peptides as they are strictly linked to peptides bioeffectiveness (Dunker et al., 2000, Ward et al., 2004, Shimizu et al., 2007).

Figure 12: flDPnn propensity of disorder in investigated peptides sequences.

Conclusion

This study investigated twenty bioeffective peptides from marine derivation. It focused on comparative exploration for amino acids-% in the different peptides, alignment, phylogeny tree, chemical and physical parameters, hydrophobic/hydrophilic amino acids, cyclic structures based on disulfide bonds, and intrinsic disorder. The study included reducing the original sequences to map the internal similarity that might be hidden in the different peptides. The achieved results can be beneficial in planning de novo or hybrids of these peptides.

Author Contributions:

N.A.E.-B. and A.A.A.: Conceptualization, Methodology, Data curation and analysis, draft writing.  N.A.E.-B. revised and finalized the work. Both authors read and agreed to the manuscript.

Conflicts of Interest:

The authors declare no conflict of interest.

References

  1. Ahmed, S, Mirzaei, H, Aschner, M, Khan, A, Al-Harrasi, A et al. (2021). Marine peptides in breast cancer: Therapeutic and mechanistic understanding. Biomedicine & Pharmacotherapy 142: 112038
    View at Publisher | View at Google Scholar
  2. Bringezu, F, Majerowicz, M, Maltseva, E, Wen, S, Brezesinski, G et al. (2007a). Penetration of the antimicrobial peptide dicynthaurin into phospholipid monolayers at the liquid-air interface. Chembiochem 8: 1038-1047
    View at Publisher | View at Google Scholar
  3. Bringezu, F, Wen, S, Dante, S, Hauss, T, Majerowicz, M et al. (2007b). The insertion of the antimicrobial peptide dicynthaurin monomer in model membranes: thermodynamics and structural characterization. Biochemistry 46: 5678-5686
    View at Publisher | View at Google Scholar
  4. Carriel-Gomes, MC, Kratz, JM, Barracco, MA, Bachére, E, Barardi, CRM et al. (2007). In vitro antiviral activity of antimicrobial peptides against herpes simplex virus 1, adenovirus, and rotavirus. Memórias do Instituto Oswaldo Cruz 102: 469-472
    View at Publisher | View at Google Scholar
  5. Cheung, RCF, Ng, TB and Wong, JH (2015). Marine Peptides: Bioactivities and Applications. Marine drugs 13: 4006-4043
    View at Publisher | View at Google Scholar
  6. Cho, J and Lee, DG (2011). The antimicrobial peptide arenicin-1 promotes generation of reactive oxygen species and induction of apoptosis. Biochimica et biophysica acta 1810: 1246-1251
    View at Publisher | View at Google Scholar
  7. De Clercq, E (2000). Novel compounds in preclinical/early clinical development for the treatment of HIV infections. Rev Med Virol 10: 255-277
    View at Publisher | View at Google Scholar
  8. Deléage, G (2012). An interactive 3D viewer of molecules compatible with the suite of ANTHEPROT programs. Journal of Biophysical Chemistry 03: 35-38
    View at Publisher | View at Google Scholar
  9. Deléage, G, Clerc, FF and Roux, B (1989). ANTHEPROT: IBM PC and Apple Macintosh versions. Bioinformatics 5: 159-160
    View at Publisher | View at Google Scholar
  10. Deléage, G, Clerc, FF, Roux, B and Gautheron, DC (1988). ANTHEPROT: a package for protein sequence analysis using a microcomputer. Bioinformatics 4: 351-356
    View at Publisher | View at Google Scholar
  11. Deléage, G, Combet, C, Blanchet, C and Geourjon, C (2001). ANTHEPROT: An integrated protein sequence analysis software with client/server capabilities. Computers in Biology and Medicine 31: 259-267
    View at Publisher | View at Google Scholar
  12. Dunker, AK, Obradovic, Z, Romero, P, Garner, EC and Brown, CJ (2000). Intrinsic protein disorder in complete genomes. . Genome Inform Ser Workshop Genome Inform 11: 161-171
    View at Publisher | View at Google Scholar
  13. Fan, L, Wei, Y, Chen, Y, Jiang, S, Xu, F, et al. (2023). Epinecidin-1, a marine antifungal peptide, inhibits Botrytis cinerea and delays gray mold in postharvest peaches. Food Chem 403: 134419
    View at Publisher | View at Google Scholar
  14. Gogineni, V, Schinazi, RF and Hamann, MT (2015). Role of Marine Natural Products in the Genesis of Antiviral Agents. Chemical reviews 115: 9655-9706
    View at Publisher | View at Google Scholar
  15. Hammer, O, Harper, DAT and Ryan, PD (2001). PAST: Paleontological Statistic package for education and data analysis. Palaeontologia Electronica 4: 9
    View at Publisher | View at Google Scholar
  16. Han, J, Jyoti, MA, Song, HY and Jang, WS (2016). Antifungal Activity and Action Mechanism of Histatin 5-Halocidin Hybrid Peptides against Candida ssp. PLoS One 11: e0150196
    View at Publisher | View at Google Scholar
  17. Han, Y, Cui, Z, Li, YH, Hsu, WH and Lee, BH (2015). In Vitro and in Vivo Anticancer Activity of Pardaxin against Proliferation and Growth of Oral Squamous Cell Carcinoma. Marine drugs 14: 2
    View at Publisher | View at Google Scholar
  18. Hsu, J-C, Lin, L-C, Tzen, JTC and Chen, J-Y (2011). Pardaxin-induced apoptosis enhances antitumor activity in HeLa cells. Peptides 32: 1110-1116
    View at Publisher | View at Google Scholar
  19. Hu, G, Katuwawala, A, Wang, K, Wu, Z, Ghadermarzi, S, et al. (2021). flDPnn: Accurate intrinsic disorder prediction with putative propensities of disorder functions. Nature communications 12: 4438-4438
    View at Publisher | View at Google Scholar
  20. Ibrahim, AH, Attia, EZ, Hajjar, D, Anany, MA, Desoukey, SY, et al. (2018). New Cytotoxic Cyclic Peptide from the Marine Sponge-Associated Nocardiopsis sp. UR67. Marine drugs 16: 290
    View at Publisher | View at Google Scholar
  21. Ireland, DC, Wang, CKL, Wilson, JA, Gustafson, KR and Craik, DJ (2008). Cyclotides as natural anti-HIV agents. Biopolymers 90: 51-60
    View at Publisher | View at Google Scholar
  22. Ishida, T and Kinoshita, K (2007). PrDOS: prediction of disordered protein regions from amino acid sequence. Nucleic acids research 35: W460-W464
    View at Publisher | View at Google Scholar
  23. Jang, WS, Kim, CH, Kim, KN, Park, SY, Lee, JH, Son, SM and Lee, IH (2003). Biological activities of synthetic analogs of halocidin, an antimicrobial peptide from the tunicate Halocynthia aurantium. Antimicrob Agents Chemother 47: 2481-2486
    View at Publisher | View at Google Scholar
  24. Jang, WS, Kim, HK, Lee, KY, Kim, SA, Han, YS et al. (2006).Antifungal activity of synthetic peptide derived from halocidin, antimicrobial peptide from the tunicate, Halo Cynthia aurantium. FEBS Lett 580: 1490-1496
    View at Publisher | View at Google Scholar
  25. Jang, WS, Kim, KN, Lee, YS, Nam, MH and Lee, IH (2002). Halocidin: a new antimicrobial peptide from hemocytes of the solitary tunicate, Halocynthia aurantium. FEBS Lett 521: 81-86
    View at Publisher | View at Google Scholar
  26. Juliano, SA, Pierce, S, deMayo, JA, Balunas, MJ and Angeles-Boza, AM (2017). Exploration of the Innate Immune System of Styela clava: Zn (2+) Binding Enhances the Antimicrobial Activity of the Tunicate Peptide Clavanin A. Biochemistry 56: 1403-1414
    View at Publisher | View at Google Scholar
  27. Kolusheva, S, Lecht, S, Derazon, Y, Jelinek, R and Lazarovici, P (2008). Pardaxin, a fish toxin peptide interaction with a biomimetic phospholipid/polydiacetylene membrane assay. Peptides 29: 1620-1625
    View at Publisher | View at Google Scholar
  28. Krenev, IA, Umnyakova, ES, Eliseev, IE, Dubrovskii, YA, Gorbunov, NP, et al. (2020). Antimicrobial Peptide Arenicin-1 Derivative Ar-1-(C/A) as Complement System Modulator. Marine drugs 18
    View at Publisher | View at Google Scholar
  29. Lee, IH, Lee, YS, Kim, CH, Kim, CR, Hong, T, et al.(2001a). Dicynthaurin: an antimicrobial peptide from hemocytes of the solitary tunicate, Halocynthia aurantium. Biochimica et biophysica acta 1527: 141-148
    View at Publisher | View at Google Scholar
  30. Lee, IH, Lee, YS, Kim, CH, Kim, CR, Hong, T, Menzel, L, et al.(2001b). Dicynthaurin: an antimicrobial peptide from hemocytes of the solitary tunicate, Halocynthia aurantium. Biochimica et biophysica acta 1527: 141-148
    View at Publisher | View at Google Scholar
  31. Lee, IH, Zhao, C, Cho, Y, Harwig, SS, Cooper, EL et al. (1997). Clavanins, alpha-helical antimicrobial peptides from tunicate hemocytes. FEBS letters 400: 158-162
    View at Publisher | View at Google Scholar
  32. Lee, Y-S, Feng, C-W, Peng, M-Y, Chen, Y-C and Chan, T-F (2022). Antiosteoporosis effects of a marine antimicrobial peptide pardaxin via regulation of the osteogenesis pathway. Peptides 148: 170686-170686
    View at Publisher | View at Google Scholar
  33. Lehrer, RI, Lee, IH, Menzel, L, Waring, A and Zhao, C (2001). Clavanins and styelins, alpha-helical antimicrobial peptides from the hemocytes of Styela clava. Advances in experimental medicine and biology 484: 71-76
    View at Publisher | View at Google Scholar
  34. Li, W, Li, Y, Sun, P, Zhang, N, Zhao, Y, et al. (2020). Antimicrobial peptide-modified silver nanoparticles for enhancing the antibacterial efficacy. RSC Adv 10: 38746-38754
    View at Publisher | View at Google Scholar
  35. Li, X, Romero, P, Rani, M, Dunker, AK and Obradovic, Z (1999). Predicting protein disorder for N-, C-, and internal regions. Genome Inform SerWorkshop Genome Inform 10: 30-40
    View at Publisher | View at Google Scholar
  36. Lin, MC, Hui, CF, Chen, JY and Wu, JL (2013). Truncated antimicrobial peptides from marine organisms retain anticancer activity and antibacterial activity against multidrug-resistant Staphylococcus aureus. Peptides 44: 139-148
    View at Publisher | View at Google Scholar
  37. Liu, CB, Shan, B, Bai, HM, Tang, J, Yan, LZ et al. (2015). Hydrophilic/hydrophobic characters of antimicrobial peptides derived from animals and their effects on multidrug resistant clinical isolates. Dongwuxue Yanjiu 36: 41-47
    View at Publisher | View at Google Scholar
  38. Majerowicz, M, Waring, AJ, Wen, S and Bringezu, F (2007). Interaction of the antimicrobial peptide dicynthaurin with membrane phospholipids at the air-liquid interface. J Phys Chem B 111: 3813-3821
    View at Publisher | View at Google Scholar
  39. Marggraf, MB, Panteleev, PV, Emelianova, AA, Sorokin, MI, Bolosov, IA, et al. (2018). Cytotoxic Potential of the Novel Horseshoe Crab Peptide Polyphemusin III. Marine drugs 16
    View at Publisher | View at Google Scholar
  40. Menzel, LP, Lee, IH, Sjostrand, B and Lehrer, RI (2002). Immunolocalization of clavanins in Styela clava hemocytes. Developmental and comparative immunology 26: 505-515
    View at Publisher | View at Google Scholar
  41. Ngo, D-H, Vo, T-S, Ngo, D-N, Wijesekara, I and Kim, S-K (2012). Biological activities and potential health benefits of bioactive peptides derived from marine organisms. International Journal of Biological Macromolecules 51: 378-383
    View at Publisher | View at Google Scholar
  42. Ohta, M, Ito, H, Masuda, K, Tanaka, S, Arakawa, Y, et al. (1992). Mechanisms of antibacterial action of tachyplesins and polyphemusins, a group of antimicrobial peptides isolated from horseshoe crab hemocytes. Antimicrob Agents Chemother 36: 1460-1465
    View at Publisher | View at Google Scholar
  43. Okonechnikov, K, Golosova, O and Fursov, M (2012). Unipro UGENE: a unified bioinformatics toolkit. Bioinformatics 28: 1166-1167
    View at Publisher | View at Google Scholar
  44. Oren, Z and Shai, Y (1996). A class of highly potent antibacterial peptides derived from pardaxin, a pore-forming peptide isolated from Moses sole fish Pardachirus marmoratus.. . European Journal of Biochemistry 237: 303–310
    View at Publisher | View at Google Scholar
  45. Ovchinnikova, TV (2019). Structure, Function, and Therapeutic Potential of Marine Bioactive Peptides. Marine drugs 17: 505
    View at Publisher | View at Google Scholar
  46. Ovchinnikova, TV, Aleshina, GM, Balandin, SV, Krasnosdembskaya, AD, Markelov, ML, et al. (2004). Purification and primary structure of two isoforms of arenicin, a novel antimicrobial peptide from marine polychaeta Arenicola marina. FEBS letters 577: 209-214
    View at Publisher | View at Google Scholar
  47. Ovchinnikova, TV, Balandin, SV, Aleshina, GM, Tagaev, AA, Leonova, YF, et al. (2006). Aurelin, a novel antimicrobial peptide from jellyfish Aurelia aurita with structural features of defensins and channel-blocking toxins. Biochemical and biophysical research communications 348: 514-523
    View at Publisher | View at Google Scholar
  48. Pangestuti, R and Kim, S-K (2017). Bioactive Peptide of Marine Origin for the Prevention and Treatment of Non-Communicable Diseases. Marine drugs 15: 67
    View at Publisher | View at Google Scholar
  49. Panteleev, PV, Myshkin, MY, Shenkarev, ZO and Ovchinnikova, TV (2017). Dimerization of the antimicrobial peptide arenicin plays a key role in the cytotoxicity but not in the antibacterial activity. Biochemical and biophysical research communications 482: 1320-1326
    View at Publisher | View at Google Scholar
  50. Panteleev, PV, Tsarev, AV, Safronova, VN, Reznikova, OV, Bolosov, IA, et al. (2020). Structure Elucidation and Functional Studies of a Novel β-hairpin Antimicrobial Peptide from the Marine Polychaeta Capitella teleta. Marine drugs 18
    View at Publisher | View at Google Scholar
  51. Park, C, Cho, J, Lee, J and Lee, DG (2011). Membranolytic antifungal activity of arenicin-1 requires the N-terminal tryptophan and the beta-turn arginine. Biotechnology letters 33: 185-189
    View at Publisher | View at Google Scholar
  52. Phyo, YZ, Ribeiro, J, Fernandes, C, Kijjoa, A and Pinto, MMM (2018). Marine Natural Peptides: Determination of Absolute Configuration Using Liquid Chromatography Methods and Evaluation of Bioactivities. Molecules (Basel, Switzerland) 23: 306
    View at Publisher | View at Google Scholar
  53. Raihan, T, Rabbee, MF, Roy, P, Choudhury, S, Baek, K-H et al. (2021). Microbial Metabolites: The Emerging Hotspot of Antiviral Compounds as Potential Candidates to Avert Viral Pandemic Alike COVID-19. Frontiers in molecular biosciences 8: 732256-732256
    View at Publisher | View at Google Scholar
  54. Rawat, DS, Joshi, MC, Joshi, P and Atheaya, H (2006). Marine Peptides and Related Compounds in Clinical Trial+. Anti-Cancer Agents in Medicinal Chemistry 6: 33-40
    View at Publisher | View at Google Scholar
  55. Riccio, G, Ruocco, N, Mutalipassi, M, Costantini, M, Zupo, V, et al. (2020). Ten-Year Research Update Review: Antiviral Activities from Marine Organisms. Biomolecules 10: 1007
    View at Publisher | View at Google Scholar
  56. Rose, R, Golosova, O, Sukhomlinov, D, Tiunov, A and Prosperi, M (2018). Flexible design of multiple metagenomics classification pipelines with UGENE. Bioinformatics 35: 1963-1965
    View at Publisher | View at Google Scholar
  57. Saint Jean, KD, Henderson, KD, Chrom, CL, Abiuso, LE, Renn, LM et al. (2017). Effects of Hydrophobic Amino Acid Substitutions on Antimicrobial Peptide Behavior. Probiotics and Antimicrobial Proteins 10: 408-419
    View at Publisher | View at Google Scholar
  58. Saito, T, Kawabata, S, Shigenaga, T, Takayenoki, Y, Cho, J, et al. (1995). A novel big defensin identified in horseshoe crab hemocytes: isolation, amino acid sequence, and antibacterial activity. J Biochem 117: 1131-1137
    View at Publisher | View at Google Scholar
  59. Serna, N, Carratalá, JV, Conchillo-Solé, O, Martínez-Torró, C, Unzueta, U, et al. (2021). Antibacterial Activity of T22, a Specific Peptidic Ligand of the Tumoral Marker CXCR4. Pharmaceutics 13
    View at Publisher | View at Google Scholar
  60. Shai, Y (1994). Pardaxin: channel formation by a shark repellant peptide from fish. Toxicology 87: 109-129
    View at Publisher | View at Google Scholar
  61. Shai, Y, Fox, J, Caratsch, C, Shih, YL, Edwards, C et al. (1988). Sequencing and synthesis of pardaxin, a polypeptide from the Red Sea Moses sole with ionophore activity. FEBS Letters 242: 161–166
    View at Publisher | View at Google Scholar
  62. Shenkarev, ZO, Panteleev, PV, Balandin, SV, Gizatullina, AK, Altukhov, DA, et al. (2012). Recombinant expression and solution structure of antimicrobial peptide aurelin from jellyfish Aurelia aurita. Biochemical and biophysical research communications 429: 63-69
    View at Publisher | View at Google Scholar
  63. Shimizu, K, Muraoka, Y, Hirose, S, Tomii, K and Noguchi, T (2007). Predicting mostly disordered proteins by using structure-unknown protein data. . BMC Bioinformatics 8: 78
    View at Publisher | View at Google Scholar
  64. Shin, SH, Lee, YS, Shin, YP, Kim, B, Kim, MH, et al. (2013). Therapeutic efficacy of halocidin-derived peptide HG1 in a mouse model of Candida albicans oral infection. Journal of Antimicrobial Chemotherapy 68: 1152-1160
    View at Publisher | View at Google Scholar
  65. Shin, YP, Park, HJ, Shin, SH, Lee, YS, Park, S, et al. (2010). Antimicrobial activity of a halocidin-derived peptide resistant to attacks by proteases. Antimicrob Agents Chemother 54: 2855-2866
    View at Publisher | View at Google Scholar
  66. Silva, ON, Fensterseifer, IC, Rodrigues, EA, Holanda, HH, Novaes, NR, et al. (2015). Clavanin A improves outcome of complications from different bacterial infections. Antimicrob Agents Chemother 59: 1620-1626
    View at Publisher | View at Google Scholar
  67. Simonetti, O, Ganzetti, G, Arzeni, D, Campanati, A, Marconi, B, et al. (2009). In vitro activity of Tachyplesin III alone and in combination with terbinafine against clinical isolates of dermatophytes. Peptides 30: 1794-1797
    View at Publisher | View at Google Scholar
  68. Tziveleka, LA, Vagias, C and Roussis, V (2003). Natural products with anti-HIV activity from marine organisms. Curr Top Med Chem 3: 1512-1535
    View at Publisher | View at Google Scholar
  69. Umnyakova, ES, Gorbunov, NP, Zhakhov, AV, Krenev, IA, Ovchinnikova, TV, et al. (2018). Modulation of Human Complement System by Antimicrobial Peptide Arenicin-1 from Arenicola marina. Marine drugs 16
    View at Publisher | View at Google Scholar
  70. Vernen, F, Craik, DJ, Lawrence, N and Troeira Henriques, S (2019a). Cyclic Analogues of Horseshoe Crab Peptide Tachyplesin I with Anticancer and Cell Penetrating Properties. ACS Chem Biol 14: 2895-2908
    View at Publisher | View at Google Scholar
  71. Vernen, F, Harvey, PJ, Dias, SA, Veiga, AS, Huang, YH, et al. (2019b). Characterization of Tachyplesin Peptides and Their Cyclized Analogues to Improve Antimicrobial and Anticancer Properties. Int J Mol Sci 20: 4184
    View at Publisher | View at Google Scholar
  72. Ward, JJ, Sodhi, JS, McGuffin, LJ, Buxton, BF and Jones, DT (2004). Prediction and functional analysis of native disorder in proteins from the three kingdoms of life. J Mol Biol 337: 635-645
    View at Publisher | View at Google Scholar
  73. Wen, S, Majerowicz, M, Waring, A and Bringezu, F (2007). Dicynthaurin (ala) monomer interaction with phospholipid bilayers studied by fluorescence leakage and isothermal titration calorimetry. J Phys Chem B 111: 6280-6287
    View at Publisher | View at Google Scholar
  74. Xia, X (2018). DAMBE7: New and Improved Tools for Data Analysis in Molecular Biology and Evolution. Molecular biology and evolution 35: 1550-1552
    View at Publisher | View at Google Scholar
  75. Xie, H, Wei, J and Qin, Q (2016. Antiviral function of Tachyplesin I against iridovirus and nodavirus. Fish Shellfish Immunol 58: 96-102
    View at Publisher | View at Google Scholar
  76. Zhang, Q-T, Liu, Z-D, Wang, Z, Wang, T, Wang, N, et al. (2021). Recent Advances in Small Peptides of Marine Origin in Cancer Therapy. Marine drugs 19: 115
    View at Publisher | View at Google Scholar
a