Antibacterial Oligomeric Polyphenols from the Green Alga Cladophora socialis.

A series of oligomeric phenols including the known natural product 3,4,3',4'-tetrahydroxy-1,1'-biphenyl (3), the previously synthesized 2,3,8,9-tetrahydroxybenzo[ c]chromen-6-one (4), and eight new related natural products, cladophorols B-I (5-12), were isolated from the Fijian green alga Cladophora socialis and identified by a combination of NMR spectroscopy, mass spectrometric analysis, and computational modeling using DFT calculations. J-resolved spectroscopy and line width reduction by picric acid addition aided in resolving the heavily overlapped aromatic signals. A panel of Gram-positive and Gram-negative pathogens used to evaluate pharmacological potential led to the determination that cladophorol C (6) exhibits potent antibiotic activity selective toward methicillin-resistant Staphylococcus aureus (MRSA) with an MIC of 1.4 μg/mL. Cladophorols B (5) and D-H (7-11) had more modest but also selective antibiotic potency. Activities of cladophorols A-I (4-12) were also assessed against the asexual blood stages of Plasmodium falciparum and revealed cladophorols A (4) and B (5) to have modest activity with EC50 values of 0.7 and 1.9 μg/mL, respectively.

Introduction

Aromatic systems are ubiquitous within natural products, reflecting the availability of their biosynthetic precursors in plants, algae, and microorganisms via the shikimate and polyketide pathways and the varied biological effects of such functional groups.[1,2] When incorporated as building blocks in vascular plants, aromatic polymers[3−6] such as lignin provide structural support that ultimately enabled plants to occupy terrestrial ecosystems leading to the evolution of trees.[7] Aromatic natural products conferring structural rigidity were unknown in the marine realm until lignin was reported in the coralline red alga Calliarthron cheilosporioides (phylum Rhodophyta).[8] This startling discovery suggested that an ancestor of green and red algae acquired the capacity to polymerize phenolic monomers but that this feature was silenced in most evolving lineages. Other red algae produce a variety of structurally simple, often halogenated phenolics of mostly unknown biological function,[9] although more complex, aromatic natural products of mixed shikimate and isoprenoid biosynthetic origin have also been observed.[10] Many members of the brown algae (phylum Ochrophyta), which evolutionarily diverged just before land plants, produce polymeric, acetogenic, phenolic compounds (phlorotannins) constituting 1–20% of dry mass.[11] These oligomers of phloroglucinol protect tissues against UV radiation and herbivory[11,12] and may contribute to structural support.[12] In contrast, significantly fewer phenolic natural products have been reported in green algae (phylum Chlorophyta). Brominated polyphenolic compounds including the feeding deterrent avrainvilleol and the condensation product of two molecules of avrainvilleol, rawsonol, are produced by green algae belonging to the genus Avrainvillea.[13−15] Another example of phenolic compounds from green algae are two vanillic acid analogues (1–2) from an Australian Cladophora socialis that have been found to disrupt insulin cell signaling.[16]

The present study describes a series of related oligomeric phenols, cladophorols A–I (4–12), from the green alga C. socialis that attracted our attention for its antibiotic activity against methicillin-resistant Staphylococcus aureus (MRSA). Analysis of mass spectrometric and NMR spectroscopic data enabled structure elucidation despite the high molecular weights and repeating structural patterns of members of this natural product family. We used J-resolved NMR spectroscopy to reduce signal overlap in the 1H and 13C NMR spectra and employed picric acid line width-narrowing to reveal the exact linkage pattern of the subunits.[17] In addition, we propose a structural revision of known vanillic acid analogues (1–2) based on a combination of spectroscopic and computational approaches.

Results and Discussion

Collection and Isolation

Offshore of Vanua Levu Island in Fiji, C. socialis, identified by morphological features and nuclear small-subunit 18S rRNA analyses (Figure S1), was collected floating at the ocean surface. The aqueous methanol extract was partitioned using a sequence of biphasic organic solvents and water. Following evaluation of MRSA activity, the ethyl acetate-soluble fraction was chromatographed on silica gel followed by repeated reversed-phase semipreparative HPLC, yielding the known compound 3,4,3′,4′-tetrahydroxy-1,1′-biphenyl (3), previously synthesized 2,3,8,9-tetrahydroxybenzo[c]chromen-6-one (4),[18] and eight new phenolic natural products (5–12) (Figure ).

Figure 1

Structures of compounds 3–12 from C. socialis.

Characterization of 4

HRESIMS was used to assign the molecular formula of 4 as C13H8O6. The 1H NMR spectrum exhibited only four singlets at δH 7.58, 7.33, 7.32, and 6.73 (Table ), attributed to two aromatic rings based on HSQC correlations with δC 115.4, 108.3, 107.4, and 104.4, respectively, and the absence of coupling between the protons. Correlations observed in HMBC from the signal at δH 6.73 (H-4) to δC 111.4 (C-10b), 144.4 (C-2), 146.0 (C-4a), and 148.6 (C-3) established the first ring. The second ring was determined from HMBC correlations from δH 7.58 (H-7) to δC 131.7 (C-10a), 147.3 (C-8), 154.8 (C-9), and 164.0 (C-6). Direct linkage between the two aromatic rings was established by HMBC correlations from δH 7.32 (H-1) to C-10a and from δH 7.33 (H-10) to C-10b. With connectivity of the carbonyl functionality and the two aromatic rings established, a third ring was necessary to account for the 10 degrees of unsaturation calculated for this compound.

Table 1

