Molecular dynamics, MM/PBSA and in vitro validation of a novel quinazoline-based EGFR tyrosine kinase inhibitor identified using structure-based in silico screening

Computational investigation of Moringa oleifera phytochemicals targeting EGFR: molecular docking, molecular dynamics simulation and density functional theory studies

Epidermal growth factor receptor (EGFR) is a prominent target for anticancer therapy due to its role in activating several cell signaling cascades. Clinically approved EGFR inhibitors are reported to show treatment resistance and toxicity, this study, therefore, investigates Moringa oleifera phytochemicals to find potent and safe anti-EGFR compounds. For that, phytochemicals were screened based on drug-likeness and molecular docking analysis followed by molecular dynamics simulation, density functional theory analysis and ADMET analysis to identify the effective inhibitors of EGFR tyrosine kinase (EGFR-TK) domain. Known EGFR-TK inhibitors (1-4 generations) were used as control. Among 146 phytochemicals, 136 compounds showed drug-likeness, of which Delta 7-Avenasterol was the most potential EGFR-TK inhibitor with a binding energy of -9.2 kcal/mol followed by 24-Methylenecholesterol (-9.1 kcal/mol), Campesterol (-9.0 kcal/mol) and Ellagic acid (-9.0 kcal/mol). In comparison, the highest binding affinity from control drugs was displayed by Rociletinib (-9.0 kcal/mol). The molecular dynamics simulation (100 ns) exhibited the structural stability of native EGFR-TK and protein-inhibitor complexes. Further, MM/PBSA computed the binding free energies of protein complex with Delta 7-Avenasterol, 24-Methylenecholesterol, Campesterol and Ellagic acid as -154.559 ± 18.591 kJ/mol, -139.176 ± 19.236 kJ/mol, -136.212 ± 17.598 kJ/mol and -139.513 ± 23.832 kJ/mol, respectively. Non-polar interactions were the major contributors to these energies. The density functional theory analysis also established the stability of these inhibitor compounds. ADMET analysis depicted acceptable outcomes for all top phytochemicals without displaying any toxicity. In conclusion, this report has identified promising EGFR-TK inhibitors to treat several cancers that can be further investigated through laboratory and clinical tests.

Structure-based virtual screening, molecular docking, molecular dynamics simulation and MM/PBSA calculations towards identification of steroidal and non-steroidal selective glucocorticoid receptor modulators

Glucocorticoids have been used in the treatment of many diseases including inflammatory and autoimmune diseases. Despite the wide therapeutic effects of synthetic glucocorticoids, the use of these compounds has been limited due to side effects such as osteoporosis, immunodeficiency, and hyperglycaemia. To this end, extensive studies have been performed to discover new glucocorticoid modulators with the aim of increasing affinity for the receptor and thus less side effects. In the present work, structure-based virtual screening was used for the identification of novel potent compounds with glucocorticoid effects. The molecules derived from ZINC database were screened on account of structural similarity with some glucocorticoid agonists as the template. Subsequently, molecular docking was performed on 200 selected compounds to obtain the best steroidal and non-steroidal conformations. Three compounds, namely ZINC_000002083318, ZINC_000253697499 and ZINC_000003845653, were selected with the binding energies of -11.5, -10.5, and -9.5 kcal/mol, respectively. Molecular dynamic simulations on superior structures were accomplished with the glucocorticoid receptor. Additionally, root mean square deviations, root mean square fluctuation, radius of gyration, hydrogen bonds, and binding-free energy analysis showed the binding stability of the proposed compounds compared to budesonide as an approved drug. The results demonstrated that all the compounds had suitable binding stability compared to budesonide, while ZINC_000002083318 showed a tighter binding energy compared to the other compounds.Communicated by Ramaswamy H. Sarma.

In Vitro and In Silico Studies of Kinase Inhibitor of MAPK3 Protein to Determine Leishmania martiniquensis Treatment

PURPOSE

Leishmaniasis is a parasitic disease transmitted by the bite of the phlebotomine female sand fly. Currently, no reported effective vaccines are available for the treatment of leishmaniasis; consequently, restricting this disease completely depends on controlling its transmission. Mitogen-activated protein (MAP) kinases have been reported to be involved in the regulation of the flagellum length and hence play an important role in disease transmission, especially the MAPK3 protein. Therefore, the current work focused on identifying approved drugs that can inhibit the MAPK3 protein.

METHODS

First, the recombinant plasmid (pET28b( +) MAPK3) was cloned into E. coli strain BL21 using the heatshock method. Afterward, E. coli was induced using IPTG, and cells were harvested for protein purification in the next step. After that, the MAPK3 protein was purified using Ni-NTA column. Then, the inhibition kinase activity of the purified MAPK3 protein was performed using an ADP-Glo^™ Kinase Assay kit. Furthermore, the cytotoxicity of Leishmania cells were detected by alamarBlue^™ Cell Viability Reagent. Finally, the binding affinity within the binding site of MAPK3 protein was performed by computational methods.

RESULTS

Purification of the MAPK3 protein was done using an Ni-NTA column and a protein band was identified at the expected 44 kDa molecular weight. Afterward, the ability of commercial drugs (afatinib and lapatinib) to inhibit the purified MAPK3 kinase activity was performed using an ADP-Glo^™ Kinase Assay kit. The half-maximal inhibitory concentrations (IC₅₀) of two drugs inhibited the MAPK3 protein within the same range of IC₅₀ values (3.27 and 2.22 µM for afatinib and lapatinib, respectively). Furthermore, the cytotoxicity assay of compounds toward the extracellular promastigote and intracellular amastigote stages was investigated using alamarBlue^™ Cell Viability Reagent. The results showed that both drugs were more efficient against extracellular promastigotes and intracellular amastigotes of both Leishmania donovani and Leishmania martiniquensis. Finally, the molecular dynamics simulation (MD) was performed to study the intermolecular interactions of both drugs with MAPK3 protein. From 100 ns molecular dynamics simulation, the structural stability of both drugs in a complex with MAPK3 was quite stable.

CONCLUSION

This work was suggesting that afatinib and lapatinib act as MAPK3 inhibitors and might be developed for leishmaniasis treatment.

Siamenflavones A-C, three undescribed biflavonoids from Selaginella siamensis Hieron. and biflavonoids from spike mosses as EGFR inhibitor

Three undescribed biflavonoids (BFVs), siamenflavones A-C along with twelve BFVs were isolated from Selaginella siamensis Hieron. and Selaginella bryopteris (L.) Baker (Selaginellaceae). The chemical structures of undescribed compounds were established through comprehensive spectroscopic techniques, chemical correlations, and X-ray crystallography. The ten isolated BFVs, siamenflavones A-C, delicaflavone, chrysocauflavone, robustaflavone, robustaflavone-4-methylether, amentoflavone, tetrahydro-amentoflavone, and sciadopitysin were evaluated for the antiproliferative effects against four human cancer cell lines A549, H1975, HepG2 and T47D. Delicaflavone and robustaflavone 4'-methylether exerted strong effects on the four human cancer cell lines. Siamenflavone B, delicaflavone and robustaflavone 4'-methylether showed potent inhibitory activities against wild-type EGFR. The inhibition of the compounds was further supported by molecular docking and predictive intermolecular interactions. Molecular dynamics simulation studies of siamenflavone B and robustaflavone-4'-methylether complexed to EGFR-TK further supported inhibition of the compounds to the ATP binding site. Finally, analysis of pharmacokinetic and electronic properties using density-functional theory and known drug index calculations suggest that the compounds are pharmaceutically compatible for drug administration.