Computational investigations of physicochemical, pharmacokinetic, toxicological properties and molecular docking of betulinic acid, a constituent of Corypha taliera (Roxb.) with Phospholipase A2 (PLA2)

Background 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. Methods 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. Results 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). Conclusions Our computed properties may assist the development of analytical method to assay BA or to develop BA derivatives with better pharmacokinetic and toxicological profile.

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) [39]. Thus, inhibiting PLA2 activity and hence regulating the production of proinflammatory 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 [43] and the high molecular weight (85 kDa) cytosolic PLA2 [44]. Convincing evidence suggests that the14 kDa PLA2 might be a potential target for the modulation of inflammatory diseases [40], 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 [45].
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 (P Skin ), plasma protein binding (PPB), Fig. 2 Interaction of (a). BR4 (green) and (b). BA (green) with PLA2 visualized in pymol. The calcium ion (Ca + 2 ) is presented as magenta color sphere. Black dash indicates H bonding. Residues involve in hydrophobic interactions are shown in gray color Fig. 3 Interaction of (a). BR4 (gray) and (b). BA (gray) with PLA2 visualized in Drug Discovery studio 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 [46] or to develop BA derivatives with better pharmacokinetic & toxicological profile and having more potent antiinflammatory activity.

Computational methods
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 [52] 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 [53].  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) [54] was used for all calculations involving the solvents. All calculations involving solvation were performed using the optimized solutionphase 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 (P Skin ), plasma protein binding (PPB) and penetration of the blood-brain barrier (BBB), interaction with Pglycoprotein (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 R work /R free , 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) [55] since the structure was solved at 2.1 Å resolution with R work /R free 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 1.7.4.4, 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 [53] 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.

Protein-ligand docking
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 [60] 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 [60] 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 v4.0.100.13345 (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.

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 [61] and 6.85 [62].   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 [54] The values are summarized in 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).

Dipole moment
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, noctanol, 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 [63]. 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 = -E HOMO Electron affinity, A = -E LUMO The electrophilicity index (ω) is calculated according to equation derived (by Parr et al., [68]) [58] as follows: 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 study
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 [69]. BA was moderately permeable since values for absorption through Caco-2 cell (P Caco-2 ) was 21.86 [70]. The skin permeability (P Skin ) is a vital parameter for assessment of drugs and chemical that might require transdermal administration [71]. BA was found to be impermeable through skin since a calculated P Skin 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 (C brain /C blood ) 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 C brain /C blood ratio all chemicals fall under three categories namely high absorption to CNS (C brain /C blood value more than 2.0), moderate absorption to CNS (C brain /C blood value within 2.0-0.1) and low absorption to CNS (C brain /C blood value less than 0.1) [72]. The ratio of C brain /C blood (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  [73]. 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 [77]. Overexpression 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 [82]. Computational classification model can be used to screen molecules and predict the likeliness to be substrate for Pgp [82]. The in silico screening revealed that BA is a dual inhibitor and substrate for Pgp like quinidine [83]. This is because due to complex modulatory interactions with the Pgp which make BA to function as combination of substrate and inhibitor [83].
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).

Toxicological study
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 [84]. The results have been in summarized in Table 8.

Protein optimization
YASARA Energy Minimization Server (http://www.yasara.org/minimizationserver.htm) [85] 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 [85]. 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 1.7.4.4, 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
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 [55]. 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 [55]. 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 [5]. Therefore, BA has the ability to interact with the catalytic residue HIS47 and the base forming residue HIS6 of human PLA2.

Binding affinity and in vitro inhibitory activity
A comparison was made between the binding affinities and in vitro inhibitory activities (IC 50 ) 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% [5] inhibition which is consistent in accordance to our findings of binding affinity and in vitro inhibitory activity.

Conclusions
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 antiinflammatory 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 [83], 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 [46] 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.