In vitro and in silico evaluation of Ononis isoflavonoids as molecules targeting the central nervous system

Abstract

Isoflavonoids with various structural elements show a promising potential effect on central nervous system activities. Despite their favorable medicinal properties, the pharmacokinetic characteristics of this thoroughly investigated group of natural phenolics have only been described to a limited extent. Regarding the lack of information about the BBB permeability of isoflavones, isoflavanones, and pterocarpans found in Ononis species, the aim of our study was to investigate their physico-chemical properties influencing their absorption and distribution. Furthermore, we aimed to characterize the possible MAO-B inhibiting features of Ononis isoflavonoids in silico. Octanol-water partitioning and BBB-PAMPA permeability of formononetin, calycosin D, onogenin, sativanone, medicarpin and maackiain were assessed for the first time in our study. The log P values ranged from 2.21 to 3.03 and log D7.4 values from 2.48 to 3.03, respectively, indicating optimal polarity for BBB permeation. The results of PAMPA-BBB expressed as log Pe values fell between -5.60 and -4.45, predicting their good permeation capability as well. The effective permeability values showed structure-dependent differences, indicating that the pterocarpan type skeleton was the most preferred type, followed by isoflavanones, then isoflavones. The methoxy or methylenedioxy substitution of the same skeleton did not influence the permeability significantly, contrary to an additional hydroxyl group. Membrane retention showed a similar structure dependent pattern to that of effective permeability, ranging from 16% to 70%. For the identification of volumes of chemical space related to particular biological activities the ChemGPS-NP framework was used. The MAO-B inhibitory potency and selectivity were also predicted and validated. Based on our results, MAO-B inhibitory potency could be predicted with good precision, but in the case of selectivity, only the direction could be concluded (favors MAO-B or MAO-A), not the magnitude. Our finding reflects that Ononis isoflavonoid aglycones show an excellent fit with the suggested parameters for BBB permeability and this is the first study to confirm the highly favorable position of these natural products for MAO-B inhibition.

About the Authors:

Nóra Gampe

Roles Conceptualization, Funding acquisition, Investigation, Methodology, Writing – original draft

Affiliation: Department of Pharmacognosy, Semmelweis University, Budapest, Hungary

ORCID logo https://orcid.org/0000-0001-7208-9372

Dominika Noémi Dávid

Roles Investigation

Affiliation: Department of Pharmacognosy, Semmelweis University, Budapest, Hungary

ORCID logo https://orcid.org/0000-0002-9305-1581

Krisztina Takács-Novák

Roles Methodology, Supervision, Writing – review & editing

Affiliation: Department of Pharmaceutical Chemistry, Semmelweis University, Budapest, Hungary

Anders Backlund

Roles Software, Supervision, Validation, Writing – review & editing

Affiliation: Department of Pharmaceutical Biosciences, Pharmacognosy, Uppsala University, Uppsala, Sweden

ORCID logo https://orcid.org/0000-0003-2043-4183

Szabolcs Béni

Roles Conceptualization, Funding acquisition, Writing – review & editing

* E-mail: [email protected]

Affiliation: Department of Pharmacognosy, Semmelweis University, Budapest, Hungary

ORCID logo https://orcid.org/0000-0001-7056-6825

Introduction

Isoflavonoids are biologically relevant plant metabolites. The most studied group of isoflavonoids is isoflavones, the main phytochemical metabolites of e.g. soy, red clover, and kudzu root. The intensity and volume of research focusing on these plants and the isoflavones they contain is understandable, as their consumption as nutrient or dietary supplements are orders of magnitudes higher, than those of any other isoflavonoid containing plants. However, in the family Fabaceae, a wide range of herbs can be found with more exotic isoflavonoid patterns [1]. For example, in the genus Ononis (restharrow), beside isoflavones and isoflavanones, pterocarpans are also present [2–6]. The isoflavones of soy and red clover (genistein, daidzein, and biochanin A) are mainly regarded as phytoestrogenic compounds [7], but other isoflavonoids with various structural elements show diverse biological effects, and they have the potential to act on the central nervous system (CNS). For example, formononetin proved to be neuroprotective in traumatic brain injuries in vivo [8, 9], by protecting cells from oxidative stress [10, 11] and neuroinflammation [12, 13]. Formononetin also showed neuroprotective effect both in vitro [12, 13] and in vivo [14] models of Alzheimer’s disease. Other Ononis isoflavonoids, such as maackiain and calycosin, showed a selective MAO-B inhibiting effect in vitro, which was a magnitude higher than that of genistein. Lee et al. found, that maackiain is a new potent (IC50 = 0.68 μM), selective (selectivity index = 126.2), and reversible MAO-B inhibitor [15]. Furthermore, Oh and co-authors found similar activity for calycosin (IC50 = 0.24 μM, selectivity index = 293.8) and medicarpin (IC50 = 0.30 μM) with even more promising results [16].

Despite the promising medicinal properties of isoflavonoids, their pharmacokinetic characteristics have been described only to a limited extent, focusing mainly on the isoflavones of soy [17–20]. Generally speaking of flavonoids and isoflavonoids, it can be concluded that the absorption rates of the glycosidic forms of these compounds are very, but they are metabolized to their aglycone forms in the alimentary tract by gut microbiome [21–23]. Information is available about the oral bioavailability of medicarpin and formononetin, indicating their fast metabolization after absorption, and the low bioavailability of their aglycone form [24, 25]. On the other hand, in vivo studies after per os administration showed positive results [14, 26, 27]. The blood-brain barrier (BBB) permeability had previously been only investigated in the case of isoflavones only. These results showed, that isoflavone aglycones can permeate the BBB and that they are not P-gp substrates [20, 28, 29].

