Structure-based inhibitory peptide design targeting peptide-substrate binding site in EGFR tyrosine kinase

Design and investigation of novel iridoid-based peptide conjugates for targeting EGFR and its mutants L858R and T790M/L858R/C797S: an in silico study

In this work, we designed novel peptide conjugates with plant-based iridoid and lichen-derived depside derivatives to target the wild-type EGFR (WT) and its mutants, L858R and T790M/L858R/C797S triple mutant. These mutations are often expressed in multiple cancers, particularly lung cancer. Specifically, the iridoids included 7-deoxyloganetic acid (7-DGA) and loganic acid (LG), while the depside derivative was sekikaic acid (SK). These compounds are known for their innate anticancer properties and were conjugated with two separate peptide sequences KLPGWSG (K) and YSIPKSS (Y). These sequences have been shown to target EGFR in previous phage display library screening, although the mechanism is unknown. Thus, we created the di-conjugates for dual targeting and investigated their interactions of the di-conjugates and that of the neat peptides with the kinase domain of EGFR (WT) and the two mutants using molecular docking, molecular dynamics (MD) simulations, and MM-GBSA analysis. Docking studies revealed that the (7-DGA)2-K showed the highest binding affinity at - 9.3 kcal/mol with the L858R mutant, while (LG)2-Y displayed the highest binding affinity at - 9.0 kcal/mol for the triple mutant receptor. Our results indicated that several of the conjugates interacted with crucial residues of the kinase domain, including ASP855 and THR854 (activation loop), MET793 and PRO794 (hinge region), ARG841 (catalytic loop), and LYS728 and LEU718 of the glycine-rich P-loop. Interestingly, strong hydrophobic interactions were also observed with the C-terminal tail residues, such as PHE997 and ALA1000 as well as with ARG999 for the YSIPKSS peptide and most of the conjugates. The hydroxyl group of the cyclopentane ring and the oxygen of the pyran ring of the (7-DGA)2-peptide conjugates contributed to binding particularly in the hinge region, while the peptide components formed an extended structure that bound well into the C-lobe. The (SK)2-Y di-conjugate and KLPGWSG peptide formed hydrogen bonds with the SER797 residue of the triple mutant. Overall, our results show that the (7-DGA)2-K, di-conjugate, the (7-DGA)2-Y di-conjugate, and the neat YSIPKSS demonstrated strong and stable binding with the L858R mutant and the highly resistant triple mutant EGFR, respectively. The novel designed conjugates demonstrate potential for further optimization for laboratory studies aimed at developing new therapeutics for targeting specific EGFR mutant expressing cells.

The prospect of substrate-based kinase inhibitors to improve target selectivity and overcome drug resistance

Human kinases are recognized as one of the most important drug targets associated with cancer. There are >80 FDA-approved kinase inhibitors to date, most of which work by inhibiting ATP binding to the kinase. However, the frequent development of single-point mutations within the kinase domain has made overcoming drug resistance a major challenge in drug discovery today. Targeting the substrate site of kinases can offer a more selective and resistance-resilient solution compared to ATP inhibition but has traditionally been challenging. However, emerging technologies for the discovery of drug leads using recombinant display and stabilization of lead compounds have increased interest in targeting the substrate site of kinases. This review discusses recent advances in the substrate-based inhibition of protein kinases and the potential of such approaches for overcoming the emergence of resistance.

A Novel Dual-labeled Peptide for Multimodal Imaging of EGFR with L858R Mutation

BACKGROUND

The development of molecular imaging agents targeting epidermal growth factor receptor (EGFR) with L858R mutation may help with the selection of non-small cell lung carcinoma (NSCLCL) patients who may benefit from EFGR tyrosine kinase inhibitor (TKI) therapy.

OBJECTIVE

In this study, we developed 99mTc STHHYYP-GHEG-ECGK-tetramethylrhodamine (STHHYYP-ECGK-TAMRA) to target EGFR with L858R mutation in NSCLC tumors and verified its probability as a molecular imaging agent.

METHODS

Fmoc solid-phase peptide synthesis was used to synthesize STHHYYP-ECGKTAMRA. 99mTc labelled STHHYYP-ECGK-TAMRA was prepared. Gamma imaging, fluorescent imaging and biodistribution were performed in murine models bearing NCI-H1975 and NCI-H1650 tumors.

