Next Article in Journal
Vibrational Biospectroscopy: An Alternative Approach to Endometrial Cancer Diagnosis and Screening
Next Article in Special Issue
Comparative Analyses of Medicinal Chemistry and Cheminformatics Filters with Accessible Implementation in Konstanz Information Miner (KNIME)
Previous Article in Journal
Inactivation of E. coli, S. aureus, and Bacteriophages in Biofilms by Humidified Air Plasma
Previous Article in Special Issue
In Silico Drug Repositioning to Target the SARS-CoV-2 Main Protease as Covalent Inhibitors Employing a Combined Structure-Based Virtual Screening Strategy of Pharmacophore Models and Covalent Docking
 
 
Font Type:
Arial Georgia Verdana
Font Size:
Aa Aa Aa
Line Spacing:
Column Width:
Background:
Article

The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors—In Silico and In Vitro Studies

by
Dominik Tarabasz
1,
Paweł Szczeblewski
2,
Tomasz Laskowski
2,
Wojciech Płaziński
3,4,
Ewa Baranowska-Wójcik
5,
Dominik Szwajgier
5,
Wirginia Kukula-Koch
1,*,† and
Henry O. Meissner
6,†
1
Department of Pharmacognosy with Medicinal Plants Garden, Medical University of Lublin, 1 Chodzki Str., 20-093 Lublin, Poland
2
Department of Pharmaceutical Technology and Biochemistry and BioTechMed Centre, Faculty of Chemistry, Gdańsk University of Technology, Gabriela Narutowicza Str. 11/12, 80-233 Gdańsk, Poland
3
Jerzy Haber Institute of Catalysis and Surface Chemistry, Polish Academy of Sciences, Niezapominajek 8, 30-239 Krakow, Poland
4
Department of Biopharmacy, Faculty of Pharmacy, Medical University of Lublin, Chodźki 4A, 20-093 Lublin, Poland
5
Department of Biotechnology, Microbiology and Human Nutrition, University of Life Sciences in Lublin, Skromna 8 Street, 20-704 Lublin, Poland
6
Therapeutic Research, TTD International Pty Ltd., 39 Leopard Ave., Elanora, Gold Coast 4221, Australia
*
Author to whom correspondence should be addressed.
These authors contributed equally to this work.
Int. J. Mol. Sci. 2022, 23(9), 4858; https://doi.org/10.3390/ijms23094858
Submission received: 19 March 2022 / Revised: 22 April 2022 / Accepted: 25 April 2022 / Published: 27 April 2022
(This article belongs to the Special Issue Application of In Silico Techniques in Drug Design)

Abstract

:
The aim of the study was to present the fingerprint of different Lepidium peruvianum tuber extracts showing glucosinolates-containing substances possibly playing an important role in preventinting dementia and other memory disorders. Different phenotypes of Lepidium peruvianum (Brassicaceae) tubers were analysed for their glucosinolate profile using a liquid chromatograph coupled with mass spectrometer (HPLC-ESI-QTOF-MS/MS platform). Qualitative analysis in 50% ethanolic extracts confirmed the presence of ten compounds: aliphatic, indolyl, and aromatic glucosinolates, with glucotropaeolin being the leading one, detected at levels between 0–1.57% depending on phenotype, size, processing, and collection site. The PCA analysis showed important variations in glucosinolate content between the samples and different ratios of the detected compounds. Applied in vitro activity tests confirmed inhibitory properties of extracts and single glucosinolates against acetylcholinesterase (AChE) (15.3–28.9% for the extracts and 55.95–57.60% for individual compounds) and butyrylcholinesterase (BuChE) (71.3–77.2% for the extracts and 36.2–39.9% for individual compounds). The molecular basis for the activity of glucosinolates was explained through molecular docking studies showing that the tested metabolites interacted with tryptophan and histidine residues of the enzymes, most likely blocking their active catalytic side. Based on the obtained results and described mechanism of action, it could be concluded that glucosinolates exhibit inhibitory properties against two cholinesterases present in the synaptic cleft, which indicates that selected phenotypes of L. peruvianum tubers cultivated under well-defined environmental and ecological conditions may present a valuable plant material to be considered for the development of therapeutic products with memory-stimulating properties.

1. Introduction