Regarding the lack of information about the BBB permeability of isoflavones, isoflavanones, and pterocarpans found in Ononis species, the aim of our study was to isolate these compounds in their aglycone form and investigate their physico-chemical properties influencing the absorption and distribution. In order to estimate their BBB permeability, octanol-water partitioning and the parallel artificial membrane permeability assay (PAMPA) methods were chosen, based on the known absorption characteristics of isoflavonoids—which is mainly passive, transcellular diffusion. The PAMPA was firstly introduced by Kansy et al. [30], as a cost-effective and high-throughput approach to predict the oral absorption of drugs. Later, Di and co-authors developed a PAMPA model which is suitable for the estimation of BBB permeability [31]. More recently, other research groups further optimized the assay to reach different requirements, e.g. faster incubation time or the evaluation of natural products [32–34] and applied it to investigate various natural products [35–43]. The PAMPA-BBB can be used as a screening method for plant extracts, highlighting the permeable compounds, which worth isolation [38, 39, 43], or it can be used for the assessment of BBB permeability of pure, isolated natural products and their semi-synthetic derivatives [35–37, 40, 42].

Ononis species contains isoflavones (calycosin D, pseudobaptigenin), isoflavanones (onogenin, sativanone), and pterocarpan (medicarpin) which are structurally similar to formononetin, calycosin, and maackiain, compounds with a probable CNS activity. The observed structural similarities raised the questions if the compounds from Ononis would also be active in these biological tests. To address this question, a charting of biologically relevant physico-chemical features using the ChemGPS-NP framework was employed. This framework has previously proven useful in identifying and defining volumes of chemical space related to particular biological activities [44, 45]. ChemGPS-NP also has the capacity of providing a reference system which allows the characterization and comparison of molecules of natural origin as well as those in routinely prescribed medicines [46, 47].

The aim of the study is assessing the physico-chemical properties of Ononis isoflavonoids influencing their permeability through the BBB for the first time by the means of octanol-water partitioning and BBB-PAMPA. Moreover, in silico experiments were carried out in order to investigate the permeability, the MAO-B inhibitor potency and selectivity of the target molecules.

Materials and methods

Solvents and chemicals

Daidzein, quinine, caffeine, naringenin, salicylic acid, and rutin standards were obtained from Sigma-Aldrich. HPLC grade acetonitrile and methanol were purchased from Merck. Ethyl-acetate, formic acid, methanol, and acetone of reagent grade were purchased from Reanal-Ker. HPLC grade water was prepared with a Millipore Direct Q5 water purification system. All aqueous eluents for HPLC were filtered through MF-Millipore membrane filters (0.45 μm, mixed cellulose esters).

Ononis isoflavonoid aglycones

The isolation and identification of isoflavonoid aglycones were carried out based on the methods described in our previous papers [4, 5, 48]. The quantity and purity of the isolated substances were as follows: maackiain (12.9 mg, 100%), medicarpin (13.4mg, 100%), onogenin (12.3mg, 95%), sativanone (23.1 mg, 98%), formononetin (7.3mg, 94%), pseudobaptigenin (14.8 mg, 93%), calycosin D (4.4 mg, 92%). The purity of the isolated compounds was evaluated by UPLC and Max Plot chromatograms, see Supplementary for chromatograms. Calycosin was obtained from Sigma-Aldrich, as it could not be isolated in appropriate quantity and purity for further testing.

Determination of log P values

The logarithm of octanol/water partition coefficient was measured by the validated protocol of traditional shake-flask method [49]. Both the true partition coefficient (log P) value and the distribution coefficient (log DpH) at pH 7.4 of the samples were determined. Britton-Robinson buffer (BR) served as the aqueous phase, while n-octanol as the organic phase. The two phases were mutually saturated with each-other. The samples were obtained in methanol. From this concentrated solution the appropriate dilution was made with the aqueous phase. Different phase ratios (R) were used for equilibration (10ml BR buffer: 0.1 mL or 0.2 mL or 0.5 mL octanol) at 25°C for 1 hour. The phases were separated by centrifugation. The decrease in concentration was determined in the aqueous phase by spectroscopy measuring the absorbance before (Ao) and after (A1) the partitioning. The log P and log DpH value was calculated using the following equation:

Four parallel experiments were conducted, and average value and standard deviations (SD) were calculated. Predicted values of log P and pKa were calculated using ChemAxon Marvin.

Parallel artificial membrane permeability assay

The applied PAMPA-BBB setup was based on the work of Könczöl et al. [32] who validated the method for the evaluation of natural products. PAMPA “sandwiches” were formed from a Stirwell (Pion Inc.) 96-well donor and acceptor plates with a polyvinylidene difluoride filter bottom. The wells were coated with 5 μL of the solution of 16.0 mg polar porcine brain lipid (Avanti Polar Lipids Inc.) and 8.0 mg cholesterol (Sigma-Aldrich) in 600 μL of dodecane (Acros Organics). The concentrations of standard compounds and isoflavonoids of initial donor samples were approximately 10 mM in DMSO, then these samples were diluted a hundredfold with phosphate buffer saline pH 7.4 to obtain donor solutions with 1% DMSO content of which 150 μL was applied. The acceptor compartment was filled with 180 μL phosphate buffer saline pH 7.4. The sandwiches were incubated in a water vapor-saturated atmosphere at 37°C for 4 hours in the Gut-Box (Pion Inc.) module with stirring to adjust the thickness of the aqueous boundary layer to 60 μm. Three to six parallel measurements were made for each sample and the assays were repeated three consecutive days. Sample concentrations in the acceptor and donor wells were determined by HPLC [33]. For this, a Waters Acquity UPLC system was furnished with a sample manager, a binary solvent manager, and a PDA detector (Waters Corporation). The samples were subsequently analyzed on the Acquity UPLC with a DEH C18 column (2.1 x 100 mm, 1.7 μm) using eluents 0.1% formic acid (A) and acetonitrile (B) at 0.3 mL/min flowrate. To determine the standard compounds, the following gradient program was used: 0 min 5% B, 5 min 100% B.

