Skip to main content

ORIGINAL RESEARCH article

Front. Mol. Biosci., 06 June 2023
Sec. Metabolomics

Discovery of anti-SARS-CoV-2 secondary metabolites from the heartwood of Pterocarpus santalinus using multi-informative molecular networking

Andreas Wasilewicz,&#x;Andreas Wasilewicz1,2Julia Zwirchmayr&#x;Julia Zwirchmayr1Benjamin Kirchweger&#x;Benjamin Kirchweger1Denisa Bojkova&#x;Denisa Bojkova3Jindrich Cinatl Jr.&#x;Jindrich Cinatl Jr.3Holger F. Rabenau&#x;Holger F. Rabenau3Judith M. Rollinger&#x;Judith M. Rollinger1Mehdi A. Beniddir&#x;Mehdi A. Beniddir4Ulrike Grienke
Ulrike Grienke1*
  • 1Division of Pharmacognosy, Department of Pharmaceutical Sciences, Faculty of Life Sciences, University of Vienna, Vienna, Austria
  • 2Vienna Doctoral School of Pharmaceutical, Nutritional, Sport Sciences, University of Vienna, Vienna, Austria
  • 3Institute of Medical Virology, University Hospital Frankfurt, Frankfurt am Main, Germany
  • 4Équipe Chimie des Substances Naturelles, BioCIS, Centre National de la Recherche Scientifique (CNRS), Université Paris-Saclay, Orsay, France

A pigment-depleted extract from the heartwood of Pterocarpus santalinus L. f. (PS-DE) showed promising anti-SARS-CoV-2 activity with an IC50 of 29.9 μg/mL in Caco-2-F03 cells. To determine the potential active constituents within the extract prior to isolation, multi-informative molecular network (MN) was applied. Therefore, the extract was separated by high-performance counter-current chromatography (HPCCC) into 11 fractions which were subsequently tested for anti-SARS-CoV-2 activity and analysed by UPLC-tandem mass spectrometry (MS2). The resulting MN combines the bioactivity data of the fractions with the MS2 data. The MN analysis led to the targeted isolation of seven compounds including one pterocarpan (7) reported for the first time as constituent of P. santalinus and four so far undescribed natural products (NPs) that belong to the compound classes of arylpropanes (9), isoflavanones (10) coumestans (16) and 3-arylcoumarins (17), respectively. In total, 15 constituents from the heartwood of P. santalinus and one synthetic isoflavonoid that is structurally related to the natural metabolites were tested for anti-SARS-CoV-2 activity. Thereby, the two pterocarpans (−)-homopterocarpin (5) and (−)-medicarpin (2), the stilbene (E)-pterostilbene (1) and the isoflavonoid 7-O-methylgenistein (11) showed a distinct antiviral activity with IC50 values of 17.2, 33.4, 34.7, and 37.9 µM, respectively, and no cytotoxic effects against Caco-2-F03 cells (CC50 > 100 µM). In addition, a structure-activity relationship (SAR) was proposed indicating structural requirements of pterocarpans for anti-SARS-CoV-2 activity. The herein presented results support the implementation of multi-informative molecular networks as powerful tool for dereplication and targeted isolation of bioactive NPs.

1 Introduction

Pterocarpus santalinus L. f. (Fabaceae), also known as red sandalwood, grows as small tree native to southern India. The red-coloured heartwood is used in Indian Ayurvedic medicine to treat various ailments such as headache, fever and chronic inflammatory disorders including chronic bronchitis (Dahat et al., 2021). In Europe, the heartwood is occasionally included in herbal detox tea preparations and as brightening agent (Wichtl, 2002).

In recent years, several in vitro studies have reported anti-inflammatory properties for extracts and constituents of P. santalinus heartwood using different cell model systems (Cho et al., 2001; Wu et al., 2011a; Wu et al., 2011b; Ham et al., 2019; Natalia et al., 2021; Zwirchmayr et al., 2023). Moreover, a heartwood extract also showed promising anti-influenza activity (IC50 of 12 μg/mL against influenza virus A in MDCK cells (CC50 of 54 μg/mL)) (Grienke et al., 2018). Two constituents from P. santalinus heartwood belonging to different compound classes were reported to have antiviral effects against coronaviruses (CoV). The lignan savinin exerted potent anti-SARS-CoV activity in Vero E6 cells with an EC50 in the low micromolar range (Wen et al., 2007). The activity of savinin was explained by the competitive inhibition of the SARS-CoV main protease (Mpro, 3CL protease) which is essential for viral replication. Recently, the stilbene pterostilbene (1) has shown anti-SARS-CoV-2 activity in Vero E6 cells and in primary human bronchial epithelial cells under air-liquid interface cell culture conditions (ter Ellen et al., 2021). To date, no studies investigated the anti-SARS-CoV-activity of multi-component mixtures from P. santalinus heartwood.

Molecular networking, introduced by (Watrous et al., 2012), has rapidly emerged as an established and efficient dereplication tool to streamline the isolation of natural products (NPs) of interest from complex mixtures (Fox Ramos et al., 2019; Beniddir et al., 2021). The principle of molecular networking is based on the alignment and visualization of tandem mass spectrometry (MS2) data of biological samples such as extracts, fractions or pure compounds in a molecular network (MN) (Aron et al., 2020). Thereby, detected compounds are presented as “nodes.” The nodes of structurally related compounds with similar MS2 spectra are connected by “edges,” leading to the formation of molecular families. By comparison of the MS2 data from the MN with reference MS2 spectra data, nodes of known compounds can be annotated. These structural annotations can be used as starting point to annotate neighbour nodes via the calculation and interpretation of mass differences as well as literature. In recent years, the integration of bioactivity data into MNs has been successfully applied to identify bioactive NPs prior to their isolation (Olivon et al., 2017; Nothias et al., 2018). Thereby, the arrangement of structural and biological data in a MN gives certain advantages for the isolation of bioactives compared to classical bioassay-guided fractionation. While the targeted isolation of unknown bioactives can be prioritized, the re-isolation of NPs with already known bioactivity can be avoided.

Due to i) the traditional use against cold-related symptoms, ii) the described anti-inflammatory and anti-influenza properties, and iii) the reported anti-SARS-CoV activity of known constituents, the heartwood of P. santalinus was considered as promising natural source for the search of NPs with so far unknown anti-SARS-CoV-2 activity. In this study, a multi-informative MN was applied for the targeted isolation of NPs with anti-SARS-CoV-2 activity from the heartwood of P. santalinus.

2 Results and Discussion

2.1 Assessment of anti-SARS-CoV-2 activity of Pterocarpus santalinus heartwood