According to available statistics, average life expectancy and longevity increased from about 45 years in 1900 to over 75 years today. This can be attributed to advances in medical, biochemical, nutritional, and environmental sciences which, in turn, constantly contribute to better medical treatment and socially implemented progress in living conditions and extended life expectancy.
The extension of human life span over the last century has translated into an increase in the incidence of diseases that disrupt the work of the Central Nervous System (CNS), especially in highly developed countries. Neurodegenerative diseases are considered a global epidemic due to their prevalence and constantly increasing incidence rates. According to data reported by the World Health Organization (WHO), in 2050, the number of people with dementia is expected to exceed 115 million worldwide [1]. Therefore, the search for new active compounds with neuroprotective activity is of the highest importance. So far, several natural products have been introduced to the treatment of neurodegenerative disorders. Galantamine, an isoquinoline alkaloid from Amaryllidaceae, is certainly a well-known example of a drug introduced to the therapeutic strategies in the treatment of the symptomatic Alzheimer’s disease (AD) as an acetylcholinesterase (AChE) inhibitor [2]. This group of medicines regulates the level of acetylcholine (ACh) in the synaptic cleft. It was proved that during the aging of the organism, the amount of ACh decreases, and additionally, an increase in the activity of esterases (acetylcholine- and butyrylcholinesterase (BuChE) that decompose ACh is observed that decreases the effectiveness of stimuli transmission. For that reason, the esterase inhibitors are first-line drugs that improve the transmission of the stimulus between neuronal cells [3]. The drugs available on the market, however, have a number of side effects or an insufficient duration of action, as discussed in comprehensive reviews, e.g., by Haake et al. 2020 [4]. Therefore, scientists continuously search for other medicines that could be more efficient in regulating the process described above.
During the last 25 years, one of the plants intensively researched for its influence on a wide spectrum of health and medical conditions contributing to improved gender-specific healthy living status and extended longevity is a tuber (or hypocotyl) of a native plant from the Brassicacae botanical family—Peruvian maca. This plant represents one of some 175 identified Lepidium species—Lepidium peruvianum synonym L. meyenii is known under the common name maca. In this work, Peruvian maca will be referred to as L. peruvianum [5]. When harvested in the high Andes, the maca crops are composed of some 13 differently coloured hypocotyls termed in this work as “phenotypes”, with yellow, red, black, and purple appearing to be dominant [6].
Maca has a long history of cultivation and consumption since the time of the Incas established in its acknowledged native habitat—the Junín plateau in the highlands of the Peruvian Andes [7]. Now it is commercially grown as an important cash crop and therapeutic commodity for domestic use and for export [8]. Blends or individual maca phenotypes, distinguished by hypocotyls’ colour, have induced different gender- and age-specific physiological responses, alleviating the severity of several known health conditions [9,10,11]. When taken orally as a dietary supplement, Peruvian maca is considered a non-hormonal, adaptogenic herb exhibiting hormone-regulating properties for a wide range of age groups of men and women [6,12,13], helping to restore metabolic harmony and induce hormonal balance specific to gender and appropriate age-stage.
Amongst the well-established multi-pharmacological functions of Peruvian maca, the confirmed antidepressant effect of this plant was found interesting to the authors as it implies the possibility for maca metabolites to display the ability to penetrate the blood–brain barrier [14].
During the last 30 years, Peruvian maca hypocotyls were the subject of extensive chemical and nutritional studies [6,7,15,16], with glucosinolates (GLSN) identified as the main active natural constituents and glucotropaeolin (GLT) and methoxyglucotropaeolin (MGT) as the leading components [16,17,18]. Additionally, unsaturated fatty acids and benzylated alkamides (macamides and macaenes) have been detected as unique to maca, helping in chemical identification of this plant [18,19].
Glucosinolates are well-defined biochemical plant components belonging to a group of metabolites characterized by the presence of nitrogen and sulphur atoms in their structure. They are spread across the representatives of the following botanical families: Brassicaceae, Tovariaceae, Capparidaceae, Stegnospermaceae, Moringaceae, Resedaceae, and other unrelated plant families, such as Limnanthaceae, Caricaceae, and Tropaeolaceae. When the structure of plant cells is physically disrupted or damaged, glucosinolates are hydrolyzed by the enzyme called myrosinase, which usually occurs during the physical processing of tubers. As a consequence, a sugar molecule is released, whereas the aglycon degrades, forming oxazolidine-2-thiones, nitriles, isothiocyanates, epithionitriles, thiocyanates, and other derivatives of interest from a pharmacological perspective [20,21]. Recent studies confirm the broad medicinal properties of glucosinolates and isothiocyanates that include anticancer and antimicrobial activity [22] and provide protection to the cardiovascular and central nervous systems, which may be important in the prevention of diabetic nephropathy and neuropathy [23].
Only a few up-to-date publications dealt with the anticholinesterase properties of this group of secondary metabolites. Previously, Burcul and colleagues determined the inhibitory properties of 11 compounds from the isothiocyanate group. Of these, 2-methoxyphenyl isothiocyanate showed the best inhibitory properties against AChE with IC50 of 0.57 mM, while in the BuChE inhibition test, 3-methoxyphenyl isocyanate showed the best result of 49.2% at 1.14 mM [24]. Additionally, Blazevic et al. [25] investigated volatile isolates from Bunias erucago, with glucosinalbin and glucoraphanin as the leading components inhibiting both AChE and BuChE. Two components of extracts from maca tubers, namely N-benzylhexadecanamide and N-acetylbenzylamine, were already reported to exhibit neuroprotective effects, presumably partly due to inhibition of the AChE and BuChE enzymes [26], but the scientific literature still lacks information on the potential of glucosinolates from maca.
This fact encouraged the authors to embark on further study to assess inhibitory properties of glucosinolates present specifically in Peruvian maca (Lepidium peruvianum syn. L. meyenii) and to explain their effect towards the AChE and BuChE enzymes. Another reason for choosing this plant for the study is the fact that Peruvian maca has been demonstrated to contain about 100 times more glucosinolates than other well-known cruciferous plants, such as cabbage or broccoli [14].
These goals are complemented by a desire to compare the content of glucosinolates and to study the fingerprints of the rich collection of five prime Peruvian maca phenotypes labelled as “Yellow”, “Black”, “Red”, “Purple”, and “Gray” that were grown in two locations of the Peruvian Andean Highlands. In this research, the HPLC-ESI-QTOF-MS/MS technique and the chemometric approach were adopted to analyse the obtained results, which were extended for measuring the in vitro inhibitory potential of individual glucosinolate fractions towards AChE and BuChE. Additionally, the explanation of the grounds for the memory-enhancing activity of glucosinolates using an in silico approach linked to modeling studies was undertaken and was specifically adopted in an attempt to explain the dynamics behind utilizing the inhibitory potential of these metabolites towards both the AChE and BuChE enzymes.

2. Results and Discussion

2.1. Qualitative Analysis of Extracts from Lepidium peruvianum by Liquid Chromatography Coupled with Mass Spectrometry (HPLC-ESI-QTOF-MS/MS)

Glucosinolates (Figure 1) are composed of a glucose structure which is combined with an O-sulfated thiohydroximate moiety and are variable in length—side chain is variable in derivatives of amino acids. The aliphatic chain is biosynthesized most often from methionine, aromatic moiety from phenylalanine or tyrosine, and the indole components from tryptophan [27,28]. Isothiocyanates are formed during the hydrolysis of glucosinolates. Using the example of glucotropaeolin (GTP), we will follow the stages of hydrolysis and the end products. During the hydrolysis of GTP with the participation of myrosinase, the glucose molecule and the hydroxysulfate ion are detached from the structure. As a result, a biologically active compound called benzyl isothiocyanate is formed, which is the major intermediate to benzylamine formation [29]. An understanding of the mechanisms of glucosinolates transformation is essential to determine their presence in plant extracts.
The applied chromatographic conditions allowed for the determination of twenty metabolites in the extracts of L. peruvianum that were visible in the negative ionization mode (also see Figures S1 and S2 in the Supplementary File), as presented in Table 1. The identification was based on high-resolution molecular weight measurements, analysis of their MS/MS spectra, retention times, open MS databases, and scientific literature and is in line with the previously published data of other researchers [14,26,30]. The total ion chromatogram (see Figure S3 in Supplementary File) of the selected extract, the extracted ion chromatograms of all glucosinolates whose presence was confirmed in the negative mode [26], and the MS/MS spectra of the tentatively identified structures are presented in Table S1 in the Supplementary File.
The obtained data clearly indicate that glucosinolates are the major constituents visible in the negative ionization mode in the majority of the tested samples. Among them, glucotropaeolin (GTP, G2)—an aromatic glucosinolate—eluted in the ninth minute was the major peak. Additionally, other representatives from different groups of sulphur-containing compounds were identified in the fingerprints of maca 50% EtOH (v/v) extracts: 1-isothiocyanato-9-methanesulfinylnonane (G6); aromatic glucosinolates, such as glucosinalbin (G1); and glucolimnanthin (G3). Indolyl glucosinolates: methoxyglucobrassicin (G4), 4-methoxyindolyl-3-hexylhydroxyglucosinolate (G5), and indolyl-3-hexyl-4-methyl-cyclohexaneglucosinolate (G9). Aliphatic glucosinolates: glucoalyssin (G8) and pent-4-enylglucosinolate (G7). Moreover, the HPLC-ESI-QTOF-MS/MS analysis revealed the presence of several organic acids, such as gluconic, malic, citric, succinic, and phenolic acids; benzoic acid and its derivatives; tannins like fucodiphlorethol; and fatty acids represented by trihydroxy-(3)-octadecadienoic acid and pinellic acid.

2.2. Chemometric Assessment

2.2.1. Data Preparation

Glucosinolates were found to constitute the leading role in the analyzed L. peruvianum extracts. Based on the qualitative determination of the extracts’ constituents presented in the Table 1, the further aim of the study was to compare the content of glucosinolates and to analyse their relative concentration in the extracts.
For that reason, the chemometric analyses of the presented dataset of L. peruviaum extracts were prepared in order to compare absolute (ABS) and relative (REL) glucosinolate compositions of the studied samples. For this purpose, the peak areas of nine compounds (G1–G9, see Table 1) were collected in triplicate from the recorded chromatograms of every extract injected into the HPLC-MS instrument. For further information on data preparation, please consult Materials and Methods (Section 3.4).

2.2.2. The “Absolute” Approach (ABS)

Principal component analysis (PCA) performed for the ABS dataset (Table S2 in the Supplementary File) resulted in the selection of the first three principal components (PCs), explaining 79.8% of the total variance (Table S3 in the Supplementary File). After the VARIMAX rotation, the following variables were associated with the resulting dimensions: (1) first dimension (Dim1)—G1, G2, G3, G8; (2) second dimension (Dim2)—G4, G9; (3) third dimension (Dim3)—G6, G7. The remaining G5 variable constituted the fourth dimension of the studied dataset, yet it was not included in the further analyses (Table S4 in the Supplementary File).
PCA and the correlation analysis (Figure 2 and Figure S4) have shown that the absolute extraction ratios of aromatic glucosinolates (G1, G2, G3) and aliphatic glucoalyssin (G8), constituting Dim1, were strongly correlated. This means that the content of G1, G2, and G3 is related to the content of G8. The higher content of the aromatic glucosinolates in the extracts of maca, the higher the amount of G8.
It is worth noting that the first dimension (Dim1), represented by a bimodal distribution, was responsible for the division of the dataset into two subsets of samples, which was also clearly demonstrated by the hierarchical clustering (HC, Figure 3). The obtained results classify the tested samples into two different groups—one of which collects juice samples obtained from fresh tubers, technologically altered gelatinized samples, maca tubers collected in the Ancash area, and yellow and purple tubers that seem to have different profiles from the remaining phenotypes—whereas the other collects large tubers of differently coloured phenotypes, mostly collected from Junin location.
Linear maps of samples (Figure 2), combined by the cluster dendrogram resulting from HC (Figure 3), have revealed that juices contained very similar amounts of glucosinolates, with BL_JU (No 7) slightly differentiating from the rest (Figure 3). These observations are in line with previous findings on the formation of glucosinolates in the plant tissues that underline that the process of oxidation is necessary to trigger the conformational changes in their structures [27]. Juices that were squeezed out from fresh tubers and immediately analysed did not contain high quantity of glucosinolates yet.
Moreover, the following samples collected in mid-October 2017: MIX_gel-AN, BL_gel-AN, BL_Bp-AN, RE_Bp-J, RE_gel-AN, RE_SM, RE_LA, and RE_LAgel contained similar amounts of glucosinolates (Figure 3). Interestingly, gelatinized samples in general contained low amounts of aromatic glucosinolates and G8 (Dim1) in comparison with non-gelatinized tubers, except for BL_gel-AN (located at the middle of the Dim1 axis). However, they all exhibited similar, medium ratios of G4, lower than unprocessed samples.
As a raw material containing a large quantity of starch, in order to improve their nutritional properties, maca tubers are subjected to a pretreatment process—gelatinization. The results presented above show that this process, which takes place at an elevated temperature, induces thermal decomposition of glucosinolates because their content in gelatinized samples is lower in relation to the untreated samples.
In the study, the composition and content of glucosinolates was found to be colour- and size- dependent. As previously mentioned, the above presented dendrogram (Figure 3) classifies the analysed extracts into two major groups. Interestingly, a majority of the extracts from black, grey, and red phenotypes tend to show a similar composition, whereas the extracts from yellow and purple phenotypes are more distant from other samples in terms of glucosinolate content. Specifically, small, grey tubers (GR_D-J, GR_Z-J) contained higher amounts of indolyl glucosinolates, i.e., G4 and G9 (Dim2) than the rest of the tubers. Large red, grey, and black tubers (A, B, X) were poor in G4 and G9, in contrast to the small ones (for which the trend was not applicable to purple and yellow species). Red tubers in general contained low amounts of indolyl glucosinolates (Dim2), while all the red species (excluding RE_Bp-J) contained smaller quantities of glucosinolates and G8, constituting the Dim1. Red and purple tubers produced lower amounts of G6 and aliphatic glucosinolate G7 (Dim3). On the contrary, black tubers excelled in ratios of G6 and G7 produced, with the exclusion of juiced samples (BL_JU) and extra small phenotypes (BL_Z-J). Generally, large tubers seem to differ in composition in comparison to small tubers, which can be seen in the above dendrogram. The glucosinolates’ profiles of yellow and purple phenotypes are alike. These results are in accordance with the study of Clement et al. [41], who calculated total glucosinolate content in their work on four phenotypes of maca. In their study, the grey phenotype was found to contain the highest amount of total glucosinolates (54.78 μmol/g d.w.), followed by the red (44.10 μmol/g d.w.), yellow (37.23 μmol/g d.w.), and purple (33.22 μmol/g d.w.) phenotypes. For that reason, the maca phenotypes tested herein are divided into two groups—one containing black, grey, and red tubers and another with yellow and purple tubers—seems to be in accordance with former observations on the total glucosinolates content.
Finally, the chicha morada (CH) species was also poor in G1, G2, G3, and G8 (Dim1), probably due to the fact that they were not dried sufficiently at the place of collection.

2.2.3. The “Relative” Approach (REL)

According to the authors’ knowledge, the relative approach to the content of glucosinolates in maca tubers was applied for the first time. This method of data analysis allows for a study on the relations between these metabolites in the maca extracts. Principal component analysis (PCA) performed for the REL dataset (Table S5 in the Supplementary File) also resulted in selection of three principal components (PCs) which explain 90.5% of the total variance (Table S6 in the Supplementary File). Correlation analysis, as well as the VARIMAX rotation, resulted in a linear map of variables in the space of three varivectors, further referred to as “dimensions” (Figure 4 and Figure S5 in the Supplementary File). This revealed that the relative amounts of indolyl glucosinolates G5 and G6 were strongly yet reversely correlated with aromatic glucosinolates G2 and G3 and indolyl glucosinolate G8, while all five compounds constituted the first dimension of the REL dataset, Dim1. These observations show that the synthesis of the individual groups of glucosinolates shall engage different biosynthetic pathways which are dependent on one another.
Additionally, aliphatic glucosinolate G7 was reversely correlated with G2, G3, and G8, yet PCA revealed its stronger connections with indolyl glucosinolate G9 and aromatic glucosinolate G1. Therefore, G1, G7, and G9 defined the second dimension of the studied dataset, Dim2. The relative ratios of indolyl glucosinolate G4 solely constituted the third dimension of the REL matrix, Dim3, which means that the biosynthesis of this compound has no relations to other glucosinolates analysed in this work. This dataset is well explained by a three-dimensional set of factors since no original variables escaped into higher dimensions (Tables S6 and S7 in the Supplementary).
Further PCA and HC studies (Figure 4 and Figure 5) have shown that the size of the tubers was correlated with the relative composition of the extracts, which additionally confirms the outcomes of the described above ABS approach.
The extracts that are located close to one another are characterized by similar relative proportions of glucosinolates.
Juices were observed to form a very separated subgroup within the REL dataset, with BL_JU being a “black sheep” of the juicy family, as it contains different proportions of the analysed glucosinolates as other samples. This differentiation was observed on the cluster dendrogram and the linear map of samples and mainly resulted from the first dimension (G2, G3, G5, G6, and G8). Gelatinization did not impact the relative composition of the glucosinolates within the black phenotype family, which explains that the eventual loss in the content of glucosinolates is the same for all compounds. However, the gelatinized MIX_gel-AN sample was a clear outlier in the REL dataset, substantially differentiating its relative composition from the rest of the extracts mainly due to its high Dim1 score (yet lower than juiced) and highest-by-far Dim2 score (G1, G7, and G9, Figure 4). A high Dim1 score indicates a marked relative concentration of G5 and G6, but low G2, G3, and G8, whereas the high Dim2 score of the sample is induced by increased relative contents of G1, G7, and G9 in comparison to other glucosinolates. Therefore, MIX_gel-AN was also the last sample to join the cluster dendrogram with a Euclidean distance above ~13 (Figure 5), which indicates that this sample exhibited extremely different relative composition in comparison to the others.
The REL dataset analysis shows that there are important correlations between certain groups of glucosinolates that may come from the similar, yet competing pathways of their biosynthesis in the plant.

2.3. Quantitative Determination of Glucotropaeolin Content in Lepidium peruvianum Extracts Using HPLC-ESI-QTOF-MS/MS

A comparative quantitative analysis of the glucotropaeolin (G2) content—the leading glucosinolate of the tested samples—was performed (Table 2) using the HPLC-MS methodology to visualize the existing differences between the analysed extracts. The analysis performed based on triplicate injections determined the content of this major glucosinolate in the 50% (v/v) alcoholic extract as ranging from 0 to 1.572%. Fresh juice samples did not contain G2, whereas its presence was confirmed in all remaining dried tuber samples.
Next, the relationship between the content of G2 and the color of tubers was investigated. For all five tested varieties of maca—black, purple, grey, yellow, and red—marked differences between the injections were noted. When analyzing the size B of all varieties, grey, red and black maca tubers were found to contain the highest quantity of this compound. The extract from black maca (BL_B-J) contained 0.965% of G2; from red maca (RE_B-J), 0.923%; whereas the extract from grey maca (GR_B-J) contained 1.010% (see Table 2). The smallest content of G2 was calculated for the purple variety: 0.239%. Interestingly, the level of this glucosinolate in the Chica morada (CH) variety was equal to 0.089%.
A drop in G2 concentration was noted for the gelatinized samples. Black gelatinized tubers (BL_gel-AN) contained 0.624% G2 (a drop from 0.965%), small gelatinized red tubers contained (RE _gel-AN) 0.103% G2 (without processing: 0.242%), and large red gelatinized maca tubers (RE_LAgel) 0.061% G2 (without processing: 0.344%).
Previous studies discuss the content of glucosinolates in maca samples of different origins. It has been observed that batches from different producers and differently coloured phenotypes significantly differed in quantity of glucosinolates. According to Clement et al. [41], the content of glucotropaeolin was the highest in lead-coloured (grey) hypocotyls compared to red, yellow, and purple phenotypes, which is in accordance with the herein presented studies. Another study by the authors [42] confirms an important role of cultivation and cultivation conditions in the biosynthesis of secondary metabolites. The extracts from the same phenotypes of maca (lead-colored, red, yellow, and purple) were examined, but they differed with collection site. The content of G2 in the lead-coloured phenotype was calculated as 44.23–45.35 µmol/g of dry matter from the Patala area and only 2.37–2.89 µmol/g of dry matter from Alpacayan area. Additionally, our studies confirmed marked differences between the samples collected in the Junín and Ancash areas, with a higher content of these secondary metabolites from the former group.

2.4. Inhibitory Activity of Individual Glucosinolates and the Total Extract towards AChE and BuChE

The obtained results confirmed the inhibitory activity of the tested total extracts and the standards.
Concerning the results of anti-enzymatic activity exhibited by the reference compounds, G4 was found to be the strongest component. Glucobrassicin derivative showed the best inhibitory properties (enzyme inhibition value: 57.6%) towards AChE among the tested compounds in relation to BuChE, for which the calculated inhibition percentage of G4 was equal to 39.9% (see, Table 3). Then, in descending order, for AChE G1 > G3 > G2, with 38.9, 37.05, and 36.2 for BuChE and 56.9, 56.5, and 55.95 for AChE, respectively. It should be underlined that G4 was the leading enzyme inhibitor, whereas the following compounds were characterized by a similar inhibitory strength.
The results obtained in the tests on reference compounds were set together with the outcomes of the anti-enzymatic assays performed for some selected total extracts. All samples were tested at the concentration of 10 mg/mL, so they can be directly compared with one another. It is important to note that, based on the results of reference compounds, the share of glucosinolates in the total enzyme inhibitory activity can be marked, especially for AChE enzymes. The AChE inhibitory activity of the total extracts ranges from approximately 15–25%, whereas the activities of single glucosinolates are within the range of 55–57%. Based on these observations, they can noticeably influence the total activity of the extract. In the case of BuChE, where the activity of the reference compounds is between 37–40%, the total inhibition rate of the extracts is higher (71–77%), so it must be influenced also by other components present in the extracts.
However, the quantitative data obtained in this study seem to correspond with the results of the enzymatic assay presented herein. Both AChE and BuChE inhibition decrease together with a decrease in size of tubers and with a decrease in the content of glucosinolates, as can be noted in the case of black maca (see Table 3).
Additionally, color-dependent variations were observed between the samples. Black, grey, and red tubers (size B) show a similar anti-enzymatic potential, whereas the yellow variety is both a weaker enzyme inhibitor and a weaker source of glucosinolates.
Rubio et al. reported that aqueous and hydroalcoholic extracts of Lepidium meyenii reduced activity of AChE in the brain by 45% after memory impairment caused by scopolamine [43]. Three different maca varieties (black, yellow, and red) were tested to evaluate the effects on memory and depression in ovariectomized mice. Of the three cultivars, the black maca variety showed the best activity [44].
L. peruvianum is rich in compounds called glucosinolates. It has been proven that glucosinolates have the ability to inhibit AChE in a dose-dependent manner. Blazevic et al. demonstrated that these secondary metabolites from Bunias erucago extracts, like gluconapin, glucoerucin, and glucoraphanin, inhibited the AChE enzyme by about 53% [45], whereas glucosinalbin did so by about 40.9% [25]. Recent studies also confirm that the products of glucosinolate hydrolysis, isothiocyanates, have AChE-inhibiting properties. Phenyl derivatives of isothiocyanates show promising AChE-inhibitory properties [24].
In the presented work, we used a standard method evaluating the AChE activity of newly characterized glucosinolates present in maca. Due to the use of unified methodology, we are able to make direct comparisons with previously characterized compounds. We must emphasize that the inhibitory activity of glucosinolates was within one order of magnitude (at a similar level) of a number of efficient cholinesterase inhibitors—terpenes and phenylpropanoids—previously studied by us in the very same manner [46]. Moreover, glucosinolates exerted much higher inhibitory activity than physodic acid isolated from Hypogymnia physodes, characterized in one of our recent works [47]. Although the cholinesterase-inhibiting activity of glucosinolates is lower than the activity of very efficient inhibitors—flavonoids (e.g., quercetin, rutin, or kaempferol) recently identified in selected fruits [48])—this activity can be judged applicable.
The results presented herein underline that next to the colour of Peruvian maca phenotype, the tuber size is an important factor in terms of glucosinolate content. This observation is crucial, as so far there is no specific method to classify tubers according to their size, and the cultivated tubers from a given plantation are collected together. The obtained results draw attention to the fact that in the context of glucosinolate content and possible preventive use of the tubers in memory impairment conditions, the important issue remains proper crop collection.

