Allosteric covalent inhibition of TOE1 as potential unexplored anti-cancer target: structure-based virtual screening and covalent molecular dynamics analysis.

Cancer remains a formidable challenge in therapeutic development owing to its complex molecular mechanisms and resistance to conventional treatments. Recent evidence suggests that TOE1 may play a role in cancer progression, making it an attractive target for therapeutic interventions, nevertheless, very limited research in literature has explored the potential of TOE1 inhibitors as anti-cancer. Herein, by exploring a library of 13,900 cysteine-targeted covalent inhibitors via a comprehensive virtual screening process, we sought to identify potential compounds that could be developed into effective cancer therapies against TOE1. The compounds were first screened based on their binding affinity, followed by their compliance with drug-like properties, and finally, by their effective covalent modeling to a reactive cysteine (Cys80). A total of 66 compounds, 28 compounds, and 3 compounds were found to have higher binding affinities, optimum drug-likeness, and higher covalent docking scores, respectively, than the reference compound. The top three screened compounds, 0462, 2204, and 7034, demonstrated favorable interaction profiles, covalent binding dynamics, free binding energetics, and per-residue energy contributions as compared to the reference compound. Notably, compound 0462 contributed to the highest free binding energy and significantly enhanced the stability and rigidity of TOE1, while restricting residue flexibility. This study provides an account of the molecular mechanics underpinning the covalent inhibition of TOE1, while providing a compelling case for further investigation and translation of the screened TOE1 inhibitors, particularly compound 0462, as novel therapeutics against cancer.

De Novo Rational Design of Peptide-Based Protein-Protein Inhibitors (Pep-PPIs) Approach by Mapping the Interaction Motifs of the PP Interface and Physicochemical Filtration: A Case on p25-Cdk5-Mediated Neurodegenerative Diseases.

Protein-protein interactions, or PPIs, are a part of every biological activity and have been linked to a number of diseases, including cancer, infectious diseases, and neurological disorders. As such, targeting PPIs is considered a strategic and vital approach in the development of new medications. Nonetheless, the wide and flat contact interface makes it difficult to find small-molecule PP inhibitors. An alternative strategy would be to use the PPI interaction motifs as building blocks for the design of peptide-based inhibitors. Herein, we designed 12-mer peptide inhibitors to target p25-inducing-cyclin-dependent kinase (Cdk5) hyperregulation, a PPI that has been shown to perpetuate neuroinflammation, which is one of the major causal implications of neurodegenerative disorders such as Alzheimer's disease, Parkinson's disease, and frontotemporal dementia. We generated a library of 5 062 500 peptide combination sequences (PCS) derived from the interaction motif of Cdk5/p25 PP interface. The 20 amino acids were differentiated into six groups, namely, hydrophobic (aliphatic), aromatic, basic, acidic, unique, and polar uncharged, on the basis of their physiochemical properties. To preserve the interaction motif necessary for ideal binding, de novo modeling of all possible peptide sequence substitutions was considered. A set of filters, backed by the Support Vector Machine (SVM) algorithm, was then used to create a shortlisted custom peptide library that met specific bioavailability, toxicity, and therapeutic relevance, leading to a refined library of 15 PCS. A greedy algorithm and coarse-grained force field were used to predict peptide structure and folding before subsequent modeling studies. Molecular docking was performed to estimate the relative binding affinities, and out of the top hits, Pep15 was subjected to molecular dynamics simulations and binding free-energy calculations in comparison to a known peptide inhibitor with experimental data (template peptide). Interestingly, the identified peptide through our protocol, Pep15, was found to show a significantly higher binding affinity than the reference template peptide (-48.10 ± 0.23 kcal/mol and -17.53 ± 0.27 kcal/mol, respectively). In comparison to the template peptide, Pep15 was found to possess a more compact and buried surface area, tighter binding landscape, and reduced conformational variability, leading to enhanced structural and kinetic stability of the Cdk5/p25 complex. Notably, both peptide inhibitors were found to have a minimal impact on the architectural integrity of the Cdk5/p25 secondary structure. Herein, we propose Pep15 as a novel and potentially disruptive peptide drug for Cdk5/p25-mediated neurodegenerative phenotypes that require further clinical investigation. The systematic protocol and findings of this report would serve as a valuable tool in the identification of critical PPI interface reactive residues, designing of analogs, and identification of more potent peptide-based PPI inhibitors.

Heat shock protein (Hsp27)-ceramide synthase (Cers1) protein-protein interactions provide a new avenue for unexplored anti-cancer mechanism and therapy.

