Identification of phytochemical as a dual inhibitor of PI3K and mTOR: a structure-based computational approach

Breast cancer is a common form of cancer that affects both men and women. One of the most common types of genomic flaws in cancer is the aberrations in the PI3K/AKT/mTOR pathway. The benefit of dual targeting PI3K as well as mTOR is that the kinase-positive feedback loops are more effectively inhibited. Therefore, in the current study, structure-based models like molecular docking, MM-GBSA, Qikprop, induced fit docking, simulated molecular dynamics (MD), and thermal MM-GBSA were used to identify the phytochemicals from the zinc 15 database, which may inhibit PI3K and mTOR. After docking the phytochemicals with PI3K (PDB 4FA6), ten ligands based on the docking score were selected, among which salvianolic acid C had the highest docking score. Hence, salvianolic acid A was also docked. All the ligands taken showed a binding energy of greater than − 30 kcal/mol. The predicted ADME showed that the ligands have druggable properties. By performing MD of the top five ligands and salvianolic acid A, it was found that ZINC000059728582, ZINC000257545754, ZINC000253532301, and salvianolic acid A form a stable complex with PI3K protein, among which ZINC000014690026 showed interaction with Val 882 for more than 89% of the time. Salvianolic acid A is already proven to suppress tumor growth in acute myeloid leukemia by inhibiting PI3K/AKT pathway, but the exact protein target is unknown. Therefore, the present study identifies new molecules and provides evidence for salvianolic acid A for dual inhibition. Further experiments must be performed both in vitro and in vivo to support the predictions of these computational tools. Graphical abstract


Graphical abstract
Molecular docking (HTVS, SP, XP precision) Top 10 molecules with Salvianolic acid A

