Understanding the mechanism of atovaquone drug resistance in Plasmodium falciparum cytochrome b mutation Y268S using computational methods

Computational exploration of compounds in Xylocarpus granatum as a potential inhibitor of Plasmodium berghei using docking, molecular dynamics, and DFT studies

Malaria remains a global health concern, with the emergence of resistance to the antimalarial drug atovaquone through cytochrome b (cyt b) being well-documented. This study was prompted by the presence of this mutation in cyt b to enable new drug candidates capable of overcoming drug resistance. Our objective was to identify potential drug candidates from compounds of Xylocarpus granatum by computationally assessing their interactions with Plasmodium berghei cyt b. Using computational methods, we modeled cyt b (GenBank: AF146076.1), identified the binding cavity, and analyzed the Ramachandran plot against cyt b. Additionally, we conducted drug-likeness and absorption, distribution, metabolism, excretion, and toxicity (ADMET) studies, along with density functional theory (DFT) analysis of the compounds. Molecular docking and molecular dynamics simulation (MDS) were used to evaluate the binding energy and stability of the cyt b-ligand complex. Notably, our investigation highlighted kaempferol as a promising compound due to its high binding energy of 7.67 kcal/mol among all X. granatum compounds, coupled with favorable pharmacological properties (ADMET) and antiprotozoal properties at Pa 0.345 > Pi 0.009 (PASS value). DFT analysis showed that kaempferol has an energy gap of 4.514 eV. MDS indicated that all tested ligands caused changes in bonding and affected the structural conformation of cyt b, as observed before MDS (0 ns) and after MDS (100 ns). The most notable differences were observed in the types of hydrogen bonds between 0 and 100 ns. Nevertheles, MDS results from a 100 ns simulation revealed consistent behavior for kaempferol across various parameters including root mean square deviation (RMSD), root mean square fluctuation (RMSF), radius of gyration (Rg), solvent-accessible surface area (SASA), molecular mechanics-Poisson Boltzmann surface area (MM-PBSA), and hydrogen bonds. The cyt b-kaempferol complex demonstrated favorable energy stability, as supported by the internal energy distribution values observed in principal component analysis (PCA), which closely resembled those of the atovaquone control. Additionally, trajectory stability analysis indicated structural stability, with a cumulative eigenvalue of 24.7 %. Dynamic cross-correlation matrix (DCCM) analysis revealed a positive correlation among catalytic cytochrome residues within the amino acid residues range 119-268. The results of our research indicate that the structure of kaempferol holds promise as a potential candidate against Plasmodium.

Systematic review on buparvaquone resistance associated with non-synonymous mutation in drug binding genes site of Theileria annulate

Theileria annulata (T. annulata) is intra-erythrocytic protozoan parasite which is more prevalent in tropical and sub-tropical countries. It has a significant economic impact on the productivity of the dairy industry, and buparvaquone is used to treat infected animals in the prevalent regions of the world. Systematically, buparvaquone targets the cyto-b gene to break the electron transport chain (ETC) and Theileria annulata peptidyl-prolyl isomerase 1 (TaPIN1) gene to destabilize transcription factor JUN (c-JUN) to inhibit proliferation of infected cells, which ultimately leads to the death of T. annulata. The reported studies on drug resistance is due to inappropriate drug application, evolutionary characteristics of the cytochrome b (cyto-b) gene and oncogenic signaling pathways gene (TaPIN1) make the parasite resistant against buparvaquone. Hence, this systematic review was designed to find out non-synonymous mutation in genes (cyto-b and TaPIN1) responsible for drug resistance reported from Tunisia, Turkey, Egypt, Sudan, Iran, Pakistan, China and Germany with reference to the T. annulata Ankara strain of cyto-b (accession no. XM_949625.1) and TaPIN1 (accession no. TA18945) wild type genes. Non-synonymous point mutations were found in cyto-b (Q01 at 130-148 and Q02 at 253-262 regions) and TaPIN1 (A53P and A53T) genes. These point mutations are responsible for developing buparvaquone resistance against T. annulata infection. These genes can be used as biomarkers for the identification of drug resistance in any endemic area. To avoid the complication of drug resistance, development of genetically resistant cattle breeds, potent vaccines and anti-theilerial drugs (Trifloxystrobin and anti-cancerous) are currently required to control proliferating economically important T. annulata parasites.