1H (800 MHz) and 13C (200 MHz) NMR Spectroscopic Data of Cladophorols A and B (4–5) in CD3OD

	4		5
position	δC	δH	δC	δH
Benzo[c]chromenone
1	107.4 (d)	7.32 s	108.2 (d)	8.37 d (2.8)
2	144.4 (s)		155.6 (s)
3	148.6 (s)		106.9 (d)	6.69 d (2.8)
4	104.4 (d)	6.73 s	146.6 (s)
4a	146.0 (s)		136.3 (s)
6	164.0 (s)		163.4 (s)
6a	112.9 (s)		114.5 (s)
7	115.4 (d)	7.58 s	123.9 (d)	7.84 d (8.5)
8	147.3 (s)		116.4 (d)	7.08 d (8.5)
9	154.8 (s)		152.7 (s)
10	108.3 (d)	7.33 s	144.3 (s)
10a	131.7 (s)		123.1 (s)
10b	111.4 (s)		120.7 (s)
Biphenyl
1A			139.6 (s)
2A			116.0 (d)	7.14 d (2.2)
3A			150.1 (s)
4A			144.0 (s)
5A			122.1 (d)	6.98 d (8.3)
6A			119.1 (d)	7.02 dd (8.3, 2.2)
1B			133.9 (s)
2B			114.9 (d)	7.04 d (2.2)
3B			146.5 (s)
4B			145.9 (s)
5B			116.6 (d)	6.82 d (8.2)
6B			119.2 (d)	6.93 dd (8.2, 2.2)

Since the 1H and 13C NMR chemical shifts of 4 were similar to those reported for 1, the structure depicted in Figure a was considered. Although it did not conflict with any of the observed NMR spectroscopic features, full consideration of all plausible structures and evaluation of these data to determine the best fit led to a more plausible alternative (Figure b). An HMBC spectrum recorded in DMSO-d6 with picric acid,[17] which reduced the line widths of the phenolic hydroxy region, revealed four exchangeable protons at δH 10.41, 9.89, 9.74, and 9.23, each correlating with three carbons (Figure S13). In the originally considered structure (Figure a), only three hydroxyl groups are within this range; the carboxylic acid group is expected to shift the fourth labile proton downfield. Additionally, HMBC signals from δH 9.89 to C-7, C-8, and C-9 indicated an alcohol group at C-8. Since none of the labile protons correlated with carbonyl C-6, we concluded that it was part of the third ring. Compound 4 was thus identified as 2,3,8,9-tetrahydroxybenzo[c]chromen-6-one and named cladophorol A. This compound was previously reported as a synthetic analogue of ellagic acid created to improve potency for DNA gyrase inhibition, but the spectroscopic data were not disclosed at that time.[18]

Figure 2

Possible structures (a and b) for cladophorol A (4) based on selected observed HMBC correlations indicated by arrows.

Characterization of 5

Compound 5 was obtained as a pale red amorphous solid with a molecular formula of C25H16O9 deduced by HRESIMS. The 1H and COSY NMR spectra exhibited resonances accounting for four aromatic rings. 1H–1H coupling constants were used to establish substitution patterns including one 1,2,3,5-tetrasubstituted ring, one 1,2,3,4-tetrasubstituted ring, and two 1,2,4-trisubstituted rings. HMBC correlations from H-1 to δC 123.1 (C-10a), 136.3 (C-4a), and 155.6 (C-2), from H-7 to C-10a, δC 152.7 (C-9), and 163.4 (C-6), and from H-8 to 114.5 (C-6a) and 144.3 (C-10) indicated the first two rings were connected in a benzo[c]chromen-6-one moiety similar to 4 (Figure S2), although the location of hydroxy groups differed with substitution at positions C-2, C-4, C-9, and C-10. The two other rings formed a biphenyl group based on HMBC correlations from H-2A and H-6A to δC 133.9 (C-1B), from H-5A to δC 139.6 (C-1A), from H-2B and H-6B to C-1A, and from H-5B to C-1B. Despite careful examination of the HMBC spectra, a connecting point between the biphenyl and benzo[c]chromen-6-one could not be identified. Whereas a 4J correlation was reported between H-5′ and C-4 for 1,[16] a similar correlation was not found in 5, despite recording HMBC spectra in DMSO-d6 or acetone-d6 and optimizing for 2 Hz long-range coupling. Adding picric acid to DMSO-d6 for the NMR acquisition led to six labile proton signals (Figure S23) which showed correlations with all carbons bearing an alcohol group and their neighbor (Figure S25). The four singlets at δH 11.04, 9.99, 9.68, and 9.57 exhibited HMBC correlations with δC 106.9 (C-3), 146.6 (C-4), C-4a, 116.4 (C-8), C-9, C-10, C-10a, 116.0 (C-2B), 150.1 (C-3B), and 144.0 (C-4B), clarifying that the benzo[c]chromen-6-one and the biphenyl were linked via an oxygen attached to C-2 and C-4B. For this reason, 5 was identified as 4,9,10-trihydroxy-2-[(3,3′,4′-trihydroxy-(1,1′-biphenyl)-4-yl)oxy]benzo[c]chromen-6-one and named cladophorol B.

Structure Revision of 1

