Literature DB >> 29949853

Cyanopeptolins with Trypsin and Chymotrypsin Inhibitory Activity from the Cyanobacterium Nostoc edaphicum CCNP1411.

Hanna Mazur-Marzec1,2, Anna Fidor3, Marta Cegłowska4, Ewa Wieczerzak5, Magdalena Kropidłowska6, Marie Goua7, Jenny Macaskill8, Christine Edwards9.   

Abstract

Cyanopeptolins (CPs) are one of the most frequently occurring cyanobacterial peptides, many of which are inhibitors of serine proteases. Some CP variants are also acutely toxic to aquatic organisms, especially small crustaceans. In this study, thirteen CPs, including twelve new variants, were detected in the cyanobacterium Nostoc edaphicum CCNP1411 isolated from the Gulf of Gdańsk (southern Baltic Sea). Structural elucidation was performed by tandem mass spectrometry with verification by NMR for CP962 and CP985. Trypsin and chymotrypsin inhibition assays confirmed the significance of the residue adjacent to 3-amino-6-hydroxy-2-piperidone (Ahp) for the activity of the peptides. Arginine-containing CPs (CPs-Arg²) inhibited trypsin at low IC50 values (0.24⁻0.26 µM) and showed mild activity against chymotrypsin (IC50 3.1⁻3.8 µM), while tyrosine-containing CPs (CPs-Tyr²) were selectively and potently active against chymotrypsin (IC50 0.26 µM). No degradation of the peptides was observed during the enzyme assays. Neither of the CPs were active against thrombin, elastase or protein phosphatase 1. Two CPs (CP962 and CP985) had no cytotoxic effects on MCF-7 breast cancer cells. Strong and selective activity of the new cyanopeptolin variants makes them potential candidates for the development of drugs against metabolic disorders and other diseases.

Entities:  

Keywords:  Nostoc; cyanobacteria; cyanopeptolins; protease inhibitors

Mesh:

Substances:

Year:  2018        PMID: 29949853      PMCID: PMC6070996          DOI: 10.3390/md16070220

Source DB:  PubMed          Journal:  Mar Drugs        ISSN: 1660-3397            Impact factor:   5.118


1. Introduction

Cyanobacterial peptides belong to the most interesting group of natural bioactive products. Initially, they were recognized as hepatotoxic and inflammatory agents. In the late 1970s, the therapeutic potential of the compounds attracted the attention of the scientific community [1,2]. Since then, numerous and structurally diverse cyanopeptides have been identified, including potent anticancer agents such as dolastatin 10, cryptophycin 52, largazol, and apratoxin [3,4,5,6,7]. Cyanopeptides have also been recognized as potent inhibitors of key metabolic enzymes, targeting mainly serine proteases and protein phosphatases [7,8,9]. Among these compounds, cyanopeptolins (CPs), a large family of cyclic depsipeptides (peptidolactones), are commonly produced by different cyanobacterial genera, including Microcystis, Planktothrix, Anabaena and Nostoc. These nonribosomal peptides (NRPs) are composed of a six-amino acid ring and a side chain with one or two residues. All CPs are characterized by the presence of 3-amino-6-hydroxy-2-piperidone (Ahp) in position 3 (Figure 1 and Figure 2). In a few CP-type peptides, the occurrence of O-methylated Ahp (Amp) in this position has been reported [10,11,12]. Ahp is also present in other cyclodepsipeptides such as aeruginopeptins, micropeptins, microcystilide, nostopeptins, and oscillapeptins. Position 1 in CP-type peptides is conserved and occupied by l-threonine with β-hydroxy group linked by an ester bond to carboxy terminus of amino acid in position 6 (i.e., Val, Ile or allo-Ile). Position 2 is most variable and occupied by residues differing in structure and polarity (e.g., Arg, Leu, Gln, Tyr, Phe, MeTyr, H4Tyr, Dhb) [13]. In position 4, Leu/Ile, Phe or Thr can be found, while in position 5, N-methylated aromatic amino acids, N-MeTyr or MePhe, or their homo-variants are usually present. In some CPs, Tyr or N-MeTyr were modified by chlorination [14]. The side chain is attached via the amino group on Thr. Usually, two major types of side chains were reported: with 1–2 amino acid residues (e.g., Asp, Glu) and aliphatic fatty acid of variable length or with glyceric acid directly linked to Thr or to an amino acid side chain [15,16,17]. Glyceric acid can be modified by sulfation or/and O-methylation [11,18]. In the side chain of aeruginopeptins, micropeptins or microcystilide, a hydroxyphenyl lactic acid (Hpla) was reported [19,20,21].
Figure 1

COSY, ROESY, and HMBC correlations in cyanopeptolin CP985.

Figure 2

Key COSY, ROESY, and HMBC correlations in cyanopeptolin CP962.

Cyanopeptolins, like other nonribosomal peptides, are synthesized on large multi-enzyme complexes with modular structure. The general organization of the gene clusters encoding the enzymes in different cyanobacteria is similar, however, some differences in the specificity of adenylation domains and in the presence of tailoring domains exist [22,23,24]. In the CP gene cluster from Microcystis and Anabaena, a halogenase gene was present [23,24], whereas in Planktothrix a glyceric acid-activating domain and sulfotransferase domain occurred [22]. The modifications in gene clusters and differences in substrate specificity of adenylation domains result in intra- and interspecies diversity of CP structures. The majority of CPs showed inhibitory activity against serine proteases, such as trypsin, chymotrypsin, thrombin, and elastase (e.g., [12,25,26,27,28,29,30]). Cyanopeptolins with one (CP S) or two sulphate groups (CP SS) also inhibited plasmin [18]. The activity of the peptides was found to be determined by the residue in position 2, however, the significance of other structural elements was also reported [12]. In ichthyopeptins, CP analogues with 2-hydroxy-3-(4′-hydroxyphenyl)lactic acid (PAA) in the side chain, strong antiviral activity against influenza A virus was observed [31]. Tests on small crustaceans revealed the harmful effects of Ahp-containing cyclic depsipeptides [16,32,33]. For CP SS, the toxicity against Daphnia magna was even higher than for microcystin-LR [10], the most widely studied cyanobacterial toxin. In cyanobacterial strains from the Nostoc genus, typical CP variants produced by Microcystis have not been reported. However, several other CP-type structures, namely nostopeptins, insulapeptolides, and nostocyclins were identified (Table 1) [34,35,36,37,38]. Nostopeptin A and B from N. minutum NIES-26, with 3-hydroxy-4-methylproline (Hmp) in position 1, showed inhibitory activity against elastase and chymotrypsin, but were inactive against papain, trypsin, thrombin, and plasmin [35]. Insulapeptolides A–D from N. insulare are characterized by the presence of Hmp in position 1 and citrulline (Cit) in the side chain. Extracts containing these peptides potently and selectively inhibited human leukocyte elastase (HLE) [38]. Nostocyclin from Nostoc sp. DUN901 has d-Hpla in the side chain and two homoserine residues (Hse): one in a ring part and one in a side chain of the molecule [34]. The peptide was not toxic in mouse bioassay, but showed weak activity against protein phosphatases [34,39].
Table 1