A crude extract (PS-E) and a pigment-depleted extract (PS-DE) were prepared from the P. santalinus heartwood. The anti-SARS-CoV-2 activity of PS-E and PS-DE was first evaluated with a pre-screening assay which is based on the measurement of caspase 3/7 activity of cells in response to SARS-CoV-2 infection (Bojkova et al., 2023). This assay represents an indirect readout method based on the measurement of the virus-induced caspase 3/7 activity of cells 48 h after SARS-CoV-2 infection. The results were then confirmed by means of a spike protein immunostaining assay in which the amount of spike protein was quantified 24 h after SARS-CoV-2 infection. For both assays, the SARS-CoV-2 strain G614 and Caco-2-F03 cells were used. The crude extract PS-E showed promising activity in the pre-screening assay but showed no activity in the spike immunostaining assay (IC50 > 100 μg/mL). The pigment-depleted extract PS-DE was determined as active in both assays showing an IC50 value of 29.9 μg/mL in the spike immunostaining assay. PS-DE was separated by high-performance counter-current chromatography (HPCCC+) and pooled into eleven fractions (F1—F11) according to TLC fingerprints (Supplementary Figure S1) (Zwirchmayr et al., 2023). All fractions were tested for anti-SARS-CoV-2 activity (Table 1). F2—F9 showed distinct activity in both assays (>60% of caspase 3/7 inhibition at 50 μg/mL and IC50 < 50 μg/mL in the spike immunostaining assay), F10 showed weak activity in the spike immunostaining assay. F1 and F11 were inactive in the pre-screening assay and were not further tested. Based on these findings, it was assumed that the P. santalinus heartwood extracts contain constituents with anti-SARS-CoV-2 activity but also metabolites that interfere with the used pre-screening assay. A false positive effect was confirmed for the crude extract PS-E. Other compounds have been already reported to give similar false positive effects in the used assay, i.e., the caspase inhibitor emricasan supressed SARS-CoV-2-induced caspase 3/7 activity but did not affect SARS-CoV-2 replication (Bojkova et al., 2023).

TABLE 1
www.frontiersin.org

TABLE 1. Anti-SARS-CoV-2 activity and Caco-2-F03 cell cytotoxicity of Pterocarpus santalinus heartwood extracts and fractions.

2.2 Molecular network generation and mapping

Each fraction (F1—F11) was subjected to UHPLC-MS2 (Supplementary Figures S2, S3). The generated MS2 data along with the bioactivity data (metadata table) were then uploaded to GNPS for MN generation (Aron et al., 2020). Figure 1 shows the entire MN. In order to visualize the distribution of the constituents over the fractions, a colour scale was used. Thereby, lipophilic fractions are coloured in yellow and green, while hydrophilic fractions are depicted in different shades of blue (Figure 2A). For the visualisation of bioactivity, the fractions were divided into three color-coded activity levels. Fractions with an IC50 < 10 μg/mL in the spike immunostaining assay were determined as active (coloured in fiery red) while fractions with an IC50 from 20 to 50 μg/mL in the spike immunostaining assay were rated as moderately active (coloured in ruby red). Fractions with IC50 > 50 μg/mL in the spike immunostaining assay or caspase 3/7 inhibition of <60% at 50 μg/mL were classified as inactive (coloured in black red). Hence, fractions F2 and F3 were categorized as active, F4—F9 were determined as moderately active and F1, F10 and F11 were classified as inactive. The size of the nodes reflects the relative quantity of a detected compound based on MS ionization (the higher the intensity, the bigger the node). The ionization rate of most of the constituents aligns well with the detected abundance in the UPLC-ELSD chromatograms (Supplementary Figures S4–S7). Structural network annotation was performed by comparison with i) previously isolated compounds from P. santalinus heartwood (orange node border) (Zwirchmayr et al., 2023), ii) GNPS spectra library (hexagon-shaped node) and iii) literature search (red node border).

FIGURE 1
www.frontiersin.org

FIGURE 1. Workflow used for the generation of the global molecular network based on UHPLC-MS2 analyses of fractions F1-11. MF, molecular family.

FIGURE 2
www.frontiersin.org

FIGURE 2. Molecular family A (MF-A) with compound 1 and molecular family B (MF-B) showing compounds 210. (A) MF-A and MF-B with fraction-based mapping indicating the distribution of the compounds over F1—F11. (B) MF-A and MF-B with bioactivity-based mapping displaying the presence of the compounds in active, moderately active and inactive fractions. The chemical structures are shown for the identified compounds.

2.3 Molecular network analysis and targeted isolation

The generated MN consists of 1,957 nodes including 23 molecular families (MF; >3 nodes). The nodes of five MFs (MF-A—MF-E) were explored in detail (Supplementary Table S1). Compound 1 [(E)-pterostilbene, m/z 257.1172 [M + H]+] with known anti-SARS-CoV-2 activity was found in MF-A. According to UPLC-ELSD, it is the major constituent of the moderately active fraction F5 (Figure 2). Notably, the lignan savinin, the second known constituent from P. santalinus heartwood with reported anti-SARS-CoV activity (Wen et al., 2007), could not be detected in any extract or fraction. Hence, it was quickly determined that the heartwood contained additional NPs with so far unknown anti-SARS-CoV-2 activity.

MF-B shows a dumbbell-shaped structure. The top part of MF-B consists of compounds of the structure class of pterocarpans such as 24. Compounds 2 [(−)-medicarpin, m/z 271.0966 [M + H]+], 3 [(−)-maackiain, m/z 285.0753 [M + H]+] and 4 [(−)-3,8-dihydroxy-9-methoxypterocarpan, m/z 287.0912 [M + H]+] were detected in F6, F7 and F10, respectively. Compounds 57, which were all connected to each other and were the largest nodes in the top part of MF-B, remained unannotated after comparison with the previously isolated compounds. Compound 5 (m/z 285.1122 [M + H]+), annotated by the GNPS spectra library as homopterocarpin, and compound 6 (m/z 299.0912 [M + H]+) were of highest interest due to their presence in the active fractions F2 and F3 as major constituents determined by UPLC-ELSD. It should be noted that compound 5 was also detected in the inactive fraction F1 as indicated by the black red part of the chart. However, its relative abundance in F1 estimated by MS is overrated, since the major constituents of this fraction, presumably triterpenes, are difficult to ionize. Compound 7 (m/z 301.1074 [M + H]+) was found as major constituent in the moderately active fraction F6. Targeted isolation followed by NMR analysis and optical rotation experiments confirmed that (−)-homopterocarpin (5), (−)-pterocarpin (6) and (−)-8-hydroxyhomopterocarpan (7) belong to the structure class of pterocarpans (Figure 3). Compounds 5 and 6 are well-known constituents from P. santalinus, while compound 7 has been previously isolated from Pterocarpus soyauxii (Bezuidenhoudt et al., 1987). In the bottom part of MF-B, the isoflavanone 8 [(−)-pterosantalin C, m/z 333.0970 [M + H]+] was the only annotated compound. The three largest nodes in the bottom part of MF-B, the compounds 9 (m/z 303.1592 [M + H]+), 10 (m/z 331.1177 [M + H]+) and m/z 347.1130 [M + H]+, were investigated in detail. Compound 9 was the major component in F4 and was selected for targeted isolation. The 1H NMR spectrum of compound 9 showed six signals of aromatic protons [δH 6.41 (dd, J = 8.2, 2.4), δH 6.44 (d, J = 2.3), δH 6.69 (d, J = 1.7), δH 6.69 (dd, J = 8.4, 1.9) δH 6.82 (d, J = 8.5) δH 7.02 (d, J = 8.2)], three signals assignable to methoxy groups [δH 3.79 (6H, s), δH 3.87 (3H,s)] and three signals of methylene groups [δH 1.85 (2H, m), δH 2.58 (4H, t, J = 7.7)]. These chemical shifts corresponded well with signals determined for NPs with a diphenylpropane scaffold (Pan et al., 2010). The substitution pattern of ring A was assigned by the HMBC correlation from H-6′ (δH 7.02) to C-1 (δC 29.40), C-2′ (158.46) and C-4′ (δC 159.17), while the substitution pattern of ring B was determined by the HMBC correlation from H-5′′ (δH 6.82) to C-1′′ (δC 134.83), C-3′′ (δC 146.38), C-4′′ (δC 143.59) and C-6′′ (δC 121.09). Eventually, compound 9 was identified as new NP belonging to the class of arylpropanes (i.e. 1,3-diphenylpropanes) and was named pterosantalin D (Table 2). Compound 10 in F8 and m/z 347.1130 in F5 were only detected in traces by UPLC-ELSD and were not further considered for isolation. However, compound 10 was co-isolated during targeted isolation of compound 16. It was identified as an isoflavanone which was previously known as synthetic compound and was named (−)-pterosantalin E (10) (Maranduba et al., 1979). Due to direct connections to compounds 8 and 10, m/z 347.1130 was supposed to be an isoflavonone substituted by three methoxy and two hydroxy groups with the molecular formula of C18H18O7. The relatively large node size of m/z 347.1130 can be explained by its higher ionization rate in comparison to the major constituent, compound 1, in F5.