Since 4 and 5 were well characterized, the validity of the structure previously reported for 1 was questioned.[16] When recorded in DMSO-d6, the carboxyl carbons of 4 and 5, which are part of a lactone ring, produced NMR signals at δC 164.5 and 160.3, respectively. The carboxyl carbon of 1 was reported at δC 161.1, a value upfield of expectations for a free carboxylic acid.[16] DFT calculations of the 1H and 13C chemical shift values were performed for the published structure 1 and for two other hypothetical structures (13 and 14) having a benzo[c]chromen-6-one unit (Figure ,).[19,20] Indeed, the mPW1PW91/6-311+G(d,p) 13C chemical shift values did not match the experimental values for the originally published structure 1 with an rmsd of 3.3 ppm (Table ). The largest discrepancies were due to C-4a and C-6 which are involved in the lactone ring. The hypothetical structures 13 and 14 showed smaller deviations for δC with rmsd of 1.7 and 1.8 ppm, respectively, but only 13 exhibited a small δH rmsd (0.08 ppm). When DP4+ probabilities were calculated, structure 13 scored 100% when considering chemical shift data for both 1H and 13C. For this reason, 1 should be revised to 13 and named 2,8,9-trihydroxy-4-((3,3′,4′-trihydroxy[1,1′-biphenyl]-4-yl)oxy)benzo[c]chromen-6-one.

Figure 3

Original and hypothetical structures proposed for 1, supporting reassignment of previously published 1 as 13.

Table 2

Experimental and Calculated NMR Chemical Shifts for 1, 13, and 14

δ values calculated at the mPW1PW91/6-311+G(d,p)//B3LYP/6-31G(d) level of theory performed with the polarizable continuum model (PCM) using DMSO. Red color highlights large differences in calculated values: |ΔδH| > 0.2, |ΔδC| > 3.0.

Root-mean-square deviation.

Maximum absolute error.

DP4+ probability.[21]

Characterization of 6

Compound 6 was isolated as a pale red solid with a molecular formula of C54H38O18 based on an [M – H]− pseudomolecular ion in the HRESIMS at m/z 973.1967. The 1H NMR spectrum was difficult to interpret because of severe signal overlap in the aromatic region. J-resolved 2D NMR spectroscopy aided the interpretation and allowed assignment of all chemical shifts and coupling constants (Table ). Using DQF-COSY and J-resolved spectroscopy, seven AMX aromatic spin systems resulting from 1,2,4-trihydroxyphenyl groups were observed as doublets at δH 6.67–6.91 and 6.42–6.57 and doublet of doublets at δH 6.27–6.47. Interestingly, within each of these ranges, a gradient of proton chemical shifts was noted whereby each set of signals from a given ring appeared in a predictable order. This feature was highlighted by simulating each set of protons using the extracted δC and J from NMR spectra (Figure ). When comparing proton signals for H-3, H-5, and H-6, the two signals corresponding to rings Y and Z appeared to be outliers. Most of the carbon signals were resolved in the 13C NMR spectrum, although HSQC and HMBC spectra showed large clusters of correlations between protons of rings D–G and carbons averaged at 157.5 (C-4D–G), δC 151.2 (C-2D–G), 139.6 (C-1D–G), 123.1 (C-6D–G), 108.8 (C-5D–G), and 106.4 (C-3D–G). This prevented a precise assignment of signals but also indicated these rings were similar. Two other spin systems unrelated to the gradient of chemical shifts showed HMBC correlations from H-2A to δC 139.9 (C-1B) and from H-2B and H-6B to δC 133.9 (C-1A) allowing assignment of these two spin systems to a biphenyl unit (Figure S2). Given all these data, two structures were hypothesized, differing by the position of the free alcohol groups on each hydroxyphenyl (Figure ).

Table 3

1H (800 MHz) and 13C (200 MHz) NMR Spectroscopic Data of Cladophorols C and D (6–7) in CD3OD

	6		7
position	δC	δH	δC	δH
Rings A–B (Biphenyl)
1A	134.0		134.0
2A	114.9	7.01 (d, 2.2)	114.9	7.01 (d, 2.2)
3A	146.5		146.5
4A	146.0		146.0
5A	116.6	6.81 (d, 8.2)	116.6	6.81 (d, 8.2)
6A	119.3	6.912 (dd, 8.2, 2.2)	119.3	6.912 (dd, 8.2, 2.3)
1B	139.9		139.9
2B	116.1	7.11 (d, 2.2)	116.1	7.11 (d, 2.2)
3B	150.2		150.2
4B	143.9		143.9
5B	122.3	6.94 (d, 8.3)	122.2	6.94 (d, 8.3)
6B	119.2	7.00 (dd, 8.3, 2.2)	119.2	7.00 (dd, 8.3, 2.2)
Ring C
1	139.8		139.87
2	151.3		151.3
3	107.0	6.57 (d, 2.9)	107.0	6.58 (d, 2.9)
4	157.0		157.0
5	109.4	6.47 (dd, 8.8, 2.9)	109.4	6.47 (dd, 8.8, 2.9)
6	123.2b	6.905 (dd, 8.8)	123.2b	6.905 (dd, 8.8)
Middle Rings
1	139.7, 139.54, 139.53, 139.50		139.76, 139.57, 139.55, 139.52 (2×)
2	151.25, 151.23, 151.21, 151.1		151.26, 151.24, 151.23, 151.21, 151.1
3	106.5, 106.36, 106.34, 106.31	6.53, 6.51, 6.503, 6.496 (all d, 2.9)	106.5 (2×), 106.37, 106.36, 106.33	6.53, 6.512, 6.51, 6.503, 6.497 (all d, 2.9)
4	157.6, 157.54, 157.51, 157.4,		157.6 (2×), 157.54, 157.51, 157.4
5	108.8, 108.77, 108.75, 108.74	6.42, 6.41, 6.40, 6.39 (all dd, 8.8, 2.9)	108.8, 108.78, 108.77 (2×), 108.75	6.423, 6.41, 6.40, 6.399, 6.388 (all dd, 8.8, 2.9),
6	123.19 (2×), 123.14, 122.9	6.88, 6.87, 6.86, 6.84 (all d, 8.8)	123.19 (2×), 123.17, 123.13, 122.9	6.88, 6.87, 6.868, 6.863, 6.84 (all d, 8.8)
Ring Ya
1	141.0		141.0
2	150.6		150.6
3	106.34b	6.48 (d, 2.9)	106.36b	6.48 (d, 2.9)
4	156.7		156.7
5	108.66b	6.34 (dd, 8.8, 2.9)	108.68b	6.34 (dd, 8.8, 2.9)
6	121.7	6.74 (d, 8.8)	121.7	6.74 (d, 8.8)
Ring Za
1	141.7		141.7
2	147.0		147.0
3	106.4b	6.42 (d, 2.9)	106.45b	6.420 (d, 2.9)
4	152.9		152.9
5	109.1	6.27 (dd, 8.6, 2.9)	109.1	6.27 (d, 8.6, 2.9)
6	116.4	6.67 (d, 8.6)	116.4	6.67 (d, 8.6)