For the quantification of isoflavonoids an isocratic system was used with 40% B. Retention (capacity) factors were calculated from the chromatographic runs detailed above, as follows:[1]

The effective permeability and the membrane retention of the compounds were calculated by the following equations [50]:[2]where Pe is the effective permeability coefficient (cm/s), A is the filter area (0.24 cm2), VD and VA are the volumes in the donor (0.15 cm3) and acceptor phases (0.18 cm3), t is the incubation time (s), τSS is the time (s) to reach steady-state (240 s), CA(t) is the concentration (mol/cm3) of the compound in the donor phase at time t, CD(0) is the concentration (mol/cm3) of the compound in the donor phase at time 0, CD(t) is the concentration (mol/cm3) of the compound in the donor phase at time t, and MR is the estimated membrane retention factor.

[3]

Prior to analyzing isoflavonoid aglycones, standard compounds were applied (daidzein, quinine, caffeine, naringenin, salicylic acid, and rutin) to verify the system suitability (see Supplementary material).

Cheminformatics

Chemical space analysis was performed using the principal component analysis (PCA)-based chemical space navigation tool ChemGPS-NP [45, 47] which is freely available online at http://www. chemgps.bmc.uu.se/. The reference sets of BBB passive diffusers, non-diffusers and P-glycoprotein substrates were based on the thesis of L. Viklund [51] titled ‘ChemGPS-NP as a tool for predicting drug distribution across the blood-brain barrier’. Monoamine oxidase A and B inhibitors were collected from the ChEMBl database searching for the human enzymes and using the cut-off value of IC50 < 10 μM. For evaluating activity, all MAO-B inhibitors were used over the cutoff value. For evaluating the selectivity, only molecules with a known MAO-A and B inhibitor activity were used. Selectivity index was calculated based on the IC50 values of the reference structures (IC50 MAO-B/IC50 MAO-A). Euclidean distances were calculated between points P = (p1, p2,…, p8) and Q = (q1, q2,…, q8) in Euclidean 8D space provided by the ChemGPS-NP coordinates using the following equation in Excel:[4]

The first three dimensions (plotted in Plotly Chart studio) of the ChemGPS-NP map, can be interpreted in such a way that the first dimension (principal component one, PC1) represents size, shape and polarizability; PC2 corresponds to aromatic and conjugation-related properties; PC3 describes lipophilicity, polarity, and H-bond capacity. For the purpose of visualization, only the first three dimensions were used but the calculations were based on all eight dimensions.

Results and discussion

Lipophilicity characterized by the log DpH / log P values

The lipophilicity of drugs is an essential feature determining their route in living systems, including their absorption, distribution, accumulation, and elimination. The octanol/water partition coefficient is a generally accepted physico-chemical parameter for the characterization of lipophilicity. The preliminary pKa predictions suggested that methoxylated isoflavones have a pKa value around 7, which was confirmed by literature sources [52, 53]. Because of the methylation in the 4’-OH position, the 7-OH group attached to the chromenone moiety has a pKa value lower than expected (Fig 1).

[Figure omitted. See PDF.]

Fig 1. The structure and numbering of the investigated isoflavonoid aglycones.

https://doi.org/10.1371/journal.pone.0265639.g001

Therefore, the log P value of isoflavone derivatives (formononetin and calycosin) was determined at pH 5. In the case of isoflavanones, the calculations predicted a pKa of 7.78 which could not be explained by the extensive conjugated electron system of the chromenone ring as the bond between C2 and C3 is saturated (see Fig 1). On the other hand, these molecules have the ability of oxo-enol tautomerism [54] (Fig 2), which could amplify the acidic character of the 7-OH group, so their log P value was determined based on the partition at pH 5.

[Figure omitted. See PDF.]

Fig 2.

The oxo-enol tautomerism and pH dependent UV spectra of sativanone (A) and onogenin (B).

https://doi.org/10.1371/journal.pone.0265639.g002

Fig 2A and 2B show the pH dependent UV spectra of these two molecules. Based on these data we concluded that the two molecules are fully unionized at pH 5. Regarding pterocarpans, ring A, which bears the ionizable 3-OH group (Fig 1), is not linked to other parts of the molecule with conjugated bonds, and thus the calculated pKa value is close to that of a phenolic hydroxyl group.

In our results (Table 1), all compounds showed favorable log P and log D7.4 values regarding drug-likeness and Lipinski’s rule of five, indicating, that the studied compounds are able to passively permeate through biological membranes. Those structures which have lower pKa values (calycosin D, formononetin, onogenin and sativanone) showed somewhat higher log DpH values at pH 5 than at pH 7.4. The difference between the log DpH values was higher in the case of isoflavanones, whereas for formononetin, the difference is measurable but not significant. The highest log D7.4 value was measured for medicarpin, which were followed by sativanone and formononetin, then onogenin and maackiain. The least lipophilic compound is calycosin D, with an order of magnitude lower lipophilicity than that of medicarpin (Table 1).