2.5. Modelling

The ligand–protein-binding energies found during docking simulations are given in Table 4 and graphically illustrated in Figure 6. Apart from compound 1-isothiocyanato-9-methanesulfinylnonane, which displays significantly distinct structure, the rest of compounds exhibit relatively similar magnitudes of binding energies. Independently of the considered protein, the strongest binding is always exhibited by compound G9, whereas the weakest is exhibited by 1-isothiocyanato-9-methanesulfinylnonane. All negative energy values clearly speak for favorable binding in all considered cases.
A series of interesting correlations is revealed when looking closer at calculated binding energies (Figure 6). Firstly, the binding energies are fairly well but negatively correlated with the molecular sizes of docked ligands (expressed here as molecular volume). This suggests that the increased size of the -R substituent (see Figure 6) positively influences the binding strength and, therefore, inhibition potency. Secondly, the binding energies obtained for a series of X-ray crystal structures of human BuChE are strongly correlated with analogous values for equine BuChE structure obtained from homology modeling. This result is expected due to high sequence identity between these two proteins. Nevertheless, it also speaks for the correctness of the homology model. Finally, a good correlation is observed between binding energies obtained for AChE and BuChE. This also is expected due to the close nature of both enzymes and structural similarities in their binding sites. It is worth noting that in the last two cases, except for a good correlation, there exist only minor deviations in absolute binding energies determined for the same ligand and a given pair of proteins (on average, 0.35 and 0.57 kcal/mol, respectively).
The results of the docking studies have been analysed with respect to the mechanistic interaction pattern that may be significant in the context of interpretation of the obtained binding energies. The summary given below relies on analyzing the ligand–protein contacts that take place if the distance between any corresponding atom pair is smaller than the arbitrarily accepted value of 0.4 nm. Figure 7 shows the superposition of all most favorable ligand poses with the most characteristic ligand-protein interactions. Let us remark that compound 1-isothiocyanato-9-methanesulfinylnonane, sharing little structural similarity with remaining compounds, was excluded from the discussion below.
All ligands prefer roughly the same binding position in the enzyme cavity, which enables them to block catalytic sites (proximity to the tryptophan and histidine residues). The relative molecular orientation is (in majority of the cases) the same among compounds interacting with the same protein.
In the case of AChE, the –R substituent (usually of either aromatic or aliphatic character) interacts with the cluster of aromatic amino acid sidechains (Phe330, Phe331, and Tyr334). The corresponding interactions are of either p-p, CH-p stacking or simply hydrophobic nature. The opposite site of the ligand molecule (glucopyranose moiety) interacts with hydrogen-bonding-acceptor-containing amino acid residues (Glu199, Ser200, Tyr130, and His440) mainly by hydrogen bonding. Additionally, the CH-p stacking involving Trp84 and hydrophobic patch in the glucose ring is observed (This type of interaction is characteristic of non-charged carbohydrate–protein binding [49]. As indicated by the dispersion of the poses (including secondary and tertiary ones), the glucopyranosidic moiety is much more conformationally flexible in comparison to the –R substituent. Finally, the charged sulfate group, located in the central part of the ligand molecule, creates some attractive, hydrogen-bonding-mediated interactions with Ser122 and Asn85 (backbone fragment). As the sulfate moiety is located close to the entrance to the binding cavity, it is also possible that multiple interactions with polar water molecules are possible (as no explicit solvent was considered during the docking procedure, this has not been definitively confirmed).
In the case of BuChE, an important reorientation of the ligand in the binding cavity may be observed in comparison to AChE. Namely, the –R substituent is shifted toward Trp110 (structural analogue of Trp84 in Ache), interacting with it through p-p stacking. P-p interactions are possible also in the case of Phe357 but are probably significantly less intense due to increased R-Phe357 distance. At the same time, the glucose moiety exhibits hydrogen bonding with backbone fragments of Ser315, Gly145, and the side chain of His466 (structural analogue of His440 in AChE). Thus, the whole ligand molecule is slightly shifted toward the more polar edge of the binding cavity. The sulfate group is again located closest to the binding cavity entry and, thus, probably interacts with water molecules. Additionally, its hydrogen binding with Thr148 and Asp98 (backbone) were identified.
The observed differences between most favorable orientations of the docked ligands interacting either with AChE or BuChE can be ascribed to the diverse composition of the corresponding binding cavities. The cavity of BuChE lacks the cluster of aromatic residues Tyr-Phe-Phe which was responsible for favorable interactions with –R in the case of AChE. Therefore, reoriented –R moiety seeks attractive interactions with Trp110, whereas the glucose group can still interact with Ser226 and His466.
Finally, let us note that in the case of some compounds and—especially in secondary, tertiary, or quaternary poses—systematic reorientations of the ligand in the binding cavity were observed. Namely, the position of the –R substituent and glucopyranosidic group are swapped. In other words, AChE-characteristic poses become close to the BuChE ones and vice versa. Such an alternative interaction pattern is reasonable to accept, as the carbohydrate moiety can exhibit attractive interactions with either Trp or Phe/Tyr moieties. The same can be said about the aromatic –R substituent. Although the present docking study classifies such binding modes as only secondary ones, the accompanying binding energies are close to the most favorable values (difference by at most 1 kcal/mol). Thus, it cannot be excluded that in real systems, when the solvent and temperature effects are present, the ligand position can be a superposition of those two alternative orientations.

3. Materials and Methods

3.1. Plant Material and Extraction Conditions

Plant material studied in this paper represents dried tubers of Maca (Lepidium peruvianum syn. L. meyenii) collected in the highlands of the Peruvian Andes in the Junín plateau (2017 and 2019 harvesting seasons) and Ancash cultivation area (2017 harvest) as summarised in Table 5.
Maca cultivated in the Junín plateau was harvested on two occasions. Sampling J-17 was conducted on 1 October 2017 from plantation in Tarma Province, Tapo District, some 30 km west of Tarma city at 4174 m above sea level (a.s.l. ±20 m). Tubers were sampled from an area of approximately 250 m radius surrounding the geolocation coordinates 11°24′12″S, 75°35′48″W according to GPS reading.
Maca tubers from the 2019 sampling in Junín, (J-19), originate from the cultivation site on the northeast-facing slopes of the Andes and were harvested on 8 October 2019 at an elevation of 4156 m a.s.l. (±25 m) in an area of approximately 300–350 m radius surrounding the GPS location point 11°06′23″S 75°58′47″W. All traditionally harvested maca tubers of mixed phenotypes harvested in the Junín plateau are from samplings J-A and J-B were transported to a drying location site in Tarma (3148 m a.s.l.) at 11°25′19″ S 75°41′47″ W, where tubers were traditionally dried, separated into individual phenotypes black, red, purple and yellow, and selected into different sizes within individual phenotypes as per the example given in Figure 8. Selected fresh maca tubers from the J-19 location were chopped and juice was squeezed out of the crushed tubers for use in analyses. Separated fresh pulp was discarded.
Maca tubers cultivated in the Ancash region (A-17) were collected in mid-October 2017 in Ancash, Huanuco in Margos, Margos District, 105 km by car southeast of Huallanca from the cultivation area at 4242 m a. s. l. (±35 m) of an approximately 250 m radius facing the northeast slope of a hill surrounding the GPS point 9°46′8″ S, 76°59′06″ W. Harvested maca tubers were traditionally sun-dried at a plantation site in Ancash, transported for selecting phenotypes and tuber sizes in Huallanca storage, delivered to a processing plant in Lima to be pulverised and packaged, then sent by courier to Lublin laboratory, the Department of Pharmacognosy with Medicinal Plants Garden. The voucher specimen of every sample is kept by the corresponding author in the Department of Pharmacognosy with Medicinal Plants Garden, Medical University of Lublin. All original specimens used in the study were authenticated by Professor Henry O. Meissner directly after collection.
All L. peruvianum extracts used in the study were obtained with 50% aqueous ethanol [37] (v/v) by ultrasound-assisted extraction technique. Three extraction cycles of 30 min each were performed on 2 g of the same plant material using fresh portions of extracting solvent at room temperature. The selected solvent to solid ratio was 1:1 (v/v). After extraction, the extracts were joined and evaporated to dryness using a vacuum evaporator at 50 °C. The obtained residue was weighed and stored in a refrigerator at 5 °C until the analysis. Juice samples were produced by cutting fresh bulbs of L. peruvianum and placing them into a slow juicer. The collected samples were later centrifuged, filtered through a syringe filter (nominal pore size of 0.2 µm), and used for chromatographic analysis.

3.2. Reagents

Ethanol at 96% concentration was purchased from Avantor Performance Materials (Gliwice, Poland) for extraction. HPLC-MS purity water, acetonitrile and formic acid were obtained from J.T. Baker (Phillipsburg, NJ, USA). The reagents used in the AChE and BuChE inhibition assays, namely acetylthiocholine iodide (ATChI), butyrylcholine iodide (BuTChI), 5,5′-Dithiobis (2-nitrobenzoic acid) (DTNB), Tris buffer, sodium chloride (NaCl), and magnesium chloride hydrate (MgCl2·6H2O), were provided by Sigma Aldrich (St. Louis, MO, USA). The standard of glucotropaeolin (purity > 98%) used to plot the calibration curve was purchased from Extrasynthese (Genay, France). For bioactivity studies, the following standards of glucosinolates were obtained at purities exceeding 95%: glucotropaeolin (G2), glucolimnanthin (G3), methoxyglucobrassicin (G4), and glucosinalbin (G1).

3.3. Qualitative and Quantitative Analysis of the Extracts by High Performance Liquid Chromatography Coupled with Mass Spectrometry (HPLC-ESI-QTOF-MS/MS)

The HPLC-ESI-QTOF-MS/MS platform from Agilent Technologies (Santa Clara, California, USA) was used to analyse the composition of L. peruvianum extracts. For the separation process, an HPLC chromatograph (1200 series) equipped with a Zorbax Eclipse Plus RP-18 chromatography column (150 mm × 2.1 mm; dp = 3.5 µm), a degasser (G1322A), a binary pump (G1312C), an autosampler (G1329B), a photodiode array detector (G1315D), and a mass spectrometer (G6530B) were used. Agilent MassHunter Workstation Software (version B.10.00) was used to acquire the MS spectra and process the data.
The temperature of the HPLC thermostat was set at 25 °C and the UV detection wavelengths at 254, 280, 320, and 365 nm. The operating wavelength range of UV-Vis/DAD detector was 190–600 nm. Chromatographic separation was performed using a 10 µL injection volume with a flow rate of 0.2 mL/min in a 45 min gradient elution program. Mobile phases consisted of eluent A (0.1% formic acid in water, v/v) and eluent B (acetonitrile solution with 0.1% formic acid added). The gradient elution was as follows: 0 min—1% of B, 4 min—20% of B, 20 min—60% B, 35–38 min—95% B, 39–45 min—1% of B.
The mass spectrometer measurements were carried out under the following conditions: gas temperature and shield gas temperature were 350 and 325 °C, respectively; gas flows were set at 12 L/min; capillary voltage at 3500 V; fragmentor voltage at 120 V; collision energies at 10 and 20 V; skimmer voltage at 65 V; nebulizer pressure at 35 psig. The collected spectra were scanned in the m/z 100–1700 Da range in the negative ionization mode that visualizes glucosinolates, as well as in the positive ionization mode. Two of the most intense signals seen in the TIC spectrum were automatically fragmented to obtain MS/MS spectra. After collecting two spectra for a given m/z value, the selected signals were excluded for 0.2 min from further fragmentation.

3.4. Statistical Analysis of Data

3.4.1. Chemometric Analyses

All the chemometric analyses and visualizations were performed using R v4.1.2 [50] programming language in RStudio [51] software with pracma [52], factoextra [53], matlib [54], and corrplot [55] packages installed.

3.4.2. Absolute (ABS) Approach vs. Relative (REL) Approach

In the ABS approach, nine variables (nine compounds) representing the absolute detector responses, i.e., peak areas for a given compound, were autoscaled in order to enable standard procedures of exploratory analysis, i.e., principal component analysis (PCA) and hierarchical clustering (HC) (Table S2 in the Supplementary File). In the REL approach, all peak areas for nine traced compounds of every given sample were summed to 100%, hence the amount of every compound was represented as a percentage share of a total detector response of a given sample (Table S5 in the Supplementary File). After such a percentage rescaling, all the variables were also autoscaled in a similar fashion to the ABS approach.

3.4.3. Principal Component Analysis (PCA)

After the standard, formal decomposition of the covariance matrices were calculated for the ABS (absolute) and REL (relative) autoscaled datasets, three principal components (PCs) were first considered relevant in both cases. After the selection of the relevant PCs, their vectors were rotated in space in order to maximize the values of correlation coefficients between the original variables and the three orthogonal factors using the VARIMAX algorithm. In every case, samples scores in the space of the resulting varivectors (dimensions) were calculated by multiplying the matrix of the original variables’ loadings in the space of varivectors by the matrix of the autoscaled ABS/REL dataset.

3.4.4. Correlation Analysis

Matrices of the coefficients of linear correlation (R) were prepared using Pearson’s definition of linear correlation.

3.4.5. Hierarchical Clustering (HC)

Hierarchical clustering (HC) for the ABS and REL approaches was performed using the Ward method based on a standard matrix of Euclidean distances between the samples.

3.5. Quantitative Analysis of Glucotropaeolin in the Analysed Extracts by HPLC-ESI-QTOF-MS/MS

Glucotropaeolin, as the most abundant glucosinolate in the tested extracts, was selected for the quantitative analysis to present the differences in its content between the tested extracts. The quantitative results were prepared for the standard compound (purity >98%, Extrasynthese, Genay, France) on the basis of its calibration curve, plotted for 7 different concentrations within the range of 0.0039–0.25 mg/mL that were obtained by diluting the stock solution with a concentration of 1 mg/mL. The same volumes (10 µL) of every solution were injected on the column and developed in the same above-described chromatographic method as for the qualitative assessment. The linearity range was observed for the concentration range of glucotropaeolin from 0.0039–0.0625 mg/mL. For these concentrations, the calculated calibration curve equation was y = 1,970,017,678x + 2,232,216 and the coefficient of determination (R2) was equal to 0.99084.
The dried extracts from different varieties of L. peruvianum were redissolved in 50% ethanol (v/v) to the final concentration of 5 mg/mL, filtered through nylon syringe filters (nominal pore size 0.22 um), and subjected to chromatographic analysis in triplicate.

3.6. AChE- and BuChE Inhibition Assay

The selected total extracts (BL_A-J, BL_B-J, BL_C-J, BL_D-J, BL_Z-J, BL_B-J, RE_B-J, GR_B-J, YE_B-J) dissolved in DMSO at a concentration of 10 mg/mL and reference compounds (GL, GTP, GB, GS) at a concentration of 1 mg/mL. These samples were tested for their BuChE- and AChE-inhibitory activity according to the previously published protocol described below.
The assay was essentially performed as described by Szwajgier et al. (2021) [56]. AChE and BuChE were at a concentration of 0.28 U/mL, ATChI and BuTChI were at a concentration of 1.5 mmol/L, and DTNB (0.3 mmol/L) contained 11.89 mg DTNB, 58.44 mg NaCl, and 40.66 mg MgCl2·6H2O, all per 100 mL of Tris buffer. The measurement of the absorbance at 405 nm was performed using a GloMax® Explorer GM3500 plate reader (Promega, Madison, WI, USA) using the Software version 3.2.3. All samples were tested in triplicate. Three blank tests were also performed together with one false positive control for every samples. In the first blank test, the sample was replaced with 10 µL DMSO, the second contained 10 µL Tris instead of the sample; in the third, the enzyme and the sample were replaced with Tris buffer and “false-positive” blank sample according to Rhee et al. [57].

3.7. Molecular Docking of Glucosinolates to AChE and BuChE

The ligand molecules were drawn by using Avogadro 1.1.1 [58] and initially optimized within the UFF force field [59] (5000 steps, steepest descent algorithm). Flexible and optimized ligand molecules were docked into the binding pocket of the protein structures either found in the PDB database or prepared by using a homology modeling procedure. In the case of AChE, the 1EEA PDB entry (AChE from electric eel, X-ray resolution: 0.45 nm) was used. In the case of BuChE, no equine BuChE structure was deposited in PDB. Therefore, seven different human BuChE structures (PDB: 6EQP, 4BDS, 1XLW, 4XII, 6ESJ, 6ESY, and 6QAC of resolutions equal to 0.235, 0.210, 0.210, 0.270, 0.298, 0.280, and 0.277 nm, respectively) and, additionally, a homology model prepared on the basis of the highest resolution 4BDS structure, were used. Note that the equine BuChE has ~90% amino acid similarity with the human sequence [60] and nearly all differences are located outside binding cavity. The homology model was prepared by using the SWISS-MODEL server [61,62] and all inherent, default procedures.
Docking simulations were carried out in the AutoDock Vina software [63]. The procedure was performed within the cubic region of dimensions of 18 × 18 × 18 Å3, which covers all the originally cocrystallized ligands present in the considered PDB structures as well as the closest amino acid residues which exhibit contact with those ligands. All the default procedures and algorithms implemented in AutoDock Vina were applied during the docking procedure. The torsional angles in ligand molecules as well as the selected amino acid sidechains within the binding cavities of both proteins were allowed to rotate. In the case of human BuChE, the predicted binding energies were averaged over all seven protein structures. The visual inspections of each pose of the docked ligand (after aligning with respect to the given protein structure) were carried out in order to assure that the partial energies contributing the final, averaged, binding energy value, correspond to the structurally analogous orientations. Only the lowest energy values, corresponding to the most energetically favorable complexes, were considered during subsequent analysis. The procedure was initially validated by docking the ligands cocrystalized with proteins and present in the 6EQP and 1E3Q (AChE from Torpedo californica). The graphical illustration of the validation results is given in Figure 7.

4. Conclusions

The results presented above provide evidence of the existence of substantial antienzymatic activity of several well-identified maca compounds against AChE and BuChE, thus confirming that L. peruvianum is a rich source of active glucosinolates, thus showing promise for possible application of selected maca phenotypes in developing therapeutic products for treating medical conditions related to memory impairment. It was proved that the potential of the plant is related both to the phenotype (color) and to the size of Peruvian maca tubers.

Supplementary Materials

The following supporting information can be downloaded at: https://www.mdpi.com/article/10.3390/ijms23094858/s1.

Author Contributions

Conceptualization, H.O.M., W.K.-K.; methodology, W.K.-K., D.S., P.S., T.L., W.P. formal analysis, D.T., W.K.-K., D.S., E.B.-W., P.S., T.L., W.P. investigation, D.T., W.K.-K., P.S., T.L., D.S.; resources, H.O.M.,W.K.-K.; data curation, all authors; writing—original draft preparation, all authors; writing—review and editing, all authors.; supervision, W.K-K., H.O.M.; project administration, W.K.-K.; funding acquisition, W.K.-K. All authors have read and agreed to the published version of the manuscript.

Funding

This research was funded by D.T. and W.K.-K. were funded by the Medical University of Lublin, Poland (project No PBmb/sd Nr 70, DS 24 and UM-CTW-10/2021.

Institutional Review Board Statement

Not applicable.

Informed Consent Statement

Not applicable.

Data Availability Statement

The supporting data are available the Supplementary File.

Acknowledgments

The assistance of Dina Guere VEGA and Hernan Juri Sotelo Aguilar from SILMAX- E.I.R.L., Tarma, Junín, Peru in collecting and sorting out maca tubers from contracted maca plantations in the Junín plateau area, followed by selecting maca phenotypes and sorting different-sized maca tubers for this study is greatly appreciated. Thanks are due to Messrs. Paulo and Sergio Paulet of Natural Peru, Lima for their assistance in processing maca tubers collected in Junín, high Peruvian Andes and preparing several maca samples for comparative study. Thanks to Carolina Huarto for the provision of, management, and collection of maca samples; for their processing and dispatch from Lima to the laboratory in Poland; and for maintaining documentation of sample preparation from the tuber harvest to courier dispatches for analysis in laboratories in Poland. Appreciation is also due to the management of NHI Manufacture (NHI-USA), subsidiary Symphony Natural Health, 2550 S Decker Lake Blvd Ste 28, West Valley City, UT 84119, USA for generously providing access to a monitored maca plantation in Ancash, Cordillera Blanca, Peruvian Andes, and use of their processing plant in Lima for preparing samples for dispatch to the laboratory in Poland. We also wish to thank Professor Jan Mazerski from the Department of Pharmaceutical Technology and Biochemistry, Faculty of Chemistry, Gdańsk University of Technology for his invaluable suggestions regarding chemometric analysis.

Conflicts of Interest

The authors declare no conflict of interest. The manufacturer (NHI-USA) of the commercial maca products who donated part of the samples used in this comparative study had no role in experimental design, analysis, interpretation of data, or in the decision to publish the results. Apart from H.O.M. financial interest in TTD International Pty Ltd., there is no conflict of intellectual or other interest for disclosure in relation to the material used, research results, or the conclusions presented in this publication, which are purely of academic nature.

References

  1. Alzheimer’s Association. 2014 Alzheimer’s disease facts and figures. Alzheimers Dement. 2014, 10, e47–e92. [Google Scholar]
  2. Perry, E.; Howes, M.J.R. Medicinal plants and dementia therapy: Herbal hopes for brain aging? CNS Neurosci. Ther. 2011, 17, 683–698. [Google Scholar] [CrossRef] [PubMed]
  3. Pohanka, M. Inhibitors of Cholinesterases in Pharmacology: The Current Trends. Mini.-Rev. Med. Chem. 2019, 20, 1532–1542. [Google Scholar] [CrossRef] [PubMed]
  4. Haake, A.; Nguyen, K.; Friedman, L.; Chakkamparambil, B.; Grossberg, G.T. An update on the utility and safety of cholinesterase inhibitors for the treatment of Alzheimer’s disease. Expert Opin. Drug Saf. 2020, 19, 147–157. [Google Scholar] [CrossRef] [PubMed]
  5. Meissner, H.O.; Mscisz, A.; Kedzia, B.; Pisulewski, P.; Piatkowska, E. Peruvian maca: Two scientific names Lepidium meyenii Walpers and Lepidium peruvianum Chacon–Are they phytochemically-synonymous? Int. J. Biomed. Sci. 2015, 11, 1–15. [Google Scholar]
  6. Meissner, H.O.; Mscisz, A.; Reich-Bilinska, H.; Mrozikiewicz, P.; Bobkiewicz-Kozlowska, T.; Kedzia, B.; Lowicka, A.; Barchia, I. Hormone-balancing effect of pregelatinized organic Maca (Lepidium peruvianum Chacon): (III) Clinical responses of early-postmenopausal women to Maca in double blind, randomized, placebo controlled, crossover configuration, outpatient study. Int. J. Biomed. Sci. 2006, 2, 375–394. [Google Scholar]
  7. Meissner, H.O.; Xu, L.; Wan, W.; Yi, F. Glucosinolates profiles in Maca phenotypes cultivated in Peru and China (Lepidium peruvianum syn. L. meyenii). Phytochem. Lett. 2019, 31, 208–216. [Google Scholar] [CrossRef]
  8. Meissner, H.O. The Unique Powers of the Maca Tuber: Scientific Facts behind Traditional Wisdom. (in German “Die Einzigartigen Krafte der Maca-Wurzel. Wissentshaftliche Facten Hinter Traditionallem Wissen”), 1st ed.; Michaels Verlag & Vertrieb GmbH: Munich, Germany, 2014; pp. 1–432. [Google Scholar]
  9. Gonzales, G.F. Ethnobiology and Ethnopharmacology of Lepidium meyenii (Maca), a Plant from the Peruvian Highlands. Evid. Based Complement. Alternat. Med. 2012, 2012, 193496. [Google Scholar] [CrossRef] [Green Version]
  10. Meissner, H.O.; Kedzia, B.; Mrozikiewicz, P.M.; Mscisz, A. Short and long-term physiological responses of male and female rats to two dietary levels of Pre-Gelatinized Maca (Lepidium peruvianum Chacon). Int. J. Biomed. Sci. 2006, 2, 13–28. [Google Scholar]
  11. Beharry, S.; Heinrich, M. Is the hype around the reproductive health claims of maca (Lepidium meyenii Walp.) justified? J. Ethnopharmacol. 2018, 211, 126–170. [Google Scholar] [CrossRef] [Green Version]
  12. Meissner, H.O.; Mscisz, A.; Reich-Bilinska, H.; Kapczynski, W.; Mrozikiewicz, P.; Bobkiewicz-Kozlowska, T.; Kedzia, B.; Lowicka, A.; Barchia, I. Hormone balancing effect of pre-gelatinized organic Maca (Lepidium peruvianum Chacon): (II) physiological and symptomatic responses of early-postmenopausal women to standardized doses of Maca in double blind, randomized, placebo-controlled, multicenter clinical study. Int. J. Biomed. Sci. 2006, 2, 360–374. [Google Scholar] [PubMed]
  13. Meissner, H.O.; Mrozikiewicz, P.; Bobkiewicz-Kozlowska, T.; Mscisz, A.; Kedzia, B.; Lowicka, A.; Reich-Bilinska, H.; Kapczynski, W.; Barchia, I. Hormoneb alancing effect of pre-gelatinized organic Maca (Lepidium peruvianum Chacon): (I) Biochemical and pharmacodynamic study on Maca using clinical laboratory model on ovariectomized rats. Int. J. Biomed. Sci. 2006, 2, 260–272. [Google Scholar] [PubMed]
  14. Wang, Y.; Wang, Y.; Mcneil, B.; Harvey, L.M. Maca: An Andean crop with multi-pharmacological functions. Food Res. Int. 2007, 40, 783–792. [Google Scholar] [CrossRef]
  15. Dini, A.; Migliuolo, G.; Rastrelli, L.; Saturnino, P.; Schettino, O. Chemical composition of Lepidium meyenii. Food Chem. 1994, 49, 347–349. [Google Scholar] [CrossRef]
  16. Piacente, S.; Carbone, V.; Plaza, A.; Zampelli, A.; Pizza, C. Investigation of the tuber constituents of maca Lepidium meyenii walp.). J. Agric. Food Chem. 2002, 50, 5621–5625. [Google Scholar] [CrossRef] [PubMed]
  17. Li, G.; Ammermann, U.; Quiros, Q.F. Glucosinolate contents in Maca (Lepidium peruvianum Chacon) Seeds, Sprouts, Mature plants and Several Derived Commercial Products. Econ. Botany 2001, 55, 255–262. [Google Scholar] [CrossRef]
  18. Muhammad, I.; Zhao, J.; Dunbar, D.C.; Khan, I.A. Constituents of Lepidium meyenii ‘maca’. Phytochemistry 2002, 59, 105–110. [Google Scholar] [CrossRef]
  19. Zhang, L.; Wang, L.; Zhao, B. Relationship between glucosinolates and colors in eight accessions of maca. J. Food Eng. Technol. 2005, 4, 45–52. [Google Scholar]
  20. Vaughn, S.F.; Berhow, M.A. Glucosinolate hydrolysis products from various plant sources: pH effects, isolation, and purification. Ind. Crops Prod. 2005, 21, 193–202. [Google Scholar] [CrossRef]
  21. Lee, S.G.; Nwumeh, R.; Jez, J.M. Structure and mechanism of isopropylmalate dehydrogenase from Arabidopsis thaliana: Insights on leucine and aliphatic glucosinolate biosynthesis. J. Biol. Chem. 2016, 291, 13421–13430. [Google Scholar] [CrossRef] [Green Version]
  22. Melrose, J. The Glucosinolates: A Sulphur Glucoside Family of Mustard Anti-Tumour and Antimicrobial Phytochemicals of Potential Therapeutic Application. Biomedicines 2019, 7, 62. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  23. Dinkova-Kostova, A.T.; Kostov, R.V. Glucosinolates and isothiocyanates in health and disease. Trends Mol. Med. 2012, 18, 337–347. [Google Scholar] [CrossRef] [PubMed]
  24. Bur Cul, A.Ã.F.; Generali, C.; Mekini, C.B.I.; Radan, M.; Rollin, P.; Bla Zevi, C.E.Ã.I. Isothiocyanates: Cholinesterase inhibiting, antioxidant, and anti-inflammatory activity. J. Enzym. Inhib. Med. Chem. 2018, 33, 577–582. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  25. Blažević, I.; Đulović, A.; Čikeš Čulić, V.; Burčul, F.; Ljubenkov, I.; Ruščić, M.; Generalić Mekinić, I. Bunias erucago L.: Glucosinolate Profile and In Vitro Biological Potential. Molecules 2019, 24, 741. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  26. Zhou, Y.; Li, P.; Brantner, A.; Wang, H.; Shu, X.; Yang, J.; Si, N.; Han, L.; Zhao, H.; Bian, B. Chemical profiling analysis of Maca using UHPLC-ESI-Orbitrap MS coupled with UHPLC-ESI-QqQ MS and the neuroprotective study on its active ingredients. Sci. Rep. 2017, 7, 1–14. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  27. Traka, M.H. Health Benefits of Glucosinolates. Adv. Bot. Res. 2016, 80, 247–279. [Google Scholar]
  28. Blažević, I.; Montaut, S.; Burčul, F.; Olsen, C.E.; Burow, M.; Rollin, P.; Agerbirk, N. Glucosinolate structural diversity, identification, chemical synthesis and metabolism in plants. Phytochemistry 2020, 169, 112100. [Google Scholar] [CrossRef]
  29. Esparza, E.; Yi, W.; Limonchi, F.; Cosio, E.G. Glucosinolate catabolism during postharvest drying determines the ratio of bioactive macamides to deaminated benzenoids in Lepidium meyenii (maca) root flour. Phytochemistry 2020, 179, 112502. [Google Scholar] [CrossRef]
  30. Wang, S.; Zhu, F. Chemical composition and health effects of maca (Lepidium meyenii). Food Chem. 2019, 288, 422–443. [Google Scholar] [CrossRef]
  31. Geng, P.; Sun, J.; Chen, P.; Brand, E.; Frame, J.; Meissner, H.; Stewart, J.; Gafner, S.; Clark, S.; Miller, J.; et al. Characterization of Maca (Lepidium meyenii/Lepidium peruvianum) Using a Mass Spectral Fingerprinting, Metabolomic Analysis, and Genetic Sequencing Approach. Planta Med. 2020, 86, 674–685. [Google Scholar] [CrossRef]
  32. Al Kadhi, O.; Melchini, A.; Mithen, R.; Saha, S. Development of a LC-MS/MS Method for the Simultaneous Detection of Tricarboxylic Acid Cycle Intermediates in a Range of Biological Matrices. J. Anal. Methods Chem. 2017, 2017, 1–12. [Google Scholar] [CrossRef] [Green Version]
  33. Perez, C.J.; Conceição, R.S.; Ifa, D.R. Chemical profiling and separation of bioactive secondary metabolites in Maca (Lepidium peruvianum) by normal and reverse phase thin layer chromatography coupled to desorption electrospray ionization-mass spectrometry. J. Mass Spectrom. 2021, 56, e4690. [Google Scholar] [CrossRef]
  34. Shirako, S.; Kojima, Y.; Tomari, N.; Nakamura, Y.; Matsumura, Y.; Ikeda, K.; Inagaki, N.; Sato, K. Pyroglutamyl leucine, a peptide in fermented foods, attenuates dysbiosis by increasing host antimicrobial peptide. NPJ. Sci. Food 2019, 3, 1–9. [Google Scholar] [CrossRef]
  35. Chen, R.; Wei, J.; Gao, Y. A review of the study of active components and their pharmacology value in Lepidium meyenii (Maca). Phyther. Res. 2021, 35, 6706–6719. [Google Scholar] [CrossRef]
  36. Andini, S. Antimicrobial Isothiocyanates from Brassicaceae Glucosinolates: Analysis, Reactivity, and Quantitative Structure-Activity Relationships. Ph.D. Thesis, Wageningen University and Research, Wageningen, The Netherlands, 2020. [Google Scholar]
  37. Meissner, H.O.; Mscisz, A.; Mrozikiewicz, M.; Baraniak, M.; Mielcarek, S.; Kedzia, B.; Piatkowska, E.; Jólkowska, J.; Pisulewski, P. Peruvian Maca (Lepidium peruvianum): (I) phytochemical and genetic differences in three Maca phenotypes. Int. J. Biomed. Sci. 2015, 11, 131–145. [Google Scholar]
  38. Isaza Martínez, J.H.; Torres Castañeda, H.G. Preparation and Chromatographic Analysis of Phlorotannins. J. Chromatogr. Sci. 2013, 51, 825–838. [Google Scholar] [CrossRef] [Green Version]
  39. Fuchs, D.; Hamberg, M.; Sköld, C.M.; Wheelock, Å.M.; Wheelock, C.E. An LC-MS/MS workflow to characterize 16 regio- and stereoisomeric trihydroxyoctadecenoic acids. J. Lipid. Res. 2018, 59, 2025–2033. [Google Scholar] [CrossRef] [Green Version]
  40. Nadeem, M.; Waseem Mumtaz, M.; Danish, M.; Rashid, U.; Mukhtar, H.; Irfan, A.; Anwar, F.; Saari, N. UHPLC-QTOF-MS/MS metabolites profiling and antioxidant/antidiabetic attributes of Cuscuta reflexa grown on Casearia tomentosa: Exploring phytochemicals role via molecular docking MS/MS metabolites profiling and antioxidant/antidiabetic attributes of Cuscuta reflexa grown on Casearia tomentosa: Exploring phytochemicals role via molecular docking. Int. J. Food Prop. 2020, 23, 918–940. [Google Scholar]
  41. Clément, C.; Diazgrados, D.A.; Avula, B.; Khan, I.A.; Mayer, A.C.; Ponce Aguirre, D.D.; Manrique, I.; Kreuzer, M. Influence of colour type and previous cultivation on secondary metabolites in hypocotyls and leaves of maca (Lepidium meyenii Walpers). J. Sci. Food Agric. 2010, 90, 861–869. [Google Scholar] [CrossRef]
  42. Clément, C.; Diaz, D.; Manrique, I.; Manrique, I.; Avula, B.; Khan, I.A.; Ponce Aguirre, D.D.; Kunz, C.; Mayer, A.C.; Kreuzer, M. Secondary Metabolites in Maca as Affected by Hypocotyl Color, Cultivation History, and Site. Agron. J. 2010, 102, 431–439. [Google Scholar] [CrossRef]
  43. Rubio, J.; Dang, H.; Gong, M.; Liu, X.; Chen, S.L.; Gonzales, G.F. Aqueous and hydroalcoholic extracts of Black Maca (Lepidium meyenii) improve scopolamine-induced memory impairment in mice. Food Chem. Toxicol. 2007, 45, 1882–1890. [Google Scholar] [CrossRef] [PubMed]
  44. Rubio, J.; Caldas, M.; Dávila, S.; Gasco, M.; Gonzales, G.F. Effect of three different cultivars of Lepidium meyenii (Maca) on learning and depression in ovariectomized mice. BMC Complement. Altern. Med. 2006, 6, 1–7. [Google Scholar] [CrossRef] [Green Version]
  45. Blazevic, I.; Burcul, F.; Ruscic, M.; Mastelic, J. Glucosinolates, volatile constituents, and acetylcholinesterase inhibitory activity of Alyssoides utriculata. Chem. Nat. Compd. 2013, 49, 374–378. [Google Scholar] [CrossRef]
  46. Szwajgier, D.; Baranowska-Wójcik, E. Terpenes and Phenylpropanoids as Acetyl- and Butyrylcholinesterase Inhibitors: A Comparative Study. Curr. Alzheimer Res. 2019, 16, 963–973. [Google Scholar] [CrossRef] [PubMed]
  47. Studzińska-Sroka, E.; Majchrzak-Celińska, A.; Zalewski, P.; Szwajgier, D.; Baranowska-Wójcik, E.; Żarowski, M.; Plech, T.; Cielecka-Piontek, J. Permeability of Hypogymnia physodes Extract Component—Physodic Acid through the Blood–Brain Barrier as an Important Argument for Its Anticancer and Neuroprotective Activity within the Central Nervous System. Cancers (Basel) 2021, 13, 1717. [Google Scholar] [CrossRef]
  48. Szwajgier, D.; Borowiec, K.; Zapp, J. Activity-guided isolation of cholinesterase inhibitors quercetin, rutin and kaempferol from Prunus persica fruit. Zeitschrift. Naturforsch.-Sect. C J. Biosci. 2020, 75, 87–96. [Google Scholar] [CrossRef]
  49. Hudson, K.L.; Bartlett, G.J.; Diehl, R.C.; Agirre, J.; Gallagher, T.; Kiessling, L.L.; Woolfson, D.N. Carbohydrate-Aromatic Interactions in Proteins. J. Am. Chem. Soc. 2015, 137, 15152–15160. [Google Scholar] [CrossRef] [Green Version]
  50. R Core Team. R: A Language and Environment for Statistical Computing; R Core Team: Vienna, Austria, 2020. [Google Scholar]
  51. RStudio Team. RStudio: Integrated Development Environment for R.; RStudio, Inc.: Boston, MA, USA, 2019. [Google Scholar]
  52. Borchers, H.W. Pracma: Practical Numerical Math Functions 2022, Version 2.3.8 from R-Forge; R Package Vignette: Madison, WI, USA, 2022. [Google Scholar]
  53. Kassambara, A.; Mundt, F. Factoextra: Extract and Visualize the Results of Multivariate Data Analyses. R Package Version 2020, 1, 337–354. [Google Scholar]
  54. Friendly, M.; Fox, J.; Chalmers, P. Matlib: Matrix Functions for Teaching and Learning Linear Algebra and Multivariate Statistics; R Package Vignette: Madison, WI, USA, 2021. [Google Scholar]
  55. Wei, T.; Simko, V. R Package “Corrplot”: Visualization of a Correlation Matrix; R Package Vignette: Madison, WI, USA, 2021. [Google Scholar]
  56. Szwajgier, D.; Baranowska-Wójcik, E.; Kukula-Koch, W.; Kowalik, K.; Polak-Berecka, M.; Waśko, A. Evolution of the anticholinesterase, antioxidant, and anti-inflammatory activity of Epilobium angustifolium L. infusion during in vitro digestion. J. Funct. Foods 2021, 85, 104645. [Google Scholar]
  57. Rhee, I.K.; van Rijn, R.M.; Verpoorte, R. Qualitative determination of false-positive effects in the acetylcholinesterase assay using thin layer chromatography. Phytochem. Anal. 2003, 14, 127–131. [Google Scholar] [CrossRef] [PubMed]
  58. Hanwell, M.D.; Curtis, D.E.; Lonie, D.C.; Vandermeerschd, T.; Zurek, E.; Hutchison, G.R. Avogadro: An advanced semantic chemical editor, visualization, and analysis platform. J. Cheminform. 2012, 4, 17. [Google Scholar] [CrossRef] [PubMed] [Green Version]
  59. Rappé, A.K.; Casewit, C.J.; Colwell, K.S.; Goddard, W.A.; Skiff, W.M. UFF, a Full Periodic Table Force Field for Molecular Mechanics and Molecular Dynamics Simulations. J. Am. Chem. Soc. 1992, 114, 10024–10035. [Google Scholar] [CrossRef]
  60. Moorad, D.R.; Luo, C.; Saxena, A.; Doctor, B.P.; Garcia, G.E. Purification and determination of the amino acid sequence of equine serum butyrylcholinesterase. Toxicol. Methods 1999, 9, 219–227. [Google Scholar] [CrossRef]
  61. Waterhouse, A.; Bertoni, M.; Bienert, S.; Studer, G.; Tauriello, G.; Gumienny, R.; Heer, F.T.; de Beer, T.A.P.; Rempfer, C.; Bordoli, L.; et al. SWISS-MODEL: Homology modelling of protein structures and complexes. Nucleic Acids Res. 2018, 46, W296–W303. [Google Scholar] [CrossRef] [Green Version]
  62. Bienert, S.; Waterhouse, A.; De Beer, T.A.P.; Tauriello, G.; Studer, G.; Bordoli, L.; Schwede. The SWISS-MODEL Repository-new features and functionality. Nucleic Acids Res. 2017, 45, D313–D319. [Google Scholar] [CrossRef] [Green Version]
  63. Trott, O.; Olson, A.J. AutoDock Vina: Improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J. Comput. Chem. 2009, 31, 455–461. [Google Scholar] [CrossRef] [Green Version]
Figure 1. General structure of glucosinolates, the group R is variable.
Figure 1. General structure of glucosinolates, the group R is variable.
Ijms 23 04858 g001
Figure 2. PCA performed for the ABS dataset. Diagonal: histograms of the first three resulting PCs subjected to VARIMAX rotation (henceforth known as dimensions, Dim). Lower triangle: linear maps of samples in the space of the first three dimensions. Upper triangle: linear maps of original variables (compounds) in the space of the first three dimensions. The numbers correspond to the numbering presented in the Section 3.1.
Figure 2. PCA performed for the ABS dataset. Diagonal: histograms of the first three resulting PCs subjected to VARIMAX rotation (henceforth known as dimensions, Dim). Lower triangle: linear maps of samples in the space of the first three dimensions. Upper triangle: linear maps of original variables (compounds) in the space of the first three dimensions. The numbers correspond to the numbering presented in the Section 3.1.
Ijms 23 04858 g002
Figure 3. Cluster dendrogram prepared on the basis of the ABS dataset.
Figure 3. Cluster dendrogram prepared on the basis of the ABS dataset.
Ijms 23 04858 g003
Figure 4. PCA performed for the REL dataset. Diagonal: histograms of the first three resulting PCs subjected to VARIMAX rotation (henceforth known as dimensions, Dim1). Lower triangle: linear maps of samples in the space of the first three dimensions. Upper triangle: linear maps of original variables (compounds) in the space of the first three dimensions.
Figure 4. PCA performed for the REL dataset. Diagonal: histograms of the first three resulting PCs subjected to VARIMAX rotation (henceforth known as dimensions, Dim1). Lower triangle: linear maps of samples in the space of the first three dimensions. Upper triangle: linear maps of original variables (compounds) in the space of the first three dimensions.
Ijms 23 04858 g004
Figure 5. Cluster dendrogram prepared on the basis of the REL dataset.
Figure 5. Cluster dendrogram prepared on the basis of the REL dataset.
Ijms 23 04858 g005
Figure 6. (A) Correlation between the average binding energies determined for X-ray crystal structure (XRD) of BuChE and the molecular volume of studied compounds. (B) Correlation between binding energies obtained for X-ray crystal structure (XRD) of human BuChE and homology model (HM) of equine BuChE. (C) Correlation of binding energies determined for AChE and BuChE and the same set of ligands. Determination coefficients are given as well.
Figure 6. (A) Correlation between the average binding energies determined for X-ray crystal structure (XRD) of BuChE and the molecular volume of studied compounds. (B) Correlation between binding energies obtained for X-ray crystal structure (XRD) of human BuChE and homology model (HM) of equine BuChE. (C) Correlation of binding energies determined for AChE and BuChE and the same set of ligands. Determination coefficients are given as well.
Ijms 23 04858 g006
Figure 7. (A,B) The graphical illustration of the validation results, i.e., the superposition of the most favorable ligand poses found during docking (colored in red) with the position of the co-crystalized ligand, present in the binding cavity of AChE (A) or BuChE (B) (colored by atom type). (C,D) The superposition of the most favorable poses of all ligands interacting with either AChE (C) or BuChE (D). The selected amino acids creating the most essential ligand–protein contacts are given separately for each of ligand moieties. (E,F) The most favorable location of compound glucotropaeolin molecule bound to either the AChE (E) or BuChE (F) structure. The ligand molecule is shown as thick sticks, whereas all the closest amino acid residues (of distances no larger than 0.4 nm) are represented by thin sticks. The description of the interaction types is given in the text.
Figure 7. (A,B) The graphical illustration of the validation results, i.e., the superposition of the most favorable ligand poses found during docking (colored in red) with the position of the co-crystalized ligand, present in the binding cavity of AChE (A) or BuChE (B) (colored by atom type). (C,D) The superposition of the most favorable poses of all ligands interacting with either AChE (C) or BuChE (D). The selected amino acids creating the most essential ligand–protein contacts are given separately for each of ligand moieties. (E,F) The most favorable location of compound glucotropaeolin molecule bound to either the AChE (E) or BuChE (F) structure. The ligand molecule is shown as thick sticks, whereas all the closest amino acid residues (of distances no larger than 0.4 nm) are represented by thin sticks. The description of the interaction types is given in the text.
Ijms 23 04858 g007
Figure 8. An example of distribution of black maca tubers allocated to size classes dependent on physical size and weight in the black phenotype of L. peruvianum synon. L. meyenii).
Figure 8. An example of distribution of black maca tubers allocated to size classes dependent on physical size and weight in the black phenotype of L. peruvianum synon. L. meyenii).
Ijms 23 04858 g008
Table 1. The composition of the 50% EtOH (v/v) extract from tubers of black A variety of Lepidium peruvianum recorded in the negative ionization mode. (RDB—rings and double bond number, delta—error of mass measurement, ion—ionization mode: positive or negative, ND—not detected). Compounds depicted in bold were subjects of chemometric considerations.
Table 1. The composition of the 50% EtOH (v/v) extract from tubers of black A variety of Lepidium peruvianum recorded in the negative ionization mode. (RDB—rings and double bond number, delta—error of mass measurement, ion—ionization mode: positive or negative, ND—not detected). Compounds depicted in bold were subjects of chemometric considerations.
No.Ion. (+/−)Rt (min)Molecular
Formula
m/z
Calculated
m/z
Experimental
Delta (ppm)RDBMS/MSProposed CompoundReferences
1 -1.88C6H12O7195.051195.052−4.971159, 129Gluconic acid[31]
2 -2.06C4H6O5133.0142133.01420.352115Malic acid[31]
3 -2.79C6H8O7191.0197191.0209−6.113155, 111Citric Acid[31]
4 -3.09C4H6O4117.0199117.0193−4.812 Succinic acid[32]
5-5.48C14H19NO10S2424.0378424.0387−2.216274, 182Glucosinalbin (G1)[33]
6-9.26C14H19NO9S2408.0428408.0430−0.386258, 195Glucotropaeolin (G2)[33]
7-12.4C15H21NO10S2438.0534438.0541−1.576196Glucolimnanthin (G3)[33]
8-13.3 C17H22N2O10S2477.0643477.0665−4.588 4-methoxyglucobrassicacin (G4)[33]
9 -15.43C11H18N2O4241.1194241.1199−2.154197Pyroglutamylleucine[34]
10-15.71C22H36ON2S551.1727551.1754−2.796389, 2354-methoxyindolyl-3-hexylhydroxyglucosinolate (G5)[35]
11-15.8C11H21NOS2246.0992246.1006−5.752210, 1301-isothiocyanato-9-methanesulfinylnonane (G6)[36]
12 -16.8C15H29O10NS2446.1154446.11169.872NDIndolyl-5-methylglucosinolate (IMG)[37]
13-17.02C12H24O7N2S2371.0941371.1000−9.92249, 121Pent-4-enylglucosinolate (G7)[33]
14-17.6C13H25NO10S3450.0568450.0570−0.482316Glucoalyssin (G8)[33]
15-17.7C22H38O9N2S2521.1985521.2016−3.675NDindolyl3-hexyl-4-methyl-cyclohexaneglucosinolate (G9)[33]
16 -20.3C11H10O6237.0405237.0406−0.587121Malic acid benzoate[16]
17 -20.35C7H602121.03121.0295−4.075NDBenzoic acid derivative[33]
18 -20.35C24H18O12497.0725497.06975.7216381, 237Fucodiphlorethol[38]
19 -25.6C18H32O5327.2177327.2195−5.493229, 171Trihydroxy-(3)-octadecadienoic acid[39]
20 -26.9C18H34O5329.2333329.2363−8.942211, 171Pinellic acid (9(S),12(S),13(S)-Trihydroxy-10(E)-octadecenoic acid)[40]
Table 2. The results of the quantitative determination of glucotropaeolin in the studied extracts (green color corresponds with the highest quantity of G2 in the extracts, wheras red colour shows the samples with the lowest calculated content).
Table 2. The results of the quantitative determination of glucotropaeolin in the studied extracts (green color corresponds with the highest quantity of G2 in the extracts, wheras red colour shows the samples with the lowest calculated content).
NoSample CodePercentage ContentSDNoSample CodePercentage ContentSD
1BL_A-J1.5140.1720GR_C-J1.5130.10
2BL_B-J0.9540.0821GR_D-J1.3010.12
3BL_C-J0.6320.0522GR_X-J0.5400.05
4BL_D-J0.8750.1023GR_Z-J1.3620.16
5BL_X-J0.2700.0224GR_JUND0.00
6BL_Z-J0.3370.0325PU_A-J0.2140.03
7BL_JUND0.0026PU_B1-J0.0780.01
8CH0.0810.0127PU_B2-J0.2390.02
9BL_gel-AN0.6240.0828PU_C-J1.5720.16
10MIX_gel-ANND0.0029PU_D-J1.1550.09
11BL_Bp-AN0.9650.1130PU_JUND0.00
12RE_Bp-J0.9230.0931YE_A-J0.7490.07
13RE_LA-J0.3440.0232YE_B-J0.8500.09
14RE_LAgel0.0610.0133YE_B20.8500.06
15RE_gel_AN0.1030.0134YE_D-J0.4740.03
16RE_SM0.2420.0235YE_JUND0.00
17GR_A-J0.5180.0436YE_Z-J1.1950.14
18GR_B1-J1.0400.0937YE_B20.8500.06
19GR_B2-J0.8040.06
Table 3. The average results of anti-enzymatic assays for the total extracts and the reference compounds, measured spectrophotometrically in the modified Ellman’s tests.
Table 3. The average results of anti-enzymatic assays for the total extracts and the reference compounds, measured spectrophotometrically in the modified Ellman’s tests.
BuChEAChE
% inhibition by the total extracts
BL_A-J77.224.1
BL_B-J77.222.5
BL_C-J77.222.5
BL_D-J73.719.3
BL_Z-J72.519.3
BL_B-J77.222.5
RE_B-J71.328.9
GR_B-J76.017.0
YE_B-J75.715.3
% inhibition by the reference compounds
GL37.0556.5
GTP36.255.95
GB39.957.6
GS38.956.9
% inhibition by a positive control
Berberine (0.5 mg/mL)33.7-
Table 4. The collection of the binding energies during docking study; the calculation results were averaged over all structures available in the same PDB record; the corresponding standard deviations are given.
Table 4. The collection of the binding energies during docking study; the calculation results were averaged over all structures available in the same PDB record; the corresponding standard deviations are given.
CompoundBinding Energy (kcal/mol)
AChEBuChE XRDBuChE Homology Model
(G2)−9.1−8.5 ± 0.1−8.0
(G3)−9.2−8.2 ± 0.3−7.5
(G1)−9.0−8.5 ± 0.3−7.9
IMG−9.6−9.4 ± 0.3−9.1
(G8)−7.9−7.2 ± 0.2−7.0
(G4)−9.5−9.3 ± 0.4−8.9
(G5)−8.9−9.0 ± 0.2−8.9
(G7)−8.2−7.3 ± 0.1−6.7
(G9)−9.6−10.2 ± 0.5−10.3
G6)−5.8−4.9 ± 0.2−4.9
Table 5. Plant material: coding, description and sampling details prior to delivery to Laboratory Lublin.
Table 5. Plant material: coding, description and sampling details prior to delivery to Laboratory Lublin.
NoSample CodeSample DescriptionSpecimen Collection LocationSampling YearTime of Collection NoSample CodeSample DescriptionSpecimen Collection LocationSampling YearTime of Collection
1BL_A-JDried black Large tubers size A (Smaller than X)Junín (J-19) Maca plantation, Peruvian Andean highlands * 8 October 2019Between 11am and 2pm local time in Peru20GR_C-JDried grey tubers size C (smaller than B)Junín (J-19)8 October 2019Between 10am and 2 pm local Peruvian time
2BL_B -JDried black tubers size B (Smaller than A)21GR_D-JDried grey tubers size D (smaller than C)
3BL_C-JDried black tubers size C (Smaller than B)22GR_X-JDried grey tubers size X- Extra Large (larger than A)
4BL_D-JDried black tubers size D (Smaller than C)23GR_Z-JDried grey tubers size E (very small–smaller than D)
5BL_X-JDried black tubers size Extra Large (X)24GR_JUSqueezed juice from fresh grey maca tubers (fresh size A)
6BL_Z-JDried black tubers size E (smaller than D)25PU_A-JDried PURPLE maca tubers size A
7BL_JUPressed juice from fresh black tubers (BL A-J and B-J) 6PU_B-JDried PURPLE maca tubers size B (smaller than A)
8CHSemi-dried dark violet tubers purchased by maca contractor from unknown sources/location origin not knownTarma (J-17) *30 September 20173 p.m.
Local time
27PU_B2-JDried PURPLE maca tubers size B-1-Replicate (smaller than A)
9BL_gel-ANOrganic gelatinized BLACK maca tubers, sampled and dried in Ancash and processed in LimaAncash (A-17)
Maca Plantation, Cordillera Blanca, Peruvian Andes.
Mid-October 2017Between 10 a.m. and 4 p.m. local Peruvian time28PU_C-JDried purple maca tubers size C (smaller than B)
10MIX_gel-ANOrganic gelatinized tubers of MIXED maca phenotypes 29PU_D-JDried purple tubers size D (smaller than C)
11BL_Bp-AN Organic BLACK tubersopen- air dried on site and pulverised in Lima processing plant30PU_X-JDried purple tubers size X Extra Large (larger than A)
12RE_Bp-JRED tubers size B dried and powdered-as per above Junín (J-17) 31PU_JUSqueezed juice from fresh purple tubers
13RE_LA-JRED maca large powdered tubers size A32YE_A-JDried yellow tubers size A
14RE_LA gelRED maca large tubers size A—gelatinised powder33YE_B-JDried yellow tubers size B (smaller than A)
15RE_gel ANorganic gelatinized RED maca powder Maca Plantation Cordillera Blanca, Peruvian AndesAncash. (A-17) Cordillera Blanca Peruvian Andes 34YE_B2-JDried yellow tubers size B2—Replicate B (smaller than A)
16RE_SM-J RED maca small tubers –powderJunín (J-17)35YE_D-JDried yellow tubers size D (smaller than C)
17GR_A-JDried grey tubers size AJunín (J-19)8 October 201936YE_JUSqueezed juice from fresh yellow tubers size A
18GR_B1-JDried grey tubers size B (smaller than A) 37YE_Z-JDried yellow tubers size E (smaller than size D) Extra small
19GR_B2-JDried grey tubers size B (smaller than A)
* Tubers from commercial drying yard in Tarma (Junín) on 30 September 2017.
Publisher’s Note: MDPI stays neutral with regard to jurisdictional claims in published maps and institutional affiliations.