Introduction
Breast cancer affects both men and women, with a peak incidence varying between people of different regions.The incidence of breast cancer in Asia is between the age of 40 and 60 years, and in western countries it is between 60 and 70 years [1].It is a leading cause of cancer deaths in this demographic globally.Albeit uncommon, breast cancer can affect even men.About 1 in 100 men (< 1%) is diagnosed with breast cancer [2,3].The crucial risk factors are family history, age, early or late menopause, and nulliparity.Other lifestyle factors like alcohol consumption, lack of physical activity, overweight, etc., also play a pivotal role in the occurrence of cancer [3,4].Based on the hormonal receptors, breast cancer can be categorized as a luminal subtype [estrogen receptor-positive (ER+) or progesterone receptor-positive (PR+)], a Her2-enriched subtype (Her2+; overexpressed Her2 receptor), and triple-negative (TNBC)/ basal-like subtype (ER−, PR−, and Her2).The pathway that is frequently activated in most cancer cases is the PI3K/Akt/mTOR pathway.It regulates cell survival, cell division, migration, cellular protein synthesis, glucose metabolism, and immune regulation [5].Studies have shown that two crucial genes of this pathway, PI3K catalytic subunit α (PIK3CA) and phosphatase and tensin homolog PTEN are the most frequently altered ones in breast cancer.Oncogenic activation in TNBC can be due to upregulation of upstream regulators like epidermal growth factor receptor (EGFR), underexpression of proline-rich inositol polyphosphatase, dysfunction or underexpression of PTEN, and activation of PI3K catalytic subunit α (PIK3CA).As a result, PI3K can also be a potential target in treating TNBC.PI3Ks have an important role in the cascade that leads to tumor cell growth.PI3Ks contain catalytic (p110) and regulatory (p85) subunits that can exist in four isoforms (α, β, δ, and γ) [6].They are classified into Class I, II, and III based on their specific different activation and downstream molecules.PTEN is a cellular antagonist for the PI3K activity with a lipid phosphatase activity, and it reduces the amount of PIP3 in the cell, hence playing the role of tumor suppressor [7].
Across different breast cancer subtypes, the most common abnormalities are aberrations in the PI3K/AKT/mTOR pathway.Therefore, in these tumors, greater than 70% harbor gene amplification or mutations in the pathway.Stimulation of this pathway is by the ligands or growth factors that are specific to receptor tyrosine kinases (RTKs) that also include EGFR, insulin, insulin-like growth factor 1 (IGF-1) receptor, fibroblastic growth factor, etc. [8].
Downstream to this pathway is mTOR, which is the complex that sets different biological functions in motion without a doubt, and is a well-studied target .[9,10].mTOR is a Ser/Thr kinase and a part of PI3K superfamily called class IV PI3K.It is present as two complexes, mTORC1 and mTORC2.mTORC1 constitutes the catalytic subunit, a regulatory associated protein of mTOR (Raptor), prolinerich AKT substrate of 40 kDa (PRAS40), and mLST8/GbL protein.mTORC2 constitutes rapamycin-insensitive of mTOR (Rictor), mammalian stress-activated protein kinase interacting protein 1 (mSIN1), and mLST8/GbL.There are well-established drugs such as rapamycin for inhibition of mTORC1.But still, mTOR is also a potential target because of the mTORC2 complex, which acts as a PDK2 that phosphorylates AKT which is required for a complete activation at Ser473 at carboxyl terminal.However, the role of PDK in tumor progression is unknown; this activation is required for tumor growth [11].Thus, an mTOR inhibitor that targets mTORC1 and mTORC2 would inhibit PI3K pathway activation more efficiently than rapamycin [10,12].
Various inhibitors of pan-PIK3 and mTOR pathways are individually approved in cancer treatment, and clinical trials for others are underway.Some of the approved inhibitors of PI3K are Idelalisib, Copanlisib, and Duvelisib.Approved mTOR3 inhibitors are rapamycin and its analogs (rapalogs) like everolimus, sirolimus, and temsirolimus as first-generation inhibitors.Currently, second-generation inhibitors are in clinical trials.However, a major setback of these inhibitors is the poor response clinically, and resistance to the pathway inhibition [8].
Targeting PI3K and mTOR simultaneously inhibits the receptor tyrosine kinase-positive feedback loops more efficiently than PI3K alone.Thus, the dual inhibitors of PI3K and mTOR can simultaneously inhibit the pathway both upstream and downstream [7].Another advantage of using a dual inhibitor is the lower chance of resistance generation [11,13].Researchers like Alejandra ortiz gonzále et al. [14], by using in silico tools, predicted the effect of Opuntia joconostle in various models of breast cancer for its antiproliferative property through PI3K/AKT/mTOR pathway, proving that computer-based prediction can be beneficial for identifying lead molecules for cancer.Therefore, the current study tries to identify a phytochemical with dual specificity for both PI3K and mTOR, which can be an effective molecule in the treatment of cancer.

Methodology
The computational studies were performed on a Linux-based system with 8 GB RAM using Maestro Schrödinger suite version 2018-3 (Schrödinger, LLC, New York).Tools like, Ligprep, Protein Preparation Wizard, GLIDE, Desmond, and WaterMap were used for the study.

Protein selection and preparation
Based on the literature, PI3K PDB 4FA6 has inbound PI3K, and mTOR dual inhibitor, 2-amino-8-cyclopentyl-4-methyl-6-(1H-pyrazol-4-yl)pyrido [2,3-d]pyrimidin-7(8H)-one, and extensive computational work has been done.Therefore, the same PDB protein was used in this study.The 4FA6 PDB was retrieved from the RCSB database and processed with Protein Preparation Wizard [15].The missing amino acid residues, hydrogen bonds, and side chains were added during protein preparation.The ionization state of the het group at pH 7.4 ± 0.5 was generated, H-bonds were assigned and optimized, and water molecules that are not essential for ligand binding (beyond 3.0 Å from het groups) were removed.In the end, the protein structure was minimized to the lowest energy level by using the Optimized Potentials for Liquid Simulations (OPLS3e) force field.