Participation of Oxidative Stress in the Activity of Compounds Isolated from Eleutherine plicata Herb

From Eleutherine plicata, naphthoquinones, isoeleutherine, and eleutherol were isolated, and previous studies have reported the antioxidant activity of these metabolites. The present work evaluated the role of oxidative changes in mice infected with Plasmodium berghei and treated with E. plicata extract, fraction, and isolated compounds, as well as to verify possible oxidative changes induced by these treatments. E. plicata extracts were prepared from powder from the bulbs, which were submitted to maceration with ethanol, yielding the extract (EEEp), which was fractionated under reflux, and the dichloromethane fraction (FDMEp) was submitted for further fractionation, leading to the isolation of isoeleutherine, eleutherine, and eleutherol. The antimalarial activity was examined using the suppressive test, evaluating the following parameters of oxidative stress: trolox equivalent antioxidant capacity (TEAC), thiobarbituric acid reactive substances (TBARS), and reduced glutathione (GSH). Furthermore, the molecular docking of naphthoquinones, eleutherol, eleutherine, and isoeleutherine interactions with antioxidant defense enzymes was investigated, which was favorable for the formation of the receptor-ligand complex, according to the re-rank score values. Eleutherine and isoeleutherine are the ones with the lowest binding energy for catalase (CAT), glutathione reductase (GR), and glutathione peroxidase (GPx1), showing themselves as possible targets of these molecules in the involvement of redox balance. Data from the present study showed that treatments with E. plicata stimulated an increase in antioxidant capacity and a reduction in oxidative stress in mice infected with P. berghei, with naphthoquinones being responsible for reducing oxidative changes and disease severity.

Malaria Molecular Surveillance in the Peruvian Amazon with a Novel Highly Multiplexed Plasmodium falciparum AmpliSeq Assay

Molecular surveillance for malaria has great potential to support national malaria control programs (NMCPs). To bridge the gap between research and implementation, several applications (use cases) have been identified to align research, technology development, and public health efforts. For implementation at NMCPs, there is an urgent need for feasible and cost-effective tools. We designed a new highly multiplexed deep sequencing assay (Pf AmpliSeq), which is compatible with benchtop sequencers, that allows high-accuracy sequencing with higher coverage and lower cost than whole-genome sequencing (WGS), targeting genomic regions of interest. The novelty of the assay is its high number of targets multiplexed into one easy workflow, combining population genetic markers with 13 nearly full-length resistance genes, which is applicable for many different use cases. We provide the first proof of principle for hrp2 and hrp3 deletion detection using amplicon sequencing. Initial sequence data processing can be performed automatically, and subsequent variant analysis requires minimal bioinformatic skills using any tabulated data analysis program. The assay was validated using a retrospective sample collection (n = 254) from the Peruvian Amazon between 2003 and 2018. By combining phenotypic markers and a within-country 28-single-nucleotide-polymorphism (SNP) barcode, we were able to distinguish different lineages with multiple resistance haplotypes (in dhfr, dhps, crt and mdr1) and hrp2 and hrp3 deletions, which have been increasing in recent years. We found no evidence to suggest the emergence of artemisinin (ART) resistance in Peru. These findings indicate a parasite population that is under drug pressure but is susceptible to current antimalarials and demonstrate the added value of a highly multiplexed molecular tool to inform malaria strategies and surveillance systems. IMPORTANCE While the power of next-generation sequencing technologies to inform and guide malaria control programs has become broadly recognized, the integration of genomic data for operational incorporation into malaria surveillance remains a challenge in most countries where malaria is endemic. The main obstacles include limited infrastructure, limited access to high-throughput sequencing facilities, and the need for local capacity to run an in-country analysis of genomes at a large-enough scale to be informative for surveillance. In addition, there is a lack of standardized laboratory protocols and automated analysis pipelines to generate reproducible and timely results useful for relevant stakeholders. With our standardized laboratory and bioinformatic workflow, malaria genetic surveillance data can be readily generated by surveillance researchers and malaria control programs in countries of endemicity, increasing ownership and ensuring timely results for informed decision- and policy-making.