Cyanopeptolin-type peptides identified in cyanobacteria from Nostoc genus.

Molecular MassPeptide NameStructureEnzyme InhibitionReferences
921Nostopeptin BN920[Thr+Leu+Ahp 1 +Phe+MeTyr+Val]Gln+Ac 2Chymotrypsin(IC50 0.11 μM)[36]
926Nostopeptin B[Hmp 3 +Leu+Ahp+Ile+MeTyr+Ile]Gln+AcElastase(IC50 11.0 μg/mL)Chymotrypsin(IC50 1.6 μg/mL)[35]
937Nostopeptin A[Hmp+Leu+Ahp+Ile+MeTyr+Ile]Gln+BA 4Elastase(IC50 1.3 μg/mL)Chymotrypsin(IC50 1.4 μg/mL)[35]
942Insulapeptolide A[Hmp+Leu+Ahp+Ile+MeTyr+Val]Cit 5 +AcHLE 10 (IC50 0.14 μM) *[38]
956Insulapeptolide B[Hmp+Leu+Ahp+Leu+MeTyr+Ile]Cit+AcHLE (IC50 0.10 μM) *
956Insulapeptolide C[Hmp+Leu+Ahp+Ile+diMeTyr 6 +Val]Cit+AcHLE (IC50 0.090 μM) *
970Insulapeptolide D[Hmp+Leu+Ahp+Ile+diMeTyr+Ile]Cit+AcHLE (IC50 0.085 μM) *
991Insulapeptolide G[Thr+Hph 7 +Ahp+Thr+MePhe+Val]Ser+Pro+BAHLE (IC50 3.5 μM) *
1005Insulapeptolide H[Thr+Hph+Ahp+Thr+MeTyr+Ile]Ser+Pro+BAHLE (IC50 2.7μM) *
1007Insulapeptolide F[Thr+Hph+Ahp+Thr+MeTyr+Val]Ser+Pro+BAHLE (IC50 1.6 μM) *
1021Insulapeptolide E[Thr+Hph+Ahp+Thr+MeTyr+Ile]Ser+Pro+BAHLE (IC50 3.2 μM) *
1116Nostocyclin[Thr+Hse 8 +Ahp+Phe+MeTyr+Val]Hse+Ile+Hpla 9PP1 (IC50 64.0 μM)[34]

* activity of extracts; 1 Ahp 3-amino-6-hydroxy-2-piperidone; 2 Ac acetic acid; 3 Hmp 3-hydroxy-4-methyl-proline; 4 BA butanoic acid; 5 Cit citrulline; 6 diMeTyr N,O-dimethyltyrosine; 7 Hph homophenylalanine; 8 Hse homoserine; 9 Hpla p-hydroxyphenyllactic acid; 10 HLE human leukocyte elastase.

Among cyanobacterial strains from the same species, significant differences in the peptide profiles are frequently reported. In our study, the structures of CPs produced by N. edaphicum CCNP1411 isolated from coastal waters of the Gulf of Gdańsk, southern Baltic Sea, were elucidated. In total, thirteen CP variants were identified. They represent structures typical of CPs from Microcystis, but different from CP-type peptides previously found in other Nostoc strains. The biological activity of the peptides against serine proteases, protein phosphatase 1, and MCF-7 breast cancer cells were assessed.

2. Results

2.1. LC-MS/MS Analysis of Cyanopeptolins

Fractionation of N. edaphicum CCNP1411 crude extract (Figure S1) resulted in isolation of thirteen CPs. Structures were identified using a quadrupole/time of flight mass spectrometer and a triple quadrupole/linear ion trap mass spectrometer (Table 2). Structural elucidation of the peptides was based on fragmentation spectra with diagnostic ions, including immonium ions and a series of other fragment ions associated with specific residues. Depending on the residue in position 2, two types of spectra were obtained. Arg2-containing CPs (CPs-Arg2), gave pseudomolecular ions [M + H]+ at m/z 1049, 1021, 1019, 1007, 979, 993, 991, and 963. The Tyr2-containing peptides (CPs-Tyr2) were detected as dehydrated protonated molecules [M + H − H2O]+ at m/z 1010, 996, 982, 968, and 952, and the Tyr-immonium ion (m/z 136) was always present in their spectra. The putative planar structures of CPs detected in N. edaphicum CCNP1411 and their fragmentation spectra are presented in Figure 1, Figure 2, Figure 3 and Figure 4 and in supplementary information (Figures S2–S12). Amino acids at positions 1, 3, 4, 6, and 7 were found to be conserved and occupied by Thr1, Ahp3, Phe4, Val6, and Asp7, respectively. The ion peak corresponding to the longest sequence of residues common to all CP variants was observed in the spectra at m/z 297 [Asp + Thr + Val + H − H2O]+ or/and at m/z 269 [Asp + Thr + Val + H − H2OCO]+. The presence of butanoic acid (BA), hexanoic acid (HA), or octanoic acid (OA) in the side chain was mainly indicated by ion peaks formed by the cleavage of the corresponding fatty acid group (FA) and the exocyclic aspartic acid (Figure 3 and Figure 4; Figures S2–S12). As this cleavage produced a stable cyclic part of the molecule, the ions [M + H − (H2O) − (FA + Asp)]+ usually belonged to the most abundant ones. The residue in position 5 (i.e., N-MePhe, N-MeTyr or N-MeHty) was identified based on immonium ion peaks at m/z 134, 150 or 164, respectively, and peaks at m/z 404, 420, and 434 corresponding to [Ahp + Phe + (MePhe/MeTyr/MeHty) + H − H2O]+. Ion peak at m/z 120, as well as peaks at 243 [Ahp + Phe + H − H2O]+ and 215 [Ahp + Phe + H − H2OCO]+ confirmed the presence of Phe in position 4.
Table 2