Hsp27 is a member of the small heat-shock proteins (sHSPs) - the known cellular line of defence against abnormal protein folding behaviors. Nevertheless, its upregulation is linked to a variety of pathological disorders, including several types of cancers. The ceramide synthases (CerS) mediate the synthesis of ceramide, a critical structural and signaling lipid. Functionally, downstream ceramide metabolites are implicated in the apoptosis process and their abnormal functionality has been linked to anticancer resistance. Studies showed that CerS1 are possibly inhibited by Hsp27 leading to biochemical anticancer effects in vitro. Nevertheless, the nature of such protein-protein interaction (PPI) has not been considerably investigated in molecular terms, hence, we present the first description of the dynamics CerS1-Hsp27 interaction landscapes using molecular dynamics simulations. Time-scale molecular dynamics simulation analysis indicated a system-wide conformational events of decreased stability, increased flexibility, reduced compactness, and decreased folding of CerS1. Analysis of binding energy showed a favorable interaction entailing 56 residues at the interface and a total stabilizing energy of -158 KJ/mol. The CerS1 catalytic domain experienced an opposite trend compared to the protein backbone. Yet, these residues adopted a highly compact conformation as per DCCM and DSSP analysis. Furthermore, conserved residues (SER 212, ASP 213, ALA 240, GLY 243, ASP 319) comprising the substrate shuttling machinery showed notable rigidity implying a restrained ceramide precursor access and assembly; hence, a possible inhibitory mechanism. Findings from this report would streamline a better molecular understanding of CerS1-Hsp27 interactions and decipher its potential avenue toward unexplored anti-cancer mechanisms and therapy.

Unveiling therapeutic frontiers: DON/DRP-104 as innovative Plasma kallikrein inhibitors against carcinoma-associated hereditary angioedema shocks - a comprehensive molecular dynamics exploration.

Human plasma kallikrein (PKa) is a member of the serine protease family and serves as a key mediator of the kallikrein-kinin system (KKS), which is known for its regulatory roles in inflammation, vasodilation, blood pressure, and coagulation. Genetic dysregulation of KKS leads to Hereditary Angioedema (HAE), which is characterized by spontaneous, painful swelling in various body regions. Importantly, HAE frequently coexists with various cancers. Despite substantial efforts towards the development of PKa inhibitors for HAE, there remains a need for bifunctional agents addressing both anti-cancer and anti-HAE aspects, especially against carcinoma-associated comorbid HAE conditions. Consequently, we investigated the therapeutic potential of the anti-glutamine prodrug, isopropyl(S)-2-((S)-2-acetamido-3-(1H-indol-3-yl)-propanamido)-6-diazo-5-oxo-hexanoate (DRP-104), and its active form, 6-Diazo-5-oxo-l-norleucine (DON), recognized for their anti-cancer properties, as novel PKa inhibitors. Utilizing structure-based in silico methods, we conducted a comparative analysis with berotralstat, a clinically approved HAE prophylactic, and sebetralstat, an investigational HAE therapeutic agent, in Phase 3 clinical trials. Inhibiting PKa with DON resulted in relatively heightened structural stability, rigidity, restricted protein folding, and solvent-accessible loop exposure, contributing to increased intra-atomic hydrogen bond formation. Conversely, PKa inhibition with DRP-104 induced restricted residue flexibility and significantly disrupted the critical SER195-HIS57 arrangement in the catalytic triad. Both DON and DRP-104, along with the reference drugs, induced strong cooperative intra-residue motion and bidirectional displacement in the PKa architecture. The results revealed favorable binding kinetics of DON/DRP-104, showing thermodynamic profiles that were either superior or comparable to those of the reference drugs. These findings support their consideration for clinical investigations into the management of carcinoma-associated HAE.

Exploring the potential of biologically active phenolic acids from marine natural products as anticancer agents targeting the epidermal growth factor receptor.

The epidermal growth factor receptor (EGFR) dimerizes upon ligand bindings to the extracellular domain that initiates the downstream signaling cascades and activates intracellular kinase domain. Thus, activation of autophosphorylation through kinase domain results in metastasis, cell proliferation, and angiogenesis. The main objective of this research is to discover more promising anti-cancer lead compound against EGRF from the phenolic acids of marine natural products using in-silico approaches. Phenolic compounds reported from marine sources are reviewed from previous literatures. Furthermore, molecular docking was carried out using the online tool CB-Dock. The molecules with good docking and binding energies scores were subjected to ADME, toxicity and drug-likeness analysis. Subsequently, molecules from the docking experiments were also evaluated using the acute toxicity and MD simulation studies. Fourteen phenolic compounds from the reported literatures were reviewed based on the findings, isolation, characterized and applications. Molecular docking studies proved that the phenolic acids have good binding fitting by forming hydrogen bonds with amino acid residues at the binding site of EGFR. Chlorogenic acid, Chicoric acid and Rosmarinic acid showed the best binding energies score and forming hydrogen bonds with amino acid residues compare to the reference drug Erlotinib. Among these compounds, Rosmarinic acid showed the good pharmacokinetics profiles as well as acute toxicity profile. The MD simulation study further revealed that the lead complex is stable and could be future drug to treat the cancer disease. Furthermore, in a wet lab environment, both in-vitro and in-vivo testing will be employed to validate the existing computational results.Communicated by Ramaswamy H. Sarma.