RESULTS

The binding affinity value (Kd) of 99mTc STHHYYP-ECGK-TAMRA was estimated to be 130.6 ± 29.2 nM in NCI-H1975 cells. The gamma camera images showed a substantial uptake of 99mTc STHHYYP-ECGK-TAMRA in the NCI-H1975 tumor. The % injected dose/gram of the NCI-H1975 tumor tissue was 2.77 ± 0.70 and 3.48 ± 1.01 at 1 and 3 h, respectively.

CONCLUSION

Specific binding of 99mTc STHHYYP-ECGK-TAMRA to L858R-mutated EGFRpositive NCI-H1975 cells and tumors was demonstrated in in vivo and in vitro studies. The results suggest that 99mTc STHHYYP-ECGK-TAMRA is a good candidate agent for dualmodality imaging targeting EGFR with L858R mutation.

Discovery of novel RARα agonists using pharmacophore-based virtual screening, molecular docking, and molecular dynamics simulation studies

Nuclear retinoic acid receptors (RARs) are ligand-dependent transcription factors involved in various biological processes, such as embryogenesis, cell proliferation, differentiation, reproduction, and apoptosis. These receptors are regulated by retinoids, i.e., retinoic acid (RA) and its analogs, as receptor agonists. RAR agonists are promising therapeutic agents for the treatment of serious dermatological disorders, including some malignant conditions. By inducing apoptosis, they are able to inhibit the proliferation of diverse cancer cell lines. Also, RAR agonists have recently been identified as therapeutic options for some neurodegenerative diseases. These features make retinoids very attractive molecules for medical purposes. Synthetic selective RAR agonists have several advantages over endogenous ones, but they suffer poor pharmacokinetic properties. These compounds are normally lipophilic acids with unfavorable drug-like features such as poor oral bioavailability. Recently, highly selective, potent, and less toxic RAR agonists with proper lipophilicity, thus, good oral bioavailability have been developed for some therapeutic applications. In the present study, ligand and structure-based virtual screening technique was exploited to introduce some novel RARα agonists. Pharmacokinetic assessment was also performed in silico to suggest those compounds which have optimized drug-like features. Finally, two compounds with the best in silico pharmacological features are proposed as lead molecules for future development of RARα agonists.

Novel in-frame duplication variant characterization in late infantile metachromatic leukodystrophy using whole-exome sequencing and molecular dynamics simulation

Metachromatic leukodystrophy (MLD) is a neurodegenerative lysosomal storage disease caused by a deficiency in the arylsulfatase A (ARSA). ARSA deficiency leads to sulfatide accumulation, which involves progressive demyelination. The profound impact of early diagnosis on MLD treatment options necessitates the development of new or updated analysis tools and approaches. In this study, to identify the genetic etiology in a proband from a consanguineous family with MLD presentation and low ARSA activity, we employed Whole-Exome Sequencing (WES) followed by co-segregation analysis using Sanger sequencing. Also, MD simulation was utilized to study how the variant alters the structural behavior and function of the ARSA protein. GROMACS was applied and the data was analyzed by RMSD, RMSF, Rg, SASA, HB, atomic distance, PCA, and FEL. Variant interpretation was done based on the American College of Medical Genetics and Genomics (ACMG) guidelines. WES results showed a novel homozygous insertion mutation, c.109_126dup (p.Asp37_Gly42dup), in the ARSA gene. This variant is located in the first exon of ARSA, fulfilling the criteria of being categorized as likely pathogenic, according to the ACMG guidelines and it was also found to be co-segregating in the family. The MD simulation analysis revealed this mutation influenced the structure and the stabilization of ARSA and led to the protein function impairment. Here, we report a useful application of WES and MD to identify the causes of a neurometabolic disorder.

Recent Developments in Free Energy Calculations for Drug Discovery

The grand challenge in structure-based drug design is achieving accurate prediction of binding free energies. Molecular dynamics (MD) simulations enable modeling of conformational changes critical to the binding process, leading to calculation of thermodynamic quantities involved in estimation of binding affinities. With recent advancements in computing capability and predictive accuracy, MD based virtual screening has progressed from the domain of theoretical attempts to real application in drug development. Approaches including the Molecular Mechanics Poisson Boltzmann Surface Area (MM-PBSA), Linear Interaction Energy (LIE), and alchemical methods have been broadly applied to model molecular recognition for drug discovery and lead optimization. Here we review the varied methodology of these approaches, developments enhancing simulation efficiency and reliability, remaining challenges hindering predictive performance, and applications to problems in the fields of medicine and biochemistry.