[Figure omitted. See PDF.]

Table 1. The calculated pKa values, experimental log DpH and log P values and wavelength used during the measurements of the six isoflavonoids.

https://doi.org/10.1371/journal.pone.0265639.t001

In vitro PAMPA-BBB permeability

In this study, a PAMPA-BBB method, previously validated for natural compounds was utilized [32] with a slight modification using stirring in wells with an unstirred water layer of 60 μm. In the first step, we used standard molecules (daidzein, quinine, caffeine, naringenin, salicylic acid, and rutin) for system validation. In our results, the measured log Pe values showed good correlation (R2 = 0.9498) with the log BB values from literature (see S1 Table and S1 Fig in S1 File), verifying the appropriate predictive power of the system.

Molecules with very high or very low log BB values on the scale could be easily defined as well permeating (BBB+) or not sufficiently permeating (BBB-). As log BB and CNS permeability are continuous data, it is hard to form only two categories, and to characterize a single value which demarcates BBB+ and BBB- molecules. In the work of Mensch et al. the condition for BBB+ category is to have a log BB value ≥ 0 and they determined the in vitro limit differently, based on the specificities of the applied PAMPA-BBB test [34]. Könczöl et al. regarded molecules as BBB+ above a log BB -0.5 value and defined the arbitrary unit of -6.0 of log Pe as an effective discriminator between BBB+ and BBB- structures. Di et al. determined an uncertain, intermediate zone between Pe = 2 x 10−6 cm/s and Pe = 4 x 10−6 cm/s (log Pe -5.69 to -5.39) [31, 32].

After analyzing the experimental log Pe values (S2 Table in S1 File), structural correlation can be observed. The isoflavones calycosin D and formononetin had the lowest log Pe values, followed by isoflavanones sativanone and onogenin, while the highest values were measured for pterocarpans maackiain and medicarpin. The methoxy or methylenedioxy substitution of the same skeleton caused no significant differences in log Pe values (onogenin vs. sativanone, medicarpin vs. maackiain), but an additional hydroxyl group contributed significantly (formononetin vs. calycosin D) (Figs 1 and 3A). Regarding the different literature cut off values in distinction of BBB+ and BBB- structures found in literature, all compounds, except calycosin D, fell into the well-permeating category, therefore they can be considered permeating through the blood-brain barrier via passive diffusion. Calycosin D could be considered as well-permeating based on the works of Könczöl et al. However, according to Mensch et al., it would be categorized as BBB-, while using the categorization of Di et al. it could be classified as uncertain with a Pe value of 2.5 x 10−6 cm/s [31, 32, 34]. Taking into consideration, that the structural isomer calycosin, which only differs slightly in the position of a methyl group, has been reported by several studies as a neuroprotective agent [55–57], it can be assumed that calycosin D can reach the cells of the brain. Comparing the permeability of structures studied, more than one order of magnitude difference could be drawn between the most and least permeable compounds (medicarpin and calycosin D, respectively). The membrane retention showed a somewhat similar pattern, with an average value of 20% for isoflavones, 33% for isoflavanones and 61% for pterocarpans, but the distinction between isoflavonoid skeletons is less pronounced (Fig 3B). Medicarpin showed the highest membrane retention, which can be correlated with its highest log D7.4 value (see Table 1).

[Figure omitted. See PDF.]

Fig 3.

Experimental log Pe values (A) and membrane retention (B) of the six isoflavonoid compounds (bright green: isoflavones, deep green: isoflavanones, bright blue: pterocarpans).

https://doi.org/10.1371/journal.pone.0265639.g003

Table 2 highlights an almost perfect fit between the isoflavonoids and the suggested BBB permeability [58]. Only calycosin D shows more than one outlier, and the polar surface area of onogenin is higher than recommended. Based on these results, the physico-chemical parameters of Ononis isoflavonoids are appropriate for further drug development aiming the CNS.

[Figure omitted. See PDF.]

Table 2. Suggested physicochemical property ranges for increasing the potential for BBB penetration and physicochemical features of the six isoflavonoid aglycones.

https://doi.org/10.1371/journal.pone.0265639.t002

Cheminformatics

Isoflavonoid BBB permeability.

In order to corroborate the in vitro results of isoflavonoids permeating the BBB, we investigated the location of the structures studied in the chemical space in the proximity of known passive diffusers and non-diffusers. Based on the plot, the isoflavonoids merge to the cluster of BBB passive diffusers and do not have common chemical space with non-diffuser molecules, confirming the transcellular passive diffusing mechanism of BBB permeability and the accuracy of the PAMPA studies (Fig 4). Most submerged in the cloud of passive diffusers are onogenin and sativanone can be observed, with relatively the highest PC1 and lowest PC2 values, which can be the result of the saturated ring C (lower aromaticity) together with the high number of O atoms (higher polarizability) (see Fig 1). In Fig 4 it can be seen that the substrates of P-gp are not clustered together with Ononis isoflavonoids, but rather with non-diffuser molecules.

[Figure omitted. See PDF.]

Fig 4. The location of BBB passive diffusers (blue), non-diffusers (red), P-glycoprotein substrates (orange) and isoflavonoids (black) in chemical space.

https://doi.org/10.1371/journal.pone.0265639.g004

Isoflavonoids and monoamine oxidases.

