Introduction

Cancer is one of the primary reasons for mortality worldwide, and its prevalence continuously increases. According to GLOBOCAN 2020, about 10 million people died from cancer in 2020, and 19.3 million new cases appeared worldwide1. Cancer in the lungs, prostate, colon, and urinary bladder are most common in men, while breast, cervical, and uterine cancer frequently occur in women2.

EGFR (epidermal growth factor receptor) is a tyrosine kinase family's 170 kDa transmembrane glycoprotein. It is mainly a cell surface receptor containing three parts: intracellular tyrosine kinase, an extracellular ligand-binding domain, and a transmembrane anchoring region3,4. It binds to epidermal growth factor, which leads to its dimerization and autophosphorylation of tyrosine residues, which activate a cascade of intracellular signal transduction events and promote cell proliferation, angiogenesis, apoptosis resistance, cellular invasion, and metastasis5.

Mutation or overexpression of the EGFR gene causes different cancers type, including breast cancer, prostate cancer, ovarian cancer, and NSCLC (non-small-cell lung cancer)6. Five amino acids (codons 746 to codon 750) are responsible for point mutation of EGFR at the L858 position (catalytic domain of EGFR) and cause autophosphorylation and unregulated cell growth. Gefitinib and erlotinib (first-generation EGFR tyrosine kinase inhibitors) are reversible inhibitors that compete with ATP (Adenosine triphosphate) to bind at the catalytic site of EGFR and cause hindrance of autophosphorylation and cell proliferation7. The need for inhibitors keeps increasing; consequently, the effectiveness of first-generation inhibitors was hindered due to resistance raised for these inhibitors, which diminished the clinical response over a year. Substitution of threonine 790 with methionine, EGFRT790M (EGFR-L858R/T790M) is one of the most frequent resistance issues against first-generation tyrosine kinase inhibitors and causes a reduction in its inhibitor's activity and prevents their interaction at the catalytic domain8. Approximately 50% of NSCLC-affected patients with mutated EGFRT790M (EGFR-T790M) were resistant to first-generation tyrosine kinase inhibitors. Thus, the alteration in EGFR at T790M established the development of resistance against first and second-generation tyrosine kinase inhibitors and provoked the finding of inhibitors of the third-generation kinase. Neratinib, Dacomitinib, and Afatinib are second-generation kinase inhibitors that interact covalently with Cys_797 residue at the binding pocket of ATP, and their clinical application is restricted because of side effects such as gastrointestinal problems and skin diseases9,10.

Third-generation kinase inhibitors showed significant efficiency in NSCLC-affected patients possessing resistance to first/second-generation inhibitors. Osimertinib is a Food and Drug Administration (FDA)-approved third-generation inhibitor, expressing effective results in NSCLC patients11. Further, the development of fourth-generation kinase inhibitors comes into the limelight with the emergence of resistance to third-generation kinase inhibitors. The third-generation EGFR tyrosine kinase inhibitors, i.e., EGFRCSTMLR (EGFR-L858R/T790M/C797S) resistance developed by converting cysteine to serine at 797, which is less reactive and inhibits the interaction between inhibitors and cysteine amino acid12. Numerous side effects are also associated with third-generation inhibitors, so an alternative drug based on natural sources, having fewer side effects and low cost, needs to be discovered13.

In this regard, Streptomyces can be a perfect alternative for discovering EGFR tyrosine kinase inhibitors due to their greater bioavailability, bountiful sources of secondary metabolites, biotic friendliness, lesser side effects, non-toxic and highly effectiveness14. Genus Streptomyces belongs to the phylum actinobacteria and is primarily ubiquitous, aerobic, Gram-positive, and filamentous soil bacteria15. They can produce many secondary metabolites due to the survival of their spores in adverse environmental conditions. Today, the genus Streptomyces is responsible for developing 80% of antibiotics, making it the most important genus in medication discovery against various diseases16. They also create complex secondary metabolites such as enzymes, protein kinase inhibitors, apoptosis inducers, and caspase-3 activators, which are helpful in cancer treatment. For example, Streptomyces verticillus synthesized the anticancer drug bleomycin17, which was FDA-approved in 1973, and has therapeutic potential in cervical carcinoma, lymphoma, head and neck cancer, and testicular cancer18.