Therapeutic Path to Triple Knockout: Investigating the Pan-inhibitory Mechanisms of AKT, CDK9, and TNKS2 by a novel 2-Phenylquinazolinone derivative in Cancer Therapy- An In-Silico Investigation.

BACKGROUND: Blocking the oncogenic Wnt//ß-catenin pathway has of late been investigated as a viable therapeutic approach in the treatment of cancer. This involves the multi-targeting of certain members of the tankyrase-kinase family; tankyrase 2 (TNKS2), protein kinase B (AKT), and cyclin-dependent kinase 9 (CDK9), which propagate the oncogenic Wnt/ß-catenin signalling pathway. METHODS: During a recent investigation, the pharmacological activity of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one was repurposed to serve as a 'triple-target' inhibitor of TNKS2, AKT and CDK9. Yet, the molecular mechanism that surrounds its multi-targeting activity remains unanswered. As such, this study aims to explore the pan-inhibitory mechanism of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one towards AKT, CDK9, and TNKS2, using in silico techniques. RESULTS: Results revealed favourable binding affinities of -34.17 kcal/mol, -28.74 kcal/mol, and -27.30 kcal/mol for 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one towards TNKS2, CDK9, and AKT, respectively. Pan-inhibitory binding of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one is illustrated by close interaction with specific residues on tankyrase-kinase. Structurally, 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one had an impact on the flexibility, solvent-accessible surface area, and stability of all three proteins, which was illustrated by numerous modifications observed in the unbound as well as the bound states of the structures, which evidenced the disruption of their biological function. Prediction of the pharmacokinetics and physicochemical properties of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one further established its inhibitory potential, evidenced by the favourable absorption, metabolism, excretion, and minimal toxicity properties. CONCLUSION: The following structural insights provide a starting point for understanding the pan-inhibitory activity of 2-(4-aminophenyl)-7-chloro-3H-quinazolin-4-one. Determining the criticality of the interactions that exist between the pyrimidine ring and catalytic residues could offer insight into the structure-based design of innovative tankyrase-kinase inhibitors with enhanced therapeutic effects.

Multi-dimensional structural footprint identification for the design of potential scaffolds targeting METTL3 in cancer treatment from natural compounds.

CONTEXT: [Formula: see text]-adenosine-methyltransferase (METTL3) is the catalytic domain of the 'writer' proteins which is involved in the post modifications of [Formula: see text]-methyladinosine ([Formula: see text]). Though its activities are essential in many biological processes, it has been implicated in several types of cancer. Thus, drug developers and researchers are relentlessly in search of small molecule inhibitors that can ameliorate the oncogenic activities of METTL3. Currently, STM2457 is a potent, highly selective inhibitor of METTL3 but is yet to be approved. METHODS: In this study, we employed structure-based virtual screening through consensus docking by using AutoDock Vina in PyRx interface and Glide virtual screening workflow of Schrodinger Glide. Thermodynamics via MM-PBSA calculations was further used to rank the compounds based on their total free binding energies. All atom molecular dynamics simulations were performed using AMBER 18 package. FF14SB force fields and Antechamber were used to parameterize the protein and compounds respectively. Post analysis of generated trajectories was analyzed with CPPTRAJ and PTRAJ modules incorporated in the AMBER package while Discovery studio and UCSF Chimera were used for visualization, and origin data tool used to plot all graphs. RESULTS: Three compounds with total free binding energies higher than STM2457 were selected for extended molecular dynamics simulations. The compounds, SANCDB0370, SANCDB0867, and SANCDB1033, exhibited stability and deeper penetration into the hydrophobic core of the protein. They engaged in relatively stronger intermolecular interactions involving hydrogen bonds with resultant increase in stability, reduced flexibility, and decrease in the surface area of the protein available for solvent interactions suggesting an induced folding of the catalytic domain. Furthermore, in silico pharmacokinetics and physicochemical analysis of the compounds revealed good properties suggesting these compounds could serve as promising MEETL3 entry inhibitors upon modifications and optimizations as presented by natural compounds. Further biochemical testing and experimentations would aid in the discovery of effective inhibitors against the berserk activities of METTL3.

