- Research article
- Open Access
- Open Peer Review
Computational investigations of physicochemical, pharmacokinetic, toxicological properties and molecular docking of betulinic acid, a constituent of Corypha taliera (Roxb.) with Phospholipase A2 (PLA2)
BMC Complementary and Alternative Medicine volume 18, Article number: 48 (2018)
Betulinic acid (BA) is a natural triterpenoid compound and exhibits a wide range of biological and medicinal properties including anti-inflammatory activity. Therefore, this theoretical investigation is performed to evaluate (a) physicochemical properties such as acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP), aqueous solubility (logS), solvation free energy, dipole moment, polarizability, hyperpolarizability and different reactivity descriptors, (b) pharmacokinetic properties like human intestinal absorption (HIA), cellular permeability, skin permeability (PSkin), plasma protein binding (PPB), penetration of the blood brain barrier (BBB), (c) toxicological properties including mutagenicity, carcinogenicity, risk of inhibition of hERG gene and (d) molecular mechanism of anti-inflammatory action which will aid the development of analytical method and the synthesis of BA derivatives.
The physicochemical properties were calculated using MarvinSketch 15.6.29 and Gaussian 09 software package. The pharmacokinetic and toxicological properties were calculated on online server PreADMET. Further, the molecular docking study was conducted on AutoDock vina in PyRx 0.8.
The aqueous solubility increased with increasing pH due to the ionization of BA leading to decrease in distribution coefficient. The solvation energies in water, dimethyl sulfoxide (DMSO), acetonitrile, n-octanol, chloroform and carbon tetrachloride were − 41.74 kJ/mol, − 53.80 kJ/mol, − 66.27 kJ/mol, − 69.64 kJ/mol, − 65.96 kJ/mol and − 60.13 kJ/mol, respectively. From the results of polarizability and softness, it was clear that BA is less stable and hence, kinetically more reactive in water. BA demonstrated good human intestinal absorption (HIA) and moderate cellular permeability. Further, BA also exhibited positive CNS activity due to high permeability through BBB. The toxicological study revealed that BA was a mutagenic compound but noncarcinogenic in mice model. Moreover, molecular docking study of BA with PLA2 revealed that BA interacts with GLY22 & GLY29 through hydrogen bond formation and LEU2, PHE5, HIS6, ALA17, ALA18, HIS47 and TYR51 through different types of hydrophobic interactions. The binding affinity of BA was − 41.00 kJ/mol which is comparable to the binding affinity of potent inhibitor 6-Phenyl-4(R)-(7-Phenyl-heptanoylamino)-hexanoic acid (BR4) (− 33.89 kJ/mol).
Our computed properties may assist the development of analytical method to assay BA or to develop BA derivatives with better pharmacokinetic and toxicological profile.
Betulinic acid (BA), (3β-hydroxy-lup-20(29)-en-28-oic acid) (Fig. 1) is a natural pentacyclic lupane type triterpenoid compound and exhibits a wide range of biological and medicinal properties such as antivenom , anti-HIV [2, 3], antibacterial , antimalarial , anti-inflammatory [5,6,7] anthelmintic , antinociceptive , anti-HSV-1 [10, 11] and anticancer activities [12,13,14]. It is abundantly distributed throughout the plant kingdom . The birch tree (Betula spp., Betulaceae) is one of the most widely reported sources of BA which can be obtained in considerable quantities [16, 17]. BA can also be isolated from various sources including Corypha taliera , Ziziphus spp. (Rhamnaceae) [19, 20], Syzygium spp. (Myrtaceae) , Diospyros spp. (Ebenaceae) [22, 23] and Paeonia spp. (Paeoniaceae) . Betulin, the reduced form of BA, was one of the first natural products to be isolated from the bark of the white birch, Betula alba .
BA was shown to exert its diverse pharmacological activities with variable median inhibitory concentrations (IC50) such as anticancer activity by inhibiting DNA Topoisomerases (Topos) II at IC50 of 56.12 μM , anti-HIV activity at IC50 of 23.65 μM , anti-malarial activity at IC50 of 56.71 μM , anti-fungal property at IC50 of 14.23 μM , anti-protozoal activity at IC50 of 50 μM .
BA also inhibited DNA polymerase beta at IC50 of 30.65 μM , protein tyrosine phosphatase 1B (PTP1B) at IC50 value of 1.5 μM , inhibited T. brucei GAPDH COX-1, COX-2 and LT formation in vitro with IC50 values of 240 μM , > 125 μM, > 125 μM and 102.2 μM , respectively. BA displayed potent anti-inflammatory activity by inhibiting Phospholipase A2 (PLA2) and showed 30% and 40% inhibition of PLA2 at concentrations of approximately 2.5 and 5 μM, respectively .
Earlier studies showed that BA exhibited anti-inflammatory activity by inhibiting TNF- alpha and increasing the production of anti-inflammatory cytokine IL-10 . Kim and colleagues demonstrated that betulinic acid exerted its anti-inflammatory activity by inhibiting the nuclear factor-kappa beta pathway, production of nitric oxide (NO), prostaglandin E2 (PGE2), tumor necrosis factor- alpha, interleukin-6 (IL-6), and interleukin-1 beta levels . Anti-inflammatory activity was also seen in encephalitogenic T cells where betulinic acid inhibited IL-17 and IFN-γ production . Arachidonoyl trifluoromethyl ketone, bromoenol lactone, varespladib, varespladib methyl, ecopladib, efipladib, giripladib, pyrrophenone, pyrroxyphene, FPL67047XX, inhibitor 22, amide 23 (GK115), 2-oxoamides, 1,3-disubstituted propan-2-ones are all synthetic phospholipase A2 inhibitors that have been used for clinical cases . Among natural compounds, extracts of curcumin, Ginkgo biloba and Centella asiatica have demonstrated to be phospholipase A2 inhibitors and have even been used to treat neurological disorders characterized by neuroinflammation . Bernard and colleagues experimentally proved that betulin and betulinic acid were potent phospholipase A2 inhibitors .
Phospholipase A2 (PLA2) hydrolyzes the membrane glycerophospholipids and releases arachidonic acid eventually leading to the production of pro-inflammatory mediators such as leukotrienes, prostaglandins, platelet activating factors (PAF) . Thus, inhibiting PLA2 activity and hence regulating the production of pro-inflammatory mediators for the development of therapeutics against inflammatory diseases [40, 41] is a plausible approach. However, the PLA2 exists in two isoforms [40, 42], which includes the low molecular weight (14 kDa) Ca+ 2 dependent extracellular PLA2 found in mammalian pancreases, several snake venoms, human platelets, human placentas, rheumatoid synovial fluids  and the high molecular weight (85 kDa) cytosolic PLA2 . Convincing evidence suggests that the14 kDa PLA2 might be a potential target for the modulation of inflammatory diseases , here we have reported the molecular docking study of BA against human secretory (14 kDa) PLA2 (1KQU) to explore the molecular basis of anti-inflammatory action of BA.
Few computational and theoretical studies of BA have been reported. BA associates with human serum albumin via hydrogen bond with PHE206 & GLU354 and hydrophobic interactions with PHE206, ARG209, ALA210, ALA213, LEU327, GLY328, LEU331, ALA350 and LYS351, in the sub-domain IIA and IIB of the large hydrophobic cavity .
In this investigation, computational studies have been carried out to evaluate (a) physical and chemical properties such as acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP), aqueous solubility (logS), solvation free energy, dipole moment, polarizability, hyperpolarizability and different reactivity descriptors (chemical hardness, softness, chemical potential, electronegativity, electrophilicity index), (b) pharmacokinetic properties like human intestinal absorption (HIA), cellular permeability using Caco-2 cell model, skin permeability (PSkin), plasma protein binding (PPB), penetration of the blood brain barrier (BBB), (c) toxicological properties including mutagenicity, carcinogenicity, risk of inhibition of human ether-a-go-go-related (hERG) gene and (d) molecular mechanism of anti-inflammatory action of BA. The purpose of this study was to investigate the physicochemical, pharmacokinetic and toxicological properties and to correlate the calculated physicochemical properties with the absorption and distribution profile of BA. These in silico investigations will provide an insight for the development of analytical method to assay BA  or to develop BA derivatives with better pharmacokinetic & toxicological profile and having more potent anti-inflammatory activity.
The acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP) and aqueous solubility (logS) of betulinic acid (BA) over the pH range of 0.0 to 14.0at 298 K were calculated using MarvinSketch 15.6.29 (ChemAxon, Hungary) (http://www.chemaxon.com). The consensus logP method was applied to calculate distribution coefficient of the molecule.
The rapid progress in computational methods such as Gaussian family of methods (G1, G2, G2MP2, G3) and complete basis set extrapolation (CBS) method enable researchers to perform highly sophisticated calculations of enthalpies and thermochemical properties with minor errors in comparison to experimental data [47,48,49,50,51]. However, these methods are computationally very expensive. An alternative to these high cost calculations is the use of density functional theory (DFT) methods. Previous report revealed that  calculation of geometries (bond length, bond angle and dihedral angle) and thermochemical properties using DFT/B3LYP level of theory showed better agreement with the experiments. So, in the current investigation, the calculation of solvation free energy, dipole moment, polarizability, hyperpolarizability and global reactivity descriptor properties such as the chemical hardness, softness, chemical potential, electronegativity, electrophilicity index were conducted in gas phase and in different solvents namely water, dimethyl sulfoxide (DMSO), acetonitrile, n-octanol, chloroform and carbon tetrachloride with the B3LYP/6-31G(d) level of theory implemented in Gaussian 09 software package . All calculations were conducted using the optimized geometry which was confirmed by the absence of imaginary frequency in the lowest energy state of the molecule. The Solvation Model on Density (SMD)  was used for all calculations involving the solvents. All calculations involving solvation were performed using the optimized solution-phase structures.
The pharmacokinetic and toxicological properties were calculated using online server PreADMET (https://preadmet.bmdrc.kr/). The pharmacokinetic properties such as human intestinal absorption (HIA), in vitro cellular permeability using Caco-2 cell model, skin permeability (PSkin), plasma protein binding (PPB) and penetration of the blood-brain barrier (BBB), interaction with P-glycoprotein (Pgp) and metabolism (both phase I and phase II) were calculated and predicted. In addition, BA was virtually screened to evaluate toxicological properties such as mutagenicity, carcinogenicity and risk of inhibition of human ether-a-go-go-related (hERG) gene.
Molecular docking study
Preparation of target protein X-ray structure
The first step of docking study is to select an appropriate X-ray crystal structure of the target protein that is already bound to its known ligand. This is because during docking, software searches complementary binding site/(s) for ligand within the search space of the target protein. That’s why ligand bound conformation of protein structure is the prerequisite to perform molecular docking. Moreover, the structure should be solved with a reasonable accuracy which is reflected in the statistics for data collection and processing of X-ray crystal structure. In addition, refinement statistics such as Rwork/Rfree, RMS deviation from ideality (bonds and angles) and validation parameters including Ramachandran outliers, rotamer outliers, bad bond count or bad angle count etc. also indicate the quality of the model structure. In our current investigation, we have selected the crystal structure of secretory PLA2 complexed with 6-Phenyl-4(R)-(7-Phenyl-heptanoylamino)-hexanoic acid (BR4) (PDB code: 1KQU)  since the structure was solved at 2.1 Å resolution with Rwork/Rfree of 0.209/0.240 and the validation parameters indicates good quality of the model structure. Water molecules and hetero atoms were then removed from the protein using PyMOL (Version 126.96.36.199, Schrödinger). Energy minimization was performed by applying YASARA force field level of theory in YASARA Energy Minimization Server (http://www.yasara.org/minimizationserver.htm).
Preparation of ligands
The structure of the BR4 and BA were drawn and optimized in Gaussian 09  with density functional theory (DFT) at the B3LYP/6-31G(d) level of theory. The optimized structures of ligands (BR4 and BA) were saved in PDB format for docking study.
Several virtual screening tools are available to perform molecular docking such as BINDSURF, METADOCK, Lead Finder, FlexScreen, AutoDock vina etc. Scoring functions implemented in these softwares calculate the protein-ligand interactions energy by computing electrostatic, Van der Waals and hydrogen bonding terms [56,57,58,59,60]. In addition to these scoring functions AutoDock vina also calculates the hydrophobic interactions  which is crucial for calculating the interaction energies of protein and hydrophobic ligand. Since BA is hydrophobic in nature, we used AutoDock vina to explore not only the electrostatic, Van der Waals and hydrogen bonding but also the hydrophobic interactions.
The docking of target protein with the ligands was conducted using AutoDock vina  in PyRx 0.8 (https://pyrx.sourceforge.io/). Docking study was performed to get a set of possible conformations and orientations for the ligand at the binding site. Using PyRx software, the PLA2 and ligands were prepared after which docking was conducted using a grid whose center was (56.5961, 34.0180, 42.4808) and the dimensions were (25.00, 25.00, 25.00) Å. During the docking analysis, the drug molecules were flexible and the macromolecule was kept rigid. Ligand displaying the lowest binding affinity and ability to bind in the binding pocket of protein was chosen as the best conformation. The interactions of different residues of protein with ligands such as hydrogen bonds, electrostatic interactions, hydrophobic interactions and bond distances were analyzed by PyMOL and Discovery Studio visualizer v188.8.131.5245 (Figs. 2 and 3).
Before docking of BA, the protocol was validated by re-docking BR4 into the binding pocket of 1KQU to get the docked pose and root mean square deviation (RMSD). The result revealed that the first pose of BR4 almost superimposed (RMSD of heavy atoms constituting the backbone of molecule is 0.5271) with the experimental crystal structure of BR4 (Fig. 4). Thus, the docking method has reasonable accuracy and reproducibility and can be used for further docking experiments.
Results and discussion
Calculation of pKa
The calculation of pKa revealed that BA has two ionized species along with the unionized form (Fig. 5). It is clear from the figure that non-ionized form (Betulinic acid_1) predominates over the pH range of 0.0 to 4.6 whereas the betulinate form (Betulinic acid_3) dominates from pH 4.8 to 14.0. However, small amount of Betulinic acid_2 was present over the pH of 0.0 to 2.4.
Distribution (logD) coefficient and partition (logP) coefficient
The logD was calculated using consensus method implemented in MarvinSketch 15.6.29. The logD vs pH of BA is presented in Fig. 6. The logD value decreases with increasing pH of the solution suggesting that the prevalence of the unionized form of BA decreases and the ionized form (betulinate) increases with increasing pH of the solution. This result is in accordance with the calculation of pKa.
The partition coefficient (logP) of BA was also computed and presented in Table 1. The logP calculated in MarvinSketch 15.6.29 was 6.64 whereas the experimental logP of BA were reported as 6.61  and 6.85 .
Aqueous solubility (logS)
The aqueous solubility of BA in terms of logS is presented in Fig. 7. The figure indicates that the solubility increases with increasing pH of the solution. However, the intrinsic solubility of BA was found − 7.34 (0.000000046 mol/L) and the aqueous solubility at pH 7.4 was − 4.79 (0.000016 mol/L) indicating that BA is practically insoluble in water.
Solvation free energy
The solvation free energies of BA were calculated with the SMD model  The values are summarized in Table 2. The solvation energies of BA in water, DMSO, acetonitrile, n-octanol, chloroform and carbontetrachloride were − 41.74 kJ/mol, − 53.80 kJ/mol, − 66.27 kJ/mol, − 69.64 kJ/mol, − 65.96 kJ/mol and − 60.13 kJ/mol, respectively. It could therefore be concluded that solvation free energy increases with decreasing polarity of polar nonprotic solvent (DMSO and acetonitrile). The highest value of solvation energy was in non-polar protic solvent (n-octanol). The free energy decreases with decreasing polarity of non-polar aprotic solvent (chloroform and carbon tetrachloride). Thus, it could be stated that BA associates with non-polar portion (hydrocarbon chain of n-octanol) via hydrophobic interactions. Electrostatic interactions might arise due to the association of the carboxyl group of BA with the polar portion (hydroxyl group of n-octanol) (Table 5).
The dipole moment of BA is found to be higher in different solvents than that of the gas phase. Table 3 presents the dipole moments computed in the gas phase and different solvents (water, DMSO, acetonitrile, n-octanol, chloroform and carbon tetrachloride) at the B3LYP level of theory with 6-31G(d) basis set using SMD solvation model. The dipole moments were 2.98D, 5.13D, 4.50D, 4.53D, 4.54D, 3.96D and 3.46D in the gas phase, water, DMSO, acetonitrile, n-octanol, chloroform and carbon tetrachloride, respectively. Therefore, increasing dielectric constant of the solvent is accompanied by a gradual increase in the dipole moment. In other words, the dipole moment increases with the increasing polarity of the solvent (Fig. 8).
Polarizability and first order hyperpolarizability
The polarizability(α) of BA was calculated using the following equation:
The quantities αxx, αyy and αzz are known as principal values of polarizability tensor.
The calculated polarizability of BA in different solvents is presented in Table 4, which showed that polarizability (α tot ) ranged from 351.90 to 448.67 a.u. The plot of polarizability vs. solvent is shown in Fig. 9. It is clear from the figure that the polarizability gradually increases when going from lower to higher dielectric constant. This suggests that the kinetic reactivity of BA increases with increasing polarity of the solvent .
The first order hyperpolarizability (β tot ) is the measure of the nonlinear optical activity and can be calculated using the following equation:
The β tot for different solvents was listed in Table 4, which displayed that the hyperpolarizability in different solvents ranged from 184.48 to 246.57 a.u. Moreover, the hyperpolarizability increases as the dielectric constant increases except for non-polar protic n-octanol where the hyperpolarizability is higher than that of acetonitrile and DMSO (Fig. 10).
Global reactivity descriptors
The highest occupied molecular orbital-lowest unoccupied molecular orbital (HOMO-LUMO) energy gap represents the stability or reactivity of molecules. Higher energy gap indicates greater stability and vice versa. The values of HOMO-LUMO energy gap in various solvents are presented in Table 5 and their trend is shown in Fig. 11. From Fig. 11, it is clear that the HOMO-LUMO energy gap is the highest in DMSO and acetonitrile, indicating greater stability of BA in these solvents, however lowest HOMO-LUMO energy gap was found in water meaning the titled molecule is less stable and hence, more reactive in water which is in agreement with calculated polarizability, chemical hardness and softness. The Table 5 and Fig. 11 suggest that the molecule is stabilized with decreasing polarity of the solvent i.e. the molecule is less likely to be kinetically reactive.
The global chemical reactivity descriptors such as softness, hardness, chemical potential and electrophilicity index can be calculated from the HOMO-LUMO energy gap of a molecule [64,65,66,67,68]. Using Koopman’s theorem for closed-shell molecules the hardness (η), chemical potential(μ), electronegativity (χ) and softness (S) are calculated according to the following equation:
Ionization potential, I = -EHOMO
Electron affinity, A = -ELUMO
The molecular properties of BA in different solvents are presented in Table 6. No systematic trend was found in the case of chemical hardness & softness, chemical potential, electronegativity and electrophilicity index. However, the highest chemical softness, electronegativity and electrophilicity index was found in water (polar protic) followed by n-octanol (non-polar aprotic) and carbon tetrachloride (non-polar aprotic).
Pharmacokinetic studies such as absorption, distribution and metabolism of BA was done by using the web based application PreADMET (https://preadmet.bmdrc.kr/). The calculated absorption, distribution and metabolism parameters are presented in Table 7.
The calculated human intestinal absorption HIA of BA (Table 7) was found to be 95.996% which suggests that BA is well absorbed through the intestinal cell . BA was moderately permeable since values for absorption through Caco-2 cell (PCaco-2) was 21.86 . The skin permeability (PSkin) is a vital parameter for assessment of drugs and chemical that might require transdermal administration . BA was found to be impermeable through skin since a calculated PSkin was − 2.11 encountered which proves BA is not permeable through the skin.
The distribution properties were assessed by measuring the brain to blood partitioning (Cbrain/Cblood) and plasma protein binding (PPB). Generally, compounds with more than 90% of PPB are classified as strongly bound chemicals whereas less than 90% are weakly bound chemicals (https://preadmet.bmdrc.kr/adme-prediction/). The calculated value of PPB for BA was 100.00% indicating BA is a strongly bound chemical and hence, the free form of BA will be less available in the systematic circulation. Based on Cbrain/Cblood ratio all chemicals fall under three categories namely high absorption to CNS (Cbrain/Cblood value more than 2.0), moderate absorption to CNS (Cbrain/Cblood value within 2.0–0.1) and low absorption to CNS (Cbrain/Cblood value less than 0.1) . The ratio of Cbrain/Cblood (8.20) suggests high absorption of BA to CNS indicating higher ability to cross blood brain barrier (BBB).
The exchangeable fraction (sum of unbound fraction and fraction dissociated from the protein) of drug determines the amount of drug that will penetrate the BBB . Besides, Pardridge and colleagues have reported the permeation of protein-bound compounds through the BBB [74,75,76]. Since BA is a strongly protein bound chemical (100%) so it could be assumed that BA will penetrate the BBB either due to the fraction dissociated from the plasma protein or as a protein-bound form or both.
P-glycoprotein (Pgp) is the product of the multi drug resistance (MDR) gene and an ATP dependent efflux transporter that affects the absorption, distribution and excretion of clinically important drugs . Over-expression of this protein, which may result in MDR is a major cause of the failure of cancer chemotherapy, and decreased efficacy of antibiotics [78, 79]. The prediction of Pgp substrates, which facilitates early identification and elimination of drug candidates of low efficacy or high potential of MDR [80, 81]. Identifying molecules that interact with Pgp transporters is important for drug discovery, but it is commonly determined through laborious in vitro and in vivo studies . Computational classification model can be used to screen molecules and predict the likeliness to be substrate for Pgp . The in silico screening revealed that BA is a dual inhibitor and substrate for Pgp like quinidine . This is because due to complex modulatory interactions with the Pgp which make BA to function as combination of substrate and inhibitor .
The computed metabolism of BA (Table 7) displayed that it is an inhibitor of CYP2C9 and a dual inhibitor and substrate for CYP3A4 due to complex modulation of CYP3A4 in phase I reaction. In phase II reaction BA is neither a substrate for UDP-glucuronosyltransferase (UGT) nor sulfotransferase (SULT).
The PreADMET server (https://preadmet.bmdrc.kr/) was used to evaluate the carcinogenicity of BA. It was found that although BA is mutagenic but it demonstrated noncarcinogenicity in mice. On the other hand, BA is unlikely to be an inhibitor of human ether-a-go-go-related (hERG) gene (low risk). Inhibition of the hERG gene has been linked to long QT syndrome . The results have been in summarized in Table 8.
YASARA Energy Minimization Server (http://www.yasara.org/minimizationserver.htm)  was used to optimize the structure of PLA2. This server utilizes a new partly knowledge-based all atom force field derived from Amber, whose parameters have been optimized to attain the protein structure as close as to its native structure with maximum accuracy . The simulation was performed in a water sphere containing ions.
The structure validation Z-scores and force field energies before and after the minimization are displayed in Table 9. The energies of the structure before and after optimization were − 68,549.6 and − 83,729.2 kJ/mol, respectively. This indicates that the protein structure was stabilized by an amount of − 15,179.6 kJ/mol. In addition, the structure validation Z-scores of before and after optimization was − 0.1 and 1.15 indicating good optimization of the target protein. The structure of PLA2 after energy minimization in a water sphere is presented in Fig. 11. PyMOL (Version 184.108.40.206, Schrödinger) was used to evaluate the root mean square deviation (RMSD) between the initial and final structure of PLA2 which was found to be 0.363.
Molecular docking of BA was carried out with human PLA2 using AutoDock Vina, to identify and understand the binding mode of BA and the intermolecular interactions with the target protein. All reported PLA2 inhibitors possess three key enzyme binding components such as Ca+ 2 binding oxygen atom, a HIS47 binding H-bond donor and a hydrophobic component that bind the active site of the enzyme . So, PLA2 inhibitors should form electrostatic interactions with Ca+ 2 ion, H bond formation with the catalytic HIS47 and different types of hydrophobic interactions with hydrophobic residues which line the active site cavity of the enzyme. In addition, the inhibitors should interact and displace the HIS6; a unique residue to human PLA2 enzyme . From our docking study, it was found that BA interacts with PLA2 at a binding affinity of − 41.00 kJ/mol whereas for BR4 it was − 33.89 kJ/mol. BA lies deep within the active site activity and makes numerous close contacts (< 5.4 Å) with the enzyme through hydrogen bond formation with GLY22 & GLY29 and hydrophobic interactions with the LEU2, PHE5, HIS6, ALA17, ALA18, TYR51 and the catalytic residue HIS47 (Figs. 2 and 3 and Table 10). However, no interactions of BA with calcium were observed as it was found in case of bovine PLA2 . Therefore, BA has the ability to interact with the catalytic residue HIS47 and the base forming residue HIS6 of human PLA2.
On the other hand, the standard compound BR4 forms hydrogen bonding with LYS62, HIS47, GLY29, GLY31, CYS28, VAL30 and LYS62. Different types of hydrophobic interactions were also displayed by BR4 with GLY29-VAL30, LEU2 and ALA1. Moreover, BR4 exhibited attractive electrostatic interactions with Ca+ 2 ions.
Binding affinity and in vitro inhibitory activity
A comparison was made between the binding affinities and in vitro inhibitory activities (IC50) of LY 311727, Manoalide, Indomethacin, Manoalogue, Aristolochic Acid and 1,1,1-Trifluoro-2-heptadecanone. The data has been presented in Table 11. Similar protocol described earlier was followed during performing docking of these compounds with PLA2 (1KQU). From the table it was clear that generally compounds with higher binding affinities also possess higher inhibitory activities. The binding affinities and reported PLA2 inhibitory concentration of BA are − 41.00 kJ/mol and ~ 2.5 μM for 30% & ~ 5 μM for 40%  inhibition which is consistent in accordance to our findings of binding affinity and in vitro inhibitory activity.
The acid dissociation constant (pKa), distribution coefficient (logD), partition coefficient (logP), aqueous solubility (logS), solvation free energy, dipole moment, polarizability, hyperpolarizability and different reactivity descriptors such as the chemical potential, electrophilicity, chemical hardness and chemical softness of betulinic acid have been calculated using MarvinSketch 15.6.29 and Gaussian 09 software. The calculated properties showed that betulinic acid has three species (one unionized and two ionized) over the pH range of 0.0 to 14.0 suggesting the unionized form is likely to exist in the acidic pH of the stomach however the betulinate form (ionized form) would be the predominant form in the alkaline pH of the intestine.
The logP, logD and logS values suggest that betulinic acid is non-polar and hydrophobic. This property would enable betulinic acid to pass through the lipid bilayer membrane of the cells. This hypothesis was corroborated by the results of HIA and permeability of betulinic acid through Caco-2 cell. It was found that BA has low skin permeability hence topical application BA as an anti-inflammatory agent can be achieved. Further, being lipophilic in nature betulinic acid would also exhibit positive CNS activity due to high permeability through BBB as evident from our in silico data.
The computational investigation of solvation free energy in different solvents explains the relative stability or reactivity of betulinic acid. Our theoretical calculation suggests that betulinic acid stabilized more in non-polar protic solvent due to hydrophobic interaction between the steroidal nucleus with the non-polar region of solvent and the electrostatic interaction between carboxyl group of betulinic acid with the polar portion of the solvent. Due to greater interaction with non-polar solvent the HOMO-LUMO gap becomes shorter as compared to the polar solvent (Table 5 and Fig. 11). It could also be inferred that it is this property which facilitates the hydrophobic interactions and formation of hydrogen bond of BA at the binding site cavity of PLA2.
BA was found to have 100% plasma protein binding and hence it could be deduced that it is unlikely to be eliminated promptly from the systemic circulation upon administration. Since BA acts as both inhibitor and substrate for P glycoprotein due to complex modulatory interactions with the Pgp , it is unlikely to be pumped out of the cells completely hence reducing the possibility of resistance mediated by efflux pump.
BA is likely to be metabolized by Cytochrome P450 3A4. However, it must be noted that it is an inhibitor for the same enzyme which could be overcome with the administration of Cytochrome P450 3A4 inducer. Since BA was also found to be an inhibitor of Cytochrome P450 2C9, it should not be used when a drug (which acts as a substrate for this enzyme) needs to be metabolized.
The toxicological study revealed that although betulinic acid is a mutagenic compound, it was noncarcinogenic in mice model. Moreover, molecular docking study revealed that betulinic acid interacts with GLY22 & GLY29 through hydrogen bond formation and LEU2, PHE5, HIS6, ALA17, ALA18, HIS47 and TYR51 through different types of hydrophobic interactions with a binding affinity of − 41.00 kJ/mol.
These computed molecular properties may assist to develop analytical method  to assay BA and the pharmacokinetic, toxicology and molecular docking study may provide a guide to synthesize betulinic acid derivatives with better pharmacokinetic & toxicological properties with potent phospholipase A2 inhibitory activity.
Highest occupied molecular orbital
Lowest unoccupied molecular orbital
Protein Data Bank
Root mean square deviation
Solvation Model on Density
Carvalho BM, Santos JD, Xavier BM, Almeida JR, Resende LM, Martins W, Marcussi S, Marangoni S, Stábeli RG, Calderon LA, Soares AM. Snake venom PLA2s inhibitors isolated from Brazilian plants: synthetic and natural molecules. Biomed Res Int. 2013;2013:1–8.
Fujoka T, Kashiwada Y, Robert EK. Anti-AIDS agents. Betulinic acid and platanic acid as anti-HIV principles from Syzygium claviflorum and the anti-HIV activity of structurally related triterpenoids. J Nat Prod. 1994;57:243–7.
Kashiwada Y, Wang HK, Nagao T, Kitanaka S, Yasuda I, Fujioka T, Yamagishi T, Cosentino LM, Kozuka M, Okabe H, Ikeshiro Y. Anti-AIDS agents. 30. Anti-HIV activity of oleanolic acid, pomolic acid, and structurally related triterpenoids 1. J Nat Prod. 1998;61(9):1090–5.
Bringmann G, Saeb W, Assi LA, Francois G, Narayanan AS, Peters K, Peters EM. Betulinic acid: isolation from Triphyophyllum peltatum and Ancistrocladus heyneanus, antimalarial activity, and crystal structure of the benzyl ester. Planta Med. 1997;63(03):255–7.
Bernard P, Scior T, Didier B, Hibert M, Berthon JY. Ethnopharmacology and bioinformatic combination for leads discovery: application to phospholipase a 2 inhibitors. Phytochemistry. 2001;58(6):865–74.
Huguet AI, del Carmen RM, Máñez S, Giner RM. Rı́os JL. Effect of triterpenoids on the inflammation induced by protein kinase C activators, neuronally acting irritants and other agents. Eur J Pharmacol. 2000;410(1):69–81.
Mukherjee PK, Saha K, Das J, Pal M, Saha BP. Studies on the anti-inflammatory activity of rhizomes of Nelumbo Nucifera. Planta Med. 1997;63(04):367–9.
Enwerem NM, Okogun JI, Wambebe CO, Okorie DA, Akah PA. Anthelmintic activity of the stem bark extracts of Berlina grandiflora and one of its active principles, Betulinic acid. Phytomedicine. 2001;8(2):112–4.
Kinoshita K, Akiba M, Saitoh M, Ye Y, Koyama K, Takahashi K, Kondo N, Yuasa H. Antinociceptive effect of triterpenes from cacti. Pharm Biol. 1998;36(1):50–5.
Ryu SY, Lee CK, Lee CO, Kim HS, Zee OP. Antiviral triterpenes from Prunella vulgaris. Arch Pharm Res. 1992;15(3):242–5.
Ryu SY, Lee CK, Ahn JW, Lee SH, Zee OP. Antiviral activity of triterpenoid derivatives. Arch Pharm Res. 1993;16(4):339–42.
Fulda S, Jeremias I, Steiner HH, Pietsch T, Debatin KM. Betulinic acid: a new cytotoxic agent against malignant brain-tumor cells. Int J Cancer. 1999;82(3):435–41.
Fulda S, Debatin KM. Betulinic acid induces apoptosis through a direct effect on mitochondria in neuroectodermal tumors. Med Pediatr Oncol. 2000;35(6):616–8.
Zuco V, Supino R, Righetti SC, Cleris L, Marchesi E, Gambacorti-Passerini C, Formelli F. Selective cytotoxicity of betulinic acid on tumor cell lines, but not on normal cells. Cancer Lett. 2002;175(1):17–25.
Cichewicz RH, Kouzi SA. Chemistry, biological activity, and chemotherapeutic potential of betulinic acid for the prevention and treatment of cancer and HIV infection. Med Res Rev. 2004;24(1):90–114.
O'Connell MM, Bentley MD, Campbell CS, Cole BJ. Betulin and lupeol in bark from four white-barked birches. Phytochemistry. 1988;27(7):2175–6.
Cole BJ, Bentley MD, Hua Y. Triterpenoid extractives in the outer bark of Betula lenta (black birch). Holzforschung. 1991;45(4):265–8.
Chowdhury A, Alam MA, Rashid RB, Al-Mansur MA, Rahman MS, Rashid MA. Steroids and triterpenoids from Corypha taliera Roxb: a critically endangered palm species of Bangladesh. Res J Med Plant. 2013;7:125–9.
Jagadeesh SG, Krupadanam GD, Srimannarayana G, Murthy SS, Kaur A, Raja SS. Tobacco caterpillar antifeedent from the gotti stem wood triterpene betulinic acid. J Agric Food Chem. 1998;46(7):2797–9.
Schühly W, Heilmann J, Calis I, Sticher O. New triterpenoids with antibacterial activity from Zizyphus joazeiro. Planta Med. 1999;65(08):740–3.
Chang CW, Wu TS, Hsieh YS, Kuo SC, Chao PD. Terpenoids of Syzygium Formosanum. J Nat Prod. 1999;62(2):327–8.
del Carmen RM, Giner RM, Manez S, Gueho J, Julien HR, Hostettmann K, Rios JL. Investigations on the steroidal anti-inflammatory activity of triterpenoids from Diospyros Leucomelas. Planta Med. 1995;61(01):9–12.
Singh P, Sharma S. Triterpenoid constituents of the seeds of Diospyros melanoxylon, Tecomella undulata and Terminalia bellirica. J Indian Chem Soc. 1997;74(6):504–5.
Lin HC, Ding HY, Wu YC. Two novel compounds from Paeonia s uffruticosa. J Nat Prod. 1998;61(3):343–6.
Krasutsky PA. Birch bark research and development. Nat Prod Rep. 2006;23(6):919–42.
Wada SI, Tanaka R. Betulinic acid and its derivatives, potent DNA topoisomerase II inhibitors, from the bark of Bischofia Javanica. Chem Biodivers. 2005;2(5):689–94.
Reutrakul V, Chanakul W, Pohmakotr M, Jaipetch T, Yoosook C, Kasisit J, Napaswat C, Santisuk T, Prabpai S, Kongsaeree P, Tuchinda P. Anti-HIV-1 constituents from leaves and twigs of Cratoxylum arborescens. Planta Med. 2006;72(15):1433–5.
Steele JP, Warhurst DC, Kirby GC, Simmonds MJ. In vitro and in vivo evaluation of betulinic acid as an antimalarial. Phytother Res. 1999;13(2):115–9.
Zhang Z, ElSohly HN, Jacob MR, Pasco DS, Walker LA, Clark AM. Natural products inhibiting Candida Albicans secreted aspartic proteases from Tovomita krukovii. Planta Med. 2002;68(01):49–54.
Domínguez-Carmona DB, Escalante-Erosa F, García-Sosa K, Ruiz-Pinell G, Gutierrez-Yapu D, Chan-Bacab MJ, Giménez-Turba A, Peña-Rodríguez LM. Antiprotozoal activity of betulinic acid derivatives. Phytomedicine. 2010;17(5):379–82.
Ma J, Starck SR, Hecht SM. DNA polymerase β inhibitors from Tetracera boiviniana. J Nat Prod. 1999;62(12):1660–3.
Choi JY, Na M, Hyun Hwang I, Ho Lee S, Young Bae E, Yeon Kim B, Seog AJ. Isolation of betulinic acid, its methyl ester and guaiane sesquiterpenoids with protein tyrosine phosphatase 1B inhibitory activity from the roots of Saussurea Lappa CB Clarke. Molecules. 2009;14(1):266–72.
Nyasse B, Nono JJ, Nganso Y, Ngantchou I, Schneider B. Uapaca genus (Euphorbiaceae), a good source of betulinic acid. Fitoterapia. 2009;80(1):32–4.
Wenzig EM, Widowitz U, Kunert O, Chrubasik S, Bucar F, Knauder E, Bauer R. Phytochemical composition and in vitro pharmacological activity of two rose hip (Rosa Canina L.) preparations. Phytomedicine. 2008;15(10):826–35.
Costa JF, Barbosa-Filho JM, de Azevedo Maia GL, Guimarães ET, Meira CS, Ribeiro-dos-Santos R, de Carvalho LC, Soares MB. Potent anti-inflammatory activity of betulinic acid treatment in a model of lethal endotoxemia. Int Immunopharmacol. 2014;23(2):469–74.
Kim KS, Lee DS, Kim DC, Yoon CS, Ko W, Oh H, Kim YC. Anti-inflammatory effects and mechanisms of action of coussaric and betulinic acids isolated from Diospyros kaki in lipopolysaccharide-stimulated RAW 264.7 macrophages, 2016. Molecules. 21(9):1206.
Blaževski J, Petković F, Momčilović M, Paschke R, Kaluđerović GN, Stojković MM, Miljković D. Betulinic acid regulates generation of neuroinflammatory mediators responsible for tissue destruction in multiple sclerosis in vitro. Acta Pharmacol Sin. 2013;34(3):424–31.
Ong WY, Farooqui T, Kokotos G, Farooqui AA. Synthetic and natural inhibitors of phospholipases A2: their importance for understanding and treatment of neurological disorders. ACS Chem Neurosci. 2015;6(6):814–31.
Irvine RF. How is the level of free arachidonic acid controlled in mammalian cells? Biochem J. 1982;204(1):3.
Mayer RJ, Marshall LA. New insights on mammalian phospholipase A2 (s); comparison of arachidonoyl-selective and-nonselective enzymes. FASEB J. 1993;7(2):339–48.
Posadas I, Terencio MC, De Rosa S, Payá M. Cavernolide: a new inhibitor of human sPLA2 sharing unusual chemical features. Life Sci. 2000;67(24):3007–14.
Friess H, Shrikhande S, Riesle E, Kashiwagi M, Baczako K, Zimmermann A, Uhl W, Büchler MW. Phospholipase A2 isoforms in acute pancreatitis. Ann Surg. 2001;233(2):204–12.
Dennis EA, Rhee SG, Billah MM, Hannun YA. Role of phospholipase in generating lipid second messengers in signal transduction. FASEB J. 1991;5(7):2068–77.
Clark JD, Milona N, Knopf JL. Purification of a 110-kilodalton cytosolic phospholipase A2 from the human monocytic cell line U937. Proc Natl Acad Sci U S A. 1990;87(19):7708–12.
Malleda C, Ahalawat N, Gokara M, Subramanyam R. Molecular dynamics simulation studies of betulinic acid with human serum albumin. J Mol Model. 2012;18(6):2589–97.
Khan MF, Rita SA, Kayser MS, Islam MS, Asad S, Rashid RB, Bari MA, Rahman MM, Al Aman DA, Setu NI, Banoo R. Theoretically guided analytical method development and validation for the estimation of rifampicin in a mixture of isoniazid and pyrazinamide by UV spectrophotometer. Front Chem. 2017;5(27):1–12.
Bach RD, Ayala PY, Schlegel HB. A reassessment of the bond dissociation energies of peroxides. An ab initio study. J Am Chem Soc. 1996;118(50):12758–65.
Glukhovtsev MN, Bach RD. Ab initio study on the thermochemistry of vinyl radical and cation. Chem Phys Lett. 1998;286(1):51–5.
Abboud JL, Castano O, Dávalos JZ, Gomperts R. The standard enthalpies of formation of 1-and 2-Adamantyl cations and radicals. An ab initio study. Chem Phys Lett. 2001;337(4):327–30.
Li ZH, Wong MW. Scaling of correlation basis set extension energies. Chem Phys Lett. 2001;337(1):209–16.
Espinosa-Garcı́a J. Theoretical C–H bond dissociation enthalpies for CH 3 OCl and CH 3 OBr. Chem Phys Lett 2000;316(5):563-568.
BecNe AD. Densityϋfunctional thermochemistry. III. The role of exact exchange. J Chem Phys. 1993;98:5648–52.
Frisch MJ, Trucks GW, Schlegel HB, Scuseria GE, Robb MA, Cheeseman JR, Scalmani G, Barone V, Mennucci B, Petersson GA, Nakatsuji H. Fox DJ Gaussian 09, revision a. 02. Wallingford, CT: Gaussian. Inc.; 2009.
Marenich AV, Cramer CJ, Truhlar DG. Universal solvation model based on solute electron density and on a continuum model of the solvent defined by the bulk dielectric constant and atomic surface tensions. J Phys Chem B. 2009;113(18):6378–96.
Hansford KA, Reid RC, Clark CI, Tyndall JD, Whitehouse MW, Guthrie T, McGeary RP, Schafer K, Martin JL, Fairlie DP. D-tyrosine as a Chiral Precusor to potent inhibitors of human nonpancreatic Secretory Phospholipase A2 (IIa) with Antiinflammatory activity. Chembiochem. 2003;4(2–3):181–5.
Sánchez-Linares I, Pérez-Sánchez H, Cecilia JM, García JM. High-throughput parallel blind virtual screening using BINDSURF. BMC Bioinf. 2012;13(14):S13.
Imbernón B, Cecilia JM, Pérez-Sánchez H, Giménez D. METADOCK: A parallel metaheuristic schema for virtual screening methods. Int J High Perform Comput Appl. 2017; p. 1–15.
Cerón-Carrasco JP, Cerezo J, Requena A, Zuñiga J, Contreras-García J, Chavan S, Manrubia-Cobo M, Pérez-Sánchez H. Labelling Herceptin with a novel oxaliplatin derivative: a computational approach towards the selective drug delivery. J Mol Model. 2014;20(9):2401.
Navarro-Fernández J, Pérez-Sánchez H, Martínez-Martínez I, Meliciani I, Guerrero JA, Vicente V, Corral J, Wenzel W. In silico discovery of a compound with nanomolar affinity to antithrombin causing partial activation and increased heparin affinity. J Med Chem. 2012;55(14):6403–12.
Trott O, Olson AJ. AutoDock Vina: improving the speed and accuracy of docking with a new scoring function, efficient optimization, and multithreading. J Comput Chem. 2010;31(2):455–61.
Elusiyan CA, Msagati TA, Shode FO, Mamba BB. Measurements of distribution coefficients and lipophilicity values for oleanolic acid and betulinic acid extracted from indigenous plants by hollow fibre supported liquid membrane. Bull Chem Soc Ethiop. 2011;25(3):321–32.
Claude B, Morin P, Lafosse M, Andre P. Evaluation of apparent formation constants of pentacyclic triterpene acids complexes with derivatized β-and γ-cyclodextrins by reversed phase liquid chromatography. J Chromatogr A. 2004;1049(1):37–42.
Suhasini M, Sailatha E, Gunasekaran S, Ramkumaar GR. Vibrational and electronic investigations, thermodynamic parameters, HOMO and LUMO analysis on Lornoxicam by density functional theory. J Mol Struct. 2015;1100:116–28.
Chattaraj PK, Maiti B, Sarkar U. Philicity: a unified treatment of chemical reactivity and selectivity. J Phys Chem A. 2003;107(25):4973–5.
Parr RG, Donnelly RA, Levy M, Palke WE. Electronegativity: the density functional viewpoint. J Chem Phys. 1978;68(8):3801–7.
Parr RG, Pearson RG. Absolute hardness: companion parameter to absolute electronegativity. J Am Chem Soc. 1983;105(26):7512–6.
Parr RG, Chattaraj PK. Principle of maximum hardness. J Am Chem Soc. 1991;113(5):1854–5.
Parr RG, Szentpaly LV, Liu S. Electrophilicity index. J Am Chem Soc. 1999;121(9):1922–4.
Yee S. In vitro permeability across Caco-2 cells (colonic) can predict in vivo (small intestinal) absorption in man—fact or myth. Pharm Res. 1997;14(6):763–6.
Yazdanian M, Glynn SL, Wright JL, Hawi A. Correlating partitioning and Caco-2 cell permeability of structurally diverse small molecular weight compounds. Pharm Res. 1998;15(9):1490–4.
Singh S, Singh J. Transdermal drug delivery by passive diffusion and iontophoresis: a review. Med Res Rev. 1993;13(5):569–621.
Chen C, Yang J. Predictive model of blood-brain barrier penetration of organic compounds. Acta Pharmacol Sin. 2005;26(4):500–12.
Tanaka H, Mizojiri K. Drug-protein binding and blood-brain barrier permeability. J Pharmacol Exp Ther. 1999;288(3):912–8.
Cornford EM, Pardridge WM, Braun LD, Oldendorf WH. Increased blood-brain barrier transport of protein-bound anticonvulsant drugs in the newborn. J Cereb Blood Flow Metab. 1983;3(3):280–6.
Pardridge WM, Sakiyama R, Fierer G. Transport of propranolol and lidocaine through the rat blood-brain barrier. Primary role of globulin-bound drug. J Clin Invest. 1983;71(4):900.
Terasaki TE, Pardridge WM, Denson DD. Differential effect of plasma protein binding of bupivacaine on its in vivo transfer into the brain and salivary gland of rats. J Pharmacol Exp Ther. 1986;239(3):724–9.
Schinkel AH. P-glycoprotein, a gatekeeper in the blood–brain barrier. Adv Drug Deliv Rev. 1999;36(2):179–94.
Cabrera MA, González I, Fernández C, Navarro C, Bermejo M. A topological substructural approach for the prediction of P-glycoprotein substrates. J Pharm Sci. 2006;95(3):589–606.
Kim RB, Fromm MF, Wandel C, Leake B, Wood AJ, Roden DM, Wilkinson GR. The drug transporter P-glycoprotein limits oral absorption and brain entry of HIV-1 protease inhibitors. J Clin Invest. 1998;101(2):289.
de Cerqueira LP, Golbraikh A, Oloff S, Xiao Y, Tropsha A. Combinatorial QSAR modeling of P-glycoprotein substrates. J Chem Inf Model. 2006;46(3):1245–54.
Wang YH, Li Y, Yang SL, Yang L. Classification of substrates and inhibitors of P-glycoprotein using unsupervised machine learning approach. J Chem Inf Model. 2005;45(3):750–7.
Joung JY, Kim HJ, Kim HM, Ahn SK, Nam KY, No KT. Prediction models of p-glycoprotein substrates using simple 2d and 3d descriptors by a recursive partitioning approach. Bull Kor Chem Soc. 2012;33(4):1123–7.
Zhou SF. Drugs behave as substrates, inhibitors and inducers of human cytochrome P450 3A4. Curr Drug Metab. 2008;9(4):310–22.
Sanguinetti MC, Tristani-Firouzi M. hERG potassium channels and cardiac arrhythmia. Nature. 2006;440(7083):463–9.
Krieger E, Joo K, Lee J, Lee J, Raman S, Thompson J, Tyka M, Baker D, Karplus K. Improving physical realism, stereochemistry, and side-chain accuracy in homology modeling: four approaches that performed well in CASP8. Proteins. 2009;77(S9):114–22.
Lister MD, Glaser KB, Ulevitch RJ, Dennis EA. Inhibition studies on the membrane-associated phospholipase A2 in vitro and prostaglandin E2 production in vivo of the macrophage-like P388D1 cell. Effects of manoalide, 7, 7-dimethyl-5, 8-eicosadienoic acid, and p-bromophenacyl bromide. J Biol Chem. 1989;264(15):8520–8.
Lobo IB, Hoult JR. Groups I, II and III extracellular phospholipases A2: selective inhibition of group II enzymes by indomethacin but not other NSAIDs. Agents Actions. 1994;41(1–2):111–3.
Santa Cruz Biotechnology. https://www.scbt.com/scbt/browse/chemicals-Protein-Interacting-Inhibitors-Activators-Substrates-protein-inhibitors-pla2-inhibitors/_/N-16mpppv. Access Date: 13 May 2017.
We would like to thank Dr. Shahidul M. Islam, Department of Chemistry, College of Science, Sultan Qaboos University for his help to conduct the research.
One of us (MAR) is thankful to the Ministry of Education, Government of the People’s Republic of Bangladesh for a research grant (No. 37.01.0000.078.02.018.13–125/10/2/2015) for the period 2014–15 to 2016–17 to carry out the work.
Availability of data and materials
Data and material are included in the article.
Ethics approval and consent to participate
Consent for publication
The authors declare that there is no competing interest.
Springer Nature remains neutral with regard to jurisdictional claims in published maps and institutional affiliations.