High-Throughput Quantitative Glycomics Enabled by 12-plex Isobaric Multiplex Labeling Reagents for Carbonyl-Containing Compound (SUGAR) Tags

Glycans, which are widely distributed on most proteins and cell surfaces, are a class of important biomolecules playing crucial roles in various biological processes such as immune response and cellular communication. Modern mass spectrometry (MS) coupled with novel chemical probes greatly facilitates routine analysis of glycans. However, the requirement of high-throughput analysis still calls for advanced tools to be developed. Recently, we devised isobaric multiplex reagents for carbonyl-containing compound (SUGAR) tags for 4-plex N-glycan analysis. To further improve the throughput, we utilized the subtle mass differences among different isotopologues and expanded the multiplexing capacity to 12 channels, a 3-fold throughput improvement for the original SUGAR tag design and achieved high-throughput N-glycan analysis in a single LC-MS/MS injection. We then applied 12-plex SUGAR tags to profile the N-glycans in four subtypes of human Immunoglobulin G (IgG) and to investigate the N-glycan changes in the endometrial cancer cells (ECC1) treated with Atovaquone, a quinone antimicrobial medication, and a dihydroorotate dehydrogenase (DHODH) inhibitor. Data are available via ProteomeXchange with the identifier PXD038501.

Parasite powerhouse: A review of the Toxoplasma gondii mitochondrion

Toxoplasma gondii is a member of the apicomplexan phylum, a group of single-celled eukaryotic parasites that cause significant human morbidity and mortality around the world. T. gondii harbors two organelles of endosymbiotic origin: a non-photosynthetic plastid, known as the apicoplast, and a single mitochondrion derived from the ancient engulfment of an α-proteobacterium. Due to excitement surrounding the novelty of the apicoplast, the T. gondii mitochondrion was, to a certain extent, overlooked for about two decades. However, recent work has illustrated that the mitochondrion is an essential hub of apicomplexan-specific biology. Development of novel techniques, such as cryo-electron microscopy, complexome profiling, and next-generation sequencing have led to a renaissance in mitochondrial studies. This review will cover what is currently known about key features of the T. gondii mitochondrion, ranging from its genome to protein import machinery and biochemical pathways. Particular focus will be given to mitochondrial features that diverge significantly from the mammalian host, along with discussion of this important organelle as a drug target.

Correlated Motions and Dynamics in Different Domains of Epidermal Growth Factor Receptor With L858R and T790M Mutations

Non-small cell lung cancer with an activating epidermal growth factor receptor (EGFR) mutation responds well to targeted drugs. In most cases, drug resistance appears after about a year. Several studies have been conducted on the kinase domain of EGFR to understand the drug resistance mechanism. Since EGFR is a multi-domain protein, mutation in the kinase domain may affect the other domains as well. In this study, we examine the complete structure of the multi-domain EGFR protein and its mutants. We performed molecular dynamics simulations for wildtype EGFR, EGFR with L858R mutation, and EGFR with L858R and T790M mutations. We applied normal mode analysis and complex network analysis to extract the correlated motions in the domains of EGFR. The normal modes are used to construct the dynamic cross-correlation map (DCCM). Simulation results show different patterns of correlated motions in each domain of EGFR mutants compared to the wildtype. In Domains 1 and 3 of the extracellular region, a small number of weak positively correlated motions are extracted. Domains 2 and 4 show large numbers of both positive and negative motions. However, the negatively correlated motions are stronger in mutant structures compared to the wildtype. In Domain 7, some residues showed a positive correlation around the main diagonal. We also identified different communities, nodes and crucial residues in the domains of the structures, which can be important for the function of EGFR. Moreover, hydrogen bond analysis is performed for the stability analysis. The mutant structures have fewer hydrogen bonds compared to the wildtype. Overall, these findings are useful for understanding the dynamics and communications in EGFR domains.

