Acetylcholinesterase inhibitors and antioxidants mining from marine fungi: bioassays, bioactivity coupled LC–MS/MS analyses and molecular networking

Marine fungi are potentially important resources for bioactive lead compounds for discovering new drugs for diseases such as Alzheimer’s disease. In this paper, the combined bioassay model of acetylcholinesterase (AChE) inhibition, 1,1-diphenyl-2-picrylhydrazyl (DPPH) free radical scavenging, and Artemia larval lethality was used to evaluate the activity and toxicity of 35 marine fungal strains from seas around China. Their bioactive constituents were revealed by thin layer chromatography (TLC) autography, bioactivity coupled LC–MS/MS and Global Natural Products Social Molecular Networking (GNPS). The results show that the extracts of five strains exhibited higher AChE inhibition ratios than the positive control compound, ‘tacrine’, for which the ratio was 89.8% at 200 μg/ml. Six strains displayed both AChE inhibition (inhibition ratios > 20% at 200 μg/ml) and DPPH scavenging activity (scavenging ratios > 30% at 200 μg/ml) together with low Artemia larval toxicity (lethal rates < 12%). TLC autography showed that the fractioned extracts of four strains contained highly diverse and different bioactive constituents, including strains Talaromyces sp. C21-1, Aspergillus terreus C23-3, Trichoderma harzianum DLEN2008005, and Penicillium corylophilum TBG1-17. From the most potent sample F-11-1-b (derived from Aspergillus terreus C23-3), five AChE inhibitors and seven antioxidants were recognized as bioactive molecules by AChE coupled ultrafiltration followed by LC–MS/MS, and LC–MS/MS coupled with DPPH incubation. Furthermore, with the aid of GNPS, the AChE inhibitors were plausibly annotated as territrem analogues including territrems A–C/D, arisugacin A and an unknown compound 4, and the seven antioxidants were assigned as butyrolactone Ι, aspernolide E, a phenolic derivative and possibly unknown compounds 8–10 and 12.