In order to find inhibitors with a strong binding affinity and ideal ADMET (Absorption, Distribution, Metabolism, Excretion, and Toxicity) profile, molecular docking and dynamics simulation studies have recently been applied to develop new drugs19,20,21,22,23. Finding new potential EGFRCSTMLR tyrosine kinase inhibitors from the StreptomeDB database is the goal of the current investigation.

Methodology

Retrieval and processing of ligand database

StreptomeDB is a curated database of Isolated or mutasynthesized Natural Products (NPs) from Streptomyces24. The StreptomeDB 3.0 is an updated database with 6524 NPs from 3302 actinomyces species25. The three-dimensional structures of these 6524 NPs were downloaded in SDF file format from http://www.pharmbioinf.uni-freiburg.de/streptomedb. All the NPs and reference molecules were converted into mae file format for easy processing by utilizing the LigPrep module of Schrödinger (Schrödinger, LLC, New York, USA, 2022). LigPrep module windows were set to OPLSe (Optimized Kanhesia for Liquid Simulations) forcefield, ion neutralization, and single isomers generation per ligand26. Reference structure CH7 (CH7233163) and co-crystal ligand Osimertinib of the target were used as the standard reference for the study27,28.

Retrieval and processing of target

A 3-D crystal structure of EGFRCSTMLR was downloaded from the RCSB-PDB (RCSB Protein Data Bank) with PDB id: 6LUD27. The receptor was already co-crystallized with a well-known anti-Tyrosine Kinase inhibitor: Osimertinib. The structure was then imported into the Protein preparation wizard of Schrödinger to rectify the problems in protein structure. The protein was further refined by optimizing hydrogen bonds and minimizing its structural energy with the OPLSe force field26,29. The receptor grid generation module in Glide-Schrödinger creates a grid box over the ligand with a Van-der-Waal scaling factor set to 1.0, and a partial charge cutoff is set to 0.25.

ADME analysis of ligand library

ADME is a test that establishes the absorption, distribution, metabolism, and excretion characteristics of a potential principal compound and helps study its disposition within an organism30. ADME was performed by incorporating the QikProp program of Schrödinger running in standard default settings (Schrödinger, LLC, New York, USA, 2022). This program helps generate the relevant descriptor (set of 46 molecular descriptors) and uses selected descriptors to carry out ADMET predictions. Drug-likeness of a lead depends on various parameters, and they were used to estimate the pharmacokinetic analysis of the StreptomeDB ligand library. Selected 26 descriptors that play a significant role in most drug-likeliness criteria were chosen on the ligand library, making them match 95% of known drug ranges.

Molecular docking simulations

Molecular docking was performed using Glide software from the Schrödinger suite with Standard-Precision protocol (SP) with ten poses included in post-docking minimization keeping all other parameters default31. Re-docking, the co-crystallized ligand with a suitable RMSD (Root Mean Square Deviation) value (less than 2 Å), enabled to verify the docking procedure. The re-docked poses gave away similar interactions compared to the active site with the native pose. The above experiment leads to the inference that the docking protocol is reliable and sustainable for predicting plausible inhibitors in the used StreptomeDB database32. For further evaluation, the compounds were selected based on the obtained glide score, binding energy, and analysis of the binding pattern of the ligand–protein complex. The PAINS (pan assay interference compounds) was also applied to the selected compound using the Knime Analytics Platform33 with PAINS, PAINS-A, PAINS-B, and PAINS-C filter criteria of RDKit Molecule catalog filter node.

Molecular dynamics simulations

