2. Materials and Methods
2.1. The analytes
The chemical structures of the investigated oleanane-type triterpenes and their glycosides of plant origin are presented in
Table 1.
2.2. Chemicals
The pharmacopoeial standards of arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid and arjunetin were purchased from Sigma Aldrich (St. Louis, MO, USA; p.a.). The organic modifiers for micellar mobile phase i.e. acetonitrile and isopropanol as well as surfactants: polyoxyethylene (23) lauryl ether (Brij35) and dodecyl sodium sulfate (SDS), were purchased from Merck (Darmstadt, Germany; p.a.). The buffer components i.e. citric acid and disodium hydrogen phosphate (Na2HPO4) were purchased from Sigma Aldrich (Sigma Aldrich, St. Louis, MO, USA; p.a.). Distilled water was obtained from the Direct-Q3 UV apparatus (Millipore, Warsaw, Poland).
2.3. Chromatographic equipment
The Shimadzu Vp liquid chromatographic system (Shimadzu, Kyoto, Japan) equipped with LC 10AT pump, SPD 10A UV-Vis detector, SCL 10A system controller, CTO-10 AS chromatographic oven and Rheodyne injector valve with a 20 µL loop was applied in the HPLC measurements.
2.4. Chromatographic conditions
The solutions of pharmacopoeial standards of the studied compounds were prepared in methanol (Merck, Darmstadt, Germany; p.a.) at a concentration of 1 mg/mL. All the oleanane-type triterpenes and their glycosides proved to be in the neutral form in solution under experimental conditions. The optimization process of the chromatographic separation was made before the experiment. The flow rate of the mobile phases was established to 1 mL/min and the temperature was set at 25oC. The tested compounds were detected with the UV light at 210 nm.
The C18 encapped column (Purosphere; 125 × 4 mm i.d., 5 µm; Merck, Darmstadt, Germany) was used as the stationary phase while the buffered solutions of both Brij35 as SDS were used as mobile phases. The mobile phases compositions were as follows: Brij35 at the concentrations: 0.04; 0.06; 0.08; 0.10 mol/dm3 (pH 7.4) with the addition of acetonitrile (10% v/v); SDS at the concentrations: 0.06; 0.08; 0.10; 0.12 mol/dm3 (pH 7.4) with an addition of isopropanol as an organic modifier (7% v/v). The buffer was prepared from the solutions of both Na2HPO4 (0.02 mol/dm3) and citric acid (0.01 mol/dm3).
The dead time values were measured from the citric acid peaks. All the reported logarithms of the retention factor were measured three times. The values of peak asymmetry factor were in the acceptable range.
2.5. Pharmacokinetic in silico studies
All the BBB-pharmacokinetic descriptors were calculated using the ACD/Percepta software (version 2012, Advanced Chemistry Development, Inc., Toronto, ON, Canada).
2.6. TLC-based bioatographic assay towards the AChE inhibitory activity
Six standards of oleanane-type triterpenes: arjunic acid, arjunolic acid, arjungenin, arjunglucoside I, sericic acid, and arjunetin - purchased in Sigma Aldrich (St. Louis, MO, USA) were prepared at the concentration of 1 mg/mL in double-distilled water : methanol (50:50 v/v) and they were applied separately at the surface of the aluminium normal phase 10 cm x 10 cm TLC plate (Silica gel 60 F254, Merck, Darmstadt, Germany) with an autosampler (Camag, Muttenz, Switzerland) as the 6 mm zones, distant from one another of 1.5 cm horizontally and 2 cm vertically. Every reference solution was applied as 4, 6, 8 and 10 µL volume bands on three TLC plates.
The TLC plates were later subjected to the enzymatic assay according to the previously published protocol [
18], with some modifications.
As the TLC plate was not developed in a TLC solvent system, but used directly in the TLC bioautographic assay, the authors modified the previously published protocol and sprayed the TLC with the substrate (2-naphtyl acetate) dissolved in distilled water at the quantity of 30 mg/20 mL. The TLC was dried in cold air and later the solution of AChE enzyme (AChE from electric eel type VI-S, Sigma Aldrich, St. Louis, CA, USA) dissolved in the aqueous solution of trizma buffer (pH 7.8) with bovine serum (500 mg/100 mL, Sigma Aldrich) at the quantity of 3 U/mL was sprayed on the TLC plate and incubated at the temperature of 37◦C for the following 20 min in a humid incubator. In the next step, the Fast Blue B solution (0.615 mg/mL) was sprayed on the plate and visualized active zones as white spots against violet background. The area of the discoloured zones was corresponding to the inhibitory strength of respective zones.
In the end the TLC plate was dried in the air, and analysed by Camag TLC visualizer at visible light. The peak areas of the discoloured zones were automatically calculated by the WinCats program (v. 1.4, Camag) and their size was compared to calculate the IC50 values that corresponded to the concentration of the standard giving half maximum inhibition of AChE enzyme.
2.7. Molecular docking procedure
The ligand molecules were obtained by the online SMILES translator (cactus.nci.nih.gov/translate) and subsequently optimized by using Avogadro 1.1.1 [
19] and the UFF force field [
20] (5000 steps, steepest descent algorithm). Flexible and optimized ligand molecules were docked into the binding pocket of the protein structure found in the PDB database (PDB:1EVE). Docking simulations were carried out in the AutoDockVina software [
21]. The procedure was performed within the cuboid region of dimensions of 22 × 30 × 34 Å
3 which covers the co-crystallized ligand present in the considered PDB record as well as the closest amino-acid residues that exhibit contact with this ligand. All the default procedures and algorithms implemented in AutoDockVina were applied during docking procedure. The rotatable torsional angles in both ligand molecules and the selected amino-acid sidechains within the binding cavity (Tyr334, Phe330, Phe75, Trp84, Glu199, Ser200, Tyr70, Tyr121, Trp279, Phe290, Phe331, Phe288, His440, Gln74, Leu282, Trp432, Asn85 and Asp285) were allowed to rotate. The visual inspections of each pose of the docked ligands were carried out in order to assure that the binding energies correspond to the structurally-analogous orientations. The procedure was validated in our previous work [
22].
2.8. Kinetics of AChE inhibition
The samples of the most potent AChE inhibitor i.e., ARG, were prepared in 12 dilutions in the concentration range of 0.00045 - 0.92 mM in dimethyl sulfoxide (DMSO≥99.7 %; Sigma Aldrich). The Ellman’s colorimetric method [
23] with some modifications [
24] was applied. Each of the tested ARG samples (15 μL) was mixed with 20 μL of the AChE solution (from electric eel, Type VI-S; Sigma Aldrich; 0.28 U/mL) and after 5 minutes completed with 35 μL of acetylthiocholine iodide (ATChI; Sigma Aldrich; 1.5 mmol/L), 175 μL of 0.3 mmol/L 5,5′-dithiobis-2-nitrobenzoic acid (DTNB containing 10 mmol/L NaCl and 2 mmol/L MgCl
2; Sigma Aldrich) and 100 μL with the Tris-HCl buffer (50 mmol/L, pH 8.0). The AChE, ATChI, and DTNB solutions were prepared in the Tris-HCl buffer. In order to eliminate the absorbance increase due to the spontaneous hydrolysis of the substrate, there were used “blank” samples composed of 15 μL of Tris-HCl buffer instead of ARG as well as of the above-mentioned compounds. The absorbance of the test samples was measured every minute for 32 minutes and it was subtracted from the absorbance of the “blank” sample. The background samples were prepared with 15 μL of each ARG solution and 330 μL of Tris-HCl buffer. The samples were incubated at room temperature for 30 minutes. The absorbance was measured at 412 nm (96-well microplate reader, Tecan Sunrise, Grödig, Austria). Each sample was analyzed in three repetitions. The linear regression analysis was conducted using the Minitab 18 Statistical Software (Minitab Inc., State College, PA, USA) and the values of the correlation coefficients, slopes, intercepts, and the standard errors were obtained.
2.9. Toxicity assay
To assess ARG toxic effect the ECOSAR (v. 1.11) free software was employed. Based on the ARG chemical structure both acute and chronic toxicity endpoints for fish, aquatic invertebrates (Daphnia), and green algae were measured.
4. Discussion
Natural products of plant origin can be interesting sources of compounds with neuroprotective properties. Among them the Traditional Chinese Medicine (TCM) herbs play an important role e.g.
Ginkgo biloba L. [
27],
Panax ginseng [
28,
29,
30,
31,
32], or
Scutellaria baicalensis [
33,
34,
35,
36,
37]. Other plants are also crucial source of a variety of compounds acting on the central nervous system (CNS) e.g.
Olea europaea L. [
37,
38,
39],
Vitis vinifera L. [
37,
38,
39],
Salvia officinalis L. [
40,
41,
42,
43,
44],
Melissa parviflora [
45],
Berberis integerrima [
46], or
Carissa edulis [
47]. The most important chemical groups of such compounds are saponins [
48,
49], tannins [
50,
51], flavonoids [
52,
53], alkaloids [
54], etc. Among the above-mentioned groups of the CNS-active compounds are triterpenes and their glycosides [
55,
56,
57,
58]. These compounds can affect the CNS including nerve cells of the brain and spinal cord which control many direct body functions and behaviour. In the context of neuroprotective properties, firstly it is important to confirm the ability of a compound to cross the BBB.
The
Terminalia arjuna accumulates bioactive triterpene glycosides (saponins) and aglycones (sapogenins), in a tissue-preferential manner [
59]. Many triterpenes demonstrate therapeutic efficacy. In most cases, they can cross the BBB and may affect the CNS including nerve cells of the brain and spinal cord which control many direct body functions and behaviours. They may also affect the autonomic nervous system which includes the regulation of internal organs, heartbeat, circulation and breathing.
Oleanane triterpenoids/saponins (derived from β-amyrin) have also been reported to have mainly cardioprotective potential [
59,
60,
61,
62]. Moreover, numerous studies have confirmed their antioxidant [
63], antimicrobial [
64,
65], anti-inflammatory [
66], anticancer [
67], precognitive [
12], hepatoprotective [
68], and others, activities.
It should be strongly emphasized that only the drug fraction unbound in media such as plasma can be transferred into body tissues.
In vitro methods including ultrafiltration or equilibrium dialysis are most often used to measure the fraction unbound value of a drug. These
in vitro obtained values are used not only for measurement of transfer rate into body tissues but also of the BBB permeability [
69]. It should be remembered that research on the penetration of compounds through the biological barriers, including the BBB one, is carried out using
in vivo methods in particular. However, for ethical and economic reasons, the need to use alternative methods to
in vivo one, including the
non-cell based-in vitro (biomimetic) and/or
in silico (computational) has been emphasized in recent years.
Both biomimetic and computational BBB-pharmacokinetic studies are commonly used in the laboratory practice at the first stages of an experiment on biologically active compounds (potential drugs) and constitute an important stage of research in the drug design process. At the stage of the
in silico studies the most important BBB-pharmacokinetic descriptors were calculated, i.e. the distribution of a substance in the blood-brain area, the rate of passive diffusion/permeability, the brain/plasma equilibration rate, the fraction unbound in plasma and the fraction unbound in brain. The blood-brain distribution (BB), frequently expressed as logBB, is defined as a ratio between the concentration in the brain and the concentration in the blood [
70,
71]. This experiment first identified 2 out of 6 tested compounds, i.e. arjunetin and arjunglucoside I, capable of crossing the BBB. However, it is commonly recognized that the most important parameter of the permeability through the BBB is the permeability – surface area product (PS) often expressed as logPS. These index is closely related to the cerebral blood flow (CBF) which is measured using various invasive as well as non-invasive techniques i.e. direct intravascular measurements, nuclear medicine, X-ray imaging, magnetic resonance imaging, ultrasound techniques, thermal diffusion, and optical methods. The most invasive methods require surgical access, arterial puncture, or catheterization while less invasive methods demand the intravenous injection of a contrast agent [
72]. The CBF is a very important parameter for brain viability and its functions because it ensures proper delivery of oxygen which is necessary for the neuronal oxidative metabolism of energy substrates. It is defined as the blood volume that flows per unit mass per unit time in brain tissue and is typically expressed in units of mL blood∕(100 g
tissue*min), or mL blood (100 mL
tissue*min) [
72] or in mL blood/(h*kg) [
73]. Taking into account PS values calculated in this experiment, arjunetin exhibited the highest BBB-permeability potential, followed by arjungenin, arjunglucoside I and sericic acid (
ex aequo) whereas both acids: arjunic and arjunolic one exhibited the lowest BBB permeability.
The scientific reports indicate that the time to reach brain equilibrium can be prolonged when the BBB permeability–surface area product (PS) or the fraction unbound in the brain decreases [
74], therefore it can be noticed that the lower values of the PS or Fb, the longer the time to reach brain equilibrium is required [
73]. In our experiment, no significant differences between Fb values were observed whereas the differences between the PS values are much greater (from 0.63 mL*h
-1*kg
-1 for both arjunic and arjunolic acids to 5 mL*h
-1*kg
-1 for arjunetin). Then, the longest time to reach brain equilibrium can be observed for the above-mentioned acids and the shortest for arjunetin. A high rate of penetration results from high BBB permeability as well as low brain tissue binding [
73].
In addition, analyzing values from
Table 2, it can be also seen that arjunetin and arjunglucoside I bind the least to blood plasma proteins (the highest value of free drug concentration, Fu) and shows the highest log BB value (0.73 and 0.12, respectively). The rest of the compounds have logBB values less than zero. Therefore, it can be presumed that among the tested compounds, arjunetin and arjunglucoside I are the substances that can penetrate the BBB to the greatest extent. However, the frequently used parameter for assessing the extent of the CNS distribution is also the ratio of brain/plasma partition coefficient, Kp,brain. This parameter – calculated for compounds that distribute solely by passive diffusion – is a function of the relative plasma and brain tissue unbound fractions at distribution equilibrium [
74]. In our case, most substances i.e. arjunic acid, arjunolic acid, arjungenin, and sericic acid, have Kp,brain values less than 1 which can result from more extensive binding to proteins in plasma than those in brain tissues. Other explanation can be significant impairment the in CNS distribution such as the efflux transport at the BBB [
74]. However, taking into account logBB values (
Table 2) it can be assumed that these compounds have simply lower CNS-distribution potential contrary to arjunetin and arjunglucoside I with Kp,brain values 5.1 and 1.25, respectively.
There exists the free drug theory that postulates that all the distribution processes of the active substance within biological barriers depend on the unbound drug concentration [
69,
75]. It must be emphasized that drug in the blood is present both in unbound and bound form to plasma proteins and erythrocytes. In our experiment, two substances i.e. arjunetin and arjunglucoside I have the highest value of the fraction unbound in plasma (0.051 and 0.050, respectively) in contrast to other compounds with Fb values in the range of 0.012 to 0.016. This could confirm earlier suppositions that arjunetin has the greatest ability among the tested compounds to cross the blood-brain barrier. Nevertheless, it is also hypothesized that drugs binding to protein can rapidly dissociate and permeate
in vivo through the BBB into the brain tissues [
69]. Therefore there may exist some differences between drug concentration obtained
in vivo in brain and that estimated
in vitro based on the free drug concentration. Nevertheless, the ability of most drugs to cross the BBB is nowadays estimated using the free drug fraction theory with reasonably acceptable results [
69].
The biomimetic studies were carried out to confirm (or not) the previously assumptions made based on the BBB-pharmacokinetic computational research. For this purpose, micellar liquid chromatography, using non-ionic Brij35(this type of chromatography is called the BMC) as well as anionic SDS surfactants, was applied. These methods are commonly used to assess the permeation of a substance through biological barriers [
76,
77,
78,
79,
80]. The concentration of a surfactant in a micellar mobile phase must be above the critical micellization concentration (cmc) whereas the commonly used stationary phase is the octadecyl-modified silica gel [
81,
82,
83]. Due to the wide application of micellar chromatography in the study of the penetration of compounds through biological barriers, it is a recognized technique in biomimetic studies on biologically active compounds.
Since the Brij35 micelle is assumed to be a kind of simple, chemical model of the biomembrane, the BMC technique can be useful in describing biological behaviours of different kinds of organic compounds. It can also mimic many biological processes such as BBB penetration, skin permeability, intestinal absorption and drug partitioning process in biological systems [
81,
82,
83], and others. In our research, the logarithms of retention factor extrapolated to pure water (log km), for both the BMC and SDS systems, have been determined. This parameter is recognized to be alternative to the logarithm of n-octanol/water partition coefficient (logPo/w) lipophilicity descriptor.
In the research, each system was previously optimized by selecting the appropriate concentrations of surfactants, selecting the organic modifier and its concentration in the mobile phase. The surfactant solutions were buffered (pH 7.4). Moreover, according to the Foley’s equation [
25], the interactions performed in the micellar systems have been characterized. Knowledge of the type of interactions between the analyte and the micelle, which in this case is a BBB model, can provide valuable information on the mechanism of interaction between a substance and a barrier. For this purpose, important physicochemical parameters such as K
MA – the analyte-micelle association constant and P
SW – the partition coefficient of an analyte between the stationary phase and water were calculated. Based on the above–mentioned parameters, one can conclude about the strength of analyte interaction with the biological membrane. Such studies can be very essential in the context of research on the biological activity of the tested compounds.
In the previous study [
84], it was proved that the logarithm of the analyte–micelle association constant (logK
MA) can characterize directly the passage of substances through the BBB comparable to the logBB pharmacokinetic parameter. Since the Foley’s model describes the retention behaviors of the analyte in the micellar system, which can be treated as a simple BBB model, therefore the parameters contained in it can characterize the biodistribution of the analyte in the BBB area. Very good linear relationships (R
2 > 0.9) between 1/k and C
M were obtained for all tested compounds (
Figure 1), confirming that the Foley’s equation correctly describes the retention of solutes in the tested BMC and SDS chromatographic systems. LogK
MA can be a useful tool for rapid assessment of the ability of a substance to cross the BBB, especially in the early stage of research. The obtained logK
MA-BMC values confirmed that both compounds: arjunetin and arjunglucoside I interact the most with Brij micelles which is recognized as a simply biological membrane model. In the SDS system, the matter is more complicated. Due to probably electrostatic interactions between the analytes and anionic micelles as well as strong retention of compounds, the intercepts for 3 out of 6 equations are negative. Unfortunately, the intercepts less than zero have no physico-chemical sense because they are equal to reciprocal of km parameter being the retention factor in the system in which the concentration of free surfactant (C
M) in the effluent is equal to zero. However, to eliminate the impact of possible electrostatic interactions, the log (km/K
MA) values calculated from the slopes of Eq. 1, were taken into account. These values have been treated as logBB values (logBB-SDS; see
Figure 2).
As shown in
Figure 2, there are no significant differences between logBB values obtained using computational and biomimetic methods. The BBB-pharmacokinetic biomimetic studies confirmed that arjunetin and arjunglucoside I can cross the BBB and have therefore the greatest BBB-penetration potential among the tested compounds while arjunic and arjunolic acids have the smallest one.
Analyzing the IC50 values obtained in the TLC-bioautography assay towards the AChE inhibition, it can be stated that among the tested compounds arjunolic acid was found to be the strongest inhibitor, whereas arjunetin was the weakest one among the tested compounds. However, the differences in the obtained IC50 values are insignificant. Thus, it can be concluded that all compounds have a very similar affinity to AChE which was later confirmed by molecular docking studies.
Taking into account the biological potential of other metabolites of plant origin, the compounds tested in the study exhibit relatively strong inhibitory potential. Previous results on triterpenoids confirmed their AChE inhibitory potential. In the study on the metabolites from
Centella asiatica, asiatic acid was found to be the strongest AChE inhibitor with the IC
50 value of 15.05 ± 0.05 µM [
85]. Also, the metabolites of
Garcinia hombroniana delivered information on a high inhibitory potential of 2β-Hydroxy-3α-O-caffeoyltaraxar-14-en28-oic acid present in the plant [
86]. In comparison to Amaryllidaceae alkaloids, like galanthamine that is registered as first line drug in the treatment of AD that was characterized by the IC
50 value of 3.520 µM [
87], the tested compounds seem to be promising.
As the ARG exhibits the lowest IC
50 value [mM] among the tested compounds, it was applied in the AChE inhibition kinetic studies using the colorimentric Ellman’s test [
23]. The Michaelis-Menten constant (Km) calculated based on the course of the curve (
Figure 6) was found to be 0.000011 mol/L. The obtained relationships show that the rate of substrate-enzyme binding is concentration-dependent and reaches the maximum velocity equal 2.2 x 10
-5.
The results of the docking study have been analyzed with respect to the mechanistic interaction pattern that may be significant in the context of interpretation of the obtained binding energies and recognizing the role of pharmacophore fragment. 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 3 (A) shows the superposition of all most favorable ligand poses whereas
Figure 3 (B) shows the most essential ligand-protein interactions. The extremely close match between superposed structures of all compounds agrees with the previous statement claiming that binding strength is determined by the pharmacophore-like, common molecular fragment of all compounds.
The detailed pattern of ligand-enzyme interactions is illustrated in
Figure 3 (B), on the example of arjunglucoside I (i.e. the compound displaying the lowest IC
50 value [mM]). However, due to similar orientations of all ligands in the binding cavity, the majority of conclusions can be transferred to remaining compounds. All ligands prefer roughly the same binding position in the enzyme cavity which enables them to block the catalytic site (the proximity to the catalytic histidine, His440, can be observed). The central fragment of ligand molecule (composed of aliphatic, cyclic moieties) interacts with aromatic cluster of sidechains, created by His440, Phe290, Trp84, Trp279, Phe288, Phe331, Tyr121, Tyr70, and Tyr334. Such contacts have a character of the CH-π interactions, supported (in some cases, e.g. His44) by hydrogen bonding with the neighboring fragments of ligand. The hydroxyl groups located at the edge of aliphatic, condensed fragment of the ligand, interact with Arg289 and Ser286. Both these contacts occur via hydrogen bonding and, surprisingly, involve backbone fragments of the protein (ligand can only be a hydrogen bonding donor). One can speculate about an analogous interaction in the case of Ile286 (also a backbone fragment) but, due to the lack of rotation around peptide bonds in the docking procedure, this was not explicitly observed. Interestingly, the ligand contacts with non-aromatic, hydrophobic sidechains are marginal and include only Ile287 and Leu127. Even in these cases, such proximitiesare rather an opportunistic consequence of much stronger interactions occurring with other, adjacent amino-acid residues.
The moiety of type and character varying between molecules (topologically-equivalent to the glucopyranose residue in the case of arjunglucoside I, illustrated in
Figure 3) is located close to a set of polar amino-acid residues, including Asn85, Ser122, Gln69, and Asp72. The dominating character of involved interactions is hydrogen bonding, where the considered fragment of ligand molecule can play a role of both donor and acceptor. In spite of the presence of tryptophan and tyrosine sidechains in the close proximity of glucopyranosidic moiety, no CH-π stacking, characteristic for carbohydrate-protein binding, was observed. This may explain why this fragment of ligand molecule (or its lack) is not particularly crucial for binding strength; the hydrogen bond donors and acceptors present in this region of cavity can equally well be saturated by water molecules, providing roughly the same balance of energy.
Taking into account toxicity predicted values, it can be stated that ARG has a low potential for chronic and acute toxicity on fish, daphnia and green algae.