Chetomin, a SARS-CoV-2 3C-like Protease (3CLpro) Inhibitor: In Silico Screening, Enzyme Docking, Molecular Dynamics and Pharmacokinetics Analysis.

The emergence of the Coronavirus Disease 2019 (COVID-19) pandemic caused by severe acute respiratory syndrome coronavirus-2 (SARS-CoV-2) has led to over 6 million deaths. The 3C-like protease (3CLpro) enzyme of the SARS-CoV-2 virus is an attractive druggable target for exploring therapeutic drug candidates to combat COVID-19 due to its key function in viral replication. Marine natural products (MNPs) have attracted considerable attention as alternative sources of antiviral drug candidates. In looking for potential 3CLpro inhibitors, the MNP database (>14,000 molecules) was virtually screened against 3CLpro with the assistance of molecular docking computations. The performance of AutoDock and OEDocking software in anticipating the ligand-3CLpro binding mode was first validated according to the available experimental data. Based on the docking scores, the most potent MNPs were further subjected to molecular dynamics (MD) simulations, and the binding affinities of those molecules were computed using the MM-GBSA approach. According to MM-GBSA//200 ns MD simulations, chetomin (UMHMNP1403367) exhibited a higher binding affinity against 3CLpro than XF7, with ΔGbinding values of −55.5 and −43.7 kcal/mol, respectively. The steadiness and tightness of chetomin with 3CLpro were evaluated, revealing the high stabilization of chetomin (UMHMNP1403367) inside the binding pocket of 3CLpro throughout 200 ns MD simulations. The physicochemical and pharmacokinetic features of chetomin were also predicted, and the oral bioavailability of chetomin was demonstrated. Furthermore, the potentiality of chetomin analogues −namely, chetomin A-D− as 3CLpro inhibitors was investigated. These results warrant further in vivo and in vitro assays of chetomin (UMHMNP1403367) as a promising anti-COVID-19 drug candidate.

Deciphering the molecular mechanisms of selective non-covalency demonstrated differentially by 9-Allylnaphtho[1,8-ef]isoindole-7,8,10(9H)-trione (C11) against fibroblast growth factor receptors 1-4.

The specific inhibition of aberrant Fibroblast Growth Factor Receptors (FGFRs) has been identified as a feasible strategy to therapeutically ameliorate their respective carcinogenic involvements. High homology among these proteins has however limited efforts towards the discovery of selective small-molecule compounds due to undesirable effects elicited by pan-FGFR inhibitors. A recent study showed the selective activity of a new compound C11 which was >52 times more potent against FGFR1 than FGFR2 and FGFR3, and 4 times than FGFR4. This C11 selective non-covalency was investigated in this study using computational methods since it has remained unresolved. Structural findings revealed that C11 enhanced structural perturbations in FGFR1 with less prominent effects in other FGFRs. High deviations also characterized the C11-bound active pocket of FGFR1 with notable fluctuations across the constituent P-loop, αC helix, hinge region, catalytic, and activation loops. These induced motions were essential for optimal C11 motion an d positioning of its phenalenone ring and prop-2-en-l-yl moiety at the FGFR1 active pocket to interact stably and strongly with A564FGFR1, L484FGFR1, Y563FGFR1, and E562FGFR1 which as well had high energy contributions. C11 exhibited highly unstable binding in F GFRs2-3 with a more steady interaction with FGFR4. Free binding energy (ΔGbind) analyses further estimated the highest interaction energy for C11-FGFR1 with favorable desolvation energy that indicated a deep hydrophobic pocket binding for C11 in FGFR1 compared to other FGFRs. We believe rational insights from this study will contribute to the structure-based design of highly specific FGFR1 inhibitors.Communicated by Ramaswamy H. Sarma.

Impact of compound mutations I1171N + F1174I and I1171N + L1198H on the structure of ALK in NSCLC pathogenesis: atomistic insights.

