1
|
Akher FB, Farrokhzadeh A, Ravenscroft N, Kuttel MM. Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation. J Chem Inf Model 2022; 62:997-1013. [PMID: 35119858 DOI: 10.1021/acs.jcim.1c01399] [Citation(s) in RCA: 3] [Impact Index Per Article: 1.5] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 11/29/2022]
Abstract
Covalent modification of the oncogenic mutant epidermal growth factor receptor (EGFR) by small molecules is an efficient strategy for achieving an enhanced and sustained pharmacological effect in the treatment of non-small-cell lung cancer. NSP-037 (18), an irreversible inhibitor of the L858R/T790M double-mutant EGFR (EGFRDM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFRDM over the wild type (EGFRWT), as compared to clinically approved osimertinib (7). Here, we employ multiple computational approaches to elucidate the mechanism underlining this improved selectivity, as well as the effect of CFA on the selectivity enhancement of inhibitor 18 over 7. We find that EGFRDM undergoes significantly larger conformational changes than EGFRWT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFRDM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFRWT and EGFRDM and 18's higher selectivity than 7. The binding free energy of the 18-bound complexes is -6.38 kcal/mol greater than that of the 7-bound complexes, explaining the difference in selectivity of these inhibitors. Further, free energy decomposition analysis indicates that the electrostatic contribution of key residues plays an important role in the 18-bound complexes. QM/MM calculations show that the most favored mechanism for the Cys797 alkylation reaction is the direct displacement mechanism through a CFA-based inhibitor, producing a reaction with the lowest energy barrier and most stable product.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Department of Computer Science, University of Cape Town, Cape Town 7700, South Africa.,Department of Chemistry, University of Cape Town, Cape Town 7700, South Africa.,Department of Biochemistry & Molecular Biology, University of Dalhousie, Halifax, NS B3H 4R2, Canada
| | | | - Neil Ravenscroft
- Department of Chemistry, University of Cape Town, Cape Town 7700, South Africa
| | - Michelle M Kuttel
- Department of Computer Science, University of Cape Town, Cape Town 7700, South Africa
| |
Collapse
|
2
|
Oluyemi WM, Samuel BB, Adewumi AT, Adekunle YA, Soliman MES, Krenn L. An Allosteric Inhibitory Potential of Triterpenes from Combretum racemosum on the Structural and Functional Dynamics of Plasmodium falciparum Lactate Dehydrogenase Binding Landscape. Chem Biodivers 2022; 19:e202100646. [PMID: 34982514 DOI: 10.1002/cbdv.202100646] [Citation(s) in RCA: 7] [Impact Index Per Article: 3.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 08/11/2021] [Accepted: 11/29/2021] [Indexed: 12/28/2022]
Abstract
Multidrug resistance is a significant drawback in malaria treatment, and mutations in the active sites of the many critical antimalarial drug targets have remained challenging. Therefore, this has necessitated the global search for new drugs with new mechanisms of action. Plasmodium falciparum lactate dehydrogenase (pfLHD), a glycolytic enzyme, has emerged as a potential target for developing new drugs due to the parasite reliance on glycolysis for energy. Strong substrate-binding is required in pfLDH enzymatic catalysis; however, there is a lack of information on small molecules' inhibitory mechanism bound to the substrate-binding pocket. Therefore, this study investigated a potential allosteric inhibition of pfLDH by targeting the substrate-binding site. The structural and functional behaviour of madecassic acid (MA), the most promising among the six triterpenes bound to pfLDH, were unravelled using molecular dynamic simulations at 300 ns to gain insights into its mechanism of binding and inhibition and chloroquine as a standard drug. The docking studies identified that the substrate site has the preferred position for the compounds even in the absence of a co-factor. The bound ligands showed comparably higher binding affinity at the substrate site than at the co-factor site. Mechanistically, a characteristic loop implicated in the enzyme catalytic activity was identified at the substrate site. This loop accommodates key interacting residues (LYS174, MET175, LEU177 and LYS179) pivotal in the MA binding and inhibitory action. The MA-bound pfLHD average RMSD (1.60 Å) relative to chloroquine-bound pfLHD RMSD (2.00 Å) showed higher stability for the substrate pocket, explaining the higher binding affinity (-33.40 kcal/mol) observed in the energy calculations, indicating that MA exhibited profound inhibitory activity. The significant pfLDH loop conformational changes and the allostery substrate-binding landscape suggested inhibiting the enzyme function, which provides an avenue for designing antimalarial compounds in the future studies of pfLDH protein as a target.
Collapse
Affiliation(s)
- Wande M Oluyemi
- Laboratory for Natural Products and Biodiscovery Research, Pharmaceutical Chemistry Department, Faculty of Pharmacy, University of Ibadan, Nigeria
| | - Babatunde B Samuel
- Laboratory for Natural Products and Biodiscovery Research, Pharmaceutical Chemistry Department, Faculty of Pharmacy, University of Ibadan, Nigeria
| | - Adeniyi T Adewumi
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Yemi A Adekunle
- Laboratory for Natural Products and Biodiscovery Research, Pharmaceutical Chemistry Department, Faculty of Pharmacy, University of Ibadan, Nigeria.,Department of Pharmaceutical Chemistry, Dora Akunyili College of Pharmacy, Igbinedion University, Okada, Benin City, Nigeria
| | - Mahmoud E S Soliman
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Liselotte Krenn
- Department of Pharmacognosy, University of Vienna, Althanstrasse 14, A-1090, Vienna, Austria
| |
Collapse
|
3
|
Adewumi AT, Elrashedy A, Soremekun OS, Ajadi MB, Soliman MES. Weak spots inhibition in the Mycobacterium tuberculosis antigen 85C target for antitubercular drug design through selective irreversible covalent inhibitor-SER124. J Biomol Struct Dyn 2020; 40:2934-2954. [PMID: 33155529 DOI: 10.1080/07391102.2020.1844061] [Citation(s) in RCA: 10] [Impact Index Per Article: 2.5] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 12/31/2022]
Abstract
Mycobacterium tuberculosis (Mtb) encoded secreted antigen 85 enzymes (Ag85A/Ag85B/Ag85C) play that critical roles in the virulence, survival and drug-resistant TB of the pathogen. Ag85 proteins are potential antitubercular drug targets because they are essential in the catalytic synthesis of trehalose moieties and mycolic acid attachment to the Mtb cell wall. Recently, experimental protocols led to the discovery of a selective covalent Ag85 inhibitor, β-isomer monocyclic enolphosphorus Cycliphostin (CyC8β) compound, which targets the Ag85 serine 124 to exhibit a promising therapeutic activity. For the first time, our study unravelled the structural features among Mtb Ag85C homologs and motions and dynamics of Ag85C when the CyC8β bound covalently and in open model conformations to the protein using bioinformatics tools and integrated Molecular dynamics simulations. Comparative Ag85C sequence analysis revealed conserved regions; 70% active site, 90% Adeniyi loop L1 and 50% loop L2, which acts as a switch between open and closed conformations. The average C-α atoms RMSD (2.05 Å) and RMSF (0.9 Å) revealed instability and high induced flexibility in the CyC8β covalent-bound compared to the apo and open model systems, which displayed more stability and lower fluctuations. DSSP showed structural transitions of α-helices to bend and loops to 310-helices in the bound systems. SASA of CyC8β covalent bound showed active site hydrophobic residues exposure to huge solvent. Therefore, these findings present the potential opportunity hotspots in Ag85C protein that would aid the structure-based design of novel chemical entities capable of resulting in potent antitubercular drugs.Communicated by Ramaswamy H. Sarma.
Collapse
Affiliation(s)
- Adeniyi T Adewumi
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Ahmed Elrashedy
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Opeyemi S Soremekun
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Mary B Ajadi
- Department of Medical Biochemistry, School of Laboratory Medicine and Medical Sciences, College of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| | - Mahmoud E S Soliman
- Molecular Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Durban, South Africa
| |
Collapse
|
4
|
Akher FB, Farrokhzadeh A, Ravenscroft N, Kuttel MM. A Mechanistic Study of a Potent and Selective Epidermal Growth Factor Receptor Inhibitor against the L858R/T790M Resistance Mutation. Biochemistry 2019; 58:4246-4259. [PMID: 31589411 DOI: 10.1021/acs.biochem.9b00710] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Abstract] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Indexed: 01/04/2023]
Abstract
Covalent targeting is a promising strategy for increasing the potency and selectivity of potential drug candidates. This therapeutic approach was recently reported for the epidermal growth factor receptor (EGFR), wherein a covalent binder, 20g [N-(3-{7-[2-methoxy-4-(4-methylpiperazin-1-yl)phenylamino]-3,4-dihydro-3-isopropyl-2,4-dioxopyrimido[4,5-d]pyrimidin-1(2H)-yl}phenyl)acrylamide], demonstrated significant selectivity and inhibitory activity toward the EGFR L858R/T790M double mutant (EGFRDM) relative to the EGFR wild-type form (EGFRWT). The enhanced therapeutic potency of 20g against EGFRDM is 263 times greater than that against EGFRWT, which necessitates a rational explanation for the underlying selective and inhibitory mechanisms. In this work, we investigate the differential binding modes of 20g with EGFRWT and EGFRDM using molecular dynamics simulations coupled with free energy calculations and further identify key residues involved in the selective targeting, binding, and inhibitory mechanisms mediated by 20g. We find that systematic orientational and conformational changes in the α-loop, p-loop, active loop, and αC-helix are responsible for the disparate binding mechanisms and inhibitory prowess of 20g with respect to EGFRWT and EGFRDM. The calculated binding free energies show good correlation with the experimental biological activity. The total binding free energy difference between EGFRWT-20g and EGFRDM-20g is -11.47 kcal/mol, implying that 20g binds more strongly to EGFRDM. This enhanced binding affinity of 20g for EGFRDM is a result of a large increase in the van der Waals and electrostatic interactions with three critical residues (Met790, Gln791, and Met793) that are chiefly responsible for the high-affinity interactions mediated by 20g with EGFRDM relative to EGFRWT.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Department of Computer Science , University of Cape Town , Cape Town 7701 , South Africa.,Department of Chemistry , University of Cape Town , Cape Town 7701 , South Africa
| | - Abdolkarim Farrokhzadeh
- School of Chemistry and Physics , University of KwaZulu-Natal , Private Bag X01 , Pietermaritzburg 3209 , South Africa
| | - Neil Ravenscroft
- Department of Chemistry , University of Cape Town , Cape Town 7701 , South Africa
| | - Michelle M Kuttel
- Department of Computer Science , University of Cape Town , Cape Town 7701 , South Africa
| |
Collapse
|
5
|
Tracing Potential Covalent Inhibitors of an E3 Ubiquitin Ligase through Target-Focused Modelling. Molecules 2019; 24:molecules24173125. [PMID: 31466292 PMCID: PMC6749425 DOI: 10.3390/molecules24173125] [Citation(s) in RCA: 4] [Impact Index Per Article: 0.8] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Download PDF] [Figures] [Journal Information] [Subscribe] [Scholar Register] [Received: 05/03/2019] [Revised: 06/13/2019] [Accepted: 06/14/2019] [Indexed: 12/20/2022] Open
Abstract
The Nedd4-1 E3 Ubiquitin ligase has been implicated in multiple disease conditions due its overexpression. Although the enzyme may be targeted both covalently and non-covalently, minimal studies provide effective inhibitors against it. Recently, research has focused on covalent inhibitors based on their characteristic, highly-selective warheads and ability to prevent drug resistance. This prompted us to screen for new covalent inhibitors of Nedd4-1 using a combination of computational approaches. However, this task proved challenging due to the limited number of electrophilic moieties available in virtual libraries. Therefore, we opted to divide an existing covalent Nedd4-1 inhibitor into two parts: a non-covalent binding group and a pre-selected α, β-unsaturated ester that forms the covalent linkage with the protein. A non-covalent pharmacophore model was built based on molecular interactions at the binding site. The pharmacophore was then subjected to virtual screening to identify structurally similar hit compounds. Multiple filtrations were implemented prior to selecting four hits, which were validated with a covalent conjugation and later assessed by molecular dynamic simulations. The results showed that, of the four hit molecules, Zinc00937975 exhibited advantageous molecular groups, allowing for favourable interactions with one of the characteristic cysteine residues. Predictive pharmacokinetic analysis further justified the compound as a potential lead molecule, prompting its recommendation for confirmatory biological evaluation. Our inhouse, refined, pharmacophore model approach serves as a robust method that will encourage screening for novel covalent inhibitors in drug discovery.
Collapse
|
6
|
Bjij I, Khan S, Ramharak P, Cherqaoui D, Soliman MES. Distinguishing the optimal binding mechanism of an E3 ubiquitin ligase: Covalent versus noncovalent inhibition. J Cell Biochem 2019; 120:12859-12869. [DOI: 10.1002/jcb.28556] [Citation(s) in RCA: 3] [Impact Index Per Article: 0.6] [Reference Citation Analysis] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 09/26/2018] [Revised: 01/07/2019] [Accepted: 01/14/2019] [Indexed: 01/18/2023]
Affiliation(s)
- Imane Bjij
- Molecular Bio‐Computation & Drug Design Lab, School of Health Sciences University of KwaZulu‐Natal Durban South Africa
- Université Cadi Ayyad, Faculté des Sciences Semlalia, Département de Chimie Marrakech Morocco
| | - Shama Khan
- Molecular Bio‐Computation & Drug Design Lab, School of Health Sciences University of KwaZulu‐Natal Durban South Africa
| | - Pritika Ramharak
- Molecular Bio‐Computation & Drug Design Lab, School of Health Sciences University of KwaZulu‐Natal Durban South Africa
| | - Driss Cherqaoui
- Université Cadi Ayyad, Faculté des Sciences Semlalia, Département de Chimie Marrakech Morocco
| | - Mahmoud E. S. Soliman
- Molecular Bio‐Computation & Drug Design Lab, School of Health Sciences University of KwaZulu‐Natal Durban South Africa
| |
Collapse
|
7
|
Akher FB, Farrokhzadeh A, Soliman MES. Covalent vs. Non-Covalent Inhibition: Tackling Drug Resistance in EGFR - A Thorough Dynamic Perspective. Chem Biodivers 2019; 16:e1800518. [PMID: 30548188 DOI: 10.1002/cbdv.201800518] [Citation(s) in RCA: 6] [Impact Index Per Article: 1.2] [Reference Citation Analysis] [Abstract] [Key Words] [Track Full Text] [Journal Information] [Subscribe] [Scholar Register] [Received: 10/02/2018] [Accepted: 12/10/2018] [Indexed: 12/20/2022]
Abstract
A persistent challenge in the treatment of non-small cell lung cancer (NSCLC) with EGFR is the emergence of drug-resistant caused by somatic mutations. The EGFR L858R/T790 M double mutant (EGFRDM ) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFRDM effectively. Recently, 11h and 45a, have been found to be potent inhibitors against EGFRDM through two distinctive mechanisms, non-covalent and covalent binding, respectively. However, the structural and dynamic implications of the two modes of inhibitions remain unexplored. Herein, two molecular dynamics simulation protocols, coupled with free-energy calculations, were applied to gain insight into the atomistic nature of each binding mode. The comparative analysis confirmed that there is a significant difference in the binding free energy between 11h and 45a (ΔΔGbind =-21.17 kcal/mol). The main binding force that governs the binding of both inhibitors is vdW, with a higher contribution for 45a. Two residues ARG841 and THR854 were found to have curtailed role in the binding of 45a to EGFRDM by stabilizing its flexible alcohol chain. The 45a binding to EGFRDM induces structural rearrangement in the active site to allow easier accessibility of 45a to target residue CYS797. The findings of this work can substantially shed light on new strategies for developing novel classes of covalent and non-covalent inhibitors with increased specificity and potency.
Collapse
Affiliation(s)
- Farideh Badichi Akher
- Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Abdolkarim Farrokhzadeh
- Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| | - Mahmoud E S Soliman
- Bio-computation and Drug Design Laboratory, School of Health Sciences, University of KwaZulu-Natal, Westville Campus, Durban, 4001, South Africa
| |
Collapse
|