Acetogenin Extracted from Annona muricata Prevented the Actions of EGF in PA-1 Ovarian Cancer Cells

BACKGROUND

In individuals with ovarian cancer, an increase in the circulating level of the epidermal growth factor (EGF) is readily apparent. Ovarian cancer cells exhibit signaling pathway of the epidermal growth factor (EGFR) and respond to the EGF. Annona muricata (AM) has been shown to decrease ovarian cell proliferation however, role of AM in regulating EGF actions is not yet to be reported.

OBJECTIVE

In this study, we proposed that the fractionated compound acetogenin can inhibit the activation of EGFR-regulated signaling cascades such as MAPK7 / PI3K-Akt / mTOR / STAT upon EGF stimulation.

METHODS

Ethanolic extract was prepared for the whole AM plant and Thin Layer Chromatography (TLC) was performed to characterize the secondary metabolites and each fraction was assessed using kedde reagent for the presence of acetogenin. The effects of acetogenins were then tested on the survival of PA-1 ovarian cancer cells under basal and EGF stimulated conditions. To delineate the role of acetogenin in EGFR signaling cascades, the in silico docking studies were conducted.

RESULTS

The fraction of acetogenin decreased the viability of EGF induced PA-1 ovarian cancer cells that indicating the EGF inhibitory effects of acetogenin. The docking studies specifically illustrated that when the acetogenin binding with tyrosine kinase (TK) and regulatory unit (RU) which subsequently resulted in a reduction in EGF induced the survival of PA-1 ovarian cancer cells.

DISCUSSION

The vital regulatory role of acetogenin reported in this study indicate significant anticancer activities of acetogenin from AM. The in silico study of the acetogenin function predicted that it binds specifically to Asp837 (phosphor-acceptor site) of EGFR, essential for phosphorylation of substrates in the TK domain and RU which promote downstream signaling.

CONCLUSION

Acetogenin isolated from AM effectively inhibited the survival of PA-1 ovarian cancer cells through impaired EGF signaling.

Rational Design of Peptide-Based Inhibitors Disrupting Protein-Protein Interactions

Protein-protein interactions (PPIs) are well-established as a class of promising drug targets for their implications in a wide range of biological processes. However, drug development toward PPIs is inevitably hampered by their flat and wide interfaces, which generally lack suitable pockets for ligand binding, rendering most PPI systems "undruggable." Here, we summarized drug design strategies for developing peptide-based PPI inhibitors. Importantly, several quintessential examples toward well-established PPI targets such as Bcl-2 family members, p53-MDM2, as well as APC-Asef are presented to illustrate the detailed schemes for peptide-based PPI inhibitor development and optimizations. This review supplies a comprehensive overview of recent progresses in drug discovery targeting PPIs through peptides or peptidomimetics, and will shed light on future therapeutic agent development toward the historically "intractable" PPI systems.

In silico investigation of critical binding pattern in SARS-CoV-2 spike protein with angiotensin-converting enzyme 2

Severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) is a newly-discovered coronavirus and responsible for the spread of coronavirus disease 2019 (COVID-19). SARS-CoV-2 infected millions of people in the world and immediately became a pandemic in March 2020. SARS-CoV-2 belongs to the beta-coronavirus genus of the large family of Coronaviridae. It is now known that its surface spike glycoprotein binds to the angiotensin-converting enzyme-2 (ACE2), which is expressed on the lung epithelial cells, mediates the fusion of the cellular and viral membranes, and facilitates the entry of viral genome to the host cell. Therefore, blocking the virus-cell interaction could be a potential target for the prevention of viral infection. The binding of SARS-CoV-2 to ACE2 is a protein-protein interaction, and so, analyzing the structure of the spike glycoprotein of SARS-CoV-2 and its underlying mechanism to bind the host cell receptor would be useful for the management and treatment of COVID-19. In this study, we performed comparative in silico studies to deeply understand the structural and functional details of the interaction between the spike glycoprotein of SARS-CoV-2 and its cognate cellular receptor ACE2. According to our results, the affinity of the ACE2 receptor for SARS-CoV-2 was higher than SARS-CoV. According to the free energy decomposition of the spike glycoprotein-ACE2 complex, we found critical points in three areas which are responsible for the increased binding affinity of SARS-CoV-2 compared with SARS-CoV. These mutations occurred at the receptor-binding domain of the spike glycoprotein that play an essential role in the increasing the affinity of coronavirus to ACE2. For instance, mutations Pro462Ala and Leu472Phe resulted in the altered binding energy from - 2 kcal mol-1 in SARS-COV to - 6 kcal mol-1 in SARS-COV-2. The results demonstrated that some mutations in the receptor-binding motif could be considered as a hot-point for designing potential drugs to inhibit the interaction between the spike glycoprotein and ACE2.