MAO inhibitors are mainly used in psychiatry for the treatment of depressive disorders, and anxiety disorders and in neurology for the treatment of Parkinson’s disease and Alzheimer’s disease. While the classical non-selective and irreversible MAO inhibitors are characterized by the risk of inducing a hypertensive crisis, the selective MAO-B inhibitor selegiline and the selective and reversible inhibitor of MAO-A, moclobemide, are free from this potential interaction [59]. This emphasizes the significance of the previous results, where the selected isoflavonoids seemed to be selective and reversible inhibitors of MAO-B [15, 16]. Having evaluated the first dataset describing the position of MAO-B inhibitors with IC50<1000 nM it can be observed that Ononis isoflavonoids can be found in the same cluster as MAO-B inhibitors (Fig 5).

[Figure omitted. See PDF.]

Fig 5. The location of MAO-B inhibitors and Ononis isoflavonoids in chemical space (blue color correlates with lower IC50 values, yellow with higher ones).

https://doi.org/10.1371/journal.pone.0265639.g005

As compounds with lower than 1 EDs are considered to be close in the chemical space [44], firstly their number were assessed around the target compounds. Calycosin and calycosin D could be regarded as identical in this system, as their ED was 0. Formononetin showed the highest number of MAO-B inhibitors in its proximity, whereas maackiain had no close neighbors. In the next step, the IC50 values of the closest 10 and 5 molecules were investigated and compared with in vitro results from literature [15, 16]. Surprisingly, maackiain with no other compounds in tight closeness, had a group of surrounding molecules with one of the most potent characteristics. The estimated IC50 values (0.41μM and 0.51 μM) fell very close to its in vitro value (0.68 μM) [15]. In the case of medicarpin, the estimation (0.60 μM and 0.21 μM) proved to be appropriate (in vitro 0.30 μM) [16]. Calycosin and calycosin D showed the highest estimated average IC50 value (1.9 and 2.1 μM) which does not correlate with the experimental data (0.24 μM). On the other hand, the compounds around calycosin showed a high deviancy in their IC50 values (45 nM– 8710 nM). Additionally, calycosin had the lowest Euclidean distance (0.06) to a standard molecule, acacetin, which has an outstandingly low IC50 value (49 nM). Other isoflavonoid compounds which had not been tested before (formononetin, onogenin, sativanone, pseudobaptigenin), showed very promising estimated results, too (see Table 3).

[Figure omitted. See PDF.]

Table 3. The number of compounds in the proximity of the target molecules, and the average MAO-B inhibitory IC50 values of the 10 and 5 closest structures in chemical space (see details in supplementary).

https://doi.org/10.1371/journal.pone.0265639.t003

Comparing the position of isoflavonoids and reference compounds with increased activity, they showed a higher value in PC3 indicating more lipophile compounds. This trend can be observed regarding the experimental log P values and the estimated activity of our molecules, too. As a consequence, increasing the apolar character of the tested molecules may increase the enzyme inhibitory activity.

In the next step, molecules from the first dataset were selected with known MAO-A inhibitory effect, then the selectivity index was calculated and plotted in chemical space (see Fig 6). Purple color shows a higher selectivity for MAO-A, while yellow is the indicator of selective MAO-B inhibition.

[Figure omitted. See PDF.]

Fig 6. The location of MAO-A and MAO-B inhibitors beside Ononis isoflavonoids in chemical space color coded by their selectivity index (purple: Inhibiting selectively MAO-A, yellow: Inhibiting selectively MAO-B).

https://doi.org/10.1371/journal.pone.0265639.g006

The same calculations were carried out, as previously, the results can be found in Table 4. The number of compounds with ED<1 is lower for each molecule, since fewer compounds were tested. Altogether, all isoflavonoids showed an estimated preference for the MAO-B enzyme, but the predicted selectivity indices are an order of magnitude lower, than the experimental ones (for maackiain, calycosin and medicarpin). These results both emphasize the strength and the limitation of the ChemGPS-NP framework, as the trends of the selectivity could be predicted, however, the exact values not. Onogenin and sativanone seems to be promising candidates for further testing based on their estimated values. More selective molecules have higher PC3 values, indicating the possible beneficial effect of increasing the lipophilicity of the tested molecules again.

[Figure omitted. See PDF.]

Table 4. The number of compounds in the proximity of the target molecules, and the average MAO-B inhibitory selectivity indices of the 10 and 5 closest structures in chemical space (see details in supplementary).

https://doi.org/10.1371/journal.pone.0265639.t004

Conclusion

In this work, the physico-chemical properties such as lipophilicity (expressed with the logarithm of octanol-water partition coefficient) and permeability (measured by BBB-PAMPA assay) of formononetin, calycosin D, onogenin, sativanone, medicarpin and maackiain were characterized for the first time. As a result, optimal log P and log D7.4 values were found for passive diffusion through the BBB. Before investigating isoflavonoid aglycones, the PAMPA-BBB system was verified using standard molecules. The obtained effective permeability values of Ononis isoflavonoids indicated excellent permeability. Structure—log Pe relationships were deduced, as major differences were found for the various aglycones, but their methoxylation or methylenedioxy-substitution showed insignificant differences. Although pterocarpan derivatives showed the most preferable log Pe values, their membrane retention was the highest, too. Plotting the target isoflavonoids in the chemical space together with passive diffusers, P-gp substrates and non-diffusers using the ChemGPS-NP framework confirmed the hypothesized transcellular passive diffusion as a route of absorption and distribution. Investigating the closest known MAO-B inhibitors in chemical space, the predicted IC50 values showed good correlation with in vitro values. When comparing the selectivity of the predicted values with the experimental ones, only the preference could be estimated, but values could not. Considering our results, formononetin, onogenin, pseudobaptigenin and sativanone are predicted inhibitors of MAO-B, making them good candidates for future in vitro and in vivo tests.

Supporting information

