Absolute Structure Determination and Kv1.5 Ion Channel Inhibition Activities of New Debromoaplysiatoxin Analogues

Potassium channel Kv1.5 has been considered a key target for new treatments of atrial tachyarrhythmias, with few side effects. Four new debromoaplysiatoxin analogues with a 6/6/12 fused ring system were isolated from marine cyanobacterium Lyngbya sp. Their planar structures were elucidated by HRESIMS, 1D and 2D NMR. The absolute configuration of oscillatoxin J (1) was determined by single-crystal X-ray diffraction, and the absolute configurations of oscillatoxin K (2), oscillatoxin L (3) and oscillatoxin M (4) were confirmed on the basis of GIAO NMR shift calculation followed by DP4 analysis. The current study confirmed the absolute configuration of the pivotal chiral positions (7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R) at traditional ATXs with 6/12/6 tricyclic ring system. Compound 1, 2 and 4 exhibited blocking activities against Kv1.5 with IC50 values of 2.61 ± 0.91 µM, 3.86 ± 1.03 µM and 3.79 ± 1.01 µM, respectively. However, compound 3 exhibited a minimum effect on Kv1.5 at 10 µM. Furthermore, all of these new debromoaplysiatoxin analogs displayed no apparent activity in a brine shrimp toxicity assay.