Rings Y and Z are the last two rings of the oligomer.

13C signal can be interchanged with one of the corresponding signals of rings D–G.

Figure 4

Gradient of NMR chemical shift along the oligomer of 1,2,4-trihydroxyphenyl units for cladophorol C (6). Experimental (black) and simulated (colored) signals are shown.

Figure 5

Hypothetical structures for cladophorol C (6a and 6b) and model structures 15a and 15b.

Two model structures, 15a and 15b, containing a biphenyl unit and three 1,2,4-trihydroxyphenyl rings (Figure ), were computationally tested to determine which arrangement was more probable. Their 1H and 13C NMR chemical shifts were predicted (Tables S26 and S27) following DFT computation at the mPW1PW91/6-311+G(d,p) level of theory and compared to the corresponding position of 6. The rmsd and MAE values calculated for both 1H and 13C were lower when 15a was considered, indicating that 6a is the correct structure. Based on these results, 6 was identified as α-hydro-ω-[3,4-dihydroxyphenyl]octa[oxy(2-hydroxyphen-4-yl)] and named cladophorol C.

Characterization of 7

The molecular formula for 7 was established as C60H42O20 based on HRESIMS, accounting for one additional dihydroxyphenyl moiety. The 1H and 13C NMR spectra almost perfectly superimposed those of 6, the only differences being three additional proton signals at δH 6.86 (d, J = 8.8 Hz), 6.50 (d, J = 2.9 Hz), and 6.40 (dd, J = 8.8 Hz) and six additional carbon signals at δC 157.6, 151.23, 139.57, 123.17, 108.77, and 106.5. As a result, 7 was identified as α-hydro-ω-[3,4-dihydroxyphenyl]nona[oxy(2-hydroxyphen-4-yl)] and named cladophorol D.

Characterization of 8

Compound 8 was isolated as a reddish solid with the molecular formula C67H44O23 based on an [M – H]− HRESIMS peak at m/z 1215.2188. The 1H and 13C NMR spectra strongly resembled those of 6 and 7, suggesting that 8 was also an oligomeric phenolic with a terminal biphenyl. A limited spread of the signals for the repeating phenol units, as illustrated by the 1H spectral comparison of H-6 (Figure ), suggested an end-capping of the oligomer with a different moiety. A carbonyl, indicated by a peak on the 13C NMR spectrum at δC 163.4, and 12 other distinct 13C NMR signals, agreed with a benzo[c]chromen-6-one unit (Table ). The substitution pattern was determined as 2,4,8,9-tetrahydroxy- from the presence of two distinct singlets at δH 7.63 and 7.32 and two doublets at δH 6.96 and 6.55 (both J = 2.7 Hz) in the 1H NMR spectrum (Table ). This assignment was corroborated by HMBC correlations from H-1 and H-3 to δC 136.2 (C-4a) and from H-1 to δC 130.7 (C-10a) (Figure S2). This benzo[c]chromen-6-one group matched revised structures 13 and 14; however, the exact point of attachment of the oligomer needed to be established to determine which configuration was more probable for 8. A comparison of 13C NMR values of 8 recorded in DMSO-d6 and those reported for 13(16) (Table S1) showed the greatest differences at C-2, C-4, and C-6, suggesting the oligomer was attached to C-2 as in 14 instead of C-4. When the calculated values for 13 and 14 were considered, a 94% DP+ value favoring 14 confirmed this hypothesis. For this reason, 8 was identified as α-[(4,8,9-trihydroxy-2-benzo[c]chromen-6-on-2-yl)oxy]-ω-[3,4-dihydroxyphenyl]octa[oxy(2-hydroxyphen-4-yl)] and named cladophorol E.

Figure 6

Chemical shift spreading for cladophorol C (6) and cladophorol E (8) for position H-6. The limited spreading of 8 suggest an end-capping of the oligomer. Simulated (blue) and experimental (gray) signals are shown.

Table 4

13C (200 MHz) NMR Spectroscopic Data of Cladophorols E–I (8–12) in CD3OD