Anaplastic lymphoma kinase (ALK) fusion genes are found in 3%-5% of non-small cell lung cancers (NSCLCs). NSCLC is the most common type of lung cancer, accounting for 84% of all lung cancer diagnoses. Available treatment options for ALK-positive NSCLCs involve the use of ALK tyrosine kinase inhibitors (ALK-TKIs) which have shown to be effective with a high response rate. Nonetheless, the emergence of multiple compound mutations such as I1171N + F1174I or I1171N + L1198H has been reported to cause resistance to all approved ALK-TKIs. However, the underlying molecular mechanisms surrounding the impact of these compound mutants remain poorly understood. Hence, we performed molecular dynamics simulations to characterize the structural effects and functional implications of these compound mutations. Findings revealed a destabilizing effect on ALK by mutants as compared to the wild-type ALK structure. Also, further insights revealed a lower root-mean-squared fluctuation, radius of gyration, and solvent-accessible surface area values of I1171N + F1174I and I1171N + L1198H ALK compound mutations suggesting that the mutants have a more compact structure and a smaller surface area than the wild-type protein. The mutants also distorted the activation loop residues (Tyr1278, Tyr1282, and Tyr1283) in the ALK structure, which further identify them as possible disruptors of phosphorylation. In contrast to wild conformation, the mutant conformations exhibited a reduced node degree in their residue interaction networks. Collectively, our findings provide deeper insights into the deleterious effects of I1171N + F1174I and I1171N + L1198H ALK compound mutations, which may contribute to NSCLC pathogenesis.Communicated by Ramaswamy H. Sarma.

Inside the cracked kernel: establishing the molecular basis of AMG510 and MRTX849 in destabilising KRASG12C mutant switch I and II in cancer treatment.

The Kirsten rat sarcoma oncoprotein (KRAS) has been punctuated by drug development failures for decades due to frequent mutations that occur mostly at codon 12 and the seemingly intractable targeting of the protein. However, with advances in covalent targeting, the oncoprotein is being expunged from the 'undruggable' list of proteins. This feat has seen some covalent drugs at different stages of clinical trials. The advancement of AMG510 and MRTX849 as inhibitors of cysteine mutated KRAS (KRASG12C) to phase-III clinical trials informed the biased selection of AMG510 and MRTX849 for this study. Despite this advance, the molecular and atomistic modus operandi of these drugs is yet to come to light. In this study, we employed computational tools to unravel the atomistic interactions and subsequent conformational effects of AMG510 and MRTX849 on the mutant KRASG12C. It was revealed that AMG510 and MRTX849 complexes presented similar total free binding energies, (ΔGbind), of -88.15 ± 5.96 kcal/mol and -88.71 ± 7.70 kcal/mol, respectively. Gly10, Lys16, Thr58, Gly60, Glu62, Glu63, Arg68, Asp69, Met72, His95, Tyr96, Gln99, Arg102 and Val103 interacted prominently with AMG510 and MRTX849. These residues interacted with the pharmacophoric moieties of AMG510 and MRTX849 via hydrogen bonds with decreasing bond lengths at various stages of the simulation. These interactions together with pi-pi stacking, pi-sigma and pi-alkyl interactions induced unfolding of switch I whiles compacting switch II, which could interrupt the binding of effector proteins to these interfaces. These insights present useful atomistic perspectives into the success of AMG510 and MRTX849 which could guide the design of more selective and potent KRAS inhibitors.Communicated by Ramaswamy H. Sarma.

Comparative Dynamic Features of Apo and Bound MDM2 Protein Reveal the Mechanism of Inhibitor Recognition for Anti-Cancer Activity.

BACKGROUND: Mouse Double Minute 2 Homolog (MDM2) oncogenic protein is the principal cellular antagonist of the p53 tumor suppressor gene. Restoration of p53 activity by inhibiting the MDM2-P53 interactions at the molecular level has become the cornerstone of cancer research due to its promising anticancer effects. Natural medicinal products possess various chemical structures and represent an essential source for drug discovery. α-Mangostin (AM) and gambogic acid (G250) are plant-derived compounds that showed inhibitory effects on MDM2-P53 interactions in vitro and in vivo. METHODS: Despite the many clinical studies which performed deeper insight about the molecular understanding of the structural mechanisms exhibited by α-Mangostin and Gambogic acid-binding to MDM2 remains critical. In this study, comparative molecular dynamics simulations were performed for each Apo and bound p53 and MDM2 proteins to shed light on the MDM2-p53 interactions and get a better understanding of the inhibition mechanisms. RESULTS: Results revealed atomistic interaction of AM and G250 within the MDM2-p53 interaction cleft. Both compounds mediate the interaction between the α-helix motifs of the p53 amino-terminal domain, which caused a significant separation between orthogonally opposed residues, specifically Lys8 and Gly47 residues of the p53 and MDM2, respectively. Contrasting changes in magnitudes were observed in per-residue fluctuation on AM and G250 (~0.04 nm and ~2.3 nm, respectively). The Radius of gyration (~0.03 nm and 0.04 nm, respectively), C-alpha deviations (~0.06 nm and 0.1 nm, respectively). The phenolic group of AM was found to establish hydrogen interactions with Glu28 and His96 residues of MDM2. The trioxahexacyclo-ring of G250 also forms hydrogen bond interactions with Lys51 and Leu26 residues of MDM2. CONCLUSION: Utilizing the information provided on the inhibitory binding mode adopted by each compound in this study may further assist in the tailored designs for cancer therapeutics.