Efforts Made to Eliminate Drug-Resistant Malaria and Its Challenges

Since 2000, a good deal of progress has been made in malaria control. However, there is still an unacceptably high burden of the disease and numerous challenges limiting advancement towards its elimination and ultimate eradication. Among the challenges is the antimalarial drug resistance, which has been documented for almost all antimalarial drugs in current use. As a result, the malaria research community is working on the modification of existing treatments as well as the discovery and development of new drugs to counter the resistance challenges. To this effect, many products are in the pipeline and expected to be marketed soon. In addition to drug and vaccine development, mass drug administration (MDA) is under scientific scrutiny as an important strategy for effective utilization of the developed products. This review discusses the challenges related to malaria elimination, ongoing approaches to tackle the impact of drug-resistant malaria, and upcoming antimalarial drugs.

Inhibition Mechanism of Antimalarial Drugs Targeting the Cytochrome bc1 Complex

Plasmodium falciparum (P. falciparum) is the main parasite known to cause malaria in humans. The antimalarial drug atovaquone is known to inhibit the Q_o-site of the cytochrome bc₁ complex of P. falciparum, which ultimately blocks ATP synthesis, leading to cell death. Through the years, mutations of the P. falciparum cytochrome bc₁ complex, causing resistance to atovaquone, have emerged. The present investigation applies molecular dynamics (MD) simulations to study how the specific mutations Y279S and L282V, known to cause atovaquone resistance in malarial parasites, affect the inhibition mechanism of two known inhibitors. Binding free energy estimates were obtained through free energy perturbation calculations but were unable to confidently resolve the effects of mutations due to the great complexity of the binding environment. Meanwhile, basic mechanistic considerations from the MD simulations provide a detailed characterization of inhibitor binding modes and indicate that the Y279S mutation weakens the natural binding of the inhibitors, while no conclusive effect of the L282V mutation could be observed.

Decoding the Molecular Effects of Atovaquone Linked Resistant Mutations on Plasmodium falciparum Cytb-ISP Complex in the Phospholipid Bilayer Membrane

Atovaquone (ATQ) is a drug used to prevent and treat malaria that functions by targeting the Plasmodium falciparum cytochrome b (PfCytb) protein. PfCytb catalyzes the transmembrane electron transfer (ET) pathway which maintains the mitochondrial membrane potential. The ubiquinol substrate binding site of the protein has heme bL, heme bH and iron-sulphur [2FE-2S] cluster cofactors that act as redox centers to aid in ET. Recent studies investigating ATQ resistance mechanisms have shown that point mutations of PfCytb confer resistance. Thus, understanding the resistance mechanisms at the molecular level via computational approaches incorporating phospholipid bilayer would help in the design of new efficacious drugs that are also capable of bypassing parasite resistance. With this knowledge gap, this article seeks to explore the effect of three drug resistant mutations Y268C, Y268N and Y268S on the PfCytb structure and function in the presence and absence of ATQ. To draw reliable conclusions, 350 ns all-atom membrane (POPC:POPE phospholipid bilayer) molecular dynamics (MD) simulations with derived metal parameters for the holo and ATQ-bound -proteins were performed. Thereafter, simulation outputs were analyzed using dynamic residue network (DRN) analysis. Across the triplicate MD runs, hydrophobic interactions, reported to be crucial in protein function were assessed. In both, the presence and absence of ATQ and a loss of key active site residue interactions were observed as a result of mutations. These active site residues included: Met 133, Trp136, Val140, Thr142, Ile258, Val259, Pro260 and Phe264. These changes to residue interactions are likely to destabilize the overall intra-protein residue communication network where the proteins' function could be implicated. Protein dynamics of the ATQ-bound mutant complexes showed that they assumed a different pose to the wild-type, resulting in diminished residue interactions in the mutant proteins. In summary, this study presents insights on the possible effect of the mutations on ATQ drug activity causing resistance and describes accurate MD simulations in the presence of the lipid bilayer prior to conducting inhibitory drug discovery for the PfCytb-iron sulphur protein (Cytb-ISP) complex.