Voltage-gated K + channels (Kv) are membrane-inserted K + selective protein complexes [15]. Kv channels are important for various physiological and pathophysiological processes [16,17]. The Shaker-related Kv1 family consisting of subtypes (Kv1.1-Kv1.8) present in most mammalian peripheral tissues such as cardiovascular, nervous and the immune system, and many of them have been identified as potential targets for a variety of marine toxins [18][19][20]. Dalazatide (ShK-186), targeting the Kv1.3 ion channel [21], has completed phase 1 clinical trials for the treatment of autoimmune diseases [22]. Gambierol, a marine polycyclic ether toxin produced by the dinoflagellate Gambierdiscus toxicus with IC 50 of 34.5 ± 1.5 nM aganist Kv1.2 channel, might be deemed as a lead compound in further studies of the treatment of pathogenic conditions [23]. Two toxins identified from the venom of Bunodosoma caissarum, BcsTx1 and BcsTx2, displayed the highest affinity for Kv1.6 with IC 50 of 1.31 ± 0.20 nM and 7.76 ± 1.90 nM, respectively [24].
The aplysiatoxins (ATXs) and their related analogues (oscillatoxins and nhatrangins) are distinct polyketide classes of marine toxins isolated from several cyanobacterial species, including Oscillatoria nigro-viridis, Schizothrix calcicola and Lyngbya majuscula [25][26][27][28]. In our previous studies, we have pre-screened many ATXs and debromoaplysiatoxin analogues (DATs) for inhibitory activity on the shaker-related subfamily of voltage-gated channels (Kv1.1-Kv1.5), and the results indicated that some analogues presented selective and strong blocking effects on potassium channel Kv1.5 [29]. The ultrarapid activating delayed rectifier K + current (IKur) carried by the Kv1.5 channel is the main repolarization current in human atria but has no effect in the ventricle [30,31]. Therefore, Kv1.5 has become a significant molecular target for the treatment of atrial tachyarrhythmias with minimum side effects. In order to find additional novel Kv1.5 inhibitors, our team has isolated four new debromoaplysiatoxin analogues (Figure 1), oscillatoxin J-M (1-4) from the cyanobacterium Lyngbya sp. The planar structures of these compounds were elucidated by analysis of MS and NMR data, and the absolute configurations were determined by single-crystal X-ray diffraction combined with gauge invariant atomic orbital (GIAO) NMR shift calculation followed by DP4 analysis. DP4 is a probability analysis method based on the errors in each 13 C or 1 H chemical shift in the GIAO NMR calculation. The bioactivity results implicated that these four new DAT analogues and DAT exhibited differential Kv1.5 blocking activities and brine shrimp toxicities in correlation with the existence of structural functionality change.
completed phase 1 clinical trials for the treatment of autoimmune diseases [22]. Gambierol, a marine polycyclic ether toxin produced by the dinoflagellate Gambierdiscus toxicus with IC50 of 34.5 ± 1.5 nM aganist Kv1.2 channel, might be deemed as a lead compound in further studies of the treatment of pathogenic conditions [23]. Two toxins identified from the venom of Bunodosoma caissarum, BcsTx1 and BcsTx2, displayed the highest affinity for Kv1.6 with IC50 of 1.31 ± 0.20 nM and 7.76 ± 1.90 nM, respectively [24].
The aplysiatoxins (ATXs) and their related analogues (oscillatoxins and nhatrangins) are distinct polyketide classes of marine toxins isolated from several cyanobacterial species, including Oscillatoria nigro-viridis, Schizothrix calcicola and Lyngbya majuscula [25][26][27][28]. In our previous studies, we have pre-screened many ATXs and debromoaplysiatoxin analogues (DATs) for inhibitory activity on the shaker-related subfamily of voltage-gated channels (Kv1.1-Kv1.5), and the results indicated that some analogues presented selective and strong blocking effects on potassium channel Kv1.5 [29]. The ultrarapid activating delayed rectifier K + current (IKur) carried by the Kv1.5 channel is the main repolarization current in human atria but has no effect in the ventricle [30,31]. Therefore, Kv1.5 has become a significant molecular target for the treatment of atrial tachyarrhythmias with minimum side effects. In order to find additional novel Kv1.5 inhibitors, our team has isolated four new debromoaplysiatoxin analogues (Figure 1), oscillatoxin J-M (1-4) from the cyanobacterium Lyngbya sp. The planar structures of these compounds were elucidated by analysis of MS and NMR data, and the absolute configurations were determined by single-crystal X-ray diffraction combined with gauge invariant atomic orbital (GIAO) NMR shift calculation followed by DP4 analysis. DP4 is a probability analysis method based on the errors in each 13 C or 1 H chemical shift in the GIAO NMR calculation. The bioactivity results implicated that these four new DAT analogues and DAT exhibited differential Kv1.5 blocking activities and brine shrimp toxicities in correlation with the existence of structural functionality change.    9 Na, 595.2883). Its 13 C and DEPT NMR spectra exhibited 32 carbon signals, attributed to one methoxy, six methyls, four methylenes, thirteen methines and eight quaternary carbons ( Table 1). The 1 H NMR spectrum and the 13 C NMR spectrum of 1 bore a close resemblance to that of oscillatoxin B1 [32], except a double bond between C-4 and C-5 instead of a hydroxyl group on C-4. The HMBC correlations of H-5 to C-3, C-4, C-6 and C-7 and H 3 -26 to C-3, C-4 and C-5 strongly supported this assignment. Hence, the planar structure of 1 was established. Furthermore, the structure of 1 was proved by X-ray analysis ( Figure 2).
Oscillatoxin  13 C and DEPT spectra, including eight quaternary carbons, twelve methines, five methylenes and seven methyls ( Figure 3). Its NMR data were similar to 2-hydorxyanhydroaplysiatoxin, except for the absence of a bromine atom on C-17 [28]. The detailed 1 H and 13 C NMR signal assignments and connectivity were determined from a combination of 1   13 C and DEPT spectra, including eight quaternary carbons, twelve methines, five methylenes and seven methyls ( Figure 3). Its NMR data were similar to 2-hydorxyanhydroaplysiatoxin, except for the absence of a bromine atom on C-17 [28]. The detailed 1 H and 13 C NMR signal assignments and connectivity were determined from a combination of 1   The relative configuration of compound 2 was ascertained by detailed Nuclear analysis Overhauser Effect Spectroscopy (NOESY) spectrum, the coupling constants, NMR analysis and biogenetically related configuration inference of compound 2. The NOE correlations between H-9 and H-8a and small couplings of H-9 to H-8b (J8b,9 =2.8 Hz) proved that H-9 showed an equatorial position on ring B. H-10 and H-11 were determined as axial orientation by the large coupling constant (J = 10.1 Hz) of H-10/H-11. On the basis of NOE correlation of H-11/H-12, Newman's projection analysis of energy equivalent isomers using nuclear coupling constant information and steric hindrance ( Figure 4) was performed. The analysis suggested a gauche conformer of H-11/H-12, which is possessed by oscillatoxin J-M (1-4). For compound 2, the three large groups (-OR1, -CH2R3 and -CH(CH3) R2) in model A2 were extremely close in space, causing a large steric hindrance; hence, model A2 was also eliminated. The 1 H-1 H coupling constant (1.8Hz) between H-11 and H-12 indicated that there was a gauche relationship between these two protons; thus, model A1 was excluded. H-11 and H-12 were oriented in the same plane, which was confirmed by remaining model A3. The NOESY spectrum correlations between H-9/H-10, H3-23/H-11/H-12 and H-10/H3-22 proved that these protons were had the same orientations. In addition, taking notice of the structural similarities of oscillatoxin J-M (1-4), these four compounds are likely to have a common biosynthetic origin [26]. The coupling constants (J = 11.9 Hz) of H-29/H-28a and (J = 2.1Hz) of H-29/H-28b, in keeping with those of aplysiatoxins, indicated the syn relationship between H-29 and H-30. The chemical shifts of H-15 and its coupling constants (J = 6.4 Hz) were similar to those of aplysiatoxins [33]. Simulated and experimental 13 C NMR chemical shifts of 2a and 2b were used for DP4 probabil- The relative configuration of compound 2 was ascertained by detailed Nuclear analysis Overhauser Effect Spectroscopy (NOESY) spectrum, the coupling constants, NMR analysis and biogenetically related configuration inference of compound 2. The NOE correlations between H-9 and H-8a and small couplings of H-9 to H-8b (J 8b,9 =2.8 Hz) proved that H-9 showed an equatorial position on ring B. H-10 and H-11 were determined as axial orientation by the large coupling constant (J = 10.1 Hz) of H-10/H-11. On the basis of NOE correlation of H-11/H-12, Newman's projection analysis of energy equivalent isomers using nuclear coupling constant information and steric hindrance ( Figure 4) was performed. The analysis suggested a gauche conformer of H-11/H-12, which is possessed by oscillatoxin J-M (1-4). For compound 2, the three large groups (-OR1, -CH 2 R3 and -CH(CH 3 ) R2) in model A2 were extremely close in space, causing a large steric hindrance; hence, model A2 was also eliminated. The 1 H-1 H coupling constant (1.8Hz) between H-11 and H-12 indicated that there was a gauche relationship between these two protons; thus, model A1 was excluded. H-11 and H-12 were oriented in the same plane, which was confirmed by remaining model A3. The NOESY spectrum correlations between H-9/H-10, H 3 -23/H-11/H-12 and H-10/H 3 -22 proved that these protons were had the same orientations. In addition, taking notice of the structural similarities of oscillatoxin J-M (1-4), these four compounds are likely to have a common biosynthetic origin [26]. The coupling constants (J = 11.9 Hz) of H-29/H-28a and (J = 2.1Hz) of H-29/H-28b, in keeping with those of aplysiatoxins, indicated the syn relationship between H-29 and H-30. The chemical shifts of H-15 and its coupling constants (J = 6.4 Hz) were similar to those of aplysiatoxins [33]. Simulated and experimental 13 C NMR chemical shifts of 2a and 2b were used for DP4 probability analysis. The calculations were performed by using the density functional theory (DFT) as carried out in the Gaussian 09 [35]. The statistical results indicated the structural equivalence of 2 to 2a (98.61% probability). In summary, the absolute configuration of 2 was established as 2R, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R.   Figure S2.3.9)). It has a molecular formula of C32H48O11, with nine degrees of unsaturation, as assigned by HR-ESI-MS data (m/z 631.3087 [M + Na] + , calculated for C32H48O11Na, 631.3094). The interpretation of the 1D and 2D NMR spectra indicated that the planar structure of 3 closely agreed with debromoaplysiatoxin except for a hydroxyl on C-4 in 3. The HMBC correlations of H3-26 to C-3, C-4 and C-5 and H-2 to C-1 and C-3 strongly supported this connection. By comparing the 1 H chemical shift and coupling constant of 3 with that of DAT 3, it was observed that they have the same stereochemical properties as DAT, except for C-3 and C-4 [36]. The DP4 analysis was again applied to the simulated 13 C NMR chemical shifts of the four possible epimers 3a-3d ( Figure 5). The results showed that the correct structure of compound 3 is the epimer 3a, with 100% probability ( Figure 5 and Table  S1.4.8., Supplementary Materials Information). Hence, the absolute configuration of 3 was determined to be 3S, 4S, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R ( Figure 5).  Figure S2.3.9)). It has a molecular formula of C 32 H 48 O 11 , with nine degrees of unsaturation, as assigned by HR-ESI-MS data (m/z 631.3087 [M + Na] + , calculated for C 32 H 48 O 11 Na, 631.3094). The interpretation of the 1D and 2D NMR spectra indicated that the planar structure of 3 closely agreed with debromoaplysiatoxin except for a hydroxyl on C-4 in 3. The HMBC correlations of H 3 -26 to C-3, C-4 and C-5 and H-2 to C-1 and C-3 strongly supported this connection. By comparing the 1 H chemical shift and coupling constant of 3 with that of DAT 3, it was observed that they have the same stereochemical properties as DAT, except for C-3 and C-4 [36]. The DP4 analysis was again applied to the simulated 13 C NMR chemical shifts of the four possible epimers 3a-3d ( Figure 5). The results showed that the correct structure of compound 3 is the epimer 3a, with 100% probability ( Figure 5 and Table S1.4.8., Supplementary Materials Information). Hence, the absolute configuration of 3 was determined to be 3S, 4S, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R ( Figure 5).  Figure S2.3.9)). It has a molecular formula of C32H48O11, with nine degrees of unsaturation, as assigned by HR-ESI-MS data (m/z 631.3087 [M + Na] + , calculated for C32H48O11Na, 631.3094). The interpretation of the 1D and 2D NMR spectra indicated that the planar structure of 3 closely agreed with debromoaplysiatoxin except for a hydroxyl on C-4 in 3. The HMBC correlations of H3-26 to C-3, C-4 and C-5 and H-2 to C-1 and C-3 strongly supported this connection. By comparing the 1 H chemical shift and coupling constant of 3 with that of DAT 3, it was observed that they have the same stereochemical properties as DAT, except for C-3 and C-4 [36]. The DP4 analysis was again applied to the simulated 13 C NMR chemical shifts of the four possible epimers 3a-3d ( Figure 5). The results showed that the correct structure of compound 3 is the epimer 3a, with 100% probability ( Figure 5 and Table  S1.4.8., Supplementary Materials Information). Hence, the absolute configuration of 3 was determined to be 3S, 4S, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R ( Figure 5).   Figure S2.4.9)). The molecular formula of C 32 H 46 O 10 with 10 degrees of unsaturation was inferred from HRESIMS data at m/z 613.2994 [M + Na] + (calcd for C 32 H 46 O 10 Na, 613.2989). The inspection of spectral data showed that the planar structure of oscillatoxin M (4) was identical to that of debromoaplysiatoxin (5), with the exception that C-3 had a methyl group instead of a hydroxyl group, and C-4 had a ketone carbonyl group instead of a methyl group. The HMBC correlations of H-5 to C-3, C-4 and C-6; H 3 -26 to C-3 and C-4; and H-2 to C-1, C-3 and C-4 strongly support this assignment. The relative configuration of H-9/H-10/H-11/H-12, H-15 and H-29/H-30 in compound 4 was determined in accordance with that of compound 2. Furthermore, the results of DP4 statistical analysis using the 13 C NMR chemical shift values proved that the correct structure for 4 is the epimer 4a ( Figure 6 and Supplementary Information Table S1.4.11.). Finally, the absolute configuration of 4 was established as 3R, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R.  Figure S2.4.9)). The molecular formula of C32H46O10 with 10 degrees of unsaturation was inferred from HRESIMS data at m/z 613.2994 [M + Na] + (calcd for C32H46O10Na, 613.2989). The inspection of spectral data showed that the planar structure of oscillatoxin M (4) was identical to that of debromoaplysiatoxin (5), with the exception that C-3 had a methyl group instead of a hydroxyl group, and C-4 had a ketone carbonyl group instead of a methyl group. The HMBC correlations of H-5 to C-3, C-4 and C-6; H3-26 to C-3 and C-4; and H-2 to C-1, C-3 and C-4 strongly support this assignment. The relative configuration of H-9/H-10/H-11/H-12, H-15 and H-29/H-30 in compound 4 was determined in accordance with that of compound 2. Furthermore, the results of DP4 statistical analysis using the 13 C NMR chemical shift values proved that the correct structure for 4 is the epimer 4a ( Figure 6 and Supplementary Information Table S1.4.11.). Finally, the absolute configuration of 4 was established as 3R, 7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R.