Investigating the Impact of Covalent and Non-covalent Binding Modes of Inhibitors on Bruton's Tyrosine Kinase in the Treatment of B Cell Malignancies - Computational Insights.

BACKGROUND: Bruton tyrosine kinase plays a key role in the survival, proliferation, activation, and differentiation of B-lineage cells and the signaling of other receptors. It is overexpressed and constitutively active in the pathogenesis of B cell malignancies and has therefore become a target for therapeutic intervention. Some success has been achieved in the discovery of small molecules, especially in the development of irreversible inhibitors. However, these inhibitors are punctuated by off target effects and have also become less effective in patients with mutations at Cys481. This motivated the search for inhibitors with improved efficacy and different binding modes. METHODS: In this study, we employed two new second generation inhibitors with different binding modes, Zanubrutinib and AS-1763, which are at various levels of clinical trials, to highlight the molecular determinants in the therapeutic inhibition of BTK through computational studies. RESULTS: This study revealed that Zanubrutinib and AS-1763 exhibited free total binding energies of -98.76 ± 4.63 kcal/mol and -51.81 ± 9.94 kcal/mol, respectively, with Zanubrutinib engaging in peculiar hydrogen bond interactions with the hinge residues Glu475 and Met477 including Asn484 and Tyr485 while AS-1763 engaged Lys430, Asp539, and Arg525. These residues contributed the most towards the free total binding energy with energies above -1.0 kcal/mol. The compounds further interacted differentially with other binding site residues through pi-alkyl, pi-cation, pianion, pi-pi-T-shaped, pi-sigma, pi-sulfur and pi-donor hydrogen bonds, and Van der Waals interactions. These interactions resulted in differential fluctuations of the residues with the consequential unfolding of the protein. CONCLUSION: Insights herein would be useful in guiding the discovery of more selective and potent small molecules.

Pharmacophore-based Identification of Potential Mutant Isocitrate Dehydrogenases I/2 Inhibitors: An Explorative Avenue in Cancer Drug Design.

BACKGROUND: Heterozygous mutations in the cytoplasmic and mitochondrial isoforms of isocitrate dehydrogenase enzymes 1 and 2 subtypes have been extensively exploited as viable druggable targets, as they decrease the affinity of isocitrate and higher affinity of D-2-hydroxyglutarate, an oncometabolite. OBJECTIVE: Vorasidenib (AG-881) has recently been reported as a promising dual inhibitor of mutant isocitrate dehydrogenase 1 and 2 with the ability to penetrate the blood-brain barrier towards the treatment of low-grade glioma. In order to combat drug resistance and toxicity levels, this compelled us to further investigate this substance as a basis for the creation of potential selective inhibitors of mutant isocitrate dehydrogenases 1 and 2. METHODS: By employing a wide range of computational techniques, binding moieties of AG-881 that contributed towards its selective binding to isocitrate dehydrogenase enzymes 1 and 2 were identified and subsequently used to generate pharmacophore models for the screening of potential inhibitor drugs that were further assessed by their pharmacokinetics and physicochemical properties. RESULTS: AG-881 was identified as the most favorable candidate for isocitrate dehydrogenase enzyme 1, exhibiting a binding free energy of -28.69 kcal/mol. ZINC93978407 was the most favorable candidatefor isocitrate dehydrogenase enzyme 2, displaying a strong binding free energy of -27.10 kcal/mol. ZINC9449923 and ZINC93978407 towards isocitrate dehydrogenase enzyme 1 and 2 showed good protein structural stability with a low radius of gyration values relative to AG-881. CONCLUSION: We investigated that ZINC9449923 of isocitrate dehydrogenase enzyme 1 and ZINC 93978407 of isocitrate dehydrogenase enzyme 2 could serve as promising candidates for the treatment of lower-grade glioma as they cross the blood-brain barrier, and present with lower toxicity levels relative to AG-881.

Exploring the composition of protein-ligand binding sites for cancerous inhibitor of PP2A (CIP2A) by inhibitor guided binding analysis: paving a new way for the Discovery of drug candidates against triple negative breast cancer (TNBC).