Introduction
Alzheimer's disease (AD) is the major form of senile dementia, which seriously threatens the physical and mental health of the elderly. Its main pathological manifestations include the formation of Aβ plaques and neurofibrillary tangles, the amyloidosis of capillaries and the abnormal damage of synapses (Goedert 1993;Haass and Selkoe 1993;Scheff et al. 1990;Verbeek et al. 1994;Wilcock and Esiri 1982). Although its exact pathological mechanism still remains undefined, the deficiency of the neurotransmitter acetylcholine (ACh) in the brain and the neuroinflammation and neuronal apoptosis induced by oxidative stress are well accepted among the main factors of AD (Crawley 1993;Kok et al. 2006;Rogers et al. 1996). Thus, inhibitors of acetylcholinesterase (AChE), the enzyme responsible for the hydrolysis of ACh, as well as antioxidant agents including those from nature, are highly valued in the treatment of AD (Butterfield et al. 2008;Mcgleenon et al. 1999;Mehta et al. 2012;Ye and Wang 2012).
Due to their high biological and chemical diversity, marine fungi are an important resource for many bioactive compounds including potent AChE inhibitors and antioxidants (Ji and Wang 2016;Jin et al. 2016;Liu et al. Edited by Chengchao Chen.

Electronic supplementary material
The online version of this article (https ://doi.org/10.1007/s4299 5-020-00065 -9) contains supplementary material, which is available to authorized users. 2019; Nong et al. 2014). However, isolating bioactive compounds from complicated microbial fermentation extracts is a time-consuming and challenging process via traditional methodology, and many known compounds are repeatedly isolated and characterized leading to a waste of resources.
To improve natural products research efficiency, there is an urgent need to quickly localize bioactive molecules in crude extracts and perform structural dereplications before the time-consuming process of traditional isolation and structure elucidation. Bioactivity-coupled liquid chromatography-(tandem) mass spectrometry (Bio-LC-MS/MS) technologies are based on the interaction between small molecules (ligands) and target molecules like receptors, enzymes or free radicals, in a high-throughput manner (Boer et al. 2007;Chen et al. 2010). For macromolecular targets like receptors and enzymes, ultrafiltration liquid chromatography-mass spectrometry (UF-LC-MS) is a suitable analytical method and utilizes the principle that the ligand-target complex can be captured on an ultrafiltration membrane and then subsequently dissociated by organic solvents (Zhao et al. 2016). For small molecular targets like free radicals and metal ions, simple co-incubation followed by LC-MS analysis is able to recognize interacting ligands via a reduction in peak height (Chen et al. 2010).
However, it is difficult to adequately describe new chemical entities by Bio-LC-MS/MS, especially when only low-resolution MS data is obtained. In recent years, the emergence of molecular networking technologies such as Global Natural Products Social Molecular Networking (GNPS; website: https ://gnps.ucsd.edu/) has provided a new perspective for the analysis of natural products at early stages of investigation (Allard et al. 2016;Wang et al. 2016). GNPS allows visualization of all of the molecular ions detected in an LC-MS/MS experiment and reveals the chemical relationships among them based on MS/MS similarity. With the accumulated MS data of natural products from global contributors, GNPS has become an important tool for the rapid and large-scale annotation of known compounds, as well as the discovery of novel compounds . Therefore, the combined use of Bio-LC-MS/MS in companion with GNPS provides a highly efficient pipeline for recognizing bioactive compounds in crude extracts and can provide the annotation of known compounds as well as the indication of structural relatives for novel ones.
In this study, the combined bioassays comprised of AChE inhibition-DPPH free radical scavenging-Artemia larval lethality was used to evaluate the anti-AD related activity and toxicity of 35 marine fungi. The most promising samples were probed for their bioactive molecules via integrated Bio-LC-MS/MS and GNPS analyses. To the best of our knowledge, this is the first report on the application of Bio-LC-MS/MS coupled with GNPS to screen for lead compounds with anti-AD potential from marine fungi.

Screening for AChE inhibition and DPPH free radical scavenging activities
Potato sucrose peptone medium (PSP) and malt extract medium (M) were used as a culture medium for 35 marine fungi and 280 fractioned samples were prepared from the 70 fermentation crude extracts by column chromatography on silica gel eluted sequentially with n-hexane, dichloromethane, ethyl acetate, and methanol. All samples were applied onto grids on TLC plates and sprayed with AChE enzyme reaction-coloring solution or DPPH solution to preliminarily demonstrate their activity by the appearance of clear zones on a colored background. The results (see details in Supplementary Table S1, Fig. S1) showed that 60 samples had AChE inhibitory activity (from 24 marine fungi) and 47 samples had DPPH radical scavenging activity (from 21 marine fungi). The samples with AChE inhibitory activity were mainly from the fractions eluted by n-hexane, dichloromethane and ethyl acetate, while the samples with DPPH scavenging activity were mainly from the eluents using dichloromethane, ethyl acetate and methanol.
The above active samples were further quantitatively determined for their potency ( Fig. 1; Table 1). At the final concentration of 200 μg/ml, seven samples (from five strains) exhibited higher AChE inhibition activity than the positive control (Tacrine: inhibition rate of 89.8% at 200 μg/ ml), and 12 samples (from nine strains) displayed DPPH radical scavenging rate > 50% (positive control vitamin C (Vc): 91.6% at 200 μg/ml). Moreover, eight extracts of six strains had both AChE inhibitory and antioxidant activities.

Toxicity evaluation of marine fungal extracts
To identify the bioactive samples with low toxicity, 26 potent samples were evaluated using the Artemia larval lethality assay. The results (Table 1) indicated that the lethal rates of strain C23-3, C3-18 and C23-18 were < 15% (the samples with AChE inhibitory activity both antioxidant activity). Among the 14 samples with AChE inhibition rates > 50%, nine samples had lethal rates < 10% (with AChE inhibition activity only). The seven samples with DPPH radical scavenging rates > 50% displayed lethal rates < 11% (only with antioxidant activity).

TLC autograph for AChE inhibitory and DPPH radical scavenging activities
To preliminarily compare the bioactive constituents in different samples, TLC autograph analysis was performed for four samples with significant activity and low toxicity. This yielded the identification of four samples with a richness of different components, including F-7-1-c (from Talaromyces sp. C21-1), F-11-1-b (from Aspergillus terreus C23-3), F-16-1-b (from Trichoderma harzianum DLEN2008005)  Table 4 in "Materials and Methods". Medium 1 and 2, respectively, represent potato medium and malt extract medium. a, b, c and d in the legends refer to the fractionated samples eluted by different solvents: (a) hexane, (b) dichloromethane, (c) ethyl acetate, and (d) methanol, respectively. This numbering scheme also applies to the other parts of this paper Table 1 Toxicity of marine fungal extracts to Artemia larvae a The sample numbering principle in this table: Fungus (F)-strain sequence number-cultural medium. The strain sequence numbers correspond to Table 4 in "Materials and Methods". Medium 1 and 2, respectively, represent potato medium and malt extract medium And the a, b, c and d in the sample numbers refer to the fractional samples eluted by different solvents hexane, dichloromethane, ethyl acetate, and methanol, respectively. This numbering rule also applies to the other sections of this paper and F-24-1-b (from Penicillium corylophilum TBG1-17). Among these, the AChE inhibitory activity and DPPH free radical scavenging activity of F-11-1-b (marked with red frame in Supplementary Fig. S2) were both significant. Thus, its active components were further analyzed by LC-MS/MS coupled with bioactivity.

Bioactivity coupled LC-MS/MS analyses and molecular networking
To localize the active molecules in sample F-11-1-b, it was incubated with AChE in an ultrafiltration-LC-MS/MS experiment and with DPPH in antioxidant coupled LC-MS/ MS experiment. The bioactive molecules were recognized by comparison of corresponding peak area in the control experiment (sample non-specifically captured by ultrafiltration membrane) vs co-incubation experiment (sample captured by AChE and by ultrafiltration membrane) or the peak area in sample control vs in co-incubation experiment (sample totally or partially consumed by DPPH). However, to avoid inaccuracies caused by directly using the peak area from the total ion chromatograms (TIC) or base peak chromatograms (BPC) due to overlapping peaks, we compared the peak areas of the extracted ion chromatograms (EIC) for the parent ions of the main peaks from the BPC.
As shown in Fig. 2 and Table 2, five parent ions at m/z 513.00-513.30 (1), m/z 497. 15-497.25 (2), m/z 527.00-527.30 (3), m/z 553.00-553.40 (4), and m/z 511.00-511.40 (5) demonstrated higher integral areas in the experimental group (sample & AChE-Ultrafiltration) than in the control group (sample-Ultrafiltration), indicating that they are specifically captured by the enzyme. Especially for compounds 4 and 5, their selective binding rates reached 27.5% and 73.7%, which suggested that they were possibly the main AChE inhibitors in this sample. The peak area of another parent ion at m/z 549.00-549.30 (the 'control peak') showed little change between the two experiments, indicating the changes of other peaks were not false positives.
To obtain more information on the structures and relationship of these bioactive molecules, the LC-MS/MS data of sample F-11-1-b were submitted to the GNPS platform for molecular networking. In the network (Fig. 4), compounds 1-5 were in the same cluster suggesting that they were structurally related. Compound 3 gave MS 2 fragment ions at m/z 509. 00, 491.01, 473.20, 373.00, 291.10, 195.07, which matched territrem B in the GNPS database with a cosine value of 0.94. Considering Aspergillus terreus was frequently reported as a producer of this AChE inhibitor (Nong et al. 2014;Peng 1995), compound 3 was annotated as territrem B. Compared with compound 3, both the parent ion and a series of fragment ions deriving from compound 1 showed a mass difference of − 14 Da (-CH 2 ). Similarly, compound 2 decreased by 30 Da (-CH 2 O) compared to compound 3, and those of compound 5 decreased by 16 Da (-CH 4 or -O). According to their MS characteristics and the MS data recorded in literature, the known territrem analogs for compounds 1, 2, and 5 were possibly territrem C (or D), arisugacin A and territrem A, respectively (Jiang et al. 2010;Lee et al. 1992). However, other isomers are also possible based on these mass spectra.
The parent ion of compound 4 showed a mass difference of + 26 Da compared with territrem B (3), and the MS 2 fragments ion were at m/z 535. 06, 509.00, 491.01, 473.20, 291.10, 195.07; thus, it was determined that compound 4 was a dehydrogenated territrem analog with two more   For the antioxidative molecules, compound 7 with a parent mass at m/z 424.97 was annotated as butyrolactone I by GNPS molecular library matching with a cosine value of 0.93, while compound 6 at m/z 423.03 in the same cluster was its dehydrogenated analog, matching the structure of aspernolide E in the Dictionary of Natural Products (https ://dnp.chemn etbas e.com) (He et al. 2013).
The potent antioxidants 8 and 9, present in the same cluster, may be long-chain lipids because 8 (at m/z 607.36) shared a mutual strong fragment at m/z 551.25-551.47 and a moderately similar MS 2 pattern (cosine value 0.74) to a phospholipid (CCMSLIB00000078672 in GNPS library) with a significantly higher parent mass (m/z 692.40).
Summarizing, compounds 4, 8-10 and 12 have a high possibility of being bioactive compounds with new structures, and this will be confirmed by ongoing investigations.
In particular, 20 bioactive fractions from 12 strains including Talaromyces sp. C21-1, Asp. terreus C23-3, Trichoderma harzianum DLEN2008005 and P. corylophilum TBG1-17 with diverse bioactive constituents also displayed low toxicity to Artemia larvae, suggesting they may be valuable in anti-AD drug discovery. The Artemia lethality model has been frequently used for screening cytotoxic and pesticidal agents due to its good correlation, low cost and convenience (Badisa et al. 2009;Cheng et al. 2003;Zhang et al. 2012Zhang et al. , 2013. In this investigation, it was used as a toxicity evaluation approach. As for AChE inhibitory samples, if they did not show significant toxicity to Artemia larvae, then they are more likely to be reversible AChE inhibitors with anti-AD potential rather than irreversible ones, the latter being preferred in pesticide discovery.
In our work, before bioactivity coupled LC-MS/MS analysis, we performed a preliminary assessment via TLC autography for surveying the diversity of bioactive components in different samples. This method is not as highresolution and informative as bioactivity coupled LC-MS/ MS; however, it has the advantage of high efficiency, lower cost, and intuition. For the discovery of bioactive natural products, it is typical to use a chemistry-guided or bioactivity-guided isolation strategy. Nevertheless, traditional chemical separation is a blind process and the conventional offline bioactivity tracing is inefficient, often leading to repeated isolation of known compounds. Hence, it is necessary to localize bioactive molecules and obtain as much structural information as possible to increase the probability of finding new compounds before the significant effort is put into compound isolation. By the strategy of combining bioactivity coupled LC-MS/MS analysis and GNPS molecular networking in the present study, we have accelerated this process. In the current study, using this integrated strategy, we were able to annotate five territrem derivatives as AChE inhibitors while two aromatic butyrolactones and five other metabolites were assigned as antioxidants. Moreover, five of the described molecules have a high possibility of being new compounds or new metabolites from Asp. terreus. This latter conclusion requires confirmation through further investigation that will require chromatography and spectroscopic experiments. Summarizing, this new integrated strategy provided a significant amount of useful information and identified key new products for further study.
The robustness of assessing novelty by this methodology is to a great extent influenced by the size of the compound MS/MS database. Up until now, GNPS has deposited the MS/MS spectra of more than 150,000 natural products (as of 2020-04-13) from contributors all over the world and from the third party MS/MS libraries like NIST, Massbank and HMDB. While a large number, this is certainly a limitation compared to the number of all known metabolites. However, the continued increase of GNPS data will improve the effectiveness of our Bio-LC-MS/MS coupled GNPS strategy in the future.
In conclusion, by a combined bioactivity-toxicity model, 20 fractionated extracts from diverse and previously less studied marine fungal taxa were assessed for their AChE inhibitory and/or antioxidant potential together with low Artemia toxicity. TLC autographies and Bio-LC-MS/MS coupled GNPS localized bioactive molecules and provided rich structural information including indications of likely new bioactive compounds for further study. This strategy could be useful for discovering new anti-Alzheimer lead compounds from the sea.

Materials
The 35 strains of marine fungi used in this study were from biological samples collected in the intertidal zone of Dalian China (Zhang et al. 2009), gill tissues from sharks captured in the East China Sea , and different corals in the Zhanjiang sea area of China, as shown in Table 4.

Preparation of samples of secondary metabolites of marine fungi
The 3-5 day agar plates of the strains were used to inoculate Erlenmeyer flasks each containing 200 ml potato sucrose peptone liquid medium (PSP, containing boiled juice from 200 g potato per liter, 20 g/L sucrose, 5 g/L peptone, and 20 g/L sea salt) or malt extract liquid medium (M, containing 15 g/L malt extract, 20 g/L sea salt), respectively. Both media were used for each strain. The static fermentation lasted for 20 days at room temperature (about 28 °C). Afterwards, the fermentation broth was extracted using ethyl acetate and the mycelia were extracted using methanol. The two extracts for each sample were combined and concentrated to dryness using a rotary evaporator at 45 °C. Then, each crude extract was mixed with silica gel, applied to a short silica gel flash column and stepwise eluted with n-hexane, dichloromethane, ethyl acetate, and methanol to yield four fractional samples.

The screening of AChE inhibitory activity and antioxidant activity
Samples were first preliminarily screened for the two activities on TLC plates by array autography. In this experiment, the samples (20 μg for each) were applied onto grids (1 cm × 1 cm) on the GF 254 TLC plates with a capillary tube. For autography of AChE inhibition, the plate was first evenly sprayed with 0.5 U/ml AChE solution, air dried to fix the enzyme, and then the plates were kept in an incubator with constant moisture at 37 °C for 20 min. Afterwards, DTNB (5 mmol/L) and ATCh (10 mmol/L) solutions were mixed (1:1, v/v) and sprayed onto the TLC plates. These were incubated for another 10 min and the results were recorded. The active spots were white on a yellow background. For autography of DPPH scavenging activity, DPPH methanol solution (1.28 mmol/L) was sprayed evenly on the TLC plates containing the samples, and the results were observed after keeping the plates in darkness for 5-8 min. The active spots had no color and contrasted to a purple background.
The positive samples from the above screening program were then quantitatively evaluated for inhibition to AChE or DPPH scavenging rates at a dose of 200 μg/ml using 96-well microplates by a previously reported method (Yang et al. 2018). Each sample was tested in duplicate and the averages were calculated. Tacrine and vitamin C were taken as positive controls for the AChE experiment and DPPH experiments, respectively.

Toxicity determination of active samples
Artemia toxicity was determined for samples with an AChE inhibition rate and/or a DPPH radical scavenging rate ≥ 50%. The final concentration of samples was 50 μg/ml and the larval number per assay was 20-30/well in 48-well microplates. CuSO 4 (50 μg/ml) was used as positive control. Each sample was tested in triplicate and the average was calculated. The experiment and lethality rate calculations were performed following a previously reported method (Kiviranta et al. 2010).

TLC-activity autograph analyses
The samples with high activity and low toxicity were compared for their diversity of active components by TLC autography. The samples were first applied onto 1 3 TLC plates, developed with an organic solvent system (chloroform:methanol = 25:1 (v/v) for n-hexane eluted fractions; chloroform:methanol = 5:1 (v/v) for other fractions), dried, observed and photographed under wavelengths of 254 and 365 nm. Then, the bioactive spots were demonstrated using the same colorimetric methods described above in the screening by autography section.

Bioactivity coupled HPLC-MS/MS analyses
Pretreatment of AChE coupled ultrafiltration-LC-MS/MS analysis: in the treatment group, the sample was first incubated with AChE (final concentration: 0.5 mg/ml for sample and 0.2 U/ml for enzyme) in 100 μl phosphate buffer solution (PBS) containing 15% MeOH in a water bath at 37 °C for 30 min. Then, the mixture was transferred into an ultrafiltration (UF) centrifugal tube (UFC510096 Millipore 100 KD) and centrifuged at 12,000 r/min for 20 min to capture the AChE-ligand complex. The UF membrane loaded with this complex was rinsed and centrifuged with PBS (100 μl, three times). Afterwards, 50% aqueous MeOH (100 μl, three times) was added to the membrane to release the small molecule ligands from the complex, followed by centrifugation. The filtrate was concentrated to dryness and finally resolved in 100 μl 50% aqueous MeOH for LC-MS/MS analysis. In the blank group, the sample was directly dissolved in 100 μl 50% aqueous MeOH. In the control group, the sample was incubated in PBS-MeOH solution without AChE, and the other procedures were the same with the treatment group. The injection volumes of the three groups were all 25 μl for the next step of LC-MS/MS analysis. Pretreatment of DPPH radical scavenging activity coupled LC-MS/MS analysis: in the treatment group, the samples (50 μl 1.0 mg/ml dissolved in methanol) reacted with DPPH methanol solution (50 μl, 0.16 mmol/L) for 30 min in darkness. In the control group, the DPPH solution was replaced by the same amount of methanol. In the blank group, 0.08 mmol/L DPPH solution was directly used as a sample for LC-MS/MS. The injection volumes of the three groups were all 25 μl in LC-MS/MS analysis.
All nominal mass resolution LC-MS/MS analyses were run on a LC-PDA-MS 2 system consisting of Thermo Finnigan Surveyor Autosampler-Plus, LC-Pump-Plus, PDA-Plus, and a Thermo Finnigan LCQ Advantage Plus mass spectrometer. The monitoring wavelength and mass scan range were set to be 190-600 nm and m/z 100-2000 in the positive mode (the unified mode for GNPS molecular networking), respectively. The column was a Phenomenex Kinetex C18 100A reverse phase column (100 × 4.60 mm, 5 μm). The flow rate was 0.6 ml/min and a gradient program of CH 3 CN/H 2 O (0.1% formic acid) was used for elution: 0-8 min for 40%, 8-11 min for 40 to 99%, 11-16 min for 99%, 16-16.2 min for 99 to 40%, 16.2-20 min for 40%.
The MS/MS spectra of sample F-11-1-b were used to establish a molecular network using the GNPS online platform and visualized by Cytoscape 3.7 software for the purpose of molecular networking and dereplication. The 'Default Molecular Networking Results views' and 'Molecular-library search functions' were used to search hits for possible known compounds. The raw data have been deposited in the MassIVE database at the GNPS Web site (https ://gnps.ucsd.edu/Prote oSAFe /stati c/gnps-splas h.jsp) and is publicly available through access number MSV000085297.