Ligand preparation
In this study, natural product molecules (170,408) obtained after applying a filter in the availability category "now" from the Zinc 15 database [16] were downloaded.After importing the molecules, they were optimized using the Maestro suite Ligprep tool to obtain appropriate three-dimensional structures optimized with the ionization state at pH 7.4 ± 0.5.The three-dimensional structures were used to determine the chirality of the molecules, and minimization was done using Optimized Potentials for Liquid Simulations force field.

Molecular docking and MM-GBSA calculation
Molecular docking was performed by using Glide module [17].First, a grid of size 20 Å was generated centroid to the inbound molecule 2-amino-8-cyclopentyl-4-methyl-6-(1H-pyrazol-4-yl)pyrido [2,3-d]pyrimidin-7(8H)-one.The receptor grid file was loaded, and the inbound ligand was redocked into the generated grid to validate the grid by calculating the root mean square deviation (RMSD) of pose generated before and after the docking.After validation, the prepared ligands of the zinc 15 database were selected from the project table to dock with the protein.Initial docking was performed in high-throughput virtual screening (HTVS) mode (123,940 molecules), standard precision (5214 molecules), and extra precision (153 molecules).The ligand poses after docking were then analyzed for interaction patterns with the protein.The ligand poses showing the required interactions with the protein were selected.The top ten hits obtained were analyzed for ligand pose and interaction with protein.After that, further evaluation was done by using molecular mechanics generalized Born surface area (MM-GBSA).
The prime MM-GBSA [18] uses the VSGB solvent model and OPLS3e force field.More negative Kcal/mol values indicate stronger binding, as MM-GBSA binding energies are approximate binding free energies.

Predicted pharmacokinetic property
The ADME properties of the selected phytochemicals were predicted by using the QikProp module.Drug-likeness of phytoconstituents was evaluated by using various descriptors and pharmacokinetic properties like predicted octanol/ water partition coefficient (QPlogPo/w), predicted aqueous solubility (QPlogPo/w), predicted apparent Caco-2 permeability (QPPCaco), predicted brain/blood partition coefficient (QPlogBB), Lipinski's rule of 5, and human oral absorption.

Induced fit docking
Induced fit docking (IFD) [19] was performed using Maestro for the selected top five ligands obtained from molecular docking and MM-GBSA.In docking studies, the amino acid interacting with ligands is rigid, which is not the case in physiological conditions.Therefore, IFD is used to permit flexibility in the protein.During IFD, the ligand and receptor van der Waals scaling was maintained at 0.50.Calculations were performed using the standard precision protocol, generating a maximum of twenty poses of the protein with ligands.All poses were analyzed individually, and the pose in which the maximum number of interactions was seen with key residues was taken forward for molecular dynamic simulation.

WaterMap analysis
To have a better understanding of water molecules' role present within the binding pocket of 4FA6, which may affect binding of ligands, the WaterMap module [20] from Maestro suite was used.The inbound ligand in the protein was defined as the region of interest for analysis, and the water molecules within 5 A° were analyzed.Before running the simulation for two nanoseconds using the OPLS3e force field, all the water molecules within the crystal structure were eliminated.After simulation, the selected ligands were evaluated for the possible overlapping of the hydration site and Gibbs energy of the water molecule.

Molecular dynamics simulation and thermal MM-GBSA
MD simulation was performed for the top five shortlisted after extra precision docking and MM-GBSA results.The Desmond module [21] of Schrödinger was used initially.Periodic boundary condition (PBC) of size of 10 × 10 × 10 Å with orthorhombic shape was prepared around the protein complex and solvated using simple point charge (SPC) as a solvent model.At this stage, the isoosmotic condition was maintained by adding sodium and chloride ions.The solvated system was minimized using the 2000 maximum iterations and 1.0 (kcal/mol/Å) convergence threshold for 100 ps.Finally, the solvated system was used in MD simulation at NPT (Normal pressure and temperature) ensemble for 100 ns at 1.01 bar pressure and 300 K temperature.A simulation interaction diagram was generated to analyze the results.
To understand the binding energy of the ligands in different frames, MD simulation and MM-GBSA were carried out, in which every tenth frame of 1002 frames generated from MD simulation was used for calculating binding energy.The binding energies of the ligands were calculated by the formula mentioned below.