Triple-negative breast cancer (TNBC) is associated with high-grade invasive carcinoma leading to a 10% to 15% death rate in younger premenopausal women. Targeting cancerous inhibitors of protein phosphatase (CIP2A) has been a highly effective approach for exploring therapeutic drug candidates. Lapatinib, a dual tyrosine kinase inhibitor, has shown promising inhibition properties by inducing apoptosis in TNBC carcinogenesis in vivo. Despite knowledge of the 3D structure of CIP2A, no reports provide insight into CIP2A ligand binding sites. To this effect, we conducted in silico site identification guided by lapatinib binding. Four of the five sites identified were cross-validated, and the stem domain revealed more excellent ligand binding affinity. The binding affinity of lapatinib in these sites was further computed using the Molecular Mechanics/Poisson-Boltzmann Surface Area (MM/PBSA) approach. According to MM/PBSA//200 ns MD simulations, lapatinib exhibited a higher binding affinity against CIP2A in site 2 with ΔG critical values of -37.1 kcal/mol. The steadiness and tightness of lapatinib with CIP2A inside the stem domain disclosed glutamic acid-318 as the culprit amino acid with the highest electrostatic energy. These results provide clear information on the CIP2A domain capable of ligand binding and validate lapatinib as a promising CIP2A inhibitor in TNBC carcinogenesis.

Adsorption of Chlormethine Anti-Cancer Drug on Pure and Aluminum-Doped Boron Nitride Nanocarriers: A Comparative DFT Study.

The efficacy of pure and aluminum (Al)-doped boron nitride nanocarriers (B12N12 and AlB11N12) in adsorbing Chlormethine (CM), an anti-cancer drug, was comparatively dissected by means of the density functional theory method. The CM∙∙∙B12N12 and ∙∙∙AlB11N12 complexes were studied within two configurations, A and B, in which the adsorption process occurred via N∙∙∙ and Cl∙∙∙B/Al interactions, respectively. The electrostatic potential affirmations confirmed the opulent ability of the studied nanocarriers to engage in delivering CM via two prominent electrophilic sites (B and Al). Furthermore, the adsorption process within the CM∙∙∙AlB11N12 complexes was noticed to be more favorable compared to that within the CM∙∙∙B12N12 analog and showed interaction and adsorption energy values up to -59.68 and -52.40 kcal/mol, respectively, for configuration A. Symmetry-adapted perturbation theory results indicated that electrostatic forces were dominant in the adsorption process. Notably, the adsorption of CM over B12N12 and AlB11N12 nanocarriers exhibited predominant changes in their electronic properties. An elemental alteration was also revealed for the softness and hardness of B12N12 and AlB11N12 nanocarriers before and following the CM adsorption. Spontaneity and exothermic nature were obviously observed for the studied complexes and confirmed by the negative values of thermodynamic parameters. In line with energetic manifestation, Gibbs free energy and enthalpy change were drastically increased by the Al doping process, with values raised to -37.15 and -50.14 kcal/mol, respectively, for configuration A of the CM∙∙∙AlB11N12 complex. Conspicuous enhancement was noticed for the adsorption process in the water phase more than that in the gas phase and confirmed by the negative values of the solvation energy up to -53.50 kcal/mol for configuration A of the CM∙∙∙AlB11N12 complex. The obtained outcomes would be the linchpin for the future utilization of boron nitride as a nanocarrier.

Probing Alterations in MDM2 Catalytic Core Structure Effect of Garcinia Mangostana Derivatives: Insight from Molecular Dynamics Simulations.

The MDM2-p53 protein-protein interaction is a promising model for researchers to design, study, and discover new anticancer drugs. The design of therapeutically active compounds that can maintain or restore the binding of MDM2 to p53 has been found to limit the oncogenic activities of both. This led to the current development of a group of xanthone-core and cis-imidazoline analogs compounds, among which γ-Mangostin (GM), α-Mangostin (AM), and Nutlin exhibited their MDM2-p53 interaction inhibitory effects. Therefore, in this study, we seek to determine the mechanisms by which these compounds elicit MDM2-p53 interaction targeting. Unique to the binding of GM, AM, and Nutlin, from our findings, they share the same three active site residues Val76, Tyr50, and Gly41, which represent the top active side residues that contribute to high electrostatic energy. Consequently, the free binding energy contributed enormously to the binding of these compounds, which culminated in the high binding affinities of GM, AM, and Nutlin with high values. Furthermore, GM, AM, and Nutlin commonly interrupted the stable and compact conformation of MDM2 coupled with its active site, where Cα deviations were relatively high. We believe that our findings would assist in the design of more potent active anticancer drugs.

Battling BTK mutants with noncovalent inhibitors that overcome Cys481 and Thr474 mutations in Waldenström macroglobulinemia therapy: structural mechanistic insights on the role of fenebrutinib.