An epigrammatic status of the 'azole'-based antimalarial drugs

The development of multidrug resistance in the malarial parasite has sabotaged majority of the eradication efforts by restraining the inhibition profile of first line as well as second line antimalarial drugs, thus necessitating the development of novel pharmaceutics constructed on appropriate scaffolds with superior potency against the drug-resistant and drug-susceptible Plasmodium parasite. Over the past decades, the infectious malarial parasite has developed resistance against most of the contemporary therapeutics, thus necessitating the rational development of novel approaches principally focused on MDR malaria. This review presents an epigrammatic collation of the epidemiology and the contemporary antimalarial therapeutics based on the 'azole' motif.

Quinoline-Based Hybrid Compounds with Antimalarial Activity

The application of quinoline-based compounds for the treatment of malaria infections is hampered by drug resistance. Drug resistance has led to the combination of quinolines with other classes of antimalarials resulting in enhanced therapeutic outcomes. However, the combination of antimalarials is limited by drug-drug interactions. In order to overcome the aforementioned factors, several researchers have reported hybrid compounds prepared by reacting quinoline-based compounds with other compounds via selected functionalities. This review will focus on the currently reported quinoline-based hybrid compounds and their preclinical studies.

Determination of glutathione redox potential and pH value in subcellular compartments of malaria parasites

The malaria parasite Plasmodium falciparum is exposed to multiple sources of oxidative challenge during its complex life cycle in the Anopheles vector and its human host. In order to further elucidate redox-based parasite host cell interactions and mechanisms of drug action, we targeted the genetically encoded glutathione redox sensor roGFP2 coupled to human glutaredoxin 1 (roGFP2-hGrx1) as well as the ratiometric pH sensor pHluorin to the apicoplast and the mitochondrion of P. falciparum. Using live cell imaging, this allowed for the first time the determination of the pH values of the apicoplast (7.12±0.40) and mitochondrion (7.37±0.09) in the intraerythrocytic asexual stages of the parasite. Based on the roGFP2-hGrx1 signals, glutathione-dependent redox potentials of -267mV and -328mV, respectively, were obtained. Employing these novel tools, initial studies on the effects of redox-active agents and clinically employed antimalarial drugs were carried out on both organelles.

Structural dynamics of Casein Kinase I (CKI) from malarial parasite Plasmodium falciparum (Isolate 3D7): Insights from theoretical modelling and molecular simulations

The protein kinases (PKs), belonging to serine/threonine kinase (STKs), are important drug targets for a wide spectrum of diseases in human. Among protein kinases, the Casein Kinases (CKs) are vastly expanded in various organisms, where, the malarial parasite Plasmodium falciparum possesses a single member i.e., PfCKI, which can phosphorylate various proteins in parasite extracts in vitro condition. But, the structure-function relationship of PfCKI and dynamics of ATP binding is yet to be understood. Henceforth, an attempt was made to study the dynamics, stability, and ATP binding mechanisms of PfCKI through computational modelling, docking, molecular dynamics (MD) simulations, and MM/PBSA binding free energy estimation. Bi-lobed catalytic domain of PfCKI shares a high degree of secondary structure topology with CKI domains of rice, human, and mouse indicating co-evolution of these kinases. Molecular docking study revealed that ATP binds to the active site where the glycine-rich ATP-binding motif (G16-X-G18-X-X-G21) along with few conserved residues plays a crucial role maintaining stability of the complex. Structural superposition of PfCKI with close structural homologs depicted that the location and length of important loops are different, indicating the dynamic properties of these loops among CKIs, which is consistent with principal component analysis (PCA). PCA displayed that the overall global motion of ATP-bound form is comparatively higher than that of apo form. The present study provides insights into the structural features of PfCKI, which could contribute towards further understanding of related protein structures, dynamics of catalysis and phosphorylation mechanism in these important STKs from malarial parasite in near future.