S1 File.

https://doi.org/10.1371/journal.pone.0265639.s001

(DOCX)

S2 File.

https://doi.org/10.1371/journal.pone.0265639.s002

(XLSX)

S3 File.

https://doi.org/10.1371/journal.pone.0265639.s003

(XLSX)

S4 File.

https://doi.org/10.1371/journal.pone.0265639.s004

(ZIP)

Acknowledgments

N. G. gratefully acknowledges the help of Ágnes Teitelbaum in carrying out the partition experiments.

Citation: Gampe N, Dávid DN, Takács-Novák K, Backlund A, Béni S (2022) In vitro and in silico evaluation of Ononis isoflavonoids as molecules targeting the central nervous system. PLoS ONE 17(3): e0265639. https://doi.org/10.1371/journal.pone.0265639

References

1. Veitch NC. Isoflavonoids of the leguminosae. Nat Prod Rep. 2013;30: 988–1027. pmid:23736284

2. Abdel-Kader MS. Preliminary pharmacological study of the pterocarpans macckian and trifolirhizin isolated from the roots of Ononis vaginalis. Pak J Pharm Sci. 2010;23: 182–187. pmid:20363697

3. Abdel-Kader MS. Two isoflavonoid glucoside derivatives from Ononis serrata growing in Egypt. Nat Prod Sci. 2004;10: 321–324.

4. Gampe N, Darcsi A, Lohner S, Béni S, Kursinszki L. Characterization and identification of isoflavonoid glycosides in the root of Spiny restharrow (Ononis spinosa L.) by HPLC-QTOF-MS, HPLC-MS/MS and NMR. J Pharm Biomed Anal. 2016;123: 74–81. pmid:26874257

5. Gampe N, Darcsi A, Nagyné Nedves A, Boldizsár I, Kursinszki L, Béni S. Phytochemical analysis of Ononis arvensis L. by liquid chromatography coupled with mass spectrometry. J Mass Spectrom. 2018; 121–133. pmid:30408845

6. Ghribi L, Waffo-Téguo P, Cluzet S, Marchal A, Marques J, Mérillon J-MM, et al. Isolation and structure elucidation of bioactive compounds from the roots of the Tunisian Ononis angustissima L. Bioorg Med Chem Lett. 2015;25: 3825–3830. pmid:26248805

7. Dixon RA, Pasinetti GM. Flavonoids and isoflavonoids: from plant biology to agriculture and neuroscience. Plant Physiol. 2010;154: 453–7. pmid:20921162

8. Li Z, Zeng G, Zheng X, Wang W, Ling Y, Tang H, et al. Neuroprotective effect of formononetin against TBI in rats via suppressing inflammatory reaction in cortical neurons. Biomed Pharmacother. 2018;106: 349–354. pmid:29966980

9. Li Z, Wang Y, Zeng G, Zheng X, Wang W, Ling Y, et al. Increased miR-155 and heme oxygenase-1 expression is involved in the protective effects of formononetin in traumatic brain injury in rats. Am J Transl Res. 2017;9: 5653–5661. pmid:29312517

10. Jia W-C, Liu G, Zhang C-D, Zhang S-P. Formononetin attenuates hydrogen peroxide (H2O2)-induced apoptosis and NF-kappaB activation in RGC-5 cells. Eur Rev Med Pharmacol Sci. 2014;18: 2191–2197. pmid:25070826

11. Wu J, Ke X, Ma N, Wang W, Fu W, Zhang H, et al. Formononetin, an active compound of Astragalus membranaceus (Fisch) Bunge, inhibits hypoxia-induced retinal neovascularization via the HIF-1α/VEGF signaling pathway. Drug Des Devel Ther. 2016;Volume 10: 3071–3081. pmid:27729769

12. Chen L, Ou S, Zhou L, Tang H, Xu J, Guo K. Formononetin attenuates Aβ25-35-induced cytotoxicity in HT22 cells via PI3K/Akt signaling and non-amyloidogenic cleavage of APP. Neurosci Lett. 2017;639: 36–42. pmid:28034780

13. Sun M, Zhou T, Zhou L, Chen Q, Yu Y, Yang H, et al. Formononetin protects neurons against hypoxia-induced cytotoxicity through upregulation of ADAM10 and sAβPPα. J Alzheimer’s Dis. 2012;28: 795–808. pmid:22085862

14. Fei H-X, Zhang Y-B, Liu T, Zhang X-J, Wu S-L. Neuroprotective effect of formononetin in ameliorating learning and memory impairment in mouse model of Alzheimer’s disease. Biosci Biotechnol Biochem. 2018;82: 57–64. pmid:29191087

15. Lee HW, Ryu HW, Kang M-GG, Park D, Oh S-RR, Kim H. Potent selective monoamine oxidase B inhibition by maackiain, a pterocarpan from the roots of Sophora flavescens. Bioorg Med Chem Lett. 2016;26: 4714–4719. pmid:27575476

16. Oh JM, Jang H-J, Kim WJ, Kang M-G, Baek SC, Lee JP, et al. Calycosin and 8-O-methylretusin isolated from Maackia amurensis as potent and selective reversible inhibitors of human monoamine oxidase-B. Int J Biol Macromol. 2020;151: 441–448. pmid:32087226

17. Zhang J, Guo Q, Wei M, Bai J, Huang J, Liu Y, et al. Metabolite Identification and Pharmacokinetic Profiling of Isoflavones from Black Soybean in Rats Using Ultrahigh-Performance Liquid Chromatography with Linear-Ion-Trap–Orbitrap and Triple-Quadrupole Tandem Mass Spectrometry. J Agric Food Chem. 2018;66: 12941–12952. pmid:30444115