Recently, the non-covalent Bruton tyrosine kinase (BTK) inhibitor fenebrutinib was presented as a therapeutic option with strong inhibitory efficacy against a single (C481S) and double (T474S/C481S) BTK variant in the treatment of Waldenström macroglobulinemia (WM). However, the molecular events surrounding its inhibition mechanism towards this variant remain unresolved. Herein, we employed in silico methods such as molecular dynamic simulation coupled with binding free energy estimations to explore the mechanistic activity of the fenebrutinib on (C481S) and (T474S/C481S) BTK variant, at a molecular level. Our investigations reveal that amino acid arginine contributed immensely to the total binding energy, this establishing the cruciality of amino acid residues, Arg132 and Arg156 in (C481S) and Arg99, Arg137, and Arg132 in (T474S/C481S) in the binding of fenebrutinib towards both BTK variants. The structural orientations of fenebrutinib within the respective hydrophobic pockets allowed favorable interactions with binding site residues, accounting for its superior binding affinity by 24.5% and relative high hydrogen bond formation towards (T474S/C481S) when compared with (C481S) BTK variants. Structurally, fenebrutinib impacted the stability, flexibility, and solvent accessible surface area of both BTK variants, characterized by various alterations observed in the bound and unbound structures, which proved enough to disrupt their biological function. Findings from this study, therefore, provide insights into the inhibitory mechanism of fenebrutinib at the atomistic level and reveal its high selectivity towards BTK variants. These insights could be key in designing and developing BTK mutants' inhibitors to treat Waldenström macroglobulinemia (WM).

Talazoparib Dual-targeting on Poly (ADP-ribose) Polymerase-1 and -16 Enzymes Offers a Promising Therapeutic Strategy in Small Cell Lung Cancer Therapy: Insight from Biophysical Computations.

In recent times, inhibition of poly (ADP-ribose) polymerase (PARP) enzymes by pharmacological drugs has attracted much attention as an anticancer therapy. As reported, PARP-16 has been discovered as a novel anticancer target for small cell lung cancer, and that the inhibition of both PARP-16 and PARP-1 by talazoparib can increase the overall effectiveness of talazoparib in the SCLC treatment. In this study, we employed computational approaches to investigate the differential inhibitory potency of Talazoparib on PARP-1 and PARP-16. Talazoparib has excellent PARP-1 and PARP-16 binding activities, as revealed by the ΔGbind (total binding energy). Pp16-tpb had binding energy of -34.85 kcal/mol, while pp1-tpb had a binding energy of -26.36 kcal/mol. The binding activity of Talazoparib on both PARP-1 and PARP-16 was significantly influenced by van der Waal and electrostatic interactions. Correspondingly, according to the findings of this study, binding residues with total binding energy greater than 1.00 kcal/mol contributed considerably to the Talazoparib's binding activities on PARP-1 and PARP-16. We believe the findings of this study will pave the way for developing dual targeting of PARP enzymes as a strategy for small-cell lung cancer treatment.

Exploring Natural Product Activity and Species Source Candidates for Hunting ABCB1 Transporter Inhibitors: An In Silico Drug Discovery Study.

The P-glycoprotein (P-gp/ABCB1) is responsible for a xenobiotic efflux pump that shackles intracellular drug accumulation. Additionally, it is included in the dud of considerable antiviral and anticancer chemotherapies because of the multidrug resistance (MDR) phenomenon. In the search for prospective anticancer drugs that inhibit the ABCB1 transporter, the Natural Product Activity and Species Source (NPASS) database, containing >35,000 molecules, was explored for identifying ABCB1 inhibitors. The performance of AutoDock4.2.6 software to anticipate ABCB1 docking score and pose was first assessed according to available experimental data. The docking scores of the NPASS molecules were predicted against the ABCB1 transporter. Molecular dynamics (MD) simulations were conducted for molecules with docking scores lower than taxol, a reference inhibitor, pursued by molecular mechanics-generalized Born surface area (MM-GBSA) binding energy estimations. On the basis of MM-GBSA calculations, five compounds revealed promising binding affinities as ABCB1 inhibitors with ΔGbinding < −105.0 kcal/mol. The binding affinity and stability of the identified inhibitors were compared to the chemotherapeutic agent. Structural and energetical analyses unveiled great steadiness of the investigated inhibitors within the ABCB1 active site throughout 100 ns MD simulations. Conclusively, these findings point out that NPC104372, NPC475164, NPC2313, NPC197736, and NPC477344 hold guarantees as potential ABCB1 drug candidates and warrant further in vitro/in vivo tests.