The enzyme–inhibitor complexes that had been docked were considered as the initial framework for future simulation runs. The Schrödinger Suite's Desmond simulation tool was used to run the MD (Molecular Dynamics) simulations (Schrödinger, LLC, New York, USA, 2022). The orthorhombic box was solvated using the TIP3P (transferable intermolecular potential with 3 points) solvation/aqueous system, and the solvent buffer was expanded around the protein by 10 Å using the OPLSe force field34. The system was then neutralized and maintained at 0.15 M salt concentration by placing the required NaCl counter ions. Before employing the simulation study, Desmond's pre-defined system relaxation algorithm was used to equilibrate the system, consisting of six phases in the protocol. The first two were set for the NVT simulation at 10 K temperature, where the first phase was carried out for 100 ps (picosecond) using the Brownian dynamics method to restrict the motion of heavy atoms present in the protein molecules, while in the second stage, the restriction was applied on the large atom of solute for next 12 ps. After that, the NPT condition was applied in the third step by restricting the atom of solute again at 10 K for 12 ps, and the constant pressure was almost one bar. The temperature was likely to be enhanced up to 310 K in the fourth stage with 10 ps time at NPT condition for relaxing the protein molecules and possessing no restriction on the system. The final simulation time was 1000 ns (nanosecond), with a 1 ns recording interval. All simulations used the OPLS4 force field parameters. To calculate the MM/GBSA (Molecular Mechanics and Generalized Bonn-Surface Area) representing the following equation, the prime module of Schrödinger was used at the last 100 ns of the MD trajectories:

$$\Delta {G}_{bind}={\Delta G}_{comp} - [{\Delta G}_{protein}+\Delta {G}_{lig}]$$

where ΔGbind represents the free binding energy, ΔGprotein is the binding energy of the target protein, and ΔGlig represents the binding energy of the ligand in the protein–ligand complex35.

Results and discussion

ADMET properties predictions