Molecular docking and selection of ligand based on binding affinity
Before screening the ligands using structure-based screening like molecular docking, the docking protocol was validated by redocking the co-crystallized ligand.The RMS deviation was found to be 0.2101, showing that not much change has occurred in the redocked ligand's pose has, as shown in Fig. 1.
After validation, the screening of 170,408 phytochemicals obtained from the Zinc 15 database was docked using different precision modes described in the methodology.After high-throughput screening (HTVS), 123,940 hits were identified, out of which ligands having docking MMGBSA dG Bind (NS) = Complex − Receptor (from optimized complex) − Ligand (from optimized complex) scores of more than − 7 kcal/mol (5214) were selected for standard precision docking.A further cutoff of − 9 kcal/ mol was applied to narrow down the number of ligands to 153, which were then docked by extra precision docking.Among them, the top ten ligands are represented in Table 1.The docking score of the co-crystallized ligand was found to be − 11.076 kcal/mol, and its aminopyrimidine formed hydrogen bonds with the hinge residue Val 88 and Asp 964.Top 10 ligands had docking scores ranging between -14.37 and − 12.67 kcal/mol.Among them, the top ligand was found to be ZINC000014690026 (Salvianolic acid C).So, a literature survey on the compound was done, and it was found that it is a phytoconstituent of a Chinese herb called Salvia miltiorrhiza, which also has salvianolic acid A and salvianolic acid B [22].Both of them halt cancer progression by affecting the cell cycle, causing apoptosis, and adjourning metastasis [23].
Both the ligands displayed activity in sensitizing cancer cells to chemotherapy [23].Therefore, we performed docking and found that Salvianolic acid A and salvianolic acid B have docking scores of − 12.266 and − 6.316 kcal/mol, respectively.MM-GBSA was further used to study the top ten ligands, and salvianolic acid A was selected based on the XP-docking scores representing the binding affinity of the molecules to the 4FA6 PDB protein.

Ligand-binding energy calculation
The binding energies for the ligand were determined using the Prime-MMGBSA, which indicates the stability of the protein-ligand complex formed after docking.All the top ten ligands showed stability when forming a complex with protein with binding energy greater than − 30 kcal/mol, tabulated in Table 1.The binding energy of the co-crystallized ligand was − 56.74 kcal/mol.The top 10 ligands had scores between − 64.52 and − 30.68 kcal/mol.Salvianolic acid A formed the most stable complex with the protein.Salvianolic acid C had the highest docking score, and less binding energy (− 56.74 kcal/mol) than Salvianolic acid A. This may be due to the ligand energy of salvianolic acid C (0.19 kcal/mol) being more than that of salvianolic acid A (− 13.898 kcal/mol).However, the complex energy for both ligands salvianolic acid C (− 36,388.9kcal/mol) and salvianolic acid A (− 36,413.1 kcal/mol) was not much different.Further, ADME analysis was performed for the top five ligands and Salvianolic acid A because they have greater than − 30 kcal/mol binding energy.

ADME analysis
The Qikprop module [24] in Maestro suite predicted the ADME properties of the top ligands by using various descriptors given in Table 2 like QPlogS, QPlogPo/w, QPPCaco, QPlogHERG, % human oral absorption, and Lipinski's rule of five.All the molecules violated more than one rule of the five, which shows drug-like properties, except for ZINC000014690026 and Salvianolic acid A. The selected ligands had good hydrophilic and hydrophobic balance and good solubility in an aqueous medium, as predicted by QPlogP0/w and QPLogS values, respectively.The QPPCaco and the % human oral absorption were very poor for the molecules.Among all, ZINC000014690026 and Salvianolic acid A had the highest values.The QlogH-ERG values were more than − 5, except for Salvianolic acid A. Therefore, most ligands have no potential to show inhibition of the HERG potassium channel.