position	8	9	10	11	12
Rings A–B (Biphenyl)
1A	134.0	134.0	134.0	134.0	134.0
2A	114.9	114.9	114.9	114.9	114.9
3A	146.5	146.5	146.5	146.5	146.5
4A	146.0	146.0	146.0	146.0	146.0
5A	116.6	116.7	116.6	116.6	116.6
6A	119.3	119.3	119.3	119.3	119.3
1B	139.93	139.91	139.92	139.92	139.93
2B	116.1	116.1	116.1	116.1	116.1
3B	150.21	150.2	150.20	150.2	150.23
4B	143.9	143.9	143.9	143.9	143.9
5B	122.2	122.2	122.2	122.2	122.3
6B	119.2	119.2	119.2	119.2	119.2
Middle Rings
1	139.87, 139.6, 139.52, 139.50, 139.49, 139.38, 139.2	139.87, 139.6, 139.5, 139.48, 139.38, 139.2	139.86, 139.6, 139.5, 139.49 (2×), 139.48, 139.37, 139.2	139.86, 139.6, 139.49 (4×), 139.47, 139.36, 139.2	139.85, 139.5, 139.48 (5×), 139.46, 139.36, 139.2
2	151.32, 151.30, 151.27, 151.25, 151.23 (3×)	151.3, 151.28, 151.26, 151.23, 151.22 (2×)	151.3, 151.29, 151.27, 151.24, 151.22 (4×)	151.32, 151.30, 151.27, 151.24, 151.23 (3×), 151.22 (2×)	151.33, 151.31, 151.28, 151.25, 151.24 (6×)
3	107.0, 106.5, 106.38 (3×), 106.37 (2×)	107.1, 106.5, 106.39, 106.38, 106.37 (2×)	107.0, 106.4, 106.37 (3×), 106.36 (3×)	107.0, 106.4, 106.37 (2×), 106.36 (5×)	107.0, 106.4, 106.3 (8×)
4	157.8, 157.7, 157.59, 157.57, 157.56, 157.50, 157.0	157.8, 157.6, 157.57, 157.55, 157.49, 157.0	157.8, 157.65, 157.58, 157.57, 157.56 (2×), 157.50, 157.0	157.8, 157.6, 157.58, 157.56 (4×), 157.50, 157.0	157.8, 157.7, 157.58 (3), 157.57 (3×), 157.52, 157.0
5	109.4, 108.85, 108.81, 108.79 (2×), 108.78 (2×)	109.4, 108.84, 108.80, 108.79, 108.78, 108.77	109.4, 108.84, 108.80, 108.78, 108.77 (4×)	109.4, 108.84, 108.80, 108.78, 108.77 (5×)	109.4, 108.8, 108.79, 108.76 (7×)
6	123.4, 123.3, 123.20 (3×), 123.18 (2×)	123.4, 123.3, 123.20 (3×), 123.19	123.4, 123.3, 123.20 (3×), 123.19 (3×)	123.4, 123.3, 123.20, 123.18 (6×)	123.41, 123.35, 123.22 (2), 123.20 (6×)
Benzo[c]chromenone
1	100.4	100.4	100.4	100.4	100.4
2	156.9	156.9	156.9	156.9	156.9
3	105.8	105.8	105.8	105.8	105.8
4	147.5	147.5	147.5	147.5	136.2
4a	136.2	136.2	136.2	136.2	147.5
6	162.7	162.8	162.8	162.7	162.7
6a	114.5	114.4	114.5	114.5	114.5
7	115.6	115.5	115.6	115.6	115.6
8	148.6	148.7	148.6	148.6	148.7
9	154.6	155.3	154.6	154.6	154.6
10	108.6	108.6	108.6	108.6	108.6
10a	130.7	130.7	130.7	130.7	130.7
10b	120.7	120.7	120.7	120.7	120.7

Table 5

1H (800 MHz) NMR Spectroscopic Data of Cladophorols E–I (8–12) in CD3OD

Characterization of 9–12

Compounds 9–12 each exhibited 1H and 13C NMR spectra similar to that of 8 (Tables and 5). Proton integrations in the three overlapping regions (δH 6.85–6.98, 6.49–6.56, and 6.39–6.45) indicated three fewer protons for 9 and 3, 6, and 9 additional protons for 10–12. This suggested 9–12 were members of a series with 8, which was further supported by high-resolution mass spectrometry with [M – H]− ions at m/z 1107.1973, 1323.2394, 1431.2617, and 1539.2797 for 9–12, respectively. Since 10 was more abundant than the other members of the series, its NMR spectra in DMSO-d6 doped with picric acid were recorded (Figure S70). Three labile protons at δH 10.52, 10.16, and 10.11 correlated with δC 153.1 (C-9), 147.4 (C-8), 146.1 (C-4), 134.4 (C-4a), 114.5 (C-7), 107.9 (C-10), and 103.7 (C-3), indicating a substituted oxygen attached to C-2 as determined previously from the comparison of chemical shifts of 8 with experimental and calculated values of 13 and 14 (Figure S2). Altogether, these data supported the assignment of these natural products as α-[(4,8,9-trihydroxy-2-benzo[c]chromen-6-on-2-yl)oxy]-ω-[3,4-dihydroxyphenyl]hepta[oxy(2-hydroxyphen-4-yl)] (9), α-[(4,8,9-trihydroxy-2-benzo[c]chromen-6-on-2-yl)oxy]-ω-[3,4-dihydroxyphenyl]nona[oxy(2-hydroxyphen-4-yl)] (10), α-[(4,8,9-trihydroxy-2-benzo[c]chromen-6-on-2-yl)oxy]-ω-[3,4-dihydroxyphenyl]deca[oxy(2-hydroxyphen-4-yl)] (11), and α-[(4,8,9-trihydroxy-2-benzo[c]chromen-6-on-2-yl)oxy]-ω-[3,4-dihydroxyphenyl]undeca[oxy(2-hydroxyphen-4-yl)] (12) and named cladophorols F–I.