It has been suggested that the compounds used as human therapeutic agents should have drug-likeliness properties or exhibit good ADMET profiles36. The ADME properties of 6524 compounds were ascertained using the QikProp web service against selected 25 molecular descriptors criteria. QikProp program predicted that the potential to penetrate the blood/brain barrier and bind with human serum albumin was portrayed by 71.8% and 85% virtual hits. Ro5 (Lipinski's rule of five) for screening drug-like compounds showed that 76% of molecules retained physicochemical properties and followed Lipinski's parameters. In addition, the Ro3 (Ghose rule of three) was applied to enhance the ADME prediction and observed that 67% of molecules were passed through it.

Further, the potassium ion (K+) channel linked to the Human ether-a-go-go-related gene (HERG) is responsible for long QT syndrome37. This channel plays a vital role in regulating heartbeats, making it a vulnerable target to the potential therapeutic drug for cardiac toxicity38. So, the value of predicted IC50 was calculated for the inhibitory action on HERG K+ channel by drugs to model drug toxicity effect, and it was found that both the standard drug lies far outside the safe threshold. 1708 (26.2%) compounds depicted a remarkable result by effectively passing through all the ADMET parameters and were selected for anticancer activity analysis by being subjected to the molecular docking study with the target protein EGFR-tyrosine kinase. The ADME properties of all 6524 compounds are presented in Fig. 1. All the compounds passed through the above criteria lie in the zero-violation category of Ro3 and Ro5 (except C_5024 with only one violation in Ro3), whereas both Osimertinib and CH7 have one violation in Ro3, 1 and 3 violations in Ro5, respectively.

Figure 1
figure 1

ADMET profile of 6524 compounds from StreptomeDB 3.0 database.

Full size image

Molecular docking study of triple mutant EGFRCSTMLR

Molecular docking of 2124 compounds at the active site of mutant EGFRCSTMLR leads to the identification of a novel potential drug for anticancer. Table 1 shows the hit compound's binding score, hydrogen, and hydrophobic bond-forming amino-acid residues. Out of 2124 compounds, five were screened out, which are C_42: Steptocarbazole A, C_4299: Limazepine H, C_4300: 3'N-Formylholyrine A, C_5024: CHEMBL1159781, and C_5702: 3′-N-acetylholyrine A, revealing a better binding score than reference inhibitor CH7: CH7233163 and co-crystallized ligand Osi: Osimertinib (well-known Tyrosine inhibitor). Duration of the interaction39. It was found that these selected molecules have also passed all criteria of the PAINS filter as well.

Table 1 Binding energy and molecular interaction analysis (Hydrogen and Hydrophobic interaction) ligands with triple mutant EGFRL858R/T790M/C797S and their production source.
Full size table

The analysis of triple mutant EGFRCSTMLR protein with interacting ligands showed that MET_793 amino acid residue forms a hydrogen bond with ligands in all five hits ligands–EGFR complex like co-crystallized ligand Osimertinib (Fig. 2). Thus, MET_793 interaction can be referred to as "vital" as it was conserved in all protein–ligands complexes40,41. In addition, all the protein–ligand complexes also shared the common hydrophobic bond-forming amino acid residues VAL_726, ALA_743, MET_790, LEU_792, PRO_794, and LEU_844.

Figure 2
figure 2

2D plot showing the network of the interaction of all compounds representing a: C_42, b: C_4299, c: C_4300, d: CH7, e: C_5024, f: C_5702 and g: Osi with active site residues of EGFR (L858R/T790M/C797S) mutant during molecular docking study.

Full size image

Out of five screened compounds from StreptomeDB 3.0 database, ligand C42 exhibited the best docking score of − 9.3 kcal/mol.

It is a derivative of staurosporine containing carbazole ring, which made hydrogen bond with MET_793 and Gln_791 in the connecting region of C and N terminal of protein. It formed another hydrogen bond via the hydroxymethyl group of the pyrrole ring with ASP_800. The drug C_42 (Steptocarbazole A or Streptocarbazole A) was investigated for cytotoxic activity by Fu et al.42 and found effective against the HL-60, A-549, P388, and Hela cell line with IC50 values of 1.4, 5.0, 18.9, and 34.5 μM, respectively. Additionally, its inhibitory effects on the cell cycle and protein kinase were also reported. The researcher also found that Steptocarbazole A stopped the cell cycle of the Hela cell line in the G2/M phase at 10 μM concentration42.

Ligands C4300 and C5702 have the same binding affinities and share a common bond interaction with the EGFR receptor because of their structural similarity. Only the presence of the formyl group in C4300 and the acetyl group in C5702 distinguishes them. Both comprise a carbazole ring bonded to MET_793 and an oxane ring that interacted with the SER_797 and ASP_800 amino acid residues of the EGFR receptor through its hydroxyl group.

Carbazole moiety in C_42, C_4300, and C_5702 is a large-size, tricyclic, strong pharmacophoric moiety that provides the best binding affinity with the receptor, which leads to the downregulation of protein kinase and eventually apoptosis of cancerous cell43.

Zhou et al.44 investigated the cytotoxic activities of C_4300 (3′-N-formylholyrine A) on the growth of PC-3 and SW-620 cell lines with the cytotoxic activity of IC50 values 2.50 and 0.73 μM, respectively.

Xiao et al.45 found C_5702 (3′-N-acetylholyrine A) compound to have cytotoxicity activity against the tumor cell lines HCT-116, K562, and Huh 7.5 as well as the normal hepatic cell line LO2 with IC50 values as > 100, 32.4, 35.2, and 41.1 μM, respectively. Wang et al.46 also tested the same compound for inhibitory action against protein kinase C enzymes θ (PKC θ) and displayed inhibitory activity with an IC50 value less than 2.5 μM.

C5024 showed good interaction with protein by forming the maximum number of hydrogen bonds (four). It contains a phosphate group that interacted with SER_797 and ASP_800 residues in the ATP binding pocket, while its unsaturated lactone group made interactions in the hinge region with MET_793, which helps in forming a strong bond with protein.

Its anticancer activity is accomplished by the presence of unsaturated lactone and phosphate ester moieties47. It is an analog of the anticancer drug Fostriecin, which is widely known for inhibiting protein phosphatase48,49.

This implies that the reported compound C5024 has both phosphatase and kinase inhibitory properties, suggesting dual action. The anti-tumor properties of C5024 (CHEMBL1159781) have been observed by McCluskey et al.49 on L1210 and HCT-8 Cell Lines, with GI50 concentrations of 0.22 and 1.5 μg/mL, respectively.

Ligand C4299 possesses amide moieties which formed two hydrogen bonds with Glu_762 and Lys_745 amino acids and improved the binding affinity. Its benzodiazepine ring is oriented in the back of the ATP binding pocket, where it interacts with MET_793 via oxo and hydroxy group, which act as donner and acceptor and lies in the hinge region of protein kinase.

Molecular docking analysis provides only a tentative result regarding protein–ligand interaction and binding score, so before concluding, molecular dynamic simulation needed to be carried out to evaluate other factors like protein stability and compactness, solvent effects, and time.

Molecular dynamics studies of the mutant protein

MD simulation was performed for ligands (C_42, C_4299, C_4300, 5024, and C_5702), reference drug CH7233163 and co-crystallized ligand Osimertinib to confirm their structural stability at the active site of triple mutant protein EGFRCSTMLR during 1000 ns MD.

RMSD analysis of the EGFRCSTMLR backbone atom confirmed that all the complexes revealed good stability (less RMSD) compared to reference inhibitor CH7233163, as shown in Fig. 3. In addition, the RMSD value for ligands C_42, C_4299, and C_5024 was comparable to co-crystallized ligand Osimertinib, while ligands C_4300 and C_5702 exhibited a lower value than Osimertinib. The average RMSD value for triple mutant EGFR + C_42, EGFR + C_4299, EGFR + C_4300, EGFR + C_5024, EGFR + C_5702, EGFR + CH7233163, EGFR + Osimertinib was 3.25 Å, 3.57 Å, 2.18 Å, 3.22 Å, 2.88 Å, 4.09 Å, 3.07 Å, respectively.

Figure 3
figure 3

RMSD profile of triple mutant EGFR (L858R/T790M/C797S) protein complex with ligands during MD simulation of 1000 ns.

Full size image

Root mean square fluctuations (RMSF) denotes the local changes in the protein chain around the ligand. The RMSF value of all the interacting residues during the simulation is presented in Fig. 4. All the complexes revealed a comparable RMSF value with the co-crystallized ligand Osimertinib. It was found that the amino acid surrounding the ligand, GLN_791, MET_793, SER_797, and ASP_800, showed significantly less RMS fluctuations and remained stable during interactions. The average RMSF value of these amino acid residues, GLN_791, MET_793, SER_797, and ASP_800 was 0.66 Å, 0.68 Å, 0.74 Å, and 0.93 Å, respectively. The average RMSF value of ligand C_42 (1.54 Å) was very similar to the co-crystallized ligand (1.54 Å) and smaller than reference inhibitor CH7233163 (1.71 Å). It proved that protein secondary structures like alpha-helix and beta sheets remain more rigid during simulation.

Figure 4
figure 4

RMSF profile of triple mutant EGFR (L858R/T790M/C797S) protein’s interacting residues with ligands during MD simulation of 1000 ns.

Full size image

The Rg (Radius of Gyration) of ligands tells the extendedness of ligands during simulation. It was assessed that ligand C_42 comprised the lowest radius of gyration value (3.87 Å) among all ligands, co-crystallized ligand, and reference drug, as shown in Fig. 5, which illustrate that C_42 remained folded and rigid state during simulation because of its aromatic structure. Rg of ligand C_4299 (4.12 Å), C_4300 (4.17 Å), C_5024 (4.69 Å), and C_5702 (4.36 Å) was very comparable to co-crystallized inhibitor Osimertinib (4.49 Å). It was also observed that all ligands exhibited less Rg value than the reference drug CH7233163 (5.88 Å) (Fig. 5).

Figure 5
figure 5

The radius of gyration of ligands during MD simulation of triple mutant EGFR(L858R/T790M/C797S) protein–receptor complex.

Full size image

Solvent Accessible Surface Area (SASA) evaluates the part of ligand surface area accessible by its surrounding water molecules or solvent. The plot of SASA of all ligands during MD simulation is presented in Fig. 6. The average SASA value of ligands C_42 was the lowest (82.9 Å2), indicating that only a tiny part of the ligand is accessible to water. Ligand C_4299 (105.79 Å2), C_4300 (153.44 Å2), and C_5702 (123.0 Å2) comprised less SASA value than reference drug CH7233163 (222.38 Å2) and co-crystallized ligand Osimertinib (178.20 Å2), while ligand C_5024 had higher value (208.60 Å2) but were comparable to the reference drug (Fig. 6).

Figure 6
figure 6

Solvent accessible surface area (SASA) of five-hit ligands, reference drug CH7233163 and co-crystallized ligand Osimertinib during MD simulation of 1000 ns.

Full size image

The number of hydrogen bonds formed by amino acid residues with all ligands during the simulation was assessed. On average, ligands C_42, C_4299, C_4300, and C_5702 formed six hydrogen bonds till 680 ns; after that, the value was enhanced up to eleven. However, ligand C_5024 shows an average of five hydrogen bonds throughout the simulation, similar to reference drug CH7233163 and co-crystallized ligand Osimertinib, as presented in Fig. 7.

Figure 7
figure 7

Hydrogen bonds formed between protein EGFRL858R/T790M/C797S and ligands during MD simulation of 1000 ns.

Full size image

MD simulation results also revealed that all the ligands form a hydrogen bond with protein residue MET_793 like reference ligand CH7233163 and co-crystallized ligand Osimertinib with 100% occupancy; in addition, ligand C_4299 remain interacted with MET_793 for the maximum period (Fig. 8).

Figure 8
figure 8

Amino acid residues contacted with ligands (A: C_42, B: C_4299, C: C_4300, D: C_5024, E: C_5702, F: CH7 and G: OSI) through molecular interaction (hydrogen-bond, hydrophobic, ionic, and water bridges) during MD simulation of 1000 ns.

Full size image

Ligand C_42 shows hydrophobic interaction with LEU_718, PHE_723, VAL_726, ALA_743, MET_790, and LEU_844 residues, which was also seen in co-crystallized ligand Osimertinib. C_42 also showed a water bridge with a SER_797, like in the co-crystallized ligand (Fig. 8).

Ligand C_42 makes hydrogen bonds with MET_793 and Gln791 with a 100% interaction percentage, which existed in reference ligand as well with 96% and 98% interaction probability, as presented in Fig. 9. Co-crystallized ligand Osimertinib interacted with MET_793 and SER_797 residue with 99% and 55% interaction percentage, respectively (Fig. 9). Ligands C_5702, C_4299, C_4300, and C_5024 also show strong interaction with MET_793, and the interaction percentage was 97% to 100%. They also interact with ASP_800 with 34%, 35%, 57%, and 98% probability, respectively. Water-mediated interaction was also revealed in ligand C_4299 with ASP_855 and LEU_718, Ligand C_5024 with MET_793, and Ligand 5720 with Ser720 (Fig. 9).

Figure 9
figure 9

2D plot represents the ligands representing a: C_42, b: C_4299, c: C_4300, d: CH7, e: C_5024, f: C_5702 and g: Osi with interacting residues of EGFRL858R/T790M/C797S and their interaction percentage during MD simulation of 1000 ns.

Full size image

Thus, based on molecular docking and molecular dynamic simulation results, the screened ligand C_42 exhibited the best interaction with triple mutant EGFRCSTMLR protein. It comprised better binding energy, the Rg, and SASA value than reference ligand CH7233163 and co-crystallized ligand Osimertinib. Its RMSD and RMSF values were comparable to the co-crystallized ligand and lower than the reference drug. In the beginning, it displayed four hydrogen bonds, but towards the end, it had increased to eleven. On the other hand, the co-crystallized ligand and reference drug revealed four hydrogen bonds throughout the simulation (Fig. 7).

Additionally, compared to the co-crystalized ligand, it also demonstrated conserved hydrogen bond interactions with the amino acids MET_793; water bridges with SER_797 residue; and hydrophobic interaction with LEU_718, PHE_723, VAL_726, ALA_743, MET_790, and LEU_844 residues. Thus, the overall analysis proved that ligand C_42 comprises better interaction and can be considered as an anticancer drug against triple mutant EGFRL858R/T790M/C797S.

MM/GBSA calculations and per-residue free energy decomposition

With the MMGBSA technique, the free-energies of protein–ligand complexes were estimated. One hundred frames were collected from the last 100 ns of MD simulation trajectories to explore the binding capabilities of ligands with the EGFRL858R/T790M/C797S. The obtained results represent that the free binding energy of C_42 (− 60.06 kcal mol−1), C_5702 (− 61.23 kcal mol−1), C_5024 (− 51.24 kcal mol−1), and C_4300 (− 59.85 kcal mol−1) are better than reference drug Osimertinib (− 50.60 kcal mol−1) and lower but comparable to CH7 (− 68.74 kcal mol−1), C4299 (− 39.16 kcal mol−1) found relatively lower to CH7 but comparable to Osimertinib. Later on, looking at the components of free energy, it was observed that coulomb energies and electrostatic energies of inhibitors are better than CH7 drugs with a remarkable difference. In contrast, electrostatic energy values of ligands found akin to Osimertinib, also lipophilic energy contribution stays equivalent with reference drugs (Table 2). Van-der-wall energy input was found to be lower than CH7 but nearly the same as Osimertinib, but in the calculation to ΔGbind, coulomb and lipo-energy contribution ratio are higher in all cases than native ligands (Table 2). Therefore, it concludes that selected ligands have a great affinity toward the active site of mutant EGFR protein.

Table 2 Free energies (kcal/mol) of complex EGFRL858R/T790M/C797S with C_42, C_4299, C_4300. C_5024, C_5702, CH7 and Osimertinib.
Full size table

The per-residue decomposition was also performed on contributing amino acids to interpretative the energy influence profile in protein–ligand interaction (Fig. 10). The residues having free energy equal to -1 kcal mol−1 or lesser were taken into account to know the depth of involvement in the protein–ligand complex stability, and it found that the residues LEU_718, PHE_723, VAL_726, ALA_743, LYS_745, MET_790, LEU_792, MET_793, GLY_796, SER_797, LEU_844, and THR_854 have played an essential role in the binding of all the seven ligands. From the analysis of Fig. 10, C_5702, C_4300, and CH7 have shown making interaction with the highest number of residues, followed by C_42, then Osimertinib, and last, C_4299 and C_5024, where all the ligands, including CH7, have shared most of the residues (8–12 residues), the ligand Osimertinib found to be gone a bit out of the league with only two residues (Fig. 10).

Figure 10
figure 10

Per-residue-energy-decomposition analysis of all ligands against EGFRL858R/T790M/C797S structure.

Full size image

Conclusions

In the current research, we have attempted to identify a novel drug against triple mutant EGFRCSTMLR as cancer therapy from streptomyces sources. In order to achieve this, the ADMET profile of all the streptomyces-derived compounds from the StreptomeDB database was assessed, and the screened hits were subjected to computational study. The docking result revealed five compounds, C_42, C_4299, C_4300, C_5024, and C_5702, with better binding affinity than both reference drug CH7233163 and co-crystallized ligand Osimertinib. As these drugs are derived from Actinomycetes so, their upscaled production can easily be achieved by the fermentation process. Further, the drugs shortlisted here were already validated experimentally by other researchers against several cancerous cell lines and protein kinases42,44,45,46,49. Based on the simulation study, ligand C_42 binds at the active site of EGFRCSTMLR with outstanding protein stability, compactness, and solvent effects, revealing excellent inhibitory potential. While Osimertinib, a well-known anti-Tyrosine inhibitor, is chemically manufactured and has several adverse effects-developed resistance. On the other hand, C_42 derived from natural sources can be considered a potential inhibitor of triple mutant EGFR in cancer treatment and an alternative to chemically synthesized.