Inhibitory Activities against Kv1.5
The Kv1.5 (ultra-fast-delay rectifier potassium channel) mediation of ultra-rapid delayed rectifier K + current (IKur) is the main current in the repolarization process of atrial action potentials. Our previous research has highlighted the capability of ATXs as ion channel blockers [34,37]. Some compounds have blocked Kv1.5 with certain selectivity [29]. In this work, a dose-response study was conducted on compounds 1-4 to evaluate their inhibitory activity against Kv1.5. Our results showed that compounds 1, 2 and 4 had significant inhibitory effects on Kv1.5, with IC50 values of 2.61 ± 0.91 µM (Figure 7), 3.86 ± 1.03 µM (Figure 8) and 3.79 ± 1.01 µM (Figure 9), respectively. Compound 1, 2 and 4 all voltage dependently inhibited the Kv1.5 current (Supplementary Materials 1.3). The inhibition was stronger at 50 mV than other tested potentials. The voltage dependence sup-  The Kv1.5 (ultra-fast-delay rectifier potassium channel) mediation of ultra-rapid delayed rectifier K + current (IKur) is the main current in the repolarization process of atrial action potentials. Our previous research has highlighted the capability of ATXs as ion channel blockers [34,37]. Some compounds have blocked Kv1.5 with certain selectivity [29]. In this work, a dose-response study was conducted on compounds 1-4 to evaluate their inhibitory activity against Kv1.5. Our results showed that compounds 1, 2 and 4 had significant inhibitory effects on Kv1.5, with IC 50 values of 2.61 ± 0.91 µM (Figure 7), 3.86 ± 1.03 µM (Figure 8) and 3.79 ± 1.01 µM (Figure 9), respectively. Compound 1, 2 and  Materials 1.3). The inhibition was stronger at 50 mV than other tested potentials. The voltage dependence supported the fact that three compounds preferentially affected the open state of the Kv1.5 channels. However, compound 3, oscillatoxin L, is structurally almost identical to DAT (IC 50 = 1.28 ± 0.08 µM), except having an additional hydroxy motif on adjacent carbon C-4, exhibiting minimum effects on the modulation of Kv1.5 at 10 µM ( Figure S1.3.1). Aplysiatoxin and its derivatives are activators of protein kinase C (PKC), which has been well researched [38]. In a previous study, we have demonstrated DAT strongly upregulating the expression of phosphor PKCδ in human hepatocellular carcinomas (HepG2) at 10 µM [29] and proposed the potential mechanism for DATs modulating the Kv channel by activating protein kinase C [39]. The results might indicate that the simple modifications of the functionalities on A ring system in DAT structure may affect either the interaction of PKC or the blocking site of the Kv1.5 ion channel. Consequently, determining the mechanism of Kv1.5 inhibition activity of DAT analogues will be our ongoing project. 4, exhibiting minimum effects on the modulation of Kv1.5 at 10 µM ( Figure S1.3.1.). Aplysiatoxin and its derivatives are activators of protein kinase C (PKC), which has been well researched [38]. In a previous study, we have demonstrated DAT strongly upregulating the expression of phosphor PKCδ in human hepatocellular carcinomas (HepG2) at 10 µM [29] and proposed the potential mechanism for DATs modulating the Kv channel by activating protein kinase C [39]. The results might indicate that the simple modifications of the functionalities on A ring system in DAT structure may affect either the interaction of PKC or the blocking site of the Kv1.5 ion channel. Consequently, determining the mechanism of Kv1.5 inhibition activity of DAT analogues will be our ongoing project.   4, exhibiting minimum effects on the modulation of Kv1.5 at 10 µM ( Figure S1.3.1.). Aplysiatoxin and its derivatives are activators of protein kinase C (PKC), which has been well researched [38]. In a previous study, we have demonstrated DAT strongly upregulating the expression of phosphor PKCδ in human hepatocellular carcinomas (HepG2) at 10 µM [29] and proposed the potential mechanism for DATs modulating the Kv channel by activating protein kinase C [39]. The results might indicate that the simple modifications of the functionalities on A ring system in DAT structure may affect either the interaction of PKC or the blocking site of the Kv1.5 ion channel. Consequently, determining the mechanism of Kv1.5 inhibition activity of DAT analogues will be our ongoing project.