Pharmacological Assays

The antimicrobial activities of all isolated natural products (3–12) were assessed against multiple drug-resistant (MDR) Gram-positive and Gram-negative human pathogens as well as amphotericin-resistant Candida albicans (Table , Figure S3). Cladophorol C (6) was found to be very potent only against MRSA with a minimal inhibitory concentration (MIC) value of 1.4 μg/mL. The bioactivity is negatively correlated by the presence of a benzo[c]chromen-6-one unit, as shown by the weaker MIC values for 5 and 8–12 (13–49 μg/mL) relative to 6 and 7 (1.4–9 μg/mL)). The total number of phenol rings is also an important factor since the most potent compounds in both series (6–7 vs 8–12) have 9–10 rings and the activity weakens as this number increases. Surprisingly, all of the other cell lines, including the other Gram-positive bacterium, VREF, were much less sensitive to this family of natural products. The bioactivities of 4–12 were also evaluated in an assay using intraerythrocytic Plasmodium falciparum.[22] Only 4 and 5 inhibited the parasite with half-maximal effective concentration (EC50) values of 0.7 and 1.9 μg/mL, respectively. In addition, the cytotoxicities of 3–12 were evaluated using immortalized human keratinocytes (HaCaT) and human kidney cells (HEK293T). Although compounds 5–11 exhibited weak toxicity (<25%) at their respective maximum concentrations against HaCaT (Figure ) and 4 and 5 had weak toxicity toward HEK293T (half maximal cytotoxicity concentration (CC50) of 9.4 and 4.5 μg/mL, respectively), no toxicity was observed at the MRSA MICs for 5–11 or the P. falciparum EC50 values for 4 and 5.

Table 6

Pharmacological Activitiesa of 3,4,3′,4′-Tetrahydroxy-1,1′-biphenyl (3) and Cladophorols A–I (4–12)

		microbial strainsc
compd	no. of ringsb	MRSA	VREF	Gram–	WTCA	ARCA	P. falciparumd
3	2	>50	>50	>50	>50	>50	NT
4	2	>50	>50	>50	>50	>50	0.7 ± 0.1
5	4	13 ± 2	>50	>50	>50	>50	1.9 ± 0.2
6	9	1.4 ± 0.3	>50	>50	46 ± 13	46 ± 9	>2
7	10	9 ± 2	>50	>50	>50	>50	>2
8e	11	25 ± 3	>50	>50	>50	>50	>2
9	10	14 ± 3	>50	>50	>50	>50	>2
10	12	31 ± 3	>50	>50	>50	>50	>2
11	13	49 ± 8	>50	>50	>50	>50	>2
12e	14	>50	>50	>50e	>50	>50	>2

The values are expressed as a 95% confidence interval.

Total number of phenyl rings.

Minimal inhibitory concentration (μg/mL). Cell lines are MRSA = methicillin-resistant Staphylococcus aureus; VREF = vancomycin-resistant Enterococcus faecium; Gram–: [Escherichia coli; multiple drug-resistant (MDR) E. coli; MDR Klebsiella pneumonia; MDR K. pneumonia; MDR Acinetobacter baumannii; A. baumannii; MDR Enterobacter aerogenes; MDR Enterobacter cloacae; wild type Pseudomonas aeruginosa]; WTCA = wild type Candida albicans; ARCA = amphotericin B-resistant C. albicans. Details for the positive controls are provided in the Experimental Section.

Half-maximal effective concentration (μg/mL).

Tested only against E. coli.

Figure 7

Cytotoxicity of 3–12 against human keratinocytes (HaCaT). The parenthetical values provide the relative cytotoxicities (%) at the maximum concentration tested (60 μg/mL for 9, 50 μg/mL for 3–5, 8, 10, 11, 45 μg/mL for 12, and 20 μg/mL for 6 and 7).

Experimental Details

General Experimental Procedures

NMR spectra (1H, HMBC, HSQC, DQF-COSY, and J-resolved) were recorded on an 18.8 T (800 MHz for 1H and 200 MHz for 13C) Bruker Avance IIIHD 800 instrument equipped with a 5 mm triple resonance broadband cryoprobe. All spectra were acquired in CD3OD or DMSO-d6, and chemical shifts were reported in ppm (δ) relative to the residual solvent peaks (δH 3.31 and δC 49.00 for CD3OD, δH 2.50 and δC 39.52 for DMSO-d6). All spectra were processed using MestReNova 11.0. HPLC separations were performed with a Waters 1525 binary pump and a Waters 2487 dual wavelength absorbance detector set at 260 nm, using two different columns: a 9.4 × 250 mm, C18 silica reversed-phase (Zorbax stable-bond, 5 μm particle size), and a 4.6 × 250 mm phenylhexyl phase (Phenomenex Luna, 5 μm particle size). High-resolution mass spectrometry was conducted on an Orbitrap spectrometer in negative ion mode. The masses of isolated compounds were estimated by qNMR using a capillary filled with benzene-d6 and calibrated against caffeine.[23]

Specimen Collection

Green alga (collection G-1240) was harvested as floating clumps on Dec 3, 2015, at the ocean surface near Titi Island, Viti Levu, Fiji (16°16′26″ S, 179°26′02″ E). This collection was identified as Cladophora socialis Kützing by morphological and 18S rRNA phylogenetic analyses. Voucher specimens were preserved in aqueous formaldehyde and stored at the University of South Pacific. The collection was stored at −80 °C until extraction.