18. Chen J, Lin H, Hu M. Absorption and metabolism of genistein and its five isoflavone analogs in the human intestinal Caco-2 model. Cancer Chemother Pharmacol. 2005;55: 159–169. pmid:15455178

19. Steensma A, Bienenmann-Ploum ME, Noteborn HPJM. Intestinal uptake of genistein and its glycoside in the rat using various isolated perfused gut segments. Environ Toxicol Pharmacol. 2004;17: 103–110. pmid:21782720

20. Tsai T-H. Concurrent measurement of unbound genistein in the blood, brain and bile of anesthetized rats using microdialysis and its pharmacokinetic application. J Chromatogr A. 2005;1073: 317–322. pmid:15909536

21. Luo L-YY, Fan M-XX, Zhao H-YY, Li M-XX, Wu X, Gao W-YY. Pharmacokinetics and Bioavailability of the Isoflavones Formononetin and Ononin and Their in Vitro Absorption in Ussing Chamber and Caco-2 Cell Models. J Agric Food Chem. 2018;66: 2917–2924. pmid:29504397

22. Setchell KD, Brown NM, Zimmer-Nechemias L, Brashear WT, Wolfe BE, Kirschner AS, et al. Evidence for lack of absorption of soy isoflavone glycosides in humans, supporting the crucial role of intestinal metabolism for bioavailability. Am J Clin Nutr. 2002;76: 447–453. pmid:12145021

23. Vitale DC, Piazza C, Melilli B, Drago F, Salomone S. Isoflavones: estrogenic activity, biological effect and bioavailability. European Journal of Drug Metabolism and Pharmacokinetics. Springer-Verlag; 2013. pp. 15–25. https://doi.org/10.1007/s13318-012-0112-y pmid:23161396

24. Taneja I, Raju KSR, Challagundla M, Raghuvanshi A, Goel A, Wahajuddin M. LC-ESI-MS/MS method for bioanalytical determination of osteogenic phytoalexin, medicarpin, and its application to preliminary pharmacokinetic studies in rats. J Chromatogr B. 2015;1001: 9–16. pmid:26241182

25. Singh SP, Wahajuddin, Tewari D, Pradhan T, Jain GK. PAMPA permeability, plasma protein binding, blood partition, pharmacokinetics and metabolism of formononetin, a methoxylated isoflavone. Food Chem Toxicol. 2011;49: 1056–1062. pmid:21266188

26. Tyagi AM, Gautam AK, Kumar A, Srivastava K, Bhargavan B, Trivedi R, et al. Medicarpin inhibits osteoclastogenesis and has nonestrogenic bone conserving effect in ovariectomized mice. Mol Cell Endocrinol. 2010;325: 101–109. pmid:20570709

27. Tyagi AM, Srivastava K, Singh AK, Kumar A, Changkija B, Pandey R, et al. Formononetin reverses established osteopenia in adult ovariectomized rats. Menopause. 2012;19: 856–863. pmid:22781783

28. Xiao BX, Feng L, Cao FR, Pan R Le, Liao YH, Liu XM, et al. Pharmacokinetic profiles of the five isoflavonoids from Pueraria lobata roots in the CSF and plasma of rats. J Ethnopharmacol. 2016;184: 22–29. pmid:26923541

29. Yang B, Du S, Lu Y, Jia S, Zhao M, Bai J, et al. Influence of paeoniflorin and menthol on puerarin transport across MDCK and MDCK-MDR1 cells as blood–brain barrier in vitro model. J Pharm Pharmacol. 2018;70: 349–360. pmid:29238980

30. Kansy M, Senner F, Gubernator K. Physicochemical high throughput screening: Parallel artificial membrane permeation assay in the description of passive absorption processes. J Med Chem. 1998;41: 1007–1010. pmid:9544199

31. Di L, Kerns EH, Fan K, McConnell OJ, Carter GT. High throughput artificial membrane permeability assay for blood–brain barrier. Eur J Med Chem. 2003;38: 223–232. pmid:12667689

32. Könczöl Á, Müller J, Földes E, Béni Z, Végh K, Kéry Á, et al. Applicability of a blood-brain barrier specific artificial membrane permeability assay at the early stage of natural product-based CNS drug discovery. J Nat Prod. 2013;76: 655–663. pmid:23565574

33. Müller J, Esso K, Dargó GG, Könczöl Á, Balogh GT, Esső K, et al. Tuning the predictive capacity of the PAMPA-BBB model. Eur J Pharm Sci. 2015;79: 53–60. pmid:26344358

34. Mensch J, Melis A, Mackie C, Verreck G, Brewster ME, Augustijns P. Evaluation of various PAMPA models to identify the most discriminating method for the prediction of BBB permeability. Eur J Pharm Biopharm. 2010;74: 495–502. pmid:20067834

35. Lemes LFN, De Andrade Ramos G, De Oliveira AS, Da Silva FMR, De Castro Couto G, Da Silva Boni M, et al. Cardanol-derived AChE inhibitors: Towards the development of dual binding derivatives for Alzheimer’s disease. Eur J Med Chem. 2016;108: 687–700. pmid:26735910

36. Al-Baghdadi OB, Prater NI, Van Der Schyf CJ, Geldenhuys WJ. Inhibition of monoamine oxidase by derivatives of piperine, an alkaloid from the pepper plant Piper nigrum, for possible use in Parkinson’s disease. Bioorganic Med Chem Lett. 2012;22: 7183–7188. pmid:23102654