Induced fit docking (IFD)-SP
After ADME analysis, it was proved that the selected molecules had druggable properties.Further, to study how ligand binds to the protein when the protein structure is flexible, induced fit docking studies were performed.IFD helps to overcome the limitations of rigid docking.In IFD, about 20 different poses of each ligand were analyzed.The comparative summary between extra precision docking pose and induced fit docking pose interactions of the ligands has been tabulated in Table 3.In the case of ZINC000014690026 (Salvianolic acid C), H-bonds with Tyr 867 and Glu 880 were observed in the XP-docking pose, which was absent in the IFD.ZINC000005004613 had an additional H-bond with Glu 880 and Thr 887 during glide XP docking, which is not seen in IFD.For Ambocin, additional H-bond with Glu 880 and Lys 807 were observed, which were not present in the IFD pose, but were seen in XP-docking pose; there is a Pi-cation interaction with Lys 833, which is an essential amino acid for inhibitory activity of the ligand.The compound ZINC000257545754 showed new H-bonds with Val 882, Ala 885, Asp 841, Tyr 867, Lys 808, and a Pi-Pi stacking interaction with Trp 812 in IFD pose.In the case of Salvianolic acid A, additional H-bonds with Asp 841, Tyr 867, Lys 833, Lys 890, Ala 885, Val 882, and Pi-Pi stacking interaction with Tyr 812 were observed in the IFD pose, which indicates strong interactions with the binding site.
For further MD studies, IFD poses in which the ligand showed maximum interactions with essential amino acids (Lys 833, Val 882) in the binding site were selected for 100 ns MD simulation, as shown in Fig. 2.

WaterMap analysis for hydration site prediction
Water map analysis involves the calculating enthalpy and entropy of every hydration site relative to bulk solvent by using thermodynamic molecular dynamics simulation [25].WaterMap analysis considers the theory of Lazaridis and Karplus, that considers a significant contribution of solvent reorganization energy and entropy to the solvationfree energy in an inhomogeneous system.WaterMap lets us determine the presence of a hydration site with an easily replaceable high-energy hydration site, represented in red spheres.The ones with less energy are represented as green spheres tightly placed around the ligand.A holistic approach was used to determine the relation between the binding affinity of the hit molecules and the computed hydration site energetics quantitatively.
Analyzing the apoprotein area (around 5Ǻ of the ligandbinding pocket) revealed the presence of 34 water molecules.The drug docking pocket contains water molecules with the hydration energy ranging from − 4.5 to 4.97 kcal/mol with the occupancy ranging from 1 to 0.28, indicating the wide range of water molecules around the pocket, as shown in Table 4.
All the ligands were found to cover the hydrophobic pockets formed by the hydration sites 37, 32, 73, 77, 72, and 64.The hydrophilic pocket of the protein comprises hydration sites 4,13,59,104, 60, which have less negative energy indicating strong binding.The predicted overlaps greater than 0.4, and the available hydration sites within 5 angstroms of the selected top ligand in the docking site of 4FA6 have been tabulated in Table 5. Compound ZINC000014690026 contains 17 water molecules which are in contact or proximity with ligand to that of relative free energy between 4.84 and − 3.59 kcal/mol.The benzyl group of ligands makes π-stacking interaction, surrounded by the hydration sites 32, 37, 73, and 77 with positive binding free energy.These hydration sites 32, 37, 73, and 77 have positive free energy indicating the possibility of easy replacement of these hydration sites by the ligand.
Similarly, 14 water molecules are present around the ligand ZINC000005004613 with occupancy in the range of 0.3-0.98 kcal/mol and ΔG with free energy from 5.29 to − 3.59 kcal/mol.Most hydration sites around this ligand have positive free energy except for sites 8, 24, and 59 (− 3.59 kcal/mol).Compound ZINC000059728582 contains 20 hydration sites with free energy ranging from 4.97 to − 3.88 kcal/mol.Compound ZINC000257545754 has 21 hydration sites with free energy ranging from 4.84 to − 4.5 kcal/mol, compound ZINC000253532301 contains 24 hydration sites with free energy ranging from 5.29 to − 4.5 kcal/mol, and Salvianolic acid A has 18 hydration sites with free energy ranging from 5.29 to − 4.5 kcal/mol.