The functional regulatory details of ERK2 in complex with RSK1: an in silico insight

Protein kinases play a significant role in cellular activation procedures by exhibiting a vivid selection in the target, as well as recognizing and phosphorylating them. Extracellular signal-regulated kinase 2 (ERK2) is one of the main kinases in the mitogen-activated protein kinase (MAPK) signaling cascade and engages in dynamically regulating the activities of signaling proteins and physiological processes, including cell proliferation, differentiation, adhesion, migration, and survival. Predicting collective dynamic and structural motions in biological macromolecules is pivotal to obtain a better understanding of the majority of biological processes. Here, through molecular dynamic simulation and normal mode analysis, we investigated ERK2 conformations, in the forms of active (phosphorylated), inactive (unphosphorylated), and in a complex with its substrate, ribosomal protein S6 kinase alpha-1 (RSK1), to determine functional characteristics. Our finding demonstrated that ERK2 plays a switch role in the regulation of pathways. In the case that this protein kinase is in the active form, all critical regions shift to be prepared to accept the substrate and catalytic action. Meanwhile, inactive ERK2 shows contrasting results in which all motions tend to close the catalytic site and cease the phosphorylation action in the MAPK cascade. These findings are in line with those from other similar studies and provide us with novel molecular target regions and recent details on how this mechanism works.

DFG, αC-helix orientation regarding the active site position and distance between K54 and Glu71 in the active and inactive states of ERK2.

Molecular Targeting of Epidermal Growth Factor Receptor (EGFR) and Vascular Endothelial Growth Factor Receptor (VEGFR)

Epidermal growth factor receptor (EGFR) and vascular endothelial growth factor receptor (VEGFR) are two extensively studied membrane-bound receptor tyrosine kinase proteins that are frequently overexpressed in many cancers. As a result, these receptor families constitute attractive targets for imaging and therapeutic applications in the detection and treatment of cancer. This review explores the dynamic structure and structure-function relationships of these two growth factor receptors and their significance as it relates to theranostics of cancer, followed by some of the common inhibition modalities frequently employed to target EGFR and VEGFR, such as tyrosine kinase inhibitors (TKIs), antibodies, nanobodies, and peptides. A summary of the recent advances in molecular imaging techniques, including positron emission tomography (PET), single-photon emission computerized tomography (SPECT), computed tomography (CT), magnetic resonance imaging (MRI), and optical imaging (OI), and in particular, near-IR fluorescence imaging using tetrapyrrolic-based fluorophores, concludes this review.

The Cytotoxicity of RNase-Derived Peptides

Bacterial ribonuclease binase exhibits a cytotoxic effect on tumor cells possessing certain oncogenes. The aim of this study was to identify the structural parts of the binase molecule that exert cytotoxicity. Out of five designed peptides, the peptides representing the binase regions 21-50 and 74-94 have the highest cytotoxic potential toward human cervical HeLa and breast BT-20 and MCF-7 cancer cells. The peptides B21-50 and B74-94 were not able to enter human lung adenocarcinoma A549 cells, unlike BT-20 cells, explaining their failure to inhibit A549 cell proliferation. The peptide B74-94 shares similarities with epidermal growth factor (EGF), suggesting the peptide's specificity for EGF receptor overexpressed in BT-20 cells. Thus, the binase-derived peptides have the potential of being further developed as tumor-targeting peptides.