Proteomics in India: the clinical aspect

Proteomics has emerged as a highly promising bioanalytical technique in various aspects of applied biological research. In Indian academia, proteomics research has grown remarkably over the last decade. It is being extensively used for both basic as well as translation research in the areas of infectious and immune disorders, reproductive disorders, cardiovascular diseases, diabetes, eye disorders, human cancers and hematological disorders. Recently, some seminal works on clinical proteomics have been reported from several laboratories across India. This review aims to shed light on the increasing use of proteomics in India in a variety of biological conditions. It also highlights that India has the expertise and infrastructure needed for pursuing proteomics research in the country and to participate in global initiatives. Research in clinical proteomics is gradually picking up pace in India and its future seems very bright.

Withanolide A offers neuroprotection, ameliorates stress resistance and prolongs the life expectancy of Caenorhabditis elegans

Withanolide A (steroidal lactone) forms the major constituent of the most popular herbal drug in Ayurvedic medicine, Ashwagandha. It has been used since ancient times as an alternative medicine for the treatment of a variety of age related disorders. Here we provide multiple lines of evidence indicating that Withanolide A improves healthspan, delays age-associated physiological changes and also extends the lifespan of Caenorhabditis elegans. We also report several neuroprotective benefits of this natural product, including its anti-amyloidogenic effects, alleviation of α-synuclein aggregation and neuroprotection through modulation of neural mediators like acetylcholine. We observed that Withanolide A mediates lifespan extension and promotes stress resistance via insulin/insulin-like growth factor signaling pathway. Such findings could be helpful to develop a therapeutic medicine from this natural product for the prevention or reversal of age-related ailments and to improve the survival of patients suffering from Alzheimer's or Parkinson's disease.

In silico study of subtilisin-like protease 1 (SUB1) from different Plasmodium species in complex with peptidyl-difluorostatones and characterization of potent pan-SUB1 inhibitors

•

Homology models of four SUB1 orthologues from P. falciparum species were produced.

•

We analyzed the binding mode of our previous difluorostatone inhibitors to six SUB1.

•

In vitro activity of our difluorostatone-based inhibitors was correctly predicted.

•

We derived a structure-based pan-SUB1 pharmacophore, and validated it in silico.

•

We confirmed that development of pan-SUB1 inhibitors is a feasible task.

Plasmodium falciparum subtilisin-like protease 1 (SUB1) is a novel target for the development of innovative antimalarials. We recently described the first potent difluorostatone-based inhibitors of the enzyme ((4S)-(N-((N-acetyl-l-lysyl)-l-isoleucyl-l-threonyl-l-alanyl)-2,2-difluoro-3-oxo-4-aminopentanoyl)glycine (1) and (4S)-(N-((N-acetyl-l-isoleucyl)-l-threonyl-l-alanylamino)-2,2-difluoro-3-oxo-4-aminopentanoyl)glycine (2)). As a continuation of our efforts towards the definition of the molecular determinants of enzyme-inhibitor interaction, we herein propose the first comprehensive computational investigation of the SUB1 catalytic core from six different Plasmodium species, using homology modeling and molecular docking approaches. Investigation of the differences in the binding sites as well as the interactions of our inhibitors 1,2 with all SUB1 orthologues, allowed us to highlight the structurally relevant regions of the enzyme that could be targeted for developing pan-SUB1 inhibitors. According to our in silico predictions, compounds 1,2 have been demonstrated to be potent inhibitors of SUB1 from all three major clinically relevant Plasmodium species (P. falciparum, P. vivax, and P. knowlesi). We next derived multiple structure-based pharmacophore models that were combined in an inclusive pan-SUB1 pharmacophore (SUB1-PHA). This latter was validated by applying in silico methods, showing that it may be useful for the future development of potent antimalarial agents.