Species Identification

The collected green alga (G-1240) was identified by comparing its morphological traits with that of previously described Cladophora species[24] and by sequence analysis of nuclear, small subunit (SSU) rRNA (18S rRNA). Genomic DNA from an ethanol-preserved algal specimen was extracted using the innuPREP plant DNA kit (Analytik Jena, Germany) according to the manufacturer’s protocol. The three overlapping 18S rRNA gene fragments from genomic DNA were amplified via the polymerase chain reaction (PCR) in three separate reactions using universal primer pairs, NS3/NS4, NS5/NS6, and NS7/NS8.[25] Each PCR amplification was performed in a 25 μL reaction volume consisting of 5–50 ng of purified genomic DNA; 200 μM of each of the dNTPs; 1 μM of each of the oligonucleotide primer; 1.0 U Taq DNA polymerase, and 1× standard PCR reaction buffer (NEB, Ipswich, MA). All PCR amplifications were performed in a GeneAmp PCR system 2700 (Applied Biosystems, Foster City, CA) thermocycler using the following temperature cycling parameters: initial denaturation at 94 °C for 5 min followed by a total of 40 cycles of amplification in which each cycle consisted of denaturation at 94 °C for 40 s, primer annealing at 50 °C for 40 s and primer extension at 72 °C for 1 min. After amplification, final extension of the incompletely synthesized DNA was carried out at 72 °C for 7 min. The PCR fragments were analyzed by agarose gel electrophoresis (1% wt/vol). The gel was stained with ethidium bromide and visualized under a UV transilluminator. All of the PCR fragments were either sequenced with forward or reverse primers, and sequences were manually edited and assembled using CAP3 Sequence Assembly Program.[26] The assembled G-1240 18S rRNA sequence (1334 bp) was submitted to GenBank (accession no. MK127549). The sequence similarity of the assembled contig of G-1240 18S rRNA to other known Cladophora spp. was determined by comparing it with the nonredundant nucleotide database (NCBI) using the blastn program.[27] The top ranked matches according to E values and maximum scores revealed high similarity to C. socialis with 96% identity.

Further, the 18S rRNA sequence of G-1240 was compared with those of other known species of the genus Cladophora and two outgroup taxa Ulva curvata and Ulothrix zonata from the same class but different orders, obtained from GenBank. Phylogenetic analysis was conducted in MEGA X[28] using the Maximum Likelihood method based on the Tamura–Nei model[29] with 1000 bootstrap iterations. The phylogenetic analysis revealed that G-1240 is closely related to C. socialis (Figure S1). Overall, morphological and phylogenetic analyses are consistent with identification of the green alga G-1240 as C. socialis.

Isolation of 3–12

Frozen material (659 g) was thawed and broken into small pieces (∼1 cm) before being agitated in 60% aqueous MeOH (1.3 L) for 3 days. The mixture was filtered with a metal wire mesh, and the residues were extracted four more times with pure MeOH (2 × 1 L) and CHCl3/MeOH 1:1 (2 × 1 L). The solutions were pooled, filtered, and evaporated under reduced pressure to yield a brown gummy extract (33.4 g). The residue was dissolved in MeOH/H2O 9:1 (500 mL) and partitioned with hexanes (4 × 250 mL, 1.56 g). The amount of water in the MeOH/H2O phase was increased to a 3:2 ratio, and the resulting solution was partitioned with CHCl3 (3 × 500 mL, 0.75 g). The methanol was removed from the MeOH/H2O phase by vacuum evaporation, and the aqueous phase was extracted with EtOAc (3 × 250 mL, 2.05 g). The EtOAc fraction was subjected to silica gel column chromatography eluting with a gradient of CH2Cl2/MeOH from 49:1 to 1:2 to yield nine fractions (3.1–3.9). Fraction 3.5 (97.4 mg) was separated by reversed-phase semipreparative HPLC eluting with H2O + 0.05% TFA and MeCN (25 to 75% gradient, in 25 min). In addition to pure 3 (tR = 3.8 min, 0.50 mg), 4 (tR = 4.9 min, 0.68 mg), and 5 (tR = 8.7 min, 3.4 mg), eight other fractions were collected (3.5.1–3.5.8). Fraction 3.5.1 was subjected to HPLC using a phenylhexyl column and eluting with H2O + 0.05% TFA and MeOH (75 to 85% gradient, in 15 min) to yield 6 (tR = 7.9 min, 0.057 mg) and 11 (tR = 12.8 min, 0.62 mg). Fraction 3.5.2 was separated by HPLC with a phenylhexyl column eluting with H2O + 0.05% TFA and MeOH (75 to 85% gradient, in 15 min) to obtain 9 (tR = 9.5 min, 0.10 mg) and 12 (tR = 12.8 min, 0.027 mg). The same conditions were used to purify fraction 3.5.3 yielding 7 (tR = 9.0 min, 0.32 mg). Fraction 3.5.4 was purified by HPLC using a phenylhexyl column and eluting with H2O + 0.05% TFA and MeOH (80 to 82.5% in 7.5 min) to obtain 8 (tR = 6.5 min, 0.60 mg). Fraction 3.5.6 was subjected to HPLC using a phenylhexyl column and eluting with H2O + 0.05% TFA and MeOH (20 to 82.5% gradient, in 7.5 min) to yield 10 (tR = 6.4 min, 2.3 mg).

Cladophorol A (4):