Toxicity of Brine Shrimp
In order to understand the toxicity effect of such compounds, brine shrimp Artemia salina (A. salina) was used as a model organism. The investigation of brine shrimp toxicity of debromoaplysiatoxin (DAT) and four DAT analogs (oscillatoxin J, oscillatoxin K, oscillatoxin L and oscillatoxin M) isolated from marine cyanobacterium Lyngbya sp., was conducted. When DAT concentration was as low as 0.1 µM, the survival of Artemia salina (A. salina) began to be affected (Table S1.3.2.1.). Compounds 1, 2, 3 and 4 had no apparent effect at 30 µM. As shown in Figure 10, compared to other tested derivatives, debromoaplysiatoxin was the most toxic compound (IC50 value = 0.34 ± 0.036 µM) ( Figure S1.3.2.1.). The current results with a previous study [37] indicated that the 3-hydroxy group at DAT seemed to play an important role in determining higher toxicity. However, the specific mechanism of toxic action of these compounds remains unclear, and further studies are needed.

General Experimental Procedure
UV was acquired on a UV/EV300 spectrometer (Thermo Scientific, Waltham, MA, USA), and IR spectra were obtained on a Nicolet iS20 instrument (Thermo Fisher Scientific, MA, USA). Optical rotation data were performed by a Jasco P-2000 polarimeter (Jasco, Hachiojishi, Tokyo, Japan). 1 H and 13 C NMR spectra were acquired on an Agilent 600 MHz spectrometer (Agilent Technologies, Santa Clara, CA, USA), with CDCl3 (δH 7.26

Toxicity of Brine Shrimp
In order to understand the toxicity effect of such compounds, brine shrimp Artemia salina (A. salina) was used as a model organism. The investigation of brine shrimp toxicity of debromoaplysiatoxin (DAT) and four DAT analogs (oscillatoxin J, oscillatoxin K, oscillatoxin L and oscillatoxin M) isolated from marine cyanobacterium Lyngbya sp., was conducted. When DAT concentration was as low as 0.1 µM, the survival of Artemia salina (A. salina) began to be affected (Table S1.3.2.1). Compounds 1, 2, 3 and 4 had no apparent effect at 30 µM. As shown in Figure 10, compared to other tested derivatives, debromoaplysiatoxin was the most toxic compound (IC 50 value = 0.34 ± 0.036 µM) (Figure S1.3.2.1). The current results with a previous study [37] indicated that the 3-hydroxy group at DAT seemed to play an important role in determining higher toxicity. However, the specific mechanism of toxic action of these compounds remains unclear, and further studies are needed. Figure 9. Kv1.5 currents were evoked by a 300 ms depolarizing pulse from -50 mV to 50 mV in 20 mV increments from a holding potential of -70 mV in the absence and presence of 0.3 µM, 1 µM, 3 µM and 10 µM oscillatoxin M (4). The current amplitudes were measured at the end of the 300 ms pulse at 50 mV. (B) Concentration-inhibition curve expressed in %. The abscissa represents the concentration, and the ordinate represents the percentage of Kv1.5 current that is blocked at different concentrations of oscillatoxin M (4). Data points represent mean ± SEM of 3 to 5 measurements, and the inhibitory effect showed an IC50 value of 3.79 ± 1.01 µM.

Toxicity of Brine Shrimp
In order to understand the toxicity effect of such compounds, brine shrimp Artemia salina (A. salina) was used as a model organism. The investigation of brine shrimp toxicity of debromoaplysiatoxin (DAT) and four DAT analogs (oscillatoxin J, oscillatoxin K, oscillatoxin L and oscillatoxin M) isolated from marine cyanobacterium Lyngbya sp., was conducted. When DAT concentration was as low as 0.1 µM, the survival of Artemia salina (A. salina) began to be affected (Table S1.3.2.1.). Compounds 1, 2, 3 and 4 had no apparent effect at 30 µM. As shown in Figure 10, compared to other tested derivatives, debromoaplysiatoxin was the most toxic compound (IC50 value = 0.34 ± 0.036 µM) ( Figure S1.3.2.1.). The current results with a previous study [37] indicated that the 3-hydroxy group at DAT seemed to play an important role in determining higher toxicity. However, the specific mechanism of toxic action of these compounds remains unclear, and further studies are needed.

Material
Cyanobacterium Lyngbya sp. was obtained from the harbor of Sanya, Hainan province, China, in November 2016. The sample was identified by Prof. Bing-Nan Han (Zhejiang Sci-Tech University, Zhejiang, China). After morphological and molecular identification, a voucher specimen (voucher number: BNH-201606; gene bank accession numbers: MH636576) has been well deposited in Zhejiang Sci-Tech University.

Electrophysiology
LTK cells cultured for at least 24 h could be used for currents recording. For Kv1.5 potassium current recording, the recording micropipettes were pulled with a P97 microelectrode puller (Sutter, CA, USA) with a resistance of about 3 MΩ when filling with internal solution containing the following: KCl 140 mM, MgCl 2 1 mM, EGTA 5 mM, HEPES 10 mM and MgATP 1 mM (pH was adjusted to 7.25 by KOH). The bath solution contained the following: NaCl 137 mM, KCl 5.4 mM, CaCl 2 1.8 mM, MgCl 2 1 mM, Glucose 10 mM and HEPES 10 mM (pH was adjusted to 7.4 by NaOH). Kv1.5 currents were recorded at room temperature (22-24 • C) by PulseMaster (Version 2.65, Heka, Lambrecht, Germany) via an EPC-10 USB amplifier (Heka, Lambrecht, Germany). In order to reduce recording errors, cells with seal resistance above 1 GΩ and series resistance that was fully compensated above 80% were used. Leak compensation was used to compensate the leak current and to subtract the capacitive artifacts (Supplementary Materials 1.3.1.2).

Brine Shrimp Toxicity Assay
Brine shrimp A. salina is an important model organism for ecosystems, and because of its high sensitivity and easy availability, it can be used in laboratory settings to study toxic effects and to provide safe results. A. salina or brine shrimp cysts were cultivated in 3.2% of saline water. After aeration with normal saline, cysts were placed at temperature for 24 h and then cultured. For toxicity screening, hatched larvae were collected and introduced in saline water. Add 0.9% brine and 30 larvae with good activity to each well to produce a 96-well test culture plate. New compounds of 0.1 µM, 1 µM, 10 µM and 30 µM and DAT were added to the experimental culture plate, respectively. Dimethyl sulfoxide (Aladdin, Shanghai, China) and dichloromethane (Aladdin, Shanghai, China) were added as blank control and positive control. The survival rate of A. salina was calculated after 24 h treatment at 25 • C.

Conclusions
In conclusion, four new debromoaplysiatoxin analogues, oscillatoxin J-M (1-4), were isolated from the cyanobacterium Lyngbya sp. The structures of the new compounds were characterized by 1D and 2D NMR and MS data. ATXs easily underwent structural rearrangement to form new structures with new functionality and new stereochemistry because of the presence of some unstable functional groups such as hemiacetal and ketal, etc. However, the current study confirmed the absolute configuration of pivotal chiral positions (7S, 9S, 10S, 11R, 12S, 15S, 29R and 30R) at traditional ATXs with 6/12/6 tricyclic ring system of compounds 1-4 via X-ray diffraction and GIAO NMR shift calculation followed by DP4 analysis. Compound 1, 2 and 4 showed inhibitory activities against Kv1.5 with IC 50 value of 2.61 ± 0.91 µM, 3.86 ± 1.03 µM and 3.79 ± 1.01 µM, respectively. Compound 3 exhibited no apparent Kv1.5 inhibition activity at 10 µM. Discovery of new DAT analogs may highlight some insightful structures and activity relationships for developing strong, effective and safe Kv1.5 inhibitors in the future. Our future work will focus on characterizing the selectivity of these compounds for Kv1.5 and how they inhibit Kv1.5 channel.