Share and Cite

MDPI and ACS Style

Tarabasz, D.; Szczeblewski, P.; Laskowski, T.; Płaziński, W.; Baranowska-Wójcik, E.; Szwajgier, D.; Kukula-Koch, W.; Meissner, H.O. The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors—In Silico and In Vitro Studies. Int. J. Mol. Sci. 2022, 23, 4858. https://doi.org/10.3390/ijms23094858

AMA Style

Tarabasz D, Szczeblewski P, Laskowski T, Płaziński W, Baranowska-Wójcik E, Szwajgier D, Kukula-Koch W, Meissner HO. The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors—In Silico and In Vitro Studies. International Journal of Molecular Sciences. 2022; 23(9):4858. https://doi.org/10.3390/ijms23094858

Chicago/Turabian Style

Tarabasz, Dominik, Paweł Szczeblewski, Tomasz Laskowski, Wojciech Płaziński, Ewa Baranowska-Wójcik, Dominik Szwajgier, Wirginia Kukula-Koch, and Henry O. Meissner. 2022. "The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors—In Silico and In Vitro Studies" International Journal of Molecular Sciences 23, no. 9: 4858. https://doi.org/10.3390/ijms23094858

APA Style

Tarabasz, D., Szczeblewski, P., Laskowski, T., Płaziński, W., Baranowska-Wójcik, E., Szwajgier, D., Kukula-Koch, W., & Meissner, H. O. (2022). The Distribution of Glucosinolates in Different Phenotypes of Lepidium peruvianum and Their Role as Acetyl- and Butyrylcholinesterase Inhibitors—In Silico and In Vitro Studies. International Journal of Molecular Sciences, 23(9), 4858. https://doi.org/10.3390/ijms23094858

Note that from the first issue of 2016, this journal uses article numbers instead of page numbers. See further details here.

Article Metrics

Back to TopTop