37. Chlebek J, De Simone A, Hoš͗álková A, Opletal L, Pérez C, Pérez DI, et al. Application of BACE1 immobilized enzyme reactor for the characterization of multifunctional alkaloids from Corydalis cava (Fumariaceae) as Alzheimer’s disease targets. Fitoterapia. 2016;109: 241–247. pmid:26779945

38. Végh K, Riethmüller E, Hosszú L, Darcsi A, Müller J, Alberti Á, et al. Three newly identified lipophilic flavonoids in Tanacetum parthenium supercritical fluid extract penetrating the Blood-Brain Barrier. J Pharm Biomed Anal. 2018;149: 488–493. pmid:29182998

39. Könczöl , Rendes K, Dékány M, Müller J, Riethmüller E, Balogh GT. Blood-brain barrier specific permeability assay reveals N-methylated tyramine derivatives in standardised leaf extracts and herbal products of Ginkgo biloba. J Pharm Biomed Anal. 2016;131: 167–174. pmid:27592255

40. Simon A, Darcsi A, Kéry Á, Riethmüller E. Blood-brain barrier permeability study of ginger constituents. J Pharm Biomed Anal. 2020;177. pmid:31476432

41. Brillatz T, Jacmin M, Queiroz EF, Marcourt L, Slacanin I, Petit C, et al. Zebrafish bioassay-guided isolation of antiseizure compounds from the Cameroonian medicinal plant Cyperus articulatus L. Phytomedicine. 2020;70. pmid:32302934

42. Ayanlowo AG, Garádi Z, Boldizsár I, Darcsi A, Nedves AN, Varjas B, et al. UHPLC-DPPH method reveals antioxidant tyramine and octopamine derivatives in Celtis occidentalis. J Pharm Biomed Anal. 2020;191. pmid:32980795

43. Zhu J, Yi X, Zhang J, Chen S, Wu Y. Rapid screening of brain-penetrable antioxidants from natural products by blood-brain barrier specific permeability assay combined with DPPH recognition. J Pharm Biomed Anal. 2018;151: 42–48. pmid:29304411

44. Buonfiglio R, Engkvist O, Várkonyi P, Henz A, Vikeved E, Backlund A, et al. Investigating Pharmacological Similarity by Charting Chemical Space. J Chem Inf Model. 2015;55: 2375–2390. pmid:26484706

45. Rosén J, Rickardson L, Backlund A, Gullbo J, Bohlin L, Larsson R, et al. ChemGPS-NP Mapping of Chemical Compounds for Prediction of Anticancer Mode of Action. QSAR Comb Sci. 2009;28: 436–446.

46. Larsson J, Gottfries J, Muresan S, Backlund A. ChemGPS-NP: Tuned for Navigation in Biologically Relevant Chemical Space. J Nat Prod. 2007;70: 789–794. pmid:17439280

47. Rosén J, Lövgren A, Kogej T, Muresan S, Gottfries J, Backlund A. ChemGPS-NPWeb: chemical space navigation online. J Comput Aided Mol Des. 2009;23: 253–259. pmid:19082743

48. Gampe N, Nagy E, Kursinszki L, Béni S. Quantitative determination of isoflavonoids in Ononis species by UPLC‐UV‐DAD. Phytochem Anal. pmid:32929804

49. Takács-Novák K, Józan M, Szász G. Lipophilicity of amphoteric molecules expressed by the true partition coefficient. Int J Pharm. 1995;113: 47–55.

50. Avdeef A, Nielsen PE, Tsinman O. PAMPA—a drug absorption in vitro model. Eur J Pharm Sci. 2004;22: 365–374. pmid:15265506

51. Viklund L. ChemGPS-NP as a tool for predicting drug distribution across the blood-brain barrier. Uppsala University. 2010.

52. Liang J, Tian YX, Fu LM, Wang TH, Li HJ, Wang P, et al. Daidzein as an antioxidant of lipid: Effects of the microenvironment in relation to chemical structure. J Agric Food Chem. 2008;56: 10376–10383. pmid:18841976

53. McLeod GS, Shepherd MJ. Determination of the ionisation constants of isoflavones by capillary electrophoresis. Phytochem Anal. 2000;11: 322–326.

54. Maul R, Schebb NH, Kulling SE. Application of LC and GC hyphenated with mass spectrometry as tool for characterization of unknown derivatives of isoflavonoids. Anal Bioanal Chem. 2008;391: 239–250. pmid:18317741

55. Guo C, Tong L, Xi M, Yang H, Dong H, Wen A. Neuroprotective effect of calycosin on cerebral ischemia and reperfusion injury in rats. J Ethnopharmacol. 2012;144: 768–774. pmid:23123262

56. Wang Y, Dong X, Li Z, Wang W, Tian J, Chen J. Downregulated RASD1 and upregulated miR-375 are involved in protective effects of calycosin on cerebral ischemia/reperfusion rats. J Neurol Sci. 2014;339: 144–148. pmid:24548484

57. Wang Y, Ren Q, Zhang X, Lu H, Chen J. Neuroprotective mechanisms of calycosin against focal cerebral ischemia and reperfusion injury in rats. Cell Physiol Biochem. 2018;45: 537–546. pmid:29402799

58. Atkinson F, Cole S, Green C, Van de Waterbeemd H. Lipophilicity and other parameters affecting brain penetration. Curr Med Chem—Cent Nerv Syst Agents. 2002;2: 229–240.

59. Yamada M, Yasuhara H. Clinical pharmacology of MAO inhibitors: Safety and future. Neurotoxicology. 2004;25: 215–221. pmid:14697896