Molecular dynamics simulation
The deviation of the complex from the reference frame which is indicated by root mean square deviation has been shown for the top five ligands with Salvianolic acid A in Fig. 3. 2D interaction diagram of the ligand after MD simulation, showing the percentage of time the amino acid had interacted with the ligand, fraction of interaction, and type of interaction, is shown in Fig. 4. The protein-ligand ZINC000014690026 complex system contains 100,299 atoms, which includes 28,567 water molecules and eight sodium molecules to neutralize the charge of the system.The RMSD plot of protein was stable throughout the simulation, with RMSD varying between 2 and 3.6 Ǻ.But, a wide deviation in the ligand was seen till 50 ns, after which it was stabilized with RMSD ranging from 4 to 8 Ǻ.The ligand stabilization after 50 ns was due to the formation of an additional bond with Lys 802.The critical residue Val 882 showed hydrogen bond interaction for 89% duration with the hydroxyl group of ZINC000014690026.Interactions with the residues Asp 841, Glu 880, Val 882, and Lys 890 observed in XP docking were retained during MD simulation, but interactions like Asp 964 and Tyr 867 were not seen.Instead, interaction with Thr 887 and Lys 833 was observed.The protein-ligand ZINC000005004613 complex system contains 100,290 atoms, of which the water molecules were 28,565.The system neutralization was done with seven sodium ions.RMSD plot showed stability after 20 ns for the entire simulation.The drift in the RMSD might be due to the disappearance of H-bond and water bridge interactions with residue Ser 806, Asp 950, hydrophobic interaction with residues Ile 831, Met 95.New interactions with residue Lys833, Glu 880, and Phe 961 were also formed after 20 ns.The important residue Val 882 showed hydrogen bond interaction for 75% duration and water bridge interaction for 68% duration with the hydroxyl group of ZINC000005004613.The interactions with residues Asp 950, Glu 880, and Thr 887, which were observed in the XP docking, were lost during MD simulation.
Ligand ZINC000059728582-protein complex consisted of 100,290 atoms, including 28,560 water molecules and seven sodium atoms to neutralize the charge.RMSD plot showed stable ligand RMSD in the range 1.5-2 Ǻ till 20 ns, followed by a steep increase ranging from 2.5 to 3.5 Ǻ till 80 ns, then a fall in the range of 2-3 Ǻ.The protein showed stable RMSD for the entire duration with the range of 2-3.5 Ǻ.The change in the RMSD after 20 ns might be due to the hydrophobic interactions' loss with the residues Tyr 867, Ile 879, and hydrogen bond interaction with Glu 880.Formation of the new interactions was observed with the residues Asn 951, Leu 969, Ser 806, Lys 808, Gly 970, and Tyr 972.The hydroxyl group of the ligand formed a hydrogen bond with the key residue like Val 882 for 85% of the simulation time.
Ligand ZINC000257545754-protein complex consisted of 100,312 atoms, including 28,570 water molecules and seven sodium atoms for neutralization of the charge.RMSD plot showed a stable ligand RMSD for the entire 100 ns with RMSD ranging from 1.2 to 2.8 Ǻ.Similarly, protein also had stable RMSD for the simulation period, showing RMSD ranging from 2.4 to 3.2 Ǻ.A stable water bridge and hydrogen bond interaction was formed for 58% of the simulation with the key residue Val 882.Loss of interactions with residues Ser806, Lys807, Lys 833, and Glu880 was observed in XP docking, and the formation of new interactions with residues Val 882, Tyr 867, Asp841, and Ala805 was seen during the simulation.
Ligand ZINC000253532301-protein complex consisted of 100,360 atoms, including 28,582 water molecules and six sodium atoms for neutralization.RMSD plot for protein shows the stability with the deviation in the range 2.8-3.6 Ǻ, and the ligand showed a deviation in the range 1-2.8 Ǻ.After 15 ns, the ligand showed a fall in RMSD from 2.8 Ǻ, which might be due to the loss of interaction with Tyr 757, Asp 836, His 948, and Asn 964 residues.However, after 20 ns duration, a stable RMSD plot was seen with the ligand.Key residue Val 882 formed water-linked hydrogen-bonding interaction for 73% duration of the simulation.Formation of new hydrogen-bonding interactions with residues His948, Ala 805, and Asp836 and loss of the interactions with residues Lys 833, Lys 890, and Ser 806 were observed in extra precision docking during MD simulation.
The protein-ligand Salvianolic acid A complex used for the simulation contained 28,560 water molecules, and eight sodium ions to neutralize the system.The protein was found