pale red amorphous solid; 1H and 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C13H7O6– 259.0248, found 259.0245.

Cladophorol B (5):

pale red amorphous solid; 1H and 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C25H15O9– 459.0722, found 459.0718.

Cladophorol C (6):

pale red amorphous solid; 1H and 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C54H37O18– 973.1985, found 973.1967.

Cladophorol D (7):

pale red amorphous solid; 1H and 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C60H41O20– 1081.2197, found 1081.2177.

Cladophorol E (8):

pale red amorphous solid; 1H NMR data, see Table ; 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C67H43O23– 1215.2201, found 1215.2188.

Cladophorol F (9):

pale red amorphous solid; 1H NMR data, see Table ; 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C61H39O21– 1107.1989, found 1107.1973.

Cladophorol G (10):

pale red amorphous solid; 1H NMR data, see Table ; 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C73H47O25– 1323.2412, found 1323.2394.

Cladophorol H (11):

pale red amorphous solid; 1H NMR data, see Table ; 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C79H51O27– 1431.2623, found 1431.2617.

Cladophorol I (12):

pale red amorphous solid; 1H NMR data, see Table ; 13C NMR data, see Table ; HRMS (ESI) m/z [M – H]− calcd for C85H55O29– 1539.2834, found 1539.2797.

Computational Methods

Model compounds (1, 13, 14, 15a, and 15b) were subjected to a Monte Carlo conformational search using the MMFF94 force field in Spartan ‘16.[30] The geometry of all conformers within 21 kJ/mol were optimized using density functional theory (DFT) in the Gaussian D09 package,[31] using the B3LYP/6-31G(d) level of theory. Solvent interactions, either DMSO (1, 13, 14) or MeOH (15a, 15b), were considered with PCM. Frequencies computation was performed to ascertain true minima were obtained (zero imaginary frequency), and the resulting free energies were extracted to calculate the Boltzmann distribution. All conformers representing more than 2% of the population were considered for the computation of the shielding tensors which was performed at the mPW1PW91/6-311+G(d,p) level of theory. The shielding tensors (GIAO) were averaged using their respective Boltzmann weight, and the chemical shifts were finally derived by scaling the shielding tensors against the experimental values.

Pharmacological Assays

Antibacterial and antifungal assays were performed as previously described[32] using methicillin-resistant Staphylococcus aureus (MRSA, ATCC 33591), vancomycin-resistant Enterococcus faecium (VREF, ATCC 700221), Escherichia coli (EC, ATCC 25922), MDR Escherichia coli (MDREC, ATCC BAA-1743), wild type Candida albicans (WTCA, ATCC 32354), and amphotericin B-resistant C. albicans (ARCA, ATCC 90873). The positive controls used were vancomycin for MRSA (MIC < 4 μg/mL), chloramphenicol for VREF (MIC < 8 μg/mL) and EC (MIC < 8 μg/mL), nitrofurantoin for MDREC (MIC = 14 ± 1 μg/mL), amphotericin B for WTCA (MIC = 16 ± 1 μg/mL), and nystatin/cycloheximide for ARCA (MIC > 500 μg/mL). Activities were expressed as a 95% confidence interval of concentration required for at least 90% inhibition of growth (MIC90).

Additional antibacterial assays were performed using MDR Klebsiella pneumonia (MDRKP1, CDC0087), MDR K. pneumonia (MDRKP2, CDC0003), MDR Acinobacter baumannii (MDRAB1, CDC0280), MDR A. baumannii (MDRAB2, CDC0033), MDR Enterobacter aerogenes (MDREA, CDC0018), MDR Enterobacter cloacae (MDRECL, CDC0050), and wild type Pseudomonas aeruginosa (WTPA, PAO1). The bacterial cultures were grown in tryptic soy broth (TSB). Minimum inhibitory concentrations were evaluated with cation-adjusted Mueller Hinton broth (CAMHB) following Clinical and Laboratory Standards Institute (CLSI) methods.[33] Overnight cultures were standardized by optical density (OD) to 5 × 105 CFU/mL, and this was confirmed by plate counts of colonies. MIC values was assigned as described;[34] this was determined by reading at an OD600 nm in a Cytation 3 multimode plate reader after 18–22 h incubation. Gentamicin was used as a positive control.

The human keratinocyte cell line (HaCaT) was maintained, and cytotoxicity of isolated compounds were assessed using LDH cytotoxicity assay as previously described.[34]

Antimalarial activity was assessed using a standard parasite proliferation assay with asexual blood-stage P. falciparum and SYBR Green detection.[22] Screening media (complete media supplemented with 0.05% Albumax II but without human serum) and fresh O+ erythrocytes (TSRI Normal Blood Donation) were used for cultures of P. falciparum strain Dd2L (gift from Dr. David Fidock, Colombia University). A Labcyte ECHO acoustic liquid handler transferred compounds into assay plates. The plates were inoculated (Multi-Flo; BioTek, VT) with parasitized erythrocytes and fresh erythrocytes prepared in screening media, resulting in a final parasitemia of 0.3% and hematocrit of 2.5%. A chamber with a low oxygen gas mixture at 37 °C was used to culture the assay plates for 72 h with daily gas exchanges. SYBR Green lysis buffer was added to the wells after incubation using a Multi-Flo liquid dispenser, and plates were incubated at room temperature for an additional 24 h to achieve optimal development of the fluorescence signal which was read by an Envision Multimode Reader (PerkinElmer, MA). The positive controls used were atovaquone (EC50 = 0.0006 μM) and artemisinin (EC50 = 0.03 μM).