Cyanopeptolins identified in Nostoc edaphicum CCNP 1411. The activities of the peptides were assessed in serine proteases (trypsin, chymotrypsin, elastase, and thrombin) and protein phosphatases 1 (PP 1) inhibition assays (- not active; */** small/medium activity; m/z of precursor ions: [M + H]+ for CPs-Arg2 and [M + H − H2O]+ for CPs-Tyr2.

Cyanopeptolin CPm/zStructureEnzyme Inhibition (IC50 [µM])
ChymotrypsinTrypsinElastaseThrombinProtein Phosphatase 1
CP 10481049[Thr+Arg+Ahp+Phe+MeHty+Val]Asp+OA***---
CP 10201021[Thr+Arg+Ahp+Phe+MeHty+Val Asp+HA3.10.25---
CP 10181019[Thr+Arg+Ahp+Phe+MePhe+Val]Asp+ OA-0.24---
CP 10061007[Thr+Arg+Ahp+Phe+MeTyr+Val]Asp+HA***---
CP 992993[Thr+Arg+Ahp+Phe+MeHty+Val]Asp+BA3.50.24---
CP 990991[Thr+Arg+Ahp+Phe+MePhe+Val]Asp+HA-**---
CP 978979[Thr+Arg+Ahp+Phe+MeTyr+Val]Asp+BA3.80.26---
CP 962963[Thr+Arg+Ahp+Phe+MePhe+Val]Asp+BA-**---
CP 10271010[Thr+Tyr+Ahp+Phe+MeHty+Val]Asp+HA0.26----
CP 1013996[Thr+Tyr+Ahp+Phe+MeTyr+Val]Asp+HA**----
CP 999982[Thr+Tyr+Ahp+Phe+MeHty+Val]Asp+BA**--.--
CP 985968[Thr+Tyr+Ahp+Phe+MeTyr+Val]Asp+BA0.26----
CP 969952[Thr+Tyr+Ahp+Phe+MePhe+Val]Asp+BA**----
Figure 3

The product ion mass spectra of CP985 [Thr+Tyr+Ahp+Phe+MeTyr+Va]Asp+BA with precursor ion [M + H − H2O] at m/z 968. The spectra were recorded with application of a hybride quadrupole/time-of-flight mass spectrometer (QTOF) (A) and a hybride triple quadrupole/linear ion trap mass spectrometer (QTRAP) (B). The mass signals were assigned to the following fragments: 950 [M + H − 2H2O]+, 869 [M + H − Val − H2O]+, 851 [M + H − Val − 2H2O]+, 773 [M + 2H − (BA+Asp) − CO]+, 765 [M + 2H − (BA + Asp) − 2H2O]+, 692 [M + H − (Val + MeTyr) − H2O]+, 674 [M+H−(Val+MeTyr)−2H2O]+, 646 [M + H − (Val + MeTyr) − 2H2O − CO]+, 432 [M + H − (Val + MeTyr + Phe + Ahp) − H2O]+, 420 [Ahp + Phe + MeTyr + H − H2O]+, 404 [M + H − (Val + MeTyr + Phe + Ahp) − H2O − CO]+, 386 [BA + Asp + Thr + Val + H]+, 308 [Phe(−N) + MeTyr + H]+, 297 [Asp + Thr + Val + H − H2O − CO]+, 243 [Ahp + Phe + H − H2O]+, 215 [Ahp + Phe + H − H2O − CO]+, 150 MeTyr immonium ion, 136 Tyr immonium ion, 120 Phe immonium ion.

Figure 4

The product ion mass spectra of CP962 [Thr+Arg+Ahp+Phe+MePhe+Val]Asp+BA with precursor ion [M + H]+ at m/z 963. The spectra were recorded with application of a hybride quadrupole/time-of-flight mass spectrometer (QTOF) (A) and a hybride triple quadrupole/linear ion trap mass spectrometer (QTRAP) (B). The mass signals were assigned to the following fragments: 945 [M + H − H2O]+, 917 [M + H − H2O − CO]+, 846 [M + H − Val − H2O]+, 828 [M + H − Val − 2H2O]+, 778 [M + 2H − (BA + Asp)]+, 760 [M + 2H − (BA + Asp) − H2O]+, 742 [M + 2H − (BA + Asp) − 2H2O]+, 685 [M + H − (Val + MePhe) − H2O]+, 425 [BA + Asp + Thr + Arg + H − H2O]+, 404 [Ahp + Phe + MePhe + H − H2O]+, 297 [Asp + Thr + Val + H − H2O]+, 243 [Ahp + Phe + H − H2O]+, 215 [Ahp + Phe + H − H2O − CO]+, 134 MePhe immonium ion, 120 Phe immonium ion, 70-Arg.

2.2. NMR Analysis

In order to confirm the structures, two CPs (i.e., CP962 and CP985) were purified in sufficient quantities for NMR spectroscopy. For both compounds the 1H-NMR spectra displayed the typical pattern of a peptide (i.e., doublet amide protons (δH 6.95–8.51 ppm) and a single amide methyl group (δH 2.71 ppm for CP985 and 2.84 ppm for CP962). The COSY, TOCSY, and HMBC experiments allowed assignment of NMR spin systems to Asp, Thr, Tyr, Ahp (3-amino-6-hydroxypiperid-2-one), Phe, MeTyr (N-methyl tyrosine), Val, and butanoic acid (BA) in the case of cyanopeptolin CP985 (Table 3, Figures S13–S18b). The presence of aromatic amino acid residues was recognized by the signals occurring in the aromatic region of the spectrum (δH 6.5–7.5 ppm). The AA’BB’ spin systems between two sets of tyrosine (Tyr-H2′/6′ and Tyr-H3′/5′, JH,H = 8.4 Hz) and N-methyl tyrosine aromatic protons (MeTyr-H2′/6′ and MeTyr-H3′/5′, JH,H = 8.4 Hz) indicated the presence of two para-di-substituted phenyl rings. The 1H-13C long range correlation from MeTyr-NH-CH3 group (δH 2.71 ppm) to the MeTyr-C2 atom (δC 61.3 ppm) revealed the presence of N-methyl tyrosine residue. Phenylalanine was found to be the third aromatic amino acid residue based on the COSY interaction between Phe-H2′/6′, Phe-H3′/5′, and Phe-H4′, and the HMBC correlation from two diastereotopic methylene protons Phe-3a (δH 2.88 ppm) and Phe-3b (δH 1.81 ppm) to the aromatic Phe-C2′/6′ carbons. The presence of Asp, Val, and Thr was confirmed by their characteristic spin systems in the COSY spectrum (Figure S17b). The macrocyclic ring closure between threonine and valine was verified by HMBC correlation between Thr-H3 (δH 5.36 ppm) and Val-C1 (δC 172.4 ppm) (Figure S15b). The presence of Ahp residue was detected by the characteristic signal of the OH proton (δH 5.99 ppm) and a broad singlet (δH 5.06 ppm) derived from H5 proton (Figure S13). The HMBC correlations from Ahp-C1 to Ahp-H5 and Ahp-H2 confirmed the cyclic nature of this residue.
Table 3

NMR Spectroscopic data (700 MHz, DMSO-d6) for cyanopeptolin CP985 [Thr+Tyr+Ahp+Phe+MeTyr+Val]Asp+BA.

UnitPositionδCδH (J in Hz)ROESYHMBC a
BA 1234172.9, C37.7, CH219.2, CH214.2, CH32.11, m1.56, m0.89, t (7.2, 7.2)Asp-NHBA-1, BA-3, BA-4BA-1, BA-2, BA-4BA-2
Asp 123a3b4NH171.7, C49.9, CH36.1, CH2172.4, C4.58, dd (8.0, 5.0)2.65, m2.46, m8.26, mThr-NHAsp-4
Thr1 1234NH168.9, C54.7, CH72.6, CH18.0, CH34.52, d (10.0)5.36, q (6.8, 6.9, 6.8)1.16, d (7.2)7.52, d (9.3)Tyr-NHTyr-NHAsp-2Thr-1, Asp-1Val-1Thr-2, Thr-3
Tyr2 123a3b1’2′/6′3′/5′4′NH169.9, C54.2, CH35.5, CH2128.5, C130.1, CH115.5, CH156.1, C4.31, m3.11, m2.55, m6.89, d (8.4)6.58, d (8.4)8.46, d (8.8)Ahp-NHThr-2, Thr-3Tyr-2, Tyr-1′Tyr-1′Tyr-4′Thr-1
Ahp3 12345NHOH169.0, C49.4, CH22.0, CH221.9, CH274.1, CH3.62, m2.41, m1.64, m5.06, brs7.06, d (8.8)5.99, d (3.1)Ahp-NHPhe-3a, Phe-3bTyr-2, Ahp-3Ahp-1, Ahp-3Ahp-2Ahp-1Tyr-1
Phe4 123a3b1′2′/6′3′/5′4′170.8, C50.7, CH35.8, CH2137.2, C129.9, CH128.2, CH126.7, CH4.76, dd (7.1, 4.6)2.88, t (12.9, 12.9)1.81, dd (10.4, 3.9)6.84, d (7.0)7.19, t (7.3, 7.3)7.14, d (7.0)Ahp-5Phe-2, Phe-1′, Phe-2′/6′Phe-3a, Phe-3b, Phe-4′Phe-1′Phe C2′/6′
MeTyr5 1231′2′/6′3′/5′4′NCH3OH169.4, C61.3, CH33.3, CH2128.2, C130.9, CH115.5, CH156.7, C33.3, CH34.89, dd (8.9, 2.5)3.10, m7.00, d (8.4)6.78, d (8.4)2.71, s9.33, sVal-NHMeTyr-1′MeTyr-1′, MeTyr-3′/5′MeTyr-2′/6′, MeTyr-4′MeTyr-2, Phe-1MeTyr-3′/5′, MeTyr-4′
Val6 12345NH172.4, C56.3, CH31.4, CH19.7, CH317.7, CH34.63, dd (4.9, 4.6)2.02, m0.84, d (6.6)0.71, d (6.6)7.39, d (9.7)MeTyr-2Val-2, Val-3, Val-5Val-2, Val-3, Val-4MeTyr-1

HMBC correlations are given from proton(s) stated to the indicated carbon atom.

The COSY, TOCSY, and HMBC data allowed identification of amino acid residues in CP CP962 as Asp, Thr, Arg, Ahp, Phe, MePhe (N-methyl phenylalanine), Val, and BA, analogously to CP985 (Table 4, Figures S19–S24b). Two aromatic residues were found: phenylalanine and N-methyl phenylalanine. The occurrence of N-methyl group was established by the 1H-13C long range correlation from MePhe-NH-CH3 group (δH 2.84 ppm) to the MePhe-C2 atom. The typical 13C chemical shift of the guanidine quaternary carbon (δC 158.7 ppm) indicated the presence of Arg whose complete spin system was assigned based on TOCSY and COSY interactions (Figure S23b). The diagnostic regions of the TOCSY, ROESY and HMBC spectra of both cyanopeptolins analyzed are presented in corresponding figures in Supplementary.
Table 4

NMR Spectroscopic data (700 MHz, DMSO-d6) for cyanopeptolin CPL962 [Thr+Arg+Ahp+Phe+MePhe+Val]Asp+BA.

UnitPositionδCδH (J in Hz)ROESYHMBC a
BA 1234172.2, C37.4, CH219.3, CH214.1, CH32.08, m1.52, m0.87, t (7.4, 7.4)Asp-NHBA-1, BA-3, BA-4BA-1, BA-2, BA-4BA-2
Asp 123a3b4NH173.7, C50.8, CH40.1, CH2175.3, C4.54, m2.51, m2.12, m8.00, d (8.0)Thr-NHAsp-2, Asp-4BA-1
Thr1 1234NH169.9, C55.0, CH72.3, CH17.9, CH34.58, d (9.0)5.30, q (6.8, 6.5, 6.8)1.15, d (6.5)7.23, d (7.6)Arg-NHArg-NHAsp-2Thr-1Thr-4, Val-1Thr-2, Thr-3
Arg2 123456NH170.4, C49.0, CH26.4, CH224.5, CH239.8, CH2158.7, C4.15, m1.88, m1.46, m2.95, m8.5, d (8.7)Ahp-NHThr-2, Thr-3Thr-1
Ahp3 12345NHOH169.6, C48.6, CH21.7, CH222.0, CH274.2, CH3.65, m2.43, m1.68, m5.03, brs6.95, d (9.7)6.01, d (2.6)Phe-3a, Phe-3bArg-2, Ahp-3 Ahp-1, Ahp-3Ahp-1Arg-1
Phe4 123a3b1′2′/6′3′/5′4′170.0, C50.6, CH35.6, CH2137.1, C129.8, CH128.2, CH126.7, CH4.74, dd (7.2, 3.5)2.86, m1.69, m6.78, d (7.1)7.17, t (7.1, 7.1)7.13, d (7.8)Ahp-5Phe-1, Phe-3a, Phe-3bPhe-1′Phe-3a, Phe-3b, Phe-4′Phe-1′Phe-3′/5′
MePhe5 1231′2′/6′3′/5′4′NCH3169.4, C60.9, CH34.4, CH2138.4, C130.0, CH129.1, CH127.2, CH35.6, CH35.02, m3.23, m7.24, d (7.6)7.41, t (7.7, 7.7)7.32, d (7.5)2.84, sVal-NHMePhe-1′MePhe-1′MePhe-4′MePhe-3′/5′MePhe-2
Val6 12345NH172.6, C56.4, CH31.6, CH19.7, CH317.8, CH34.71, dd (5.2, 4.2)2.05, m0.88, d (6.9)0.76, d (6.9)7.41, d (7.7)MePhe-2Val-3Val-2, Val-3, Val-5Val-2, Val-3, Val-4MePhe-1

HMBC correlations are given from proton(s) stated to the indicated carbon atom.

2.3. Bioassays

Preparative chromatography resulted in separation of several fractions containing thirteen pure cyanopeptolin variants (Table 2). These peptides were evaluated for inhibition against four serine proteases and protein phosphatase 1. In addition, their effect on MCF-7 breast cancer cells was tested. However, only six CPs were isolated from N. edaphicum CCNP1411 in sufficient amounts to obtain quantitative results of the assays. None of the peptides were active against thrombin, elastase, and protein phosphatase 1, even at the highest concentration used in the study (45.4 μg/mL). However, all CPs with Arg in position 2 significantly reduced the activity of trypsin. The IC50 values of trypsin inhibitors were comparable and in the range of 0.24–0.26 μM. The CPs-Arg2 with N-MeTyr or N-MeHty in position 5 were also active against chymotrypsin, but the IC50 values were lower (IC50 = 3.1–3.8 μM). Chymotrypsin inhibition activity of CPs-Arg2 was not observed when position 5 was occupied by MePhe. All CPs-Tyr2 reduced activity of chymotrypsin and were inactive against other enzymes. For CP1027 and CP985, the IC50 was 0.26 μM. Following enzyme inhibition assays, all samples with Arg2-containing CP962 and Tyr2-containing CP985 (at 4.54 μg/mL) were analyzed by LC-MS/MS. The recovery was based on extracted mass chromatogram of parent ions. The enzymes did not cause any significant loss of the peptides; their contents were in the range from 94.6% to 97.7% of that in samples without the enzyme. Due to limited amounts of pure peptides isolated from N. edaphicum CCNP1411, only Arg2-containing CP962 and Tyr2-containing CP985 were used in MTT (3-(4,5-dimethylthiazole-2-yl)-2,5-diphenyltetrazolium bromide) assay. After 24-h exposure, no cytotoxic effects on the MCF-7 cells were observed (Figure 5).
Figure 5

Viability of breast cancer cells MCF-7 treated for 24 h with Arg2-containing CP962 (white bars) and Tyr2-containing CP985 (black bars), isolated from N. edaphicum CCNP1411. Data are presented as percentage of the control, mean ± s.e.m, n = 3.

3. Discussion

Cyanobacteria from the Nostoc genus, especially the symbiotic strains and those living in waters from tropical regions, are important producers of bioactive compounds with potential biotechnological or pharmaceutical application [37,40,41,42,43]. One of the most prominent examples are cryptophycins, cyclic depsipeptides isolated from Nostoc sp. ATCC 53789 and GSV 224. They are promising candidates for anti-cancer drug development [44,45]. Another important Nostoc metabolite is cyanovirin-N, a small cyanobacterial lectin which blocks the entry of the enveloped viruses such as HIV, influenza, and Ebola [46]. In this study, the structures of CPs produced by N. edaphicum CCNP1411 from the Baltic Sea were elucidated based on the fragmentation spectra of their pseudomolecular ions (CPs-Arg2) (Figure 4; Figures S2–S8) or the dehydrated forms of the ions (CPs-Tyr2) (Figure 3; Figures S9–S12). The occurrence of the dehydrated pseudomolecular ion as a precursor ion in mass fragmentation spectra of Tyr2-(and Ile2)-containing aeruginopeptins has been previously reported [20,47]. It was suggested that the cleavage of an ester bond at Thr1 and subsequent dehydration resulted in the generation of linear peptides. For the two CPs, Arg2-containing CP962 and Tyr2-containing CP985, the composition and sequences of amino acids were confirmed by NMR spectroscopy. The results of the analyses were consistent with structure elucidation performed by tandem mass spectrometry. The core structures of these peptides were more similar to those identified in Microcystis [13,14,28,32], compared to those reported from other Nostoc isolates (Table 1). For example, instead of Arg2 or Tyr2 present in position 2 of the CPs identified in this work, nostopeptin and insulapeptolides from N. minutum and N. insulare possess Leu2 or Hph2 (homophenylalanine) [35,38]. Of the CPs identified in this study, only CP1006A (m/z 1007) was previously reported. This peptide, along with its chlorinated derivative (CP1040A) were found in Microcystis bloom and culture samples [14,48]. The structural diversity of CPs, specifically the residue in position 2, was found to have significant effect on activity of the peptides against serine proteases: trypsin and chymotrypsin [28,29,32]. The CP-type peptides with potent inhibitory activity against trypsin were characterized by the presence of basic amino acid (Arg or Lys), whereas in peptides active against chymotrypsin, position 2 was occupied by hydrophobic residues (Tyr, Phe, Hty or Leu). The same structure-activity relationship was observed in this study. Eight CPs-Arg2 inhibited the activity of trypsin with IC50 values of 0.24–0.26 µM, and five CPs-Tyr2 inhibited chymotrypsin with similar potency. In addition, the inhibitory activity of CPs produced by CCNP1411 seemed to be affected by the residue in position 5. Only those CPs-Arg2 which had N-MeTyr5 or N-MeHty5 in this position were active against chymotrypsin. However, to unequivocally prove the significance of this structure-activity relationship, more CPs should be tested. The fact that some CPs-Arg2, apart from strong inhibition of trypsin, are also active against chymotrypsin was previously documented by other authors [12,29,49]. The trypsin inhibition activity of CP-type peptides was suggested to be enhanced by the presence of isoleucine in position 6, instead of valine. Other modifications in the structure of Ahp-containing cyclic depsipeptides, such as the presence of chloride or sulfide groups, may also have an effect on enzyme inhibition activity [12]. Cyanobacteria produce many other Ahp-containing cyclic depsipeptides with inhibitory activity against proteases [28,29]. The majority of the peptides were active at micromolar concentrations, but some had even lower IC50 values. Symplocamide A, with citrulline in position 2, showed potent activity against chymotrypsin with IC50 of 0.38 µM, and was 200-times less active against trypsin (IC50 80.2 µM) [29]. Symplocamide A also had cytotoxic activity to NCI H460 lung cancer cells and neuro-2A neuroblastoma cells [29]. Chymotrypsin was most potently inhibited by the glyceric acid 3′-O-phosphate-containing micropeptin T20 from M. aeruginosa, characterized by the presence of the Thr-Phe-Ahp sequence (IC50 2.5 nM) [50]. Picomolar inhibition of trypsin was documented for CP 1020, which was also active against chymotrypsin, plasmin, human kallikrein, and factor XIa [32]. The structure of CP1020 differs from CP978 isolated in our study only in the presence of Glu in a side chain, instead of Asp. In enzymatic assays, CP978 and other CPs isolated from N. edaphicum CCNP1411 were less active than CP1020, but still belong to the most potent protease inhibitors among this class of compounds [12,28,29,32]. Ahp-containing depsipeptides were suggested to block the active center of trypsin or/and chymotrypsin, so the enzymes cannot cleave the peptide bonds at the carboxyl side of Arg or Tyr, respectively [25,28]. In this study, the Tyr2-containing CP985 and Arg2-containing CP962 exposed to proteases inhibited the activity of chymotrypsin or/and trypsin, but their concentrations remained almost unchanged. These results are in line with the hypothesis by Yamaki et al. [28] and confirm the blockage of the active centers of the enzymes by the peptides. Trypsin and chymotrypsin, the two enzymes inhibited by CPs from CCNP1411, are essential for food digestion. Their deregulation can also lead to a number of human diseases such as cancer, cardiovascular and inflammatory diseases. The molecules that modify activity of these proteases, and especially those that act selectively, are widely explored as agents of significant biotechnological and pharmaceutical potential [51]. As the new CP variants identified in N. edaphicum CCNP1411 inhibited the activity of trypsin or chymotrypsin at low concentrations and were inactive against the other tested enzymes and MCF-7 breast cancer cells, their possible use as therapeutic agents should be further explored. Besides its therapeutic potential, this class of protease inhibitors was proven to be an important group of defense agents protecting cyanobacteria from grazers. Oscillapeptin J, CP SS and CP 1020 induced acute effects in small crustaceans [18,32,33]. However, CP SS was not toxic to the isolated rat hepatocytes [18], and oscillapeptin J did not induce any harmful effects in mice when administered intraperitoneally at concentrations up to 1000 µg/kg b.w. [52]. In the case of CP1020, induction of anti-inflammatory effects in human hepatoma cell line Huh7 was observed [53]. In addition, exposure of zebrafish eleuthero-embryos to CP1020 led to transcriptional alterations of genes involved in many important processes, including DNA damage recognition and repair, and circadian rhythm [54]. These findings, along with the acute toxicity of CP-type compounds observed in crustaceans suggest that their activity is not only related to the inhibition of digestive enzymes.

4. Materials and Methods

4.1. Culture Conditions

Nostoc cf. edaphicum CCNP 1411 (GenBank accession number KJ161445) was isolated in 2010 from the Gulf of Gdańsk, southern Baltic Sea, and established as monospecies culture by Dr. Justyna Kobos. Purification of the strain was carried out by multiple transfer to a liquid or/and solid (1% bacterial agar) Z8S medium [55]. To obtain a higher biomass of the cyanobacterium, the culture was grown in 5-L bottles at 22 °C ± 1 °C, at continuous light of 10–30 µE/m2/s. The collected biomass was lyophilized and kept at −20 °C until used.

4.2. Extraction and Isolation

Freeze-dried biomass of N. edaphicum (36 g) was extracted twice with 75% methanol (MeOH) (1 L) by vortexing for 30 min. Following centrifugation at 4,000 g, the extracts were pooled and diluted in water purified using ELGA PURELAB® flex (Veoilia, London, UK) to adjust the concentration of MeOH to <10%. Then, the sample was loaded onto the preconditioned 120-g SNAP KP-C18-HS cartridge (Biotage Uppsala, Sweden). Sample components were eluted using an Isolera flash chromatography system (Biotage Uppsala, Sweden), with a step gradient (10–100% MeOH in water) with 3 column volumes of eluent at each step. Absorbance was monitored at 210 nm and 280 nm. The flow rate was 40 mL/min and 60-mL fractions were collected. Fraction composition was analyzed by UPLC-MS/MS. Fractions containing the same peptides were pooled, diluted with water to <10% MeOH, and concentrated on YMC C18 cartridges (20 mm ID × 2 cm; YMC, Gmbh, Dinslaken, Germany). Six pooled fractions were further purified using preparative HPLC (Biotage Parallex Flex, Cardiff, UK) and Flex V3 software for instrument control and data acquisition. The separation was performed on XBridge Prep C18 column (5 µm CBD, 19 mm ID × 250 mm long; Waters, Elstree, UK) using a 30-min linear gradient from 15% to 80% acetonitrile in MilliQ water with 0.1% formic acid. Absorbance was monitored at 210 nm and 280 nm. The flow rate was 20 mL/min and 4-mL fractions were collected.

4.3. LC-MS/MS Analyses

At each step of the extraction and isolation procedure, the content of the collected fractions was determined by UPLC-MS/MS. The system comprised a Waters Acquity Ultra performance LC coupled to a photodiode array detector (PDA) and a Xevo quadrupole time of flight mass detector (Waters, Elstree, UK). Samples were separated on an ethylene-bridged hybride BEH C18 column (2.1 mm ID × 100 mm; 1.7 µm, Waters) maintained at 40 °C. The mobile phase was Milli-Q water and acetonitrile (solvent B), both containing 0.1% formic acid. Separation was performed using gradient elution (0.3 mL/min) from 20% to 70% B over 10 min, followed by a 100% B wash step and re-equilibration. Data was acquired in positive ion electrospray scanning from m/z 50 to 2000 with a scan time of 2 s and inter-scan delay of 0.1 s. The capillary and cone voltages were set at 0.7 kV and 25 V, respectively. The desolvation gas was maintained at 400 L/h at a temperature of 300 °C. The cone gas was set at 50 L/h with a source temperature of 80 °C. Instrument control, data acquisition and processing were achieved using MassLynx v4.1 (Waters, Milford, MA, USA). Structures of cyanopeptolins were additionally characterized using Agilent 1200 HPLC (Agilent Technologies, Waldbronn, Germany) coupled to a hybrid triple quadrupole/linear ion trap mass spectrometer QTRAP5500 (Applied Biosystems MDS Sciex, Concord, ON, Canada). Peptides were separated on Zorbax Eclipse XDB-C18 column (4.6 mm ID × 150 mm, 5 μm; Agilent Technologies, Santa Clara, CA, USA) column. A gradient elution (0.6 mL/min) was applied with mobile phase composed of 5% acetonitrile in MilliQ water and acetonitrile (solvent B), both containing 0.1% formic acid. The gradient started at 15% B and went to 50% B within 5 min. The content of phase B was then increased to 100% within the next 3 min and kept at that level for 10 min before returning to the starting conditions. The QTRAP MS/MS system was operated in the positive mode, with turbo ion source voltage set at 5.5 kV and temperature at 550 °C. For ions within the m/z range 500–1250 and signal intensity above the threshold of 500,000 cps, fragmentation spectra were acquired within a range 50–1000 Da, at collision energy of 60 V and declustering potential set at 80 eV. Data acquisition and processing were accomplished with the Analyst® Sofware (version 1.5.1, Applied Biosystems, Concord, ON, Canada).

4.4. NMR Analyses

1D 1H-NMR and 2D homo- and heteronuclear 2D NMR (COSY, TOCSY, ROESY, HSQC, and HMBC) were acquired on a Bruker Avance III spectrometers, 500 MHz and 700 MHz. Spectra were recorded in DMSO-d6. NMR data were processed and analyzed by TopSpin (Bruker, Billerica, MA, USA) and SPARKY software (3.114, Goddard and Kneller, freeware https://www.cgl.ucsf.edu/home/sparky).

4.5. Enzyme Inhibition Assay

The chymotrypsin and trypsin inhibition assays were performed following the procedures of Ploutno and Carmeli [41]. The α-chymotrypsin from bovine pancreas (C4129), trypsin from porcine pancreas (T0303), aprotinin (1.5–200 μg/mL) as enzyme inhibitor, N-Suc-Gly-Gly-p-nitroanilide and N-α-benzoyl-dl-arginine-p-nitroanilide hydrochloride (BAPNA) as chymotrypsin and trypsin substrates, respectively were used. Enzymes (0.1 mg/mL) and substrate (2 mM) were dissolved in the same buffer (50 mM Tris-HCl, 100 mM NaCl, 1 mM CaCl2, pH 7.5). The thrombin inhibition assay was performed according to Ocampo and Bennet [56]. Thrombin (T4648) was dissolved in buffer (0.5 mg/mL; 0.2 M Tris-HCl; pH 8.0); substrate (N-p-tosyl-Gly-l-Pro-l-Lys-p-nitroanilide acetate salt, 0.5 mg/mL) and inhibitor (4-(2-aminoethyl)benzenesulfonyl fluoride hydrochloride (AEBSF; 60–2400 μg/mL)) were dissolved in MilliQ water. The elastase inhibition assay was performed according to Kwan et al. [57]. The enzyme (75 μg/mL; E0258) from porcine pancreas, substrate (2 mM; N-Suc-Ala-Ala-Ala-p-nitroanilide) and inhibitor (elastatinal; 5–125 μg/mL) were dissolved in buffer (0.2 M Tris-HCl; pH 8.0). All enzymes, substrates and inhibitors used for proteases inhibition assays were from Sigma-Aldrich (St. Louis, MO, USA). Protein phosphatase 1 inhibition assay was performed according to the procedure described by Rapala et al. [58]. PP1 was from England Biolabs, Hitchin, UK (754S) and the substrate, p-nitrophenyl phosphate disodium salt hexahydrate (p-NPP), was from Sigma-Aldrich (Irvine, UK). Microcystin MC-LR (0.125–4.0 ng/mL) from Enzo Life Sciences, Lausen, Switzerland was used as inhibitor. The enzyme was dissolved in buffer solution A (50 mM Tris at pH 7.4, 1 mg/mL bovine serum albumin (BSA; Sigma -Aldrich, St. Louis, MO, USA), 1 mM MnCl2, 2 mM dithiotreitol (DTT; Sigma-Aldrich, St. Louis, MO, USA)). p-NPP (5.5 mg/mL) was dissolved in buffer solution B (50 mM Tris, pH 8.1, 0.5 mg/mL BSA, 20 mM MgCl2 × 6H2O, 200 mM MnCl2 × 4H2O). All enzyme inhibition assays were performed in a 96 multi-well plate at 37 °C, in triplicate. The absorbance of the reaction mixtures was measured at 405 nm using a microplate reader (Molecular Devices, Sunnyvale, CA, USA). After trypsin and chymotrypsin inhibition assays, the samples with and without the enzymes and containing the highest concentration of CP962 and CP985 (45.4 µg/mL) were analyzed by LC-MS/MS. The content of the peptides was determined based of the peak area of the extracted ions.

4.6. Cytotoxicity Assay

For the test, two CPs were selected: Arg2-containing CP962 and Tyr2-containing CP985. MCF-7 breast cancer cells were seeded at 7.5 × 103 cells/100 μL in a 96-well plates and incubated at 37 °C, 5% CO2 for 24 h. The cells were then treated for a further 24 h with CP962 and CP985 (0 to 500 μg/mL). After 24 h, sterile-filtered 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide solution (MTT; 1 mg/mL) was added to each well. After 4 h incubation at 37 °C in the dark, the MTT solution was removed and formazan crystals solubilized in DMSO. The plates were shaken for 20 min, in the dark, at room temperature and absorbance was measured at 560 nm (Synergy/HT, BIOTEK, Wnooski, VT, USA). For each CP, three independent experiments were carried out and each treatment consisted of six replicates per plate. Bar charts were used to represent the viability of MCF-7 cells treated with CP962 and CP985, compared to the control (i.e., untreated cells) that represented 100% cell viability.

5. Conclusions

Nostoc edaphicum CCNP1411 isolated from the Gulf of Gdańsk (southern Baltic) produces at least thirteen CPs, including twelve variants reported here for the first time. The structures of the peptides are different from other Ahp-containing cyclic depsipeptides previously found in Nostoc. The activity of the peptides was mainly determined by the presence of Arg2 or Tyr2 in Ahp-adjacent position. The fact that trypsin and chymotrypsin did not degrade the tested CPs constitutes an additional evidence for enzyme inactivation by the peptides. Neither of the CP inhibited thrombin, elastase, and protein phosphatase 1; CP962 and CP985 also showed no cytotoxic effects on MCF-7 breast cancer cells. The CPs produced by N. edaphicum CCNP1411, as peptidic structures with selective and potent proteases inhibiting activity, are potential lead compounds in drug discovery process.
  39 in total

1.  Structural analysis of a non-ribosomal halogenated cyclic peptide and its putative operon from Microcystis: implications for evolution of cyanopeptolins.

Authors:  Ave Tooming-Klunderud; Thomas Rohrlack; Kamran Shalchian-Tabrizi; Tom Kristensen; Kjetill S Jakobsen
Journal:  Microbiology       Date:  2007-05       Impact factor: 2.777

Review 2.  A review of pharmacological and toxicological potentials of marine cyanobacterial metabolites.

Authors:  M Nagarajan; V Maruthanayagam; M Sundararaman
Journal:  J Appl Toxicol       Date:  2011-09-09       Impact factor: 3.446

3.  Molecular effects of the cyanobacterial toxin cyanopeptolin (CP1020) occurring in algal blooms: global transcriptome analysis in zebrafish embryos.

Authors:  Susanne Faltermann; Sara Zucchi; Esther Kohler; Judith F Blom; Jakob Pernthaler; Karl Fent
Journal:  Aquat Toxicol       Date:  2014-02-02       Impact factor: 4.964

4.  Two new chymotrypsin inhibitors isolated from the Cyanobacterium Microcystis aeruginosa NIES-88.

Authors:  Hiroshi Yamaki; Namthip Sitachitta; Tomoharu Sano; Kunimitsu Kaya
Journal:  J Nat Prod       Date:  2005-01       Impact factor: 4.050

5.  Cyanopeptoline CB071: a cyclic depsipeptide isolated from the freshwater cyanobacterium Aphanocapsa sp.

Authors:  Hyukjae Choi; Sun Kwan Oh; Wonho Yih; Jungwook Chin; Heonjoong Kang; Jung-Rae Rho
Journal:  Chem Pharm Bull (Tokyo)       Date:  2008-08       Impact factor: 1.645

Review 6.  Cyanobacterial cyclopeptides as lead compounds to novel targeted cancer drugs.

Authors:  Ioannis Sainis; Demosthenes Fokas; Katerina Vareli; Andreas G Tzakos; Valentinos Kounnis; Evangelos Briasoulis
Journal:  Mar Drugs       Date:  2010-03-15       Impact factor: 5.118

7.  New peptolides from the cyanobacterium Nostoc insulare as selective and potent inhibitors of human leukocyte elastase.

Authors:  Christian Mehner; Daniela Müller; Stefan Kehraus; Stephanie Hautmann; Michael Gütschow; Gabriele M König
Journal:  Chembiochem       Date:  2008-11-03       Impact factor: 3.164

8.  New prenylated aeruginosin, microphycin, anabaenopeptin and micropeptin analogues from a Microcystis bloom material collected in Kibbutz Kfar Blum, Israel.

Authors:  Shira Elkobi-Peer; Shmuel Carmeli
Journal:  Mar Drugs       Date:  2015-04-15       Impact factor: 5.118

Review 9.  Plant Protease Inhibitors in Therapeutics-Focus on Cancer Therapy.

Authors:  Sandhya Srikanth; Zhong Chen
Journal:  Front Pharmacol       Date:  2016-12-08       Impact factor: 5.810

10.  Halogenase genes in nonribosomal peptide synthetase gene clusters of Microcystis (cyanobacteria): sporadic distribution and evolution.

Authors:  Sabrina Cadel-Six; Catherine Dauga; Anne Marie Castets; Rosmarie Rippka; Christiane Bouchier; Nicole Tandeau de Marsac; Martin Welker
Journal:  Mol Biol Evol       Date:  2008-07-08       Impact factor: 16.240

View more
  5 in total

Review 1.  Algal Toxic Compounds and Their Aeroterrestrial, Airborne and other Extremophilic Producers with Attention to Soil and Plant Contamination: A Review.

Authors:  Georg Gӓrtner; Maya Stoyneva-Gӓrtner; Blagoy Uzunov
Journal:  Toxins (Basel)       Date:  2021-04-29       Impact factor: 4.546

Review 2.  Bioactive Peptides Produced by Cyanobacteria of the Genus Nostoc: A Review.

Authors:  Anna Fidor; Robert Konkel; Hanna Mazur-Marzec
Journal:  Mar Drugs       Date:  2019-09-29       Impact factor: 5.118

3.  Nostocyclopeptides as New Inhibitors of 20S Proteasome.

Authors:  Anna Fidor; Katarzyna Cekała; Ewa Wieczerzak; Marta Cegłowska; Franciszek Kasprzykowski; Christine Edwards; Hanna Mazur-Marzec
Journal:  Biomolecules       Date:  2021-10-08

4.  Anabaenopeptins from Nostoc edaphicum CCNP1411.

Authors:  Robert Konkel; Michał Grabski; Marta Cegłowska; Ewa Wieczerzak; Grzegorz Węgrzyn; Hanna Mazur-Marzec
Journal:  Int J Environ Res Public Health       Date:  2022-09-28       Impact factor: 4.614

5.  Nostoc edaphicum CCNP1411 from the Baltic Sea-A New Producer of Nostocyclopeptides.

Authors:  Anna Fidor; Michał Grabski; Jan Gawor; Robert Gromadka; Grzegorz Węgrzyn; Hanna Mazur-Marzec
Journal:  Mar Drugs       Date:  2020-08-26       Impact factor: 5.118

  5 in total

北京卡尤迪生物科技股份有限公司 © 2022-2023.