Material. The following reagents were used for the solid dispersion production: berberine (90%, Sigma-Aldrich), Kolliphor® Poloxamer 407 (encapsulant, 12,000 g/mol, Sigma-Aldrich), Tween 80 (Dinâmica) and ethanol (99.8%, Neon). All chemicals were of analytical grade and purchased from common sources unless otherwise mentioned. The following reagents were used for the AChE activity assay: acetylcholinesterase (AChE; E.C. 3.1.1.7 from electric eel, Sigma-Aldrich)), 5’,5-dithiobis(2-nitrobenzoic acid) (DTNB; 98%, Sigma-Aldrich), acetylthiocholine iodide (ASCh; 99%, Sigma-Aldrich), and potassium phosphate buffer (TFK; pH 7.5, Neon). Potassium bromide (spectroscopic standard, Sigma-Aldrich) was used in the infrared spectra analyses. Trichloroacetic acid (TCA), 6-hydroxy-2,5,7,8-tetramethylchroman-2-carboxylic acid (Trolox) and 2,2′-azobis(2-methylpropionamidine) dihydrochloride (AAPH) (from Sigma-Aldrich) were used in the oxidative hemolysis inhibition assay (OxHLIA) and thiobarbituric acid reactive substances (TBARS) assays. Acetic acid, sulforhodamine B (SRB), ellipticine, dexamethasone, trypan blue, lipopolysaccharide (LPS), tris-(hydroxymethyl) aminomethane (TRIS) (from Sigma-Aldrich), dimethyl sulfoxide (DMSO) (from Fisher Scientific), Dulbecco's Modified Eagle's Medium (DMEM) and RPMI-1640 medium, fetal bovine serum (FBS), Hank's balanced salt solution (HBSS), L-glutamine, nonessential amino acid solution (2 mmol/L), penicillin/streptomycin solution (100 U.mL− 1 and 100 mg mL− 1, respectively), trypsin and EDTA (from Hyclone) were used in the cytotoxic and anti-inflammatory activity assays. Ethanol, acetic acid, and distilled water were used in the cytotoxic and genotoxic analyses. Methanesulfonate (MMS) was used as positive control and onion bulbs (Allium cepa L, variety beta crystal) was selected from an organic garden. The purified acetylcholinesterase was purchased directly from Sigma-Aldrich, not involving manipulation of animals or parts of animals, tissues, or other derivatives. The blood used for the OxLHIA assay was obtained randomly from the normal blood analysis of the animals and not directly for the assay.
Solid dispersion production and physicochemical characterization. The berberine-loaded solid dispersion was produced according to Sá et al.43 with minor modifications. Poloxamer 407 (1.200 g), a PEO-PPO-PEO block copolymer, and Tween 80 (0.012 g), a PEO sorbitan monooleate, were added to ethanol (50 mL) and mixed under gentle stirring for 10 min. Then, berberine (0.120 g) was added and mixed for 5 min. The dispersion was sonicated (Fisher Scientific, 120W, 1/8’ probe) for 3 min under a pulse regime (30 s on and 10 s off) in an ice bath. Ethanol was evaporated in a forced air circulation oven at 40°C for 24 h and the resulting powder was stored at -10°C protected from light.
The interaction between Poloxamer and berberine was quantitatively evaluated by UV-Vis spectroscopy (OceanOptics, Red Tide USB 650 UV) as reported by Karavas et al. 60. Ethanol was used as a solvent and different amounts of the two components were added (mass proportions from 1:1 to 20:1 Poloxamer:berberine). Solutions were obtained in triplicate and absorbance at 350 nm was used to calculate the normalized interaction coefficient (Eq. 1), where A and A0 were the maximum absorbance at 350 nm for the solutions and an ethanolic solution of berberine alone.
\(F=\frac{A-{A}_{0}}{{A}_{0}}\) | (1) |
The thermal properties of the solid dispersion were analyzed by Differential Scanning Calorimetry (DSC, Perkin Elmer 4000) and samples were heated in aluminum pans (0°C to 350°C at 10°C.min-1) under nitrogen flow (100 mL.min-1). Fourier Transform Infrared spectra (FTIR; Frontier Perkin Elmer) was performed in potassium bromide pellets with a resolution of 2 cm-1 from 4500 to 425 cm-1 with 32 cumulative scans. Transmission electron microscopy (TEM; JEOL model JEM 2100, 200 kV) was performed to observe the morphology of the nanoparticles in parlodium-covered copper grids (300 mesh). In addition, a mixture of berberine and Poloxamer 407 was obtained by simply mixing them in a laboratory mortar in the same mass proportion found in the solid dispersion. The objective was to compare this physical mixture with the berberine-loaded solid dispersion43,35.
In vitro cytotoxicity, anti-inflammatory and antioxidant activity. The antihemolytic activity (OxHLIA) was assessed using the method described by Takebayashi et al.61 as fully described in previous work43. Briefly, sheep blood samples were collected from healthy animals and centrifuged at 1,000 g for 5 min at 10 ºC. Plasma and buffy coats were discarded and erythrocytes were firstly washed once with NaCl solution (150 mmol/L) and three times with phosphate-buffered saline solution (PBS, pH 7.4) 62. The erythrocyte pellet was then resuspended in PBS at 2.8% (v/v). Using a flat-bottom 48-well microplate, 200 µL of erythrocyte solution was mixed with 400 µL of either PBS solution (control) and was dispersed in PBS, or water (for complete hemolysis). After pre-incubation at 37 ºC for 10 min with shaking, AAPH (200 µL, 160 mmol/L in PBS) was added and the optical density was measured in a microplate reader (Bio-Tek Instruments, ELX800) at 690 nm. After that, the microplate was incubated under the same conditions and the optical density was measured every 10 min at the same wavelength for approximately 300 min. The percentage of the erythrocyte population that remained intact (P) was calculated by Eq. 2 (St and S0 correspond to the optical density of the sample at t and 0 min, respectively, and CH0 is the optical density of the complete hemolysis at 0 min).
$$P\%=100 \left(\frac{{S}_{t}-{CH}_{0}}{{S}_{0}-{CH}_{0}}\right)$$
2
Results were expressed as the delayed time of hemolysis (Δt), calculated by Eq. 3, where Ht50 is the 50% hemolytic time (min) graphically obtained from the hemolysis curve of each sample concentration. The Δt values were then correlated to the different sample concentrations 61 and the concentration able to promote a Δt hemolysis delay was calculated. The results were expressed as IC50 values (µg.mL− 1) at Δt 60 and 120 min, i.e. the sample concentration required to keep 50% of the erythrocyte population intact for 60 and 120 min. Trolox was used as a positive control.
$${\Delta }t\left(\text{min}\right)=\frac{{\text{H}\text{t}}_{50}\left(\text{s}\text{a}\text{m}\text{p}\text{l}\text{e}\right)}{{\text{H}\text{t}}_{50}\left(\text{c}\text{o}\text{n}\text{t}\text{r}\text{o}\text{l}\right)}$$
3
The capacity of the sample to inhibit the formation of thiobarbituric acid reactive substances (TBARS), such as malondialdehyde generated from the ex vivo decomposition of lipid peroxidation products, was evaluated using porcine brain cell homogenates, following the method described previously63. Trolox was used as the positive control. The results were expressed as IC50 values (µg.mL− 1), i.e. the sample concentration providing 50% of antioxidant activity.
To assess the cytotoxicity of the sample, the sulforhodamine (SRB) assay was performed according to a procedure previously established by Abreu et al.64) in triplicate. Berberine was dissolved in water:DMSO (1:1 vol) and berberine-loaded solid dispersion was dispersed in water, both at the same concentration of berberine (8 mg.mL1), and this stock solution was used to prepare successive dilutions. CaCo cell line, MCF-7 (breast adenocarcinoma), NCIH460 (non-small cell lung carcinoma), and VERO cells from DSMZ (Leibniz-Institute DSMZ - German Collection of Microorganisms and Cell Cultures) were selected as human tumour cell lines. Porcine liver cells (PLP2), a primary cell culture, were prepared according to the procedure described by Abreu et al. 64. These cells were treated for 48 h with the different sample solutions and the SRB assay was followed65. Ellipticine was used as a positive control. The results were expressed as GI50 values (concentration that inhibited 50% of the net cell growth).
For the anti-inflammatory activity determination66, the lipopolysaccharide (LPS)-induced nitric oxide (NO) production by a murine macrophage (RAW 264.7) cell line was quantified as nitrite concentration in the culture medium. The effect of the tested compounds in the absence of LPS was also evaluated, to observe if they induced changes in NO basal levels. In negative controls, no LPS was added. For the NO determination, a Griess Reagent System kit containing sulfanilamide, N1naphthyl ethylenediamine dihydrochloride (NED), and nitrite solutions were used. Dexamethasone was used as a positive control. The results were expressed as IC50 values (µg.mL-1), i.e. compound concentration providing 50% of NO production inhibition.
Cytotoxic and genotoxic analysis of berberine inAllium cepaL. The experiments were carried out as previously described in details by Fiskesjo67 and Sales et al.68. For the assessment of cytotoxicity and genotoxicity of berberine and the berberine-loaded solid dispersion, the onion bulbs were placed in vials with distilled water, constantly aerated, to obtain roots of 2.0 cm in length. For analysis of berberine and Poloxamer concentrations (treatments), an experimental group with five onion bulbs was set up. Before putting the roots in contact with their respective treatments, some roots were collected and fixed to serve as a control of the bulb itself, which was identified as the time of analysis of 0 h or control of the bulb itself (Co − 0h). Then, the other roots were put in their respective treatments for 24 and 48 h, procedures called exposure times 24 and 48 h, where roots were collected every 24 h. Positive control was prepared with methyl methanesulfonate (MMS), a known cytotoxic and genotoxic substance to the A. cepa test system, at the concentration 4x10− 4 mol.L− 1. All roots collected during the experiment were fixed in 3:1 Carnoy solution (ethanol: acetic acid) for up to 24 h. Glass slides were prepared according to the protocol proposed by Herrero et al.69 and analyzed under an optical microscope with a 40x objective lens. For each bulb, 1,000 cells were analyzed, totaling 5,000 cells for each control group (0 hour), each group exposure time 24 hour and each group exposure time 48 hour, totaling 15,000 cells analyzed for each concentration of tretments. For the MMS group 5,000 cells were analyzed in five independent experiments.
For estimates of the mitotic index, cells in interphase, prophase, metaphase, anaphase, and telophase were counted to determine the cytotoxic potential. The mitotic index (MI) or cell division index was calculated by Eq. 4.
$$\text{M}\text{I}=100 \text{x} \frac{\text{t}\text{o}\text{t}\text{a}\text{l} \text{n}\text{u}\text{m}\text{b}\text{e}\text{r} \text{o}\text{f} \text{d}\text{i}\text{v}\text{i}\text{d}\text{i}\text{n}\text{g} \text{c}\text{e}\text{l}\text{l}\text{s}}{\text{T}\text{o}\text{t}\text{a}\text{l} \text{n}\text{u}\text{m}\text{b}\text{e}\text{r} \text{o}\text{f} \text{c}\text{e}\text{l}\text{l}\text{s} \text{a}\text{n}\text{a}\text{l}\text{y}\text{z}\text{e}\text{d}}$$
4
Genotoxic potential (chromosomal alterations index, Eq. 5) was assessed by frequency of cell alterations such as micronuclei, colchicine metaphases, anaphase and telophase bridges, cells with adhesions, nuclear buds and multipolar anaphases.
$$\text{C}\text{h}\text{r}\text{o}\text{m}\text{o}\text{s}\text{o}\text{m}\text{a}\text{l} \text{a}\text{l}\text{t}\text{e}\text{r}\text{a}\text{t}\text{i}\text{o}\text{n}\text{s} \text{i}\text{n}\text{d}\text{e}\text{x}=100 \text{x} \frac{\text{t}\text{o}\text{t}\text{a}\text{l} \text{n}\text{u}\text{m}\text{b}\text{e}\text{r} \text{o}\text{f} \text{c}\text{e}\text{l}\text{l} \text{a}\text{l}\text{t}\text{e}\text{r}\text{a}\text{t}\text{i}\text{o}\text{n}\text{s} }{\text{T}\text{o}\text{t}\text{a}\text{l} \text{n}\text{u}\text{m}\text{b}\text{e}\text{r} \text{o}\text{f} \text{c}\text{e}\text{l}\text{l}\text{s} \text{a}\text{n}\text{a}\text{l}\text{y}\text{z}\text{e}\text{d}}$$
5
Enzyme Activity Assays
AChE activity assay and reaction kinetics. The determination of the acetylcholinesterase (AChE) activity was performed according to modifications in the methods of Ellman et al. 57 and Pereira et al. 70. In 96-well plate, the following reagents were pre-incubated for 10 min at 25°C: TFK (90 µL; 50 mmol/L), water (45 or 55 µL), AChE (15 µL; 0.09 U/mL) and different dilutions of (1) berberine (10 µL; 0.03, 0.1, 0.3, 1 and 3 µmol/L) or different dilutions of (2) berberine-loaded solid dispersion (10 µL; 0.03, 0.1, 0.3, 1 and 3 µmol/L). The control group was incubated with water. After incubation, DTNB (20 µL; 0.2 mmol/L) and acetylthiocholine iodide (ASCh; 20 µL; 800 µmol/L) were added. Absorbance was read at 412 nm for 4 min every 60 s using a Vis spectrophotometer (Thermoplate TP-Reader). The concentration that inhibited AChE activity in 50% (IC50) compared with the control (H2O curve) was determined by nonlinear regression analysis.
Enzyme kinetic analysis was performed according to modifications in the methods described by Miranda et al.71. As in the AChE activity assays, the following reagents were pre-incubated for 10 min at 25°C: TFK, ultrapure water, enzyme (AChE or S1), and different dilutions of (1) berberine or (2) berberine-loaded nanoparticles. The control group was incubated with water. After incubation, DTNB and different concentrations of ASCh (0, 25, 50, 100, 400 and 800 µmol/L) or BSCh (0, 25, 50, 100, 400 and 800 µmol/L) were added. Absorbance was read at 412 nm for 4 min every 60 s using a Vis spectrophotometer (Thermoplate TP-Reader).
Reaction kinetics modeling. The kinetic constants of the enzyme inhibition reactions were determined by different strategies: i) Docking simulations; ii) The Lineweaver-Burk reciprocal plot; iii) Non-linear parameter estimation using the Particle Swarm Optimization coupled with the Gauss-Newton algorithm.
i) Docking simulations. 2D and 3D optimizations of Berberine structure were performed using Marvin Sketch 16.4 software (www.chemaxon.com). AutoDockTools 1.5.2 (ADT) was then used to convert the 3D structure of Berberine to the PDBQT file format 72. The X-ray crystal structure of AChE (PDB: 6G1V) was obtained from the Protein Data Bank (PDB) (http://www.rcsb.org). The co-crystallized ligand was extracted from each PDB file, and ADT was used to assign polar hydrogens, add Gasteiger charges, and save the AChE structures in the required PDBQT file format 71. Autodock Vina 1.272 was then used to perform molecular docking. An XYZ grid size of 30 by 30 by 30 Å was used for all structures with an exhaustiveness parameter of 16. The XYZ coordinates used for each structure were the following: 6ESY (3.6, -12.7, -12.3) and 6G1V (3.7; -4.5; 20.9). The docking conformations obtained were analyzed using the ProteinsPlus platform using the Pose View tool for 2D73. The predicted Ki (inhibition equilibrium constant) was calculated as follows: Ki = exp((ΔG.1000)/(Rcal.TK)) where ΔG is the predicted binding energy (cal/mol), Rcal is 1.98719 cal/(mol.K), and TK is 298.15 K. Structure representations were prepared using PyMOL (The PyMOL Molecular Graphics System, Version 1.3, Schrödinger, LLC).
ii) The Lineweaver-Burk reciprocal plot. The Lineweaver-Burk method was applied using the classical equations described by Bisswanger 74. This method applies a linear transformation to the inhibition kinetic models and is largely applied to enzymatic kinetic data. The following linearized models were considered for evaluation: a) Competitive inhibition (Eq. 6); b) Uncompetitive inhibition (Eq. 7); c) Mixed Non-competitive inhibition (Eq. 8); d) Pure Non-competitive inhibition (Eq. 9), where Vmax is the maximum reaction rate (µmol/L/min/U); Km is the Michaelis constant (µmol/L); [S] is the substrate concentration (µmo/L); [I] is the inhibitor (pure berberine or nanoencapsulated berberine) concentration (µmol/L); while \({K}_{ic}\) (competitive), \({K}_{iu}\) (uncompetitive) and \({K}_{i}\) are inhibition equilibrium constants (µmol/L).
$$\frac{1}{v}=\frac{1}{Vmax}+ \frac{Km . \left(1+ \frac{\left[I\right]}{{K}_{ic}}\right) }{Vmax . \left[S\right]}$$
6
$$\frac{1}{v}=\frac{\left(1+ \frac{\left[I\right]}{{K}_{iu}}\right)}{Vmax}+ \frac{Km}{Vmax . \left[S\right]}$$
7
$$\frac{1}{v}=\frac{\left(1+ \frac{\left[I\right]}{{K}_{iu}}\right)}{Vmax}+ \frac{Km . \left(1+ \frac{\left[I\right]}{{K}_{ic}}\right) }{Vmax . \left[S\right]}$$
8
$$\frac{1}{v}=\frac{\left(1+ \frac{\left[I\right]}{{K}_{i}}\right)}{Vmax}+ \frac{Km . \left(1+ \frac{\left[I\right]}{{K}_{i}}\right) }{Vmax . \left[S\right]}$$
9
iii) Non-linear parameter estimation. Non-linear parameter estimation was performed using a hybrid optimization method, which combines the particle swarm optimization (PSO) algorithm (a heuristic optimization method, based on empirical evolutionary rules that frequently mimic successful optimization strategies found in nature), and Gauss-Newton algorithms15. This approach was used for the computation of likelihood parameter confidence regions. The same kinetic models described in the previous section were considered in their non-linear form for parameter estimation, Equations (10) for Michaelis-Menten, (11) for the competitive model, (12) for the uncompetitive model, (13) for the non-competitive mixed model, and (14) for the non-competitive pure model, except for the inclusion of the Michaelis-Menten equation for analysis o kinetic data without the addition of an inhibitor. Furthermore, in this case, the parameters “a” and “b” that represent the factor with which Km and Vmax are multiplied to calculate the apparent constants in the presence of the inhibitor75.
$$v= \frac{{V}_{max}.\left[S\right]}{Km+\left[S\right]}$$
10
$$v= \frac{{V}_{max} . \left[S\right]}{b.Km . \left(1+\frac{\left[I\right]}{kic}\right) + \left[S\right]}$$
11
$$v= \frac{a.{V}_{max} . \left[S\right]}{b.Km + \left[S\right] . \left(1+\frac{\left[I\right]}{kiu}\right) }$$
12
$$v= \frac{a. {V}_{max} . \left[S\right]}{Km . \left(1+\frac{\left[I\right]}{kic}\right) + \left[S\right] . \left(1+\frac{\left[I\right]}{kiu}\right) }$$
13
$$v= \frac{a.{V}_{max} . \left[S\right]}{Km . \left(1+\frac{\left[I\right]}{ki}\right) + \left[S\right] . \left(1+\frac{\left[I\right]}{ki}\right) }$$
14
The absorbance data measured in the kinetic assay was firstly used to calculate the reaction rate (v, µM/min/U) for each time interval (1, 2, 3, and, 4 min) with Eq. (15) where: ΔAbs is the absorbance variation during the evaluated time interval (Δt, min); VW is the volume of the reaction media in the well (µL); l represents the optic path (cm); ε- represents the molar extinction coefficient of DTNB (L/mol.cm); and VE is the volume of enzyme solution in the reaction media (µL). After that, the data was submitted to the interpolation procedure in Matlab (R2021a) by the piecewise cubic Hermite interpolating polynomials (PCHIP), which was done using MATLAB's pchip function 76 to determine the initial reaction rate at approximately 0.2 min. These results were used in the PSO procedure.
$$v= \frac{\varDelta Abs .{ V}_{w}. l}{\varDelta t . \epsilon { . V}_{E}}$$
15
The interpolated reaction rate results were submitted to the hybrid optimization method using Spyder (4.2.5) with Python (3.8). This computational code is based on ESTIMA, developed by Schwaab et al.77 originally for FORTRAN, which employs a hybrid estimation method that combines the particle swarm optimization method with a Gauss-Newton procedure, and also performs the statistical analysis of results. The following conditions were applied for parameter estimation: confidence level of 95%, 100 particles, and a maximum number of interactions equal to 100.
The weighting least squares function was considered the objective function (Fobj) to be minimized (maximum likelihood method), as described in Eq. (16):
$${F_{obj}}=\sum\limits_{{i=1}}^{{NE}} {\sum\limits_{{j=1}}^{{NY}} {\frac{{\left( {{\mathbf{y}}_{{ij}}^{e} - {\mathbf{y}}_{{ij}}^{m}} \right)\left( {{\mathbf{x}}_{i}^{m},{\mathbf{\theta }}} \right)}}{{\sigma _{{ij}}^{2}}}} }$$
16
where NE is the number of experiments, NY is the number of output variables, \({{y}}_{ij}^{e}\) is the vector of experimental values for the output variable, \({{y}}_{ij}^{m}\) is the vector of the predicted values for the output variables, \({{x}}_{i}^{m}\) is the vector of input variables, θ is the parameter cluster, and \({\sigma }_{ij}^{2}\) is the experimental variance.
For the model discrimination, the Chi2 test was applied and the models were considered proper when the objective function was found within the Chi2 range. In addition, the quality of the estimated model parameters was evaluated concerning the Likelihood confidence region and graphics comparing experimental and predicted results.
Statistical analyses of the experimental results. For the cytotoxic and genotoxic analysis in Allium cepa L., data represented in the graphs are expressed as the mean ± standard deviation of three independent experiments, and means were compared by the Scott-Knott test at 0.05 significance. In the enzymatic experiments, data represented in the graphs are expressed as the mean ± standard deviation of at least three independent experiments. Data were analyzed using one- or two-way analysis of variance (ANOVA) followed by Tukey’s post hoc test. Values of p ≤ 0.05 were considered statistically significant. The statistical analysis was performed using Prism GraphPad 5.0 software or Statistica 7.0 software.