FIGURE 3
www.frontiersin.org

FIGURE 3. Chemical structures of compounds 123.

TABLE 2
www.frontiersin.org

TABLE 2. 1H (500 MHz) and 13C (125 MHz) NMR Data of Pterosantalin D (9) (CDCl3, δ in ppm).

MF-C mostly contains isoflavones such as compounds 1115 (Figure 4) which have been previously isolated from P. santalinus heartwood (Zwirchmayr et al., 2023). Compounds 11 [7-O-methylgenistein, m/z 285.0761 [M + H]+] and 12 [7,3′-di-O-methylorobol, m/z 315.0666 [M + H]+] were found in F7 and were the major constituents of this fraction according to UPLC-ELSD analysis. Compound 13 [santal, m/z 301.0706 [M + H]+] was detected as the most abundant compound in fraction F9. Compounds 14 [khrinone C, m/z 331.0813 [M + H]+] and 15 [pterosantalin B, m/z 361.0918 [M + H]+] were mostly found in F10. Four additional nodes of minor compounds (m/z 269.0811 [M + H]+, 299.0915 [M + H]+ (RT: 6.30 min), 299.0916 [M + H]+ (RT: 5.82 min), 345.0975 [M + H]+ were putatively annotated as isoflavones which differ in the number and position of hydroxy and methoxy groups. Compounds 16 (m/z 297.0758 [M + H]+) and 17 (m/z 299.0916 [M + H]+) were not dereplicated prior to isolation. Due to their occurrence in active and moderately active fractions, compounds 16 (in F3 and F4) and 17 (in F8) were both selected for targeted isolation. Since neither compound was directly connected to the annotated isoflavonoids, it was expected that compound 16 and 17 belong to structurally related compound classes. Eventually, isolation and subsequent NMR analysis identified compound 16 as a coumestan named 3,9-O-dimethylcoumestrol and 17 as a 3-arylcoumarin named pterosantalin F. This is the first report of those previously synthesized compounds as genuine compounds (Harper et al., 1969; Naik et al., 2017).

FIGURE 4
www.frontiersin.org

FIGURE 4. Molecular family C (MF-C) including compounds 1117. (A) MF-C with fraction-based mapping indicating the distribution of the compounds over F1—F11. (B) MF-C with bioactivity-based mapping displaying the presence of the compounds in active, moderately active and inactive fractions. The chemical structures are shown for the identified compounds.

In addition, MF-D that only consists of metabolites from the inactive F11 was roughly characterised (Figure 5). Thereby, four major pigments reported for P. santalinus heartwood, compound 18 (santalin A, m/z 583.1602 [M + H]+), 19 (santalin B, m/z 597.1761 [M + H]+), 20 (santarubin A, m/z 611.1911 [M + H]+), and 21 (santarubin B, m/z 597.1755 [M + H]+) were tentatively annotated. F11 was the fraction with by far the highest determined yield (Supplementary Table S2). Hence, the lack of bioactivity of the crude extract PS-E in the spike assay could likely be due to the high content of non-active pigments occurring as bulk constituents.

FIGURE 5
www.frontiersin.org

FIGURE 5. Molecular family D (MF-D) with compounds 1821 and molecular family E (MF-E). (A) MF-D and MF-E with fraction-based mapping indicating the distribution of the compounds over F1—F11. (B) MF-D and MF-E with bioactivity-based mapping displaying the presence of the compounds in active, moderately active and inactive fractions. The chemical structures of identified compounds are shown.

The major constituent of F8 according to UPLC-ELSD was found in MF-E. The relatively small node size with m/z 541.2222 can be explained by the poor ionization rate. Since all nodes from MF-E remained unannotated, only vague assumptions concerning the potential compound class were made. If m/z 541.2222 is assumed to be an [M + H]+ ion, and considering the high abundance of isoflavonoids and related compound classes in the heartwood, this compound could possibly be a bi-isoflavonoid with the molecular formula C30H20O10. Bi-isoflavonoids have been previously found in other plants belonging to the Fabaceae family such as in Podocarpium and Dalbergia species (Yahara et al., 1985; Ma et al., 2013). Targeted isolation of m/z 541.2222 was conducted but could not be achieved due to the presence of a potential isomer. After unsuccessful separation attempts using orthogonal techniques such as HPCCC+ and size-exclusion chromatography (SEC), the isolation was not further pursued.

2.4 Anti-SARS-CoV-2 activity of pure compounds

In total, 15 compounds were subjected to the anti-SARS-CoV-2 spike immunostaining assay (Table 3) (17, 9, 1114, 16, 22, 23). The minor compounds 8 and 17 as well as compound 10 were not tested since only low yields were isolated which were insufficient for testing. Based on the network analysis, compounds 5 and 6 were expected to show distinct bioactivity, compounds 1, 2, 3, 7, 9, 11, 12, 13 and 16 were predicted to be moderately active, and compounds 4 and 14 were supposed to be inactive. The activity of compound 16 could not be assessed since it was only a minor compound in F8. The isoflavane 22 ((−)-isosativan) was not present in the MN but was found in F4 and was selected for testing due to its structural resemblance to compound 5. In addition, the synthetic compound 23 ((+)-7,4′-di-O-methyldihydrodaidzein) which is the isoflavanone structural analogue of compound 5 was included.

TABLE 3
www.frontiersin.org

TABLE 3. Anti-SARS-CoV-2 activity and Caco-2-F03 cell cytotoxicity of pure compounds.

Eventually, compound 5 was determined as the most active constituent with an IC50 of 17.2 µM. Moreover, compounds 1 (IC50 of 34.7 µM), 2 (IC50 of 33.4 µM), 11 (IC50 of 37.9 µM), 22 (IC50 of 40.8 µM) and 23 (IC50 of 53.3 µM) showed moderate anti-SARS-CoV-2 activity. The determined activity of compound 1 in this study is in line with the previously reported activity (ter Ellen et al., 2021). Compound 13 (IC50 of 74.2 µM) showed weak anti-SARS-CoV-2 activity. As expected, compounds 4 and 14 were determined as inactive as well as compounds 3, 6, 7, 9, 12 and 16 (IC50 > 100 µM). In summary, the bioactivities of F2 and F3 (containing compound 5), F4 (containing compound 22), F5 (containing compound 1), F6 (containing compound 2) and F7 (containing compound 11) were pinpointed to specific constituents, while the activity of F8 remained elusive. The activity of the mother extract PS-DE was mainly attributed to compounds 1 and 5, as they were determined to be the major constituents in the extract along with the inactive compounds 6, 13 and the red-coloured pigments.

To evaluate the potential molecular target involved in the phenotypic activity, compound 5 with the strongest anti-SARS-CoV-2 activity was tested for SARS-CoV-2 main protease (Mpro) inhibition using an enzyme-based FRET assay (Wasilewicz et al., 2023). However, compound 5 did not show any inhibitory activity against Mpro at 50 µM when compared to the DMSO control (data not shown, method in Supplementary Material). Hence, Mpro can be excluded as the molecular target for compound 5.

2.5 Structure-activity relationship of pterocarpanoids and anti-SARS-CoV-2 activity

By comparing the bioactivity from the spike assay of the most active compound 5 with compounds 24, 6, 7, 22 and 23, a structure-activity relationship (SAR) for pterocarpan compounds was proposed (Figure 6). Thereby, three structural requirements were determined. The substitution of position C-3 by a methoxy group resulted in a 2.2-fold increase of bioactivity compared to the substitution with a hydroxyl moiety as seen in compound 2. Furthermore, any type of position C-8 substitution led to a complete abolishment of bioactivity as it was the case for compounds 3, 4, 6 and 7. Finally, the ether group in position C-11 plays an essential role for bioactivity. Compound 11 possesses a 3.0-fold higher activity than its corresponding isoflavone 23 with a keto group and is 2.2-fold more active than the isoflavane 22 whose hydroxy group aligns well with the ether group according to 3D structural alignment (Supplementary Figure S8).

FIGURE 6
www.frontiersin.org

FIGURE 6. Structure-activity relationship chart of compound 5 and related compounds 24, 6, 7, 22 and 23 for anti-SARS-CoV-2 activity. Structural changes compared to compound 5 at position C-3, C-8 and C-11 are depicted by red, green and blue spheres, respectively.

In the present study, bioactivity-informed molecular networking was successfully applied to identify NPs with anti-SARS-CoV-2 activity from P. santalinus heartwood prior to isolation. Instead of laborious iterative fractionation and testing cycles needed for bioassay-guided isolation approaches, seven compounds were isolated in a target-oriented manner. Molecular networking not only revealed compounds with potential bioactivity but also facilitated the dereplication of known major constituents from the heartwood. Based on the structural annotations, MF-B and MF-C were well characterised revealing the heartwood as a rich source of isoflavonoids and pterocarpan compounds. However, MFs without structurally annotated compounds such as MF-E are difficult to explore. Since there are several MFs without any structural annotations, the heartwood of P. santalinus remains a promising source to find so far undescribed NPs.

The findings support the ethnopharmacological use of P. santalinus heartwood in Ayurvedic medicine against cold-related symptoms. Compound 5 (IC50 of 17.2 µM) was determined as the most active constituent and is assumed to be mainly responsible for the bioactivity shown by the investigated pigment-depleted extract PS-DE (IC50 of 29.9 μg/mL) alongside compound 1 (IC50 of 34.7 µM). In addition, compounds 2, 11, 22 and 23 have shown moderate activity ranging from 33.4 to 53.3 µM. All bioactive NPs were well tolerated by Caco-2-F03 cells (CC50 > 100 µM). For the evaluation of the anti-SARS-CoV-2 activity, an in vitro phenotypic assay was performed. Compound 5 had no inhibitory activity against the commonly tested SARS-CoV-2 target Mpro. Recent results from time-of-addition assays with compound 1 suggested that the anti-SARS-CoV-2 activity occurs in post-infection conditions within the cells (ter Ellen et al., 2021). The exact mechanism of action and potential molecular targets of compounds discovered in this study are still unknown and require further investigations. Nevertheless, the chemical scaffold of pterocarpan compounds presents a promising starting point for the development of more potent compounds with anti-SARS-CoV-2 activity.

3 Material and methods

3.1 General experimental procedures

The UPLC-ELSD analysis was performed on a Waters Acquity UPLC H-class system comprising a quaternary solvent manager, sample manager, column manager, PDA detector and ELSD. The instrument was controlled by the software Empower 3. An HSS T3 column (1.7 µm, 2.1 × 100 mm, Waters) was used at 40°C as stationary phase, while the mobile phase consisted of water (A) and acetonitrile (B). The flow rate was set to 0.25 mL/min. The following gradient was used: 30% B isocratic for 0.5 min, 30%–45% B in 1.5 min, 45%–50% B in 1.5 min, 50%–55% B in 8.5 min, 55%–98% B in 0.1 min, 98% B for 5.9 min.

Flash chromatography (FC) was performed on PuriFlash 4250 from Interchim equipped with a PDA detector and ELSD and was controlled by the Interchim software. High-performance counter-current chromatography (HPCCC) was performed using a HPCCC+ device (Zwirchmayr et al., 2023). Therefore, the PuriFlash instrument was hyphenated to a Spectrum HPCCC instrument from Dynamic Extraction Ltd., whereby the HPCCC device replaces a solid stationary column. An Accel 500 LC chiller from Thermo Scientific was used for cooling.

TLC was carried out using Merck silica gel 60 PF254 plates as stationary phase. The mobile phase consisted of chloroform:ethyl acetate:formic acid (8:1:1). Detection was conducted under visible light, UV254 and UV366 before and after derivatisation with vanillin/sulfuric acid (5% in methanol).

SEC was performed with Sephadex LH-20 as stationary phase and methanol as mobile phase (column size: 100 cm × 2 cm).

NMR experiments were conducted on a Bruker Advance 500 NMR spectrometer (UltraShield) equipped with a TCI Prodigy CryoProbe (5 mm, triple resonance inverse detection probe head). NMR spectra were analysed by using TopSpin 4.1.4.

The optical rotation was measured at 25°C in methanol or chloroform using an MCP 100 polarimeter.

3.2 Plant material

The P. santalinus heartwood was purchased from Kottas Pharma GmbH, Vienna, Austria (batch no. P16301836). A voucher specimen (JR-20190315-A1) is deposited at the Division of Pharmacognosy, Department of Pharmaceutical Sciences, University of Vienna, Austria.

3.3 Generation of PS-E and PS-DE

First, a heartwood crude extract (PS-E) was generated as previously described (Zwirchmayr et al., 2023). 1 kg of the dried, pulverized heartwood was defatted with 2 L of n-hexane for 3 days. Then, the remaining material was extracted three times with 4 L of dichloromethane for 3 days. The extraction procedure was repeated with methanol. The dichloromethane and methanol extracts were combined and dried under vacuum to obtain the extract PS-E. For the pigment depletion, liquid-liquid extraction (LLE) was performed using dichloromethane and 50% aqueous methanol (1:1, v/v) as biphasic solvent system. The LLE process of pigment depletion was chosen because it is considered as the most efficient method for the removal of bulk constituents (Picker et al., 2014). The solvent system was selected as the most suitable of six systems tested in preliminary trials monitored by TLC. Approximately 85 g of PS-E were separated by 16 cycles of LLE (500 mL of each phase). The combined dichloromethane phases were evaporated to dryness yielding 67.3 g of the pigment-depleted extract (PS-DE).

3.4 Generation of F1—F11

PS-DE fractionation was performed using two semi-preparative runs of HPCCC+ in normal phase mode at a flow rate of 6 mL/min and rotation speed of 1,600 rpm. About 800 mg of PS-DE was fractionated using the HEMWat biphasic solvent system, which consists of several mixtures with different proportions of hexane, ethyl acetate, methanol and water (Supplementary Table S3). The lower layer (LL) of HEMWat system 19 (300 mL) was used as the stationary phase, while the upper layers (UL) of HEMWat systems 19, 18, 17, 16 and 14 (250 mL of each layer) were successively used as mobile phase. For a stepwise gradient elution, the following procedure was used. When 200 mL of a mobile phase was consumed, the residual 50 mL were blended with 50 mL of the next mobile phase. After the consumption of 50 mL of this mixture, the remaining volume of the next mobile phase was added. 6 mL per fraction were collected. After the appearance of the first red coloured fractions, elution extrusion was performed using 50% aqueous methanol at a flow rate of 10 mL/min and a rotation speed of 200 rpm. The elution extrusion was collected as a whole. In total, 392 fractions were collected and analysed by TLC resulting in 11 pooled fractions (F1—F11), with the final fraction being the elution extrusion.

3.5 Targeted isolation of Pterocarpus santalinus constituents and pure compounds

A fractionation tree is provided in Supplementary Figure S9. Approximately 5 mg of F6 were subjected to SEC to yield three fractions (PS-DE_F6a-c). PS-DE_F6a was identified as compound 7 (0.73 mg; (−)-8-hydroxyhomopterocarpan) (Bezuidenhoudt et al., 1987).

PS-DE was fractionated using FC as described previousely (Zwirchmayr et al., 2023). Briefly, about 10 g of PS-DE were applied via dry load on a Silica HC 200 G column (stationary phase). Petroleum ether (A), dichloromethane (B) and methanol (C) were selected as mobile phase using the following gradient: 80% A/20% B isocratic for 13 min, 80% A/20% B—100% B in 3 min, 100% B isocratic for 16 min, 100% B—50% B/50% C in 45 min, 50% B/50% C—30% B/70% C in 1 min, 30% B/70% C isocratic for 15 min. The flow rate was set to 100 mL/min. The collected fractions were analysed by TLC resulting in four combined fractions PS-DE_A-D.

About 350 mg of PS-DE_A were fractionated by HPCCC+ in normal phase mode (flow rate of 6 mL/min, 1,600 rpm) using HEMWat system 23 (300 mL of LL, 1.8 L of UL) to yield three final fractions (PS-DE_A1-3). PS-DE_A2 was identified as compound 5 (235.89 mg; (−)-homopterocarpin (Bezuidenhoudt et al., 1987).

Approximately 250 mg of PS-DE_B were separated by HPCCC+ in normal phase mode (flow rate of 6 mL/min, 1,600 rpm) using HEMWat systems 23 (300 mL of LL, 400 mL of UL) and 19 (700 mL of UL) to yield eight fractions (PS-DE_B1-8). PS-DE_B3 and PS-DE_B5 were identified as compound 6 (0.8 mg (−)-pterocarpin) and the new compound 9 (2.07 mg; pterosantalin D (B-son Bredenberg and Shoolery, 1961). PS-DE_B4 was further subjected to SEC to yield two final fractions (PS-DE_B4a-b). PS-DE_B4b was identified as compound 16 (3.87 mg, 3,9-di-O-methylcoumestrol) (Harper et al., 1969).

About 1 g of PS-DE_C was fractionated by four runs of FC in reversed phase mode using a 15 µm C18 HQ 35 G column as stationary phase and water (A) and acetonitrile (B) as mobile phase. The following gradient was applied: 5% B isocratic for 10 min, 5%–30% B in 10 min, 30% B isocratic for 30 min, 30%–50% B in 30 min, 50% B isocratic for 15 min, 50%–98% B in 10 min, 98% B isocratic for 30 min. The flow rate was set to 15 mL/min. The collected fractions were analysed by TLC yielding in 10 fractions (PS-DE_C1-10). Approximately 250 mg of PS-DE_C5 were further separated by HPCCC+ in normal phase mode (flow rate of 6 mL/min, 1,600 rpm) using HEMWat systems 19 (300 mL of LL, 300 of UL), 18 (600 mL of UL), 17 (600 mL of UL) and 13 (120 mL of UL) resulting in six fractions (PS-DE_C5a-f). PS-DE_C5b was subjected to SEC and led to the isolation of compounds 10 (0.71 mg (−)-pterosantalin E) and 17 (0.52 mg, pterosantalin F) (Maranduba et al., 1979; Naik et al., 2017). Identity of all isolated compounds was assessed by 1D and 2D NMR, HR-MS and optical rotation.

Compounds 14, 1115 and 22 have been previously isolated from the heartwood of P. santalinus (Zwirchmayr et al., 2023). Compound 23 was obtained from Specs.

Pterosantalin D (9): brownish, amorphous powder; 1H (CDCl3, 500 MHz) and 13C NMR (CDCl3, 125 MHz), see Table 2 and Supplementary Figures S10–S13; HRESIMS m/z 303.1592 [M + H]+ (calcd for C18H23O4+, 303.1518).

3.6 UHPLC-HRMS2 data generation

LC-ESI-HRMS2 analyses were performed on an Agilent 1290 Infinity II UHPLC coupled to a hybrid quadrupole time of flight (QTOF) mass spectrometer Agilent 6546 (Agilent Technologies) equipped with an ESI source, operating in positive ion mode. A BEH Waters Acquity C18 UPLC column (2.1 × 150 mm; 1.7 µm) was used, while water + 0.1% formic acid (A) and acetonitrile + 0.1% (B) were used as mobile phase. The flow rate was set to 0.5 mL/min. A linear gradient from 5% to 100% B in 12 min followed by 100% B for 4 min was used. Source parameters were set as follows: capillary temperature at 320°C, source voltage at 3,500 V, sheath gas flow rate at 11 L/min. MS scans were operated in full-scan mode from m/z 100 to 1,200 (0.1 s scan time) with a mass resolution of 67.000 at m/z = 922. A MS1 scan was followed by MS2 scans of the four most intense ions above an absolute threshold of 3,000 counts. Selected parent ions were fragmented at a collision energy fixed at 45 eV and an isolation window of 1.3 amu. The purine [M + H]+ ion (C5H5N4+, m/z = 121.0509) and the hexakis (1H,1H, 3H-tetrafluoropropoxy)phosphazene [M + H]+ ion (C18H19F24N3O6P3+, m/z = 922.0098) were used as internal lock masses. A permanent MS/MS exclusion list criterion was set to prevent oversampling of the internal calibrant. LC-UV and MS data acquisition and processing were performed using MassHunter Workstation software (Agilent Technologies).

3.7 MZmine 3 processing for feature-based molecular networking

The MS2 data files were converted from the .d format to the .mzZML format using MSConvert from ProteoWizard (Chambers et al., 2012). The .mzXML files were then processed using MZmine 3 (version 3.1.0-beta) (Schmid et al., 2023). For mass detection, a noise level of 1 × 104 and 1 × 101 was specified for MS1 and MS2, respectively. The ADAP chromatogram builder was used with the following settings: minimum group size of scans = 4, group intensity threshold = 1 × 104, minimum highest intensity = 1 × 104, m/z tolerance = 10 ppm (Schmid et al., 2023). Next, the ADAP chromatogram deconvolution module was employed with the following settings: S/N threshold = 23, minimum feature height = 1 × 104, coefficient/area threshold = 10, peak duration length = 0.8 min, RT wavelet range = 0.0—0.08 (Myers et al., 2017). Isotopes were grouped by the 13C isotope filter, setting an m/z tolerance of 10 ppm, maximum charge of 1 and RT tolerance of 0.3 min. Peak alignment was performed by the join aligner module, selecting a m/z tolerance of 10 ppm, weight for m/z and RT of 50 and RT tolerance of 0.12 min. The peak list was gap-filled using the same RT and m/z range gap filler with a m/z tolerance of 10 ppm. The feature list row filter was applied which only kept peaks with m/z from 100.00 to 1,200.00 and RT from 0.04 to 16.00 min. The data from the resulting peak list was exported in .mgf and .csv format.

3.8 Generation of the metadata table

A metadata table was created using Notepad++ to incorporate the bioactivity data into the MN. Each .mzXML file was given the corresponding fraction name (F1—F11). The bioactivity was classified into three groups: IC50 < 10 μg/mL in spike assay = active, IC50 from 10 to 50 μg/mL in spike assay = moderately active, IC50 > 50 μg/mL in spike assay or % of caspase inhibition <50% at 50 μg/mL = inactive. Fractions F2 and F3 were attributed as active, F4—F9 were determined as moderately active and F1, F10 and F11 were classified as inactive. The metadata table was saved as .txt file.

3.9 Molecular network generation parameters

A feature-based MN was generated on the GNPS website with the following settings: parent and fragment ion mass tolerance = 0.02 Da, cosine score = 0.6, network topK = 10, minimum matched fragment ions = 4, maximum connected component size = 100, maximum shift between precursors = 500 Da (Wang et al., 2016). The network was searched against the GNPS spectra library. For a spectra library annotation, six peaks must match between a network and a library spectrum with a cosine score above 0.7. For quantification, row sum normalization and the summarization of all files within a group were selected. The obtained MN can be accessed by the following address: https://gnps.ucsd.edu/ProteoSAFe/status.jsp?task=104b80cc65e44c8c86e6a86b255b2045. The MN was visualized with Cytoscape 3.8.2 (Shannon et al., 2003).

3.10 Cells and virus

Caco-2-F03 cells from the Resistant Cancer Cell Line Colletion were grown at 37°C in minimal essential medium (MEM) supplemented with 10% fetal bovine serum (FBS), 100 IU/mL penicillin and 100 μg/mL streptomycin. All culture reagents were purchased from Sigma-Aldrich. The SARS-CoV-2 variant (SARS-CoV-2/FFM7, MT358643) used in the experiment was isolated from Caco-F03 cells. SARS-CoV-2 stocks were cultivated for a maximum of three passages and stored at −80°C.

3.11 Pre-screening assay

A phenotypic screening platform using the measurement of caspase 3/7 activation to monitor virus infection was utilized as previously described (Bojkova et al., 2023). Briefly, Caco-2-F03 cells were seeded into 96-well plates (50,000 cells/well) and incubated at 37°C for 4 days. After the cells reached confluence, the supernatant was replaced by 25 μL/well of fresh medium, 25 μL/well of medium containing tested samples (50 μg/mL) in singlets, and 50 μL/well SARS-CoV-2 suspension (MOI 0.01). Remdesivir (10 μM) was used as positive control. After infection of Caco-2-F03 cells and incubation for 48 h, the caspase 3/7 inhibition was measured using the Caspase-Glo assay kit (Promega), according to the manufacturer’s instructions. Briefly, 100 µL of Caspase-Glo reagent were added to each well, mixed, and incubated at room temperature for 30 min. Luminescence intensity was measured using an Infinite M200 microplate reader (Tecan). The raw values were subtracted by the cell control (no drug, no virus). The final results are expressed as percentage of inhibition relative to the virus control which received no treatment.

3.12 Spike immunostaining assay

Caco-2-F03 cells were seeded into 96-well plates (50,000 cells/well) and incubated at 37°C for 4 days. After the cells reached confluence, the supernatant was replaced by 25 μL/well of fresh medium, 25 μL/well of medium containing serially diluted fractions or pure compounds in triplicates, and 50 μL/well SARS-CoV-2 suspension (MOI 0.01). Remdesivir was used as positive control. 24 h post infection, cells were fixed with acetone:methanol (40:60) solution and immunostaining was performed using a monoclonal antibody directed against the spike protein of SARS-CoV-2 (1:1,500, Sinobiological), which was detected with a peroxidase-conjugated anti-rabbit secondary antibody (1:1,000, Dianova), followed by addition of AEC substrate. The S positive area was scanned and quantified by the Bioreader 7000-F-Z-I microplate reader (Biosys). The results are expressed as percentage of inhibition relative to virus control which received no drug.

3.13 Cytotoxicity

The cytotoxicity of tested fractions and pure compounds was determined in parallel to the spike immunostaining assay. The supernatant of confluent Caco-2-F03 was replaced by 75 μL/well of fresh medium and 25 μL/well of medium containing serially diluted fractions and pure compounds in triplicates and incubated for 24 h. Cell viability was measured by 3-(4,5-dimethylthiazol-2-yl)-2,5-diphenyltetrazolium bromide (MTT) dye reduction assay. 25 μL of MTT solution (2 mg/mL in PBS) were added per well, and the plates were incubated at 37°C for 4 h. Afterwards, the cells were lysed using 100 μL of a buffer containing 20% SDS and 50% N,N-dimethylformamide with the pH adjusted to 4.7 at 37°C for 4 h. Absorbance was determined at 560 nm (reference wavelength 620 nm) using a Tecan infinite M200 microplate reader (TECAN).

Data availability statement

The original contributions presented in the study are included in the article/Supplementary Material, further inquiries can be directed to the corresponding author.

Author contributions

UG, MB, and JR designed and supervised the study. AW wrote the manuscript including the generation of all figures and tables. DB investigated the antiviral effects and interpreted the results together with JC and HR. MB performed the UPLC-MS/MS analysis and helped with setting up the molecular networking analysis. AW performed the processing of MS data including the generation and interpretation of the molecular network. JZ, BK, and AW performed the extraction, fractionation, and structure elucidation of sandalwood constituents. UG and JR revised the manuscript. All authors contributed to the article and approved the submitted version.

Funding

This work was funded by the Austrian Science Fund (FWF, P34028 and P35115).

Acknowledgments

The authors thank M. Loidolt, R. Blažević and P. Fischer-Windsteig (Division of Pharmacognosy, Department of Pharmaceutical Sciences, University of Vienna, Austria) for technical support, and J. Orts and M. J. Abi Saad (Division of Pharmaceutical Chemistry, Department of Pharmaceutical Sciences, University of Vienna, Austria) for the production and purification of SARS-CoV-2 Mpro.

Conflict of interest

The authors declare that the research was conducted in the absence of any commercial or financial relationships that could be construed as a potential conflict of interest.

Publisher’s note

All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article, or claim that may be made by its manufacturer, is not guaranteed or endorsed by the publisher.

Supplementary material

The Supplementary Material for this article can be found online at: https://www.frontiersin.org/articles/10.3389/fmolb.2023.1202394/full#supplementary-material

References

Aron, A. T., Gentry, E. C., McPhail, K. L., Nothias, L. F., Nothias-Esposito, M., Bouslimani, A., et al. (2020). Reproducible molecular networking of untargeted mass spectrometry data using GNPS. Nat. Protoc. 15 (6), 1954–1991. doi:10.1038/s41596-020-0317-5

PubMed Abstract | CrossRef Full Text | Google Scholar

B-son Bredenberg, J., and Shoolery, J. N. (1961). A revised structure for pterocarpin. Tetrahedron Lett. 2 (8), 285–288. doi:10.1016/S0040-4039(01)99242-1

CrossRef Full Text | Google Scholar

Beniddir, M. A., Kang, K. B., Genta-Jouve, G., Huber, F., Rogers, S., and van der Hooft, J. J. J. (2021). Advances in decomposing complex metabolite mixtures using substructure- and network-based computational metabolomics approaches. Nat. Prod. Rep. 38 (11), 1967–1993. doi:10.1039/d1np00023c

PubMed Abstract | CrossRef Full Text | Google Scholar

Bezuidenhoudt, B. C. B., Brandt, E. V., and Ferreira, D. (1987). Flavonoid analogues from Pterocarpus species. Phytochemistry 26 (2), 531–535. doi:10.1016/S0031-9422(00)81448-X

CrossRef Full Text | Google Scholar

Bojkova, D., Reus, P., Panosch, L., Bechtel, M., Rothenburger, T., Kandler, J. D., et al. (2023). Identification of novel antiviral drug candidates using an optimized SARS-CoV-2 phenotypic screening platform. iScience 26 (2), 105944. doi:10.1016/j.isci.2023.105944

PubMed Abstract | CrossRef Full Text | Google Scholar

Chambers, M. C., Maclean, B., Burke, R., Amodei, D., Ruderman, D. L., Neumann, S., et al. (2012). A cross-platform toolkit for mass spectrometry and proteomics. Nat. Biotechnol. 30 (10), 918–920. doi:10.1038/nbt.2377

PubMed Abstract | CrossRef Full Text | Google Scholar

Cho, J. Y., Park, J., Kim, P. S., Yoo, E. S., Baik, K. U., and Park, M. H. (2001). Savinin, a lignan from Pterocarpus santalinus inhibits tumor necrosis factor-alpha production and T cell proliferation. Biol. Pharm. Bull. 24 (2), 167–171. doi:10.1248/bpb.24.167

PubMed Abstract | CrossRef Full Text | Google Scholar

Dahat, Y., Saha, P., Mathew, J. T., Chaudhary, S. K., Srivastava, A. K., and Kumar, D. (2021). Traditional uses, phytochemistry and pharmacological attributes of Pterocarpus santalinus and future directions: A review. J. Ethnopharmacol. 276, 114127. doi:10.1016/j.jep.2021.114127

PubMed Abstract | CrossRef Full Text | Google Scholar

Fox Ramos, A. E., Evanno, L., Poupon, E., Champy, P., and Beniddir, M. A. (2019). Natural products targeting strategies involving molecular networking: Different manners, one goal. Nat. Prod. Rep. 36 (7), 960–980. doi:10.1039/c9np00006b

PubMed Abstract | CrossRef Full Text | Google Scholar

Grienke, U., Mair, C. E., Kirchmair, J., Schmidtke, M., and Rollinger, J. M. (2018). Discovery of bioactive natural products for the treatment of acute respiratory infections – An integrated approach. Planta Med. 84 (09/10), 684–695. doi:10.1055/a-0590-5153

PubMed Abstract | CrossRef Full Text | Google Scholar

Ham, B., Kim, M., Son, Y. J., Chang, S., JungNho, S. H. C. W., Kwon, M. J., et al. (2019). Inhibitory effects of Pterocarpus santalinus extract against IgE/antigen-sensitized mast cells and atopic dermatitis-like skin lesions. Planta Med. 85 (7), 599–607. doi:10.1055/a-0824-1282

PubMed Abstract | CrossRef Full Text | Google Scholar

Harper, S. H., Kemp, A. D., Underwood, W. G. E., and Campbell, R. V. M. (1969). Pterocarpanoid constituents of the heartwoods of Pericopsis angolensis and Swartzia madagascariensis. J. Chem. Soc. C 8, 1109–1116. doi:10.1039/J39690001109

CrossRef Full Text | Google Scholar

Ma, X. Q., Zheng, C. J., Zhang, Y., Hu, C. L., Lin, B., Fu, X. Y., et al. (2013). Antiosteoporotic flavonoids from Podocarpium podocarpum. Phytochem. Lett. 6 (1), 118–122. doi:10.1016/j.phytol.2012.12.004

CrossRef Full Text | Google Scholar

Maranduba, A., De Oliveira, A. B., De Oliveira, G. G., de Reis, J. E., and Gottlieb, O. R. (1979). Isoflavonoids from Myroxylon peruiferum. Phytochemistry 18 (5), 815–817. doi:10.1016/0031-9422(79)80020-5

CrossRef Full Text | Google Scholar

Myers, O. D., Sumner, S. J., Li, S., Barnes, S., and Du, X. (2017). One step forward for reducing false positive and false negative compound identifications from mass spectrometry metabolomics data: New algorithms for constructing extracted ion chromatograms and detecting chromatographic peaks. Anal. Chem. 89 (17), 8696–8703. doi:10.1021/acs.analchem.7b00947

PubMed Abstract | CrossRef Full Text | Google Scholar

Naik, M. M., Kamat, V. P., and Tilve, S. G. (2017). Copper-mediated synthesis of coumestans via C(sp2)-H functionalization: Protective group free route to coumestrol and 4′-O-methylcoumestrol. Tetrahedron 73 (37), 5528–5536. doi:10.1016/j.tet.2017.07.057

CrossRef Full Text | Google Scholar

Natalia, P., Zwirchmayr, J., Rudžionytė, I., Pulsinger, A., Breuss, J. M., Uhrin, P., et al. (2021). Pterocarpus santalinus selectively inhibits a subset of pro-inflammatory genes in interleukin-1 stimulated endothelial cells. Front. Pharmacol. 12, 802153. doi:10.3389/fphar.2021.802153

PubMed Abstract | CrossRef Full Text | Google Scholar

Nothias, L. F., Nothias-Esposito, M., da Silva, R., Wang, M., Protsyuk, I., Zhang, Z., et al. (2018). Bioactivity-based molecular networking for the discovery of drug leads in natural product bioassay-guided fractionation. J. Nat. Prod. 81 (4), 758–767. doi:10.1021/acs.jnatprod.7b00737

PubMed Abstract | CrossRef Full Text | Google Scholar

Olivon, F., Allard, P. M., Koval, A., Righi, D., Genta-Jouve, G., Neyts, J., et al. (2017). Bioactive natural products prioritization using massive multi-informational molecular networks. ACS Chem. Biol. 12 (10), 2644–2651. doi:10.1021/acschembio.7b00413

PubMed Abstract | CrossRef Full Text | Google Scholar

Pan, E., Harinantenaina, L., Brodie, P. J., Miller, J. S., Callmander, M. W., Rakotonandrasana, S., et al. (2010). Four diphenylpropanes and a cycloheptadibenzofuran from Bussea sakalava from the Madagascar dry forest. J. Nat. Prod. 73 (11), 1792–1795. doi:10.1021/np100411d

PubMed Abstract | CrossRef Full Text | Google Scholar

Picker, P., Vogl, S., McKinnon, R., Mihaly-Bison, J., Binder, M., Atanasov, A. G., et al. (2014). Plant extracts in cell-based anti-inflammatory assays - pitfalls and considerations related to removal of activity masking bulk components. Phytochem. Lett. 10, 41–47. doi:10.1016/j.phytol.2014.04.001

CrossRef Full Text | Google Scholar

Schmid, R., Heuckeroth, S., Korf, A., Smirnov, A., Myers, O., Dyrlund, T. S., et al. (2023). Integrative analysis of multimodal mass spectrometry data in MZmine 3. Nat. Biotechnol. 41, 447–449. doi:10.1038/s41587-023-01690-2

PubMed Abstract | CrossRef Full Text | Google Scholar

Shannon, P., Markiel, A., Ozier, O., Baliga, N. S., Wang, J. T., Ramage, D., et al. (2003). Cytoscape: A software environment for integrated models of biomolecular interaction networks. Genome Res. 13 (11), 2498–2504. doi:10.1101/gr.1239303

PubMed Abstract | CrossRef Full Text | Google Scholar

ter Ellen, B. M., Dinesh Kumar, N., Bouma, E. M., Troost, B., van de Pol Dpi, , van der Ende-Metselaar, H. H., et al. (2021). Resveratrol and pterostilbene inhibit SARS-CoV-2 replication in air–liquid interface cultured human primary bronchial epithelial cells. Viruses 13 (7), 1335. doi:10.3390/v13071335

PubMed Abstract | CrossRef Full Text | Google Scholar

Wang, M., Carver, J. J., Phelan, V. V., Sanchez, L. M., Garg, N., Peng, Y., et al. (2016). Sharing and community curation of mass spectrometry data with global natural products social molecular networking. Nat. Biotechnol. 34 (8), 828–837. doi:10.1038/nbt.3597

PubMed Abstract | CrossRef Full Text | Google Scholar

Wasilewicz, A., Kirchweger, B., Bojkova, D., Abi Saad, M. J., Langeder, J., Bütikofer, M., et al. (2023). Identification of natural products inhibiting SARS-CoV-2 by targeting viral proteases: A combined in silico and in vitro approach. J. Nat. Prod. 86 (2), 264–275. doi:10.1021/acs.jnatprod.2c00843

PubMed Abstract | CrossRef Full Text | Google Scholar

Watrous, J., Roach, P., Alexandrov, T., Heath, B. S., Yang, J. Y., Kersten, R. D., et al. (2012). Mass spectral molecular networking of living microbial colonies. Proc. Natl. Acad. Sci. U. S. A. 109 (26), E1743–E1752. doi:10.1073/pnas.1203689109

PubMed Abstract | CrossRef Full Text | Google Scholar

Wen, C. C., Kuo, Y. H., Jan, J. T., Liang, P. H., Wang, S. Y., Liu, H. G., et al. (2007). Specific plant terpenoids and lignoids possess potent antiviral activities against severe acute respiratory syndrome coronavirus. J. Med. Chem. 50 (17), 4087–4095. doi:10.1021/jm070295s

PubMed Abstract | CrossRef Full Text | Google Scholar

Wichtl, M. (2002). Teedrogen und Phytopharmaka. Stuttgart: Wissenschaftliche Verlagsgesellschaft mbH, 551–552.

Google Scholar

Wu, S. F., Chang, F. R., Wang, S. Y., Hwang, T. L., Lee, C. L., Chen, S. L., et al. (2011). Anti-inflammatory and cytotoxic neoflavonoids and benzofurans from Pterocarpus santalinus. J. Nat. Prod. 74 (5), 989–996. doi:10.1021/np100871g

PubMed Abstract | CrossRef Full Text | Google Scholar

Wu, S. F., Hwang, T. L., Chen, S. L., Wu, C. C., Ohkoshi, E., Lee, K. H., et al. (2011). Bioactive components from the heartwood of Pterocarpus santalinus. Bioorg. Med. Chem. Lett. 21 (18), 5630–5632. doi:10.1016/j.bmcl.2011.06.036

PubMed Abstract | CrossRef Full Text | Google Scholar

Yahara, S., Saijo, R., Nohara, T., Konishi, R., Yamahara, J., Kawasaki, T., et al. (1985). Novel bi-isoflavonoids from Dalbergia odorifera. Chem. Pharm. Bull. 33 (11), 5130–5133. doi:10.1248/cpb.33.5130

CrossRef Full Text | Google Scholar

Zwirchmayr, J., Schachner, D., Grienke, U., Rudžionytė, E., de Martin, R., Dirsch, V. M., et al. (2023). Biochemometry identifies suppressors of pro-inflammatory gene expression in Pterocarpus santalinus heartwood. Phytochemistry 2023, 113709. doi:10.1016/j.phytochem.2023.113709

PubMed Abstract | CrossRef Full Text | Google Scholar

Keywords: molecular networking, Pterocarpus santalinus, sandalwood, analysis, natural product, SARS-CoV-2, high-performance counter-current chromatography, isolation

Citation: Wasilewicz A, Zwirchmayr J, Kirchweger B, Bojkova D, Cinatl J, Rabenau HF, Rollinger JM, Beniddir MA and Grienke U (2023) Discovery of anti-SARS-CoV-2 secondary metabolites from the heartwood of Pterocarpus santalinus using multi-informative molecular networking. Front. Mol. Biosci. 10:1202394. doi: 10.3389/fmolb.2023.1202394

Received: 08 April 2023; Accepted: 24 May 2023;
Published: 06 June 2023.

Edited by:

Evagelos Gikas, National and Kapodistrian University of Athens, Greece

Reviewed by:

Seung Hyun Kim, Yonsei University, Republic of Korea
Roland Hellinger, Medical University of Vienna, Austria

Copyright © 2023 Wasilewicz, Zwirchmayr, Kirchweger, Bojkova, Cinatl, Rabenau, Rollinger, Beniddir and Grienke. This is an open-access article distributed under the terms of the Creative Commons Attribution License (CC BY). The use, distribution or reproduction in other forums is permitted, provided the original author(s) and the copyright owner(s) are credited and that the original publication in this journal is cited, in accordance with accepted academic practice. No use, distribution or reproduction is permitted which does not comply with these terms.

*Correspondence: Ulrike Grienke, ulrike.grienke@univie.ac.at

ORCID: Andreas Wasilewicz, orcid.org/0000-0003-3396-5071; Julia Zwirchmayr, orcid.org/0000-0003-2721-3762; Benjamin Kirchweger, orcid.org/0000-0001-8948-1012; Denisa Bojkova, orcid.org/0000-0001-7914-9710; Jindrich Cinatl Jr., orcid.org/0000-0002-6744-2087; Holger F. Rabenau, orcid.org/0000-0001-9394-2176; Judith M. Rollinger, orcid.org/0000-0001-6581-0774; Mehdi A. Beniddir, //orcid.org/0000-0003-2153-4290; Ulrike Grienke, orcid.org/0000-0003-0305-9270

Disclaimer: All claims expressed in this article are solely those of the authors and do not necessarily represent those of their affiliated organizations, or those of the publisher, the editors and the reviewers. Any product that may be evaluated in this article or claim that may be made by its manufacturer is not guaranteed or endorsed by the publisher.