The Stapled AKAP Disruptor Peptide STAD-2 Displays Antimalarial Activity through a PKA-Independent Mechanism

Drug resistance poses a significant threat to ongoing malaria control efforts. Coupled with lack of a malaria vaccine, there is an urgent need for the development of new antimalarials with novel mechanisms of action and low susceptibility to parasite drug resistance. Protein Kinase A (PKA) has been implicated as a critical regulator of pathogenesis in malaria. Therefore, we sought to investigate the effects of disrupted PKA signaling as a possible strategy for inhibition of parasite replication. Host PKA activity is partly regulated by a class of proteins called A Kinase Anchoring Proteins (AKAPs), and interaction between HsPKA and AKAP can be inhibited by the stapled peptide Stapled AKAP Disruptor 2 (STAD-2). STAD-2 was tested for permeability to and activity against Plasmodium falciparum blood stage parasites in vitro. The compound was selectively permeable only to infected red blood cells (iRBC) and demonstrated rapid antiplasmodial activity, possibly via iRBC lysis (IC₅₀ ≈ 1 μM). STAD-2 localized within the parasite almost immediately post-treatment but showed no evidence of direct association with PKA, indicating that STAD-2 acts via a PKA-independent mechanism. Furosemide-insensitive parasite permeability pathways in the iRBC were largely responsible for uptake of STAD-2. Further, peptide import was highly specific to STAD-2 as evidenced by low permeability of control stapled peptides. Selective uptake and antiplasmodial activity of STAD-2 provides important groundwork for the development of stapled peptides as potential antimalarials. Such peptides may also offer an alternative strategy for studying protein-protein interactions critical to parasite development and pathogenesis.

Contribution of EGFR and ErbB-3 Heterodimerization to the EGFR Mutation-Induced Gefitinib- and Erlotinib-Resistance in Non-Small-Cell Lung Carcinoma Treatments

EGFR mutation-induced drug resistance has become a major threat to the treatment of non-small-cell lung carcinoma. Essentially, the resistance mechanism involves modifications of the intracellular signaling pathways. In our work, we separately investigated the EGFR and ErbB-3 heterodimerization, regarded as the origin of intracellular signaling pathways. On one hand, we combined the molecular interaction in EGFR heterodimerization with that between the EGFR tyrosine kinase and its inhibitor. For 168 clinical subjects, we characterized their corresponding EGFR mutations using molecular interactions, with three potential dimerization partners (ErbB-2, IGF-1R and c-Met) of EGFR and two of its small molecule inhibitors (gefitinib and erlotinib). Based on molecular dynamics simulations and structural analysis, we modeled these mutant-partner or mutant-inhibitor interactions using binding free energy and its components. As a consequence, the mutant-partner interactions are amplified for mutants L858R and L858R_T790M, compared to the wild type EGFR. Mutant delL747_P753insS represents the largest difference between the mutant-IGF-1R interaction and the mutant-inhibitor interaction, which explains the shorter progression-free survival of an inhibitor to this mutant type. Besides, feature sets including different energy components were constructed, and efficient regression trees were applied to map these features to the progression-free survival of an inhibitor. On the other hand, we comparably examined the interactions between ErbB-3 and its partners (EGFR mutants, IGF-1R, ErbB-2 and c-Met). Compared to others, c-Met shows a remarkably-strong binding with ErbB-3, implying its significant role in regulating ErbB-3 signaling. Moreover, EGFR mutants corresponding to poor clinical outcomes, such as L858R_T790M, possess lower binding affinities with ErbB-3 than c-Met does. This may promote the communication between ErbB-3 and c-Met in these cancer cells. The analysis verified the important contribution of IGF-1R or c-Met in the drug resistance mechanism developed in lung cancer treatments, which may bring many benefits to specialized therapy design and innovative drug discovery.