Thermal MM-GBSA
The trajectory generated after MD simulation was used for Thermal MM-GSBA.The average binding energy was more than − 60 kcal/mol for all the ligands.The most stable complex among all was Salvianolic acid A with − 73.87 ± 9.56 kcal/mol, which also showed a stable trajectory during MD simulation.The stability was due to interactions with Ala 805, Ala 885, Lys 833, Asp 964, Tyr 867, and Val 882.The binding energies of the top 5 ligands and salvianolic acid A are tabulated in Table 6, and variation in the binding energy of structure obtained from the various frames of MD simulation is shown in Fig. 5.

Conclusion
Targeting PI3K and mTOR can synchronously inhibit both upstream and downstream signaling proteins in the pathway.A single PI3K/mTOR dual inhibitor has an advantage of causing less drug resistance.Therefore, in this study, in silico studies were opted for by using stratified approaches like HTVS, SP, and XP docking of phytoconstituents obtained from curated databases like Zinc 15.
After the docking, MM-GBSA, ADME prediction, Water-Map, and molecular dynamics prediction were performed.By molecular docking, it was found that Salvianolic acid C had the highest docking score, and therefore, we carried out a literature survey and found that Salvianolic acid A has a similar structure to it and it was already proven to suppress tumor growth in acute myeloid leukemia by inhibiting PI3K/AKT pathway.However, the molecular mechanism by which it acts has not been established.Therefore, molecular dynamics was performed to predict if the top five ligands and salvianolic acid A form a stable complex with the PI3K protein or not.In 100 ns MD, ZINC000059728582, ZINC000257545754, ZINC000253532301, and Salvianolic acid A showed very stable complex, which can be seen in their RMSD values throughout the simulation.ZINC000014690026 showed interaction with Val 882 for more than 89% of the time, proving that it occupies the same pocket as that of the known inbound ligand of 4FA6.From the current study, it could be concluded that Salvianolic acid A can be explored for molecular mechanism by which it showed tumor suppression effect in acute myeloid leukemia.Other ligands mentioned can further be evaluated by various in vitro and in vivo studies in various types of cancer which involves PI3K/mTOR pathway.

Fig. 3 Fig. 4
Fig. 3 RMSD plots of protein and ligand after MD simulation for 100 ns

Table 4
Predicted hydration site, occupancy, and thermodynamic properties of the available hydration sites within five angstroms selected ligands

Table 6
Binding energy of molecules post MD simulation; the value represented here is the average of energy calculated from every 10th frame of 1000 frames obtained from MD simulation