Metabolites profiling and cheminformatics bioprospection of selected medicinal plants against the main protease and RNA-dependent RNA polymerase of SARS-CoV-2

Despite the existence of some vaccines, SARS-CoV-2 (S-2) infections persist for various reasons relating to vaccine reluctance, rapid mutation rate, and an absence of specific treatments targeted to the infection. Due to their availability, low cost and low toxicity, research into potentially repurposing phytometabolites as therapeutic alternatives has gained attention. Therefore, this study explored the antiviral potential of metabolites of some medicinal plants [Spondias mombin, Macaranga barteri and Dicerocaryum eriocarpum (Sesame plant)] identified using liquid chromatography-mass spectrometry (LCMS) as possible inhibitory agents against the S-2 main protease (S-2 MP) and RNA-dependent RNA polymerase (RP) using computational approaches. Molecular docking was used to identify the compounds with the best affinities for the selected therapeutics targets. Afterwards, compounds with poor physicochemical characteristics, pharmacokinetics, and drug-likeness were screened out. The top-ranked compounds were further subjected to a 120-ns molecular dynamics (MD) simulation. Only quercetin 3-O-rhamnoside (-48.77 kcal/mol) had higher binding free energy than the reference standard (zafirlukast) (-44.99 kcal/mol) against S-2 MP. Conversely, all the top-ranked compounds (ellagic acid hexoside, spiraeoside, apigenin-4'-glucoside and chrysoeriol 7-glucuronide) except gnetin L (-24.24 kcal/mol) had higher binding free energy (-55.19 kcal/mol, -52.75 kcal/mol, -47.22 kcal/mol and -43.35 kcal/mol) respectively, against S-2 RP relative to the reference standard (-34.79 kcal/mol). The MD simulations study further revealed that the investigated inhibitors are thermodynamically stable and form structurally compatible complexes that impede the regular operation of the respective S-2 therapeutic targets. Although, these S-2 therapeutic candidates are promising, further in vitro and in vivo evaluation is required and highly recommended.Communicated by Ramaswamy H. Sarma.

Probing Dual Covalent Irreversible Inhibition of EGFR/FGFR4 by Electrophilic-Based Natural Compounds to Overcome Resistance and Enhance Combination Therapeutic Potentials and Management of Hepatocellular Carcinoma (HCC)

Hepatocellular carcinoma (HCC) is one of the most prevalent cancer types in the world and accounts for the majority of cases of primary liver cancer. A crucial part of the carcinogenesis of HCC involves aberrant stimulation of the FGF19-FGFR4 signaling pathway. Therefore, FGFR4 inhibition has become a strategic therapeutic approach for the treatment of HCC. However, the clinical treatment procedure is significantly hampered by the prevalence of kinase inhibitors resistance. It was recently established that the activation of EGFR signaling was found to be one of the primary mechanisms mediating the acquired resistance to FGFR4 inhibitors, moreover, sensitivity to FGFR4 inhibitors was effectively restored by inhibiting EGFR. These results provide compelling evidence that dual inhibition of EGFR and FGFR4 could represent a viable therapeutic approach to overcome resistance, hence enhanced management of HCC. To this end, we proposed a dual irreversible inhibition strategy through covalent binding by naturally occurring electrophilic warhead-bearing compounds (curcumin, deoxyelephantopin, eupalmerin acetate, syringolin A and andrographolide) to covalently target both EGFR and FGFR4 through cysteine residues, Cys797 and Cys552, respectively. Covalent docking and covalent molecular dynamics (MM/MDcov) simulations combined with thermodynamic binding free energy calculations were performed, and the results were compared against known potent and selective covalent EGFR and FGFR4 inhibitors with available X-ray crystal structures, Afatinib and BLU9931, respectively. Curcumin, deoxyelephantopin, eupalmerin acetate, syringolin A, and andrographolide showed relative binding free energies of -22.85, -17.14, -12.98, -21.81, and - 19.00 kcal/mol against EGFR and - 41.06, -29.45, -24.76, -40.11, and - 37.55 kcal/mol against FGFR4, respectively. The mechanisms of binding were emphasized by hydrogen bonding and binding forces analysis as well as active site physicochemical profiling. The findings of this study identified that curcumin, syringolin A and andrographolide-but not eupalmerin acetate or deoxyelephantopin -could be viable dual EGFR and FGFR4 covalent irreversible inhibitors and could be implemented in HCC combination therapy protocols alone or in conjunction with other chemotherapeutic agents. Investigations of this study conclusively indicate dual blockade of EGFR and FGFR4 may be a promising future therapeutic strategy for enhanced management of HCC.

Cheminformatics identification of modulators of key carbohydrate-metabolizing enzymes from C. cujete for type-2 diabetes mellitus intervention

Purpose

The therapeutic use of oral hypoglycaemic agents in the management of type-2 diabetes mellitus (T2DM) is without adverse effects; thus, calls for alternative and novel candidates from natural products in medicinal plants.

Method

The study explored molecular docking and molecular dynamics (MD) simulation approaches to identify key antidiabetic metabolites from Crescentia cujete.

Results

Molecular docking results identified four and/or five best compounds against each target enzyme (alpha-glucosidase, dipeptidyl peptidase-IV, aldose reductase, and protein tyrosine phosphatase-1B (PTP-1B)) implicated in diabetes. The resulting complexes (except against PTP-1B) had higher docking scores above respective standards (acarbose, Diprotin A, ranirestat). The MD simulation results revealed compounds such as benzoic acid (-48.414 kcal/mol) and phytol (-45.112 kcal/mol) as well as chlorogenic acid (-42.978 kcal/mol) and naringenin (-31.292 kcal/mol) had higher binding affinities than the standards [acarbose (-28.248 kcal/mol), ranirestat (-21.042 kcal/mol)] against alpha-glucosidase and aldose reductase, respectively while Diprotin A (-45.112 kcal/mol) and ursolic acid (-18.740 kcal/mol) presented superior binding affinities than the compounds [luteolin (-41.957 kcal/mol and naringenin (-16.518 kcal/mol)] against DPP-IV and PTP-1B respectively.

Conclusion

While isoflavone (alpha-glucosidase), xylocaine (DPP-IV), luteolin (aldose reductase,) and chlorogenic acid (PTP-1B) were affirmed as the best inhibitors of respective enzyme targets, luteolin, and chlorogenic acid may be suggested and proposed as probable candidates against T2DM and related retinopathy complication based on their structural stability, compactness and affinity for three (DPP-IV, aldose reductase, and PTP-1B) of the four targets investigated. Further studies are warranted in vitro and in vivo on the antihyperglycaemic effects of these drug candidates.

Supplementary Information

The online version contains supplementary material available at 10.1007/s40200-023-01249-7.

Mechanistic insights into MARK4 inhibition by galantamine toward therapeutic targeting of Alzheimer's disease

Introduction: Hyperphosphorylation of tau is an important event in Alzheimer's disease (AD) pathogenesis, leading to the generation of "neurofibrillary tangles," a histopathological hallmark associated with the onset of AD and related tauopathies. Microtubule-affinity regulating kinase 4 (MARK4) is an evolutionarily conserved Ser-Thr (S/T) kinase that phosphorylates tau and microtubule-associated proteins, thus playing a critical role in AD pathology. The uncontrolled neuronal migration is attributed to overexpressed MARK4, leading to disruption in microtubule dynamics. Inhibiting MARK4 is an attractive strategy in AD therapeutics. Methods: Molecular docking was performed to see the interactions between MARK4 and galantamine (GLT). Furthermore, 250 ns molecular dynamic studies were performed to investigate the stability and conformational dynamics of the MARK4-GLT complex. We performed fluorescence binding and isothermal titration calorimetry studies to measure the binding affinity between GLT and MARK4. Finally, an enzyme inhibition assay was performed to measure the MARK4 activity in the presence and absence of GLT. Results: We showed that GLT, an acetylcholinesterase inhibitor, binds to the active site cavity of MARK4 with an appreciable binding affinity. Molecular dynamic simulation for 250 ns demonstrated the stability and conformational dynamics of the MARK4-GLT complex. Fluorescence binding and isothermal titration calorimetry studies suggested a strong binding affinity. We further show that GLT inhibits the kinase activity of MARK4 significantly (IC50 = 5.87 µM). Conclusion: These results suggest that GLT is a potential inhibitor of MARK4 and could be a promising therapeutic target for AD. GLT's inhibition of MARK4 provides newer insights into the mechanism of GLT's action, which is already used to improve cognition in AD patients.

In silico exploration and molecular dynamics of deleterious SNPs on the human TERF1 protein triggering male infertility

By limiting chromosome erosion and end-to-end fusions, telomere integrity is critical for chromosome stability and cell survival. During mitotic cycles or due to environmental stresses, telomeres become progressively shorter and dysfunctional, thus triggering cellular senescence, genomic instability and cell death. To avoid such consequences, the telomerase action, as well as the Shelterin and CST complexes, assure the telomere's protection. Telomeric repeat binding factor 1 (TERF1), which is one of the primary components of the Shelterin complex, binds directly to the telomere and controls its length and function by regulating the telomerase activity. Several reports about TERF1 gene variations have been associated with different diseases, and some of them have linked these variations to male infertility. Hence, this paper can be advantageous to investigate the association between the missense variants of the TERF1 gene and the susceptibility to male infertility. The stepwise prediction of SNPs pathogenicity followed in this study was based on stability and conservation analysis, post-translational modification, secondary structure, functional interaction prediction, binding energy evaluation and finally molecular dynamic simulation. Prediction matching among the tools revealed that out of 18 SNPs, only four (rs1486407144, rs1259659354, rs1257022048 and rs1320180267) were predicted as the most damaging and highly deleterious SNPs affecting the TERF1 protein and its molecular dynamics when interacting with the TERB1 protein by influencing the function, structural stability, flexibility and compaction of the overall complex. Interestingly, these polymorphisms should be considered during genetic screening so they can be used effectively as genetic biomarkers for male infertility diagnosis.Communicated by Ramaswamy H. Sarma.

Cheminformatics Identification and Validation of Dipeptidyl Peptidase-IV Modulators from Shikimate Pathway-Derived Phenolic Acids towards Interventive Type-2 Diabetes Therapy

Recently, dipeptidyl peptidase-IV (DPP-IV) has become an effective target in the management of type-2 diabetes mellitus (T2D). The study aimed to determine the efficacy of shikimate pathway-derived phenolic acids as potential DPP-IV modulators in the management of T2D. The study explored in silico (molecular docking and dynamics simulations) and in vitro (DPP-IV inhibitory and kinetics assays) approaches. Molecular docking findings revealed chlorogenic acid (CA) among the examined 22 phenolic acids with the highest negative binding energy (−9.0 kcal/mol) showing a greater affinity for DPP-IV relative to the standard, Diprotin A (−6.6 kcal/mol). The result was corroborated by MD simulation where it had a higher affinity (−27.58 kcal/mol) forming a more stable complex with DPP-IV than Diprotin A (−12.68 kcal/mol). These findings were consistent with in vitro investigation where it uncompetitively inhibited DPP-IV having a lower IC₅₀ (0.3 mg/mL) compared to Diprotin A (0.5 mg/mL). While CA showed promising results as a DPP-IV inhibitor, the findings from the study highlighted the significance of medicinal plants particularly shikimate-derived phenolic compounds as potential alternatives to synthetic drugs in the effective management of T2DM. Further studies, such as derivatisation for enhanced activity and in vivo evaluation are suggested to realize its full potential in T2D therapy.

Triazole Derivatives Target 14α-Demethylase (LDM) Enzyme in Candida albicans Causing Ergosterol Biosynthesis Inhibition

Candida albicans is the most dominant and prevalent cause of fungal infections in humans. Azoles are considered as first-line drugs for the treatment of these infections. However, their prolonged and insistent use has led to multidrug resistance and treatment failures. To overcome this, modification or derivatization of the azole ring has led to the development of new and effective antifungal molecules. In a previous study, we reported on the development of new triazole-based molecules as potential antifungal agents against Candida auris. In this study, the most potent molecules from the previous study were docked and simulated with lanosterol 14-alpha demethylase enzyme. These compounds were further evaluated for in vitro susceptibility testing against C. albicans. In silico results revealed favorable structural dynamics of the compounds, implying that the compounds would be able to effectively bind to the target enzyme, which was further manifested by the strong interaction of the test compounds with the amino acid residues of the target enzyme. In vitro studies targeting quantification of ergosterol content revealed that pta1 was the most active compound and inhibited ergosterol production by >90% in both drug-susceptible and resistant C. albicans isolates. Furthermore, RT-qPCR results revealed downregulation of ERG11 gene when C. albicans cells were treated with the test compound, which aligns with the decreased ergosterol content. In addition, the active triazole derivatives were also found to be potent inhibitors of biofilm formation. Both in silico and in vitro results indicate that these triazole derivatives have the potential to be taken to the next level of antifungal drug development.

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 (KRAS^G12C) 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 KRAS^G12C. It was revealed that AMG510 and MRTX849 complexes presented similar total free binding energies, (ΔG_bind), 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.

Structure-based identification of potential SARS-CoV-2 main protease inhibitors

To address coronavirus disease (COVID-19), currently, no effective drug or vaccine is available. In this regard, molecular modeling approaches are highly useful to discover potential inhibitors of the main protease (Mpro) enzyme of SARS-CoV-2. Since, the Mpro enzyme plays key roles in mediating viral replication and transcription; therefore, it is considered as an attractive drug target to control SARS-CoV-2 infection. By using structure-based drug design, pharmacophore modeling, and virtual high throughput drug screening combined with docking and all-atom molecular dynamics simulation approach, we have identified five potential inhibitors of SARS-CoV-2 Mpro. MD simulation studies revealed that compound 54035018 binds to the Mpro with high affinity (ΔGbind -37.40 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. We have identified promising leads to fight COVID-19 infection as these compounds fulfill all drug-likeness properties. However, experimental and clinical validations are required for COVID-19 therapy.Communicated by Ramaswamy H. Sarma.

Deciphering the Mechanism of Binding Selectivity of Chlorofluoroacetamide-Based Covalent Inhibitors toward L858R/T790M Resistance Mutation

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 (EGFR^DM) using α-chlorofluoroacetamide (CFA) as a novel warhead, has seven times the inhibition selectivity for EGFR^DM over the wild type (EGFR^WT), 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 EGFR^DM undergoes significantly larger conformational changes than EGFR^WT upon binding to 18. The conformational stability of the diamine side chain and the CFA motif of 18 in the orthosteric site of EGFR^DM is identified as key for the disparate binding mechanism and inhibitory prowess of 18 with respect to EGFR^WT and EGFR^DM 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.

Identification of Flavonoid C-Glycosides as Promising Antidiabetics Targeting Protein Tyrosine Phosphatase 1B

Protein tyrosine phosphatase 1B (PTP1B), a negative regulator of the insulin signaling pathway, has gained attention as a validated druggable target in the management of type 2 diabetes mellitus (T2DM). The lack of clinically approved PTP1B inhibitors has continued to prompt research in plant-derived therapeutics possibly due to their relatively lesser toxicity profiles. Flavonoid C-glycosides are one of the plant-derived metabolites gaining increased relevance as antidiabetic agents, but their possible mechanism of action remains largely unknown. This study investigates the antidiabetic potential of flavonoid C-glycosides against PTP1B in silico and in vitro. Of the seven flavonoid C-glycosides docked against the enzyme, three compounds (apigenin, vitexin, and orientin) had the best affinity for the enzyme with a binding score of -7.3 kcal/mol each, relative to -7.4 kcal/mol for the reference standard, ursolic acid. A further probe (in terms of stability, flexibility, and compactness) of the complexes over a molecular dynamics time study of 100 ns for the three compounds suggested orientin as the most outstanding inhibitor of PTP1B owing to its overall -34.47 kcal/mol binding energy score compared to ursolic acid (-19.24 kcal/mol). This observation was in accordance with the in vitro evaluation result, where orientin had a half maximal inhibitory concentration (IC₅₀) of 0.18 mg/ml relative to 0.13 mg/ml for the reference standard. The kinetics of inhibition of PTP1B by orientin was mixed-type with V _max and K _m values of 0.004 μM/s and 0.515 μM. Put together, the results suggest orientin as a potential PTP1B inhibitor and could therefore be further explored in the management T2DM as a promising therapeutic agent.

Myricetin inhibits breast and lung cancer cells proliferation via inhibiting MARK4

Identifying novel molecules as potential kinase inhibitors are gaining significant attention globally. The present study suggests Myricetin as a potential inhibitor of microtubule-affinity regulating kinase (MARK4), adding another molecule to the existing list of anticancer therapeutics. MARK4 regulates initial cell division steps and is a potent druggable target for various cancers. Structure-based docking with 100 ns molecular dynamics simulation depicted activity of Myricetin in the active site pocket of MARK4 and the formation of a stable complex. The fluorescence-based assay showed excellent affinity of Myricetin to MARK4 guided by static and dynamic quenching. Moreover, the assessment of enthalpy change (∆H) and entropy change (∆S) delineated electrostatic interactions as a dominant force in the MARK4-myricetin interaction. Isothermal titration calorimetric measurements revealed spontaneous binding of Myricetin with MARK4. Further, the kinase assay depicted significant inhibition of MARK4 by Myricetin with IC₅₀  = 3.11 µM. Additionally, cell proliferation studies established that Myricetin significantly inhibited the cancer cells (MCF-7 and A549) proliferation, and inducing apoptosis. This study provides a solid rationale for developing Myricetin as a promising anticancer molecule in the MARK4 mediated malignancies.

Investigating single amino acid substitutions in PIM1 kinase: A structural genomics approach

PIM1, is a serine/threonine proto-oncogene kinase, involved in many biological functions, including cell survival, proliferation, and differentiation, thus play a key role in oncogenesis. It plays a crucial role in the onset and progression of various hematopoietic and non-hematopoietic malignancies, including acute myeloid leukemia and prostate cancer. Mutations in PIM1, especially in its kinase domain, can induce abnormal structural changes and thus alter functionalities that can lead to disease progression and other complexities. Herein, we have performed an extensive analysis of the PIM1 mutations at sequence and structure level while utilizing state-of-the-art computational approaches. Based on the impact on PIM1, numerous pathogenic and destabilizing mutations were identified and subsequently analyzed in detail. Finally, two amino acid substitutions (W109C and F147C) in the kinase domain of PIM1 were selected to explore their impact on the PIM1 structure in a time evolution manner using all-atom molecular dynamics (MD) simulations for 200 ns. MD results indicate significant conformational altercations in the structure of PIM1, especially upon F147C mutation. This study provides a significant insight into the PIM1 dysfunction upon single amino acid substitutions, which can be utilized to get insights into the molecular basis of PIM1-associated disease progression.

MAP/Microtubule Affinity Regulating Kinase 4 Inhibitory Potential of Irisin: A New Therapeutic Strategy to Combat Cancer and Alzheimer's Disease

Irisin is a clinically significant protein playing a valuable role in regulating various diseases. Irisin attenuates synaptic and memory dysfunction, highlighting its importance in Alzheimer's disease. On the other hand, Microtubule Affinity Regulating Kinase 4 (MARK4) is associated with various cancer types, uncontrolled neuronal migrations, and disrupted microtubule dynamics. In addition, MARK4 has been explored as a potential drug target for cancer and Alzheimer's disease therapy. Here, we studied the binding and subsequent inhibition of MARK4 by irisin. Irisin binds to MARK4 with an admirable affinity (K = 0.8 × 107 M-1), subsequently inhibiting its activity (IC50 = 2.71 µm). In vitro studies were further validated by docking and simulations. Molecular docking revealed several hydrogen bonds between irisin and MARK4, including critical residues, Lys38, Val40, and Ser134. Furthermore, the molecular dynamic simulation showed that the binding of irisin resulted in enhanced stability of MARK4. This study provides a rationale to use irisin as a therapeutic agent to treat MARK4-associated diseases.

Impact of single amino acid substitution on the structure and function of TANK-binding kinase-1

Tank-binding kinase 1 (TBK1) is a serine/threonine protein kinase involved in various signaling pathways and subsequently regulates cell proliferation, apoptosis, autophagy, antiviral and antitumor immunity. Dysfunction of TBK1 can cause many complex diseases, including autoimmunity, neurodegeneration, and cancer. This dysfunction of TBK1 may result from single amino acid substitutions and subsequent structural alterations. This study analyzed the effect of substituting amino acids on TBK1 structure, function, and subsequent disease using advanced computational methods and various tools. In the initial assessment, a total of 467 mutations were found to be deleterious. After that, in detailed structural and sequential analyses, 13 mutations were found to be pathogenic. Finally, based on the functional importance, two variants (K38D and S172A) of the TBK1 kinase domain were selected and studied in detail by utilizing all-atom molecular dynamics (MD) simulation for 200 ns. MD simulation, including correlation matrix and principal component analysis, helps to get deeper insights into the TBK1 structure at the atomic level. We observed a substantial change in variants' conformation, which may be possible for structural alteration and subsequent TBK1 dysfunction. However, substitution S172A shows a significant conformational change in TBK1 structure as compared to K38D. Thus, this study provides a structural basis to understand the effect of mutations on the kinase domain of TBK1 and its function associated with disease progression.

Structure-based investigation of MARK4 inhibitory potential of Naringenin for therapeutic management of cancer and neurodegenerative diseases

MAP/microtubule affinity-regulating kinase 4 (MARK4) is a member of serine/threonine kinase family and considered an attractive drug target for many diseases. Screening of Indian Medicinal Plants, Phytochemistry, and Therapeutics (IMPPAT) using virtual high-throughput screening coupled with enzyme assay suggested that Naringenin (NAG) could be a potent inhibitor of MARK4. Structure-based molecular docking analysis showed that NAG binds to the critical residues found in the active site pocket of MARK4. Furthermore, molecular dynamics (MD) simulation studies for 100 ns have delineated the binding mechanism of NAG to MARK4. Results of MD simulation suggested that binding of NAG further stabilizes the structure of MARK4 by forming a stable complex. In addition, no significant conformational change in the MARK4 structure was observed. Fluorescence binding and isothermal titration calorimetric measurements revealed an excellent binding affinity of NAG to MARK4 with a binding constant (K) = 0.13 × 10⁶ M^-1 obtained from fluorescence binding studies. Further, enzyme inhibition studies showed that NAG has an admirable IC₅₀ value of 4.11 µM for MARK4. Together, these findings suggest that NAG could be an effective MARK4 inhibitor that can potentially be used to treat cancer and neurodegenerative diseases.

Impact of non-synonymous mutations on the structure and function of telomeric repeat binding factor 1

Telomeric repeat binding factor 1 (TRF1) is one of the major components of the shelterin complex. It directly binds to the telomere and controls its function by regulating the telomerase acting on it. Several variations are reported in the TRF1 gene; some are associated with variety of diseases. Here, we have studied the structural and functional significance of these variations in the TRFH domain of TRF1. We have used cutting-edge computational methods such as SIFT, PolyPhen-2, PROVEAN, Mutation Assessor, mCSM, SDM, STRUM, MAESTRO, and DUET to predict the effects of 124 mutations in the TRFH domain of TRF1. Out of 124 mutations, we have identified 12 deleterious mutations with high confidence based on their prediction. To see the impact of the finally selected mutations on the structure and stability of TRF1, all-atom molecular dynamics (MD) simulations on TRF1-Wild type (WT), L79R and P150R mutants for 200 ns were carried out. A significant conformational change in the structure of the P150R mutant was observed. Our integrated computational study provides a comprehensive understanding of structural changes in TRF1 incurred due to the mutations and subsequent function, leading to the progression of many diseases.Communicated by Ramaswamy H. Sarma.

Distinguishing the optimal binding mechanism through reversible and irreversible inhibition analysis of HSP72 protein in cancer therapy

Over the past two decades, covalent inhibitors have gained much interest and are living up to their reputation as a powerful tool in drug discovery. Covalent inhibitors possess several significant advantages, including increased biochemical efficiency, prolonged duration and the ability to target shallow, solvent-exposed substrate-binding domains. One of the enzymes that have been both covalently and non-covalently targeted is the heat shock protein 72 (HSP72). This elevated enzyme expression in cancer cells may be responsible for tumorigenesis and tumor progression by providing chemotherapy resistance. A critical gap remains in the molecular understanding of the structural mechanism's covalent and non-covalent binding to HSP72. In this study, we explore the most optimal binding mechanism in the inhibition of the HSP72. Based on the molecular dynamic analyses, it was evident that the non-covalent complex showed more stability than the covalent complex. The covalent ligand, however, was more able to induce and stabilize the sealed conformation of the HSP72-NBD ATP binding domain throughout the. Also, the non-covalent ligand does not induce any significant conformational change as it remained close to the shape of the unbound complex; and the affinity is only dependent on the multiple hydrogen bonds in contrast to the covalent ligand. This is supported by the secondary structure elements and principal component analysis that was more dominant in the covalently inhibited complex. Covalent bond induced the α-helices sealed conformation of the HSP72-NBD; based on our findings, this will prevent other small molecules from interacting at the ATP binding site domain. Moreover, inhibition of the ATP binding domain can directly affect the ATPs protein folding mechanism of the HSP72 enzyme. The essential dynamic analysis presented in this report compliments the binding mechanism of HSP72, establishing covalent inhibition as the preferred method of inhibiting the HSP72 protein. The findings from this study may assist in the design of more target-specific HSP72 covalent inhibitors exploring the surface-exposed lysine residues.

Investigation of sphingosine kinase 1 inhibitory potential of cinchonine and colcemid targeting anticancer therapy

Sphingosine kinase 1 (SphK1) and sphingosine-1-phosphate (S1P) signaling regulates numerous diseases such as cancer, diabetes, and inflammation-related ailments, rheumatoid arthritis, atherosclerosis, and multiple sclerosis. The importance of SphK1 in chemo-resistance has been extensively explored in breast, lung, colon, and hepatocellular carcinomas. SphK1 is considered an attractive drug target for the development of anticancer therapy. New drug molecules targeting the S1P signaling are required owing to its pleiotropic nature and association with multiple downstream targets. Here, we have investigated the binding affinity and SphK1 inhibitory potential of cinchonine and colcemid using a combined molecular docking and simulation studies followed by experimental analysis. These compounds bind to SphK1 with a significantly high affinity and subsequently inhibit kinase activity (IC₅₀ 7-9 μM). Further, MD simulation studies revealed that both cinchonine and colcemid bind to the residues at the active site pocket of SphK1 with several non-covalent interactions, which may be responsible for inhibiting its kinase activity. Besides, the binding of cinchonine and colcemid causes substantial conformational changes in the structure of SphK1. Taken together, cinchonine and colcemid may be implicated in designing potential drug molecules with improved affinity and specificity for SphK1 targeting anticancer therapy.Communicated by Ramaswamy H. Sarma.

Targeting ebola virus VP40 protein through novel inhibitors: exploring the structural and dynamic perspectives on molecular landscapes

Ebola filovirus (EBOV) is one of the deadliest known infectious agents, and a cause of Western African epidemics from 2013 to 2016. The virus has infected nearly 3000 humans and almost 1900 have died. In the past few years, various small molecules have been discovered to display efficiency against EBOV and some of them have progressed towards clinical trials. Even though continuous attempts have been made to find antiEBOV therapeutics, no potential drugs are yet approved against this viral infection. The development of small antiviral inhibitors has gained tremendous attention in the attempt to overcome EVD. With this background, we seek to offer molecular insights into EBOV VP40 protein inhibition, using all atom molecular mechanics methodology and binding free energy calculations. We have selected five novel reported inhibitors against VP40 protein, namely Comp1, Comp2, Comp3, Comp4, and Comp5, and explored their binding against the same target. It was evident from the analysis that all the inhibitors displayed stability in complex with VP40 protein; however, Comp1 exhibited enhanced stability and compactness. Comp1 unveiled favorable binding, which accounted for positive correlation motions in the active site residues. Likewise, Comp1 revealed the most promising binding (ΔG_bind - 40.3504 kcal/mol) as compared to the other four inhibitors, which disclosed relatively less favorable ΔG_bind. The highest binding energy of Comp1 to VP40 protein can be primarily endorsed to the upsurge in van der Waals energy by ΔE_vdW - 37.1609 kcal/mol and Coulomb energy by ΔE_ele - 52.7332 kcal/mol. Also, the hydrogen bond network is robust in Comp1-VP40 complex, with four hydrogen bonds, whilst it is less in other inhibitors. The outcomes from this report may assist in the advancement of novel VP40 inhibitors with high selectivity and potency for EVD therapeutics.

ABBV-744 as a potential inhibitor of SARS-CoV-2 main protease enzyme against COVID-19

A new pathogen severe acute respiratory syndrome coronavirus 2 (SARS-CoV-2) has spread worldwide and become pandemic with thousands new deaths and infected cases globally. To address coronavirus disease (COVID-19), currently no effective drug or vaccine is available. This necessity motivated us to explore potential lead compounds by considering drug repurposing approach targeting main protease (M^pro) enzyme of SARS-CoV-2. This enzyme considered to be an attractive drug target as it contributes significantly in mediating viral replication and transcription. Herein, comprehensive computational investigations were performed to identify potential inhibitors of SARS-CoV-2 M^pro enzyme. The structure-based pharmacophore modeling was developed based on the co-crystallized structure of the enzyme with its biological active inhibitor. The generated hypotheses were applied for virtual screening based PhaseScore. Docking based virtual screening workflow was used to generate hit compounds using HTVS, SP and XP based Glide GScore. The pharmacological and physicochemical properties of the selected lead compounds were characterized using ADMET. Molecular dynamics simulations were performed to explore the binding affinities of the considered lead compounds. Binding energies revealed that compound ABBV-744 binds to the M^pro with strong affinity (ΔG_bind -45.43 kcal/mol), and the complex is more stable in comparison with other protein-ligand complexes. Our study classified three best compounds which could be considered as promising inhibitors against main protease SARS-CoV-2 virus.

Weak spots inhibition in the Mycobacterium tuberculosis antigen 85C target for antitubercular drug design through selective irreversible covalent inhibitor-SER124

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 (CyC_8β) 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 CyC_8β 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 CyC_8β 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 3₁₀-helices in the bound systems. SASA of CyC_8β 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.

Imatinib exhibit synergistic pleiotropy in the prevention of colorectal cancer by suppressing proinflammatory, cell survival and angiogenic signaling

Recent global incidences and mortality rates have placed colorectal cancer (CRC) at third and second positions, respectively, among both sexes of all ages. Resistance during chemotherapy is a big problem in the treatment and disease-free survival of CRC patients. Discovery of new anticancer drug(s) is a time taking process and therefore, invites studies for repurposing the known therapeutics. The present study was conceived to analyze the anticancer role of Imatinib in experimental CRC at early stages. Different experimental procedures e.g. tumor incidences or histoarchitectural changes, gene and protein expression analysis, estimations of intracellular calcium, ROS, mitochondrial membrane potential, apoptotic index and molecular docking was performed to support the hypothesis. It was observed that Imatinib could function as an immunomodulator by breaking the feed-back loop between the proinflammatory cytokines (IL-1β and TNF-α) and transcription factors (NF-κB, Jak3/Stat3) knowingly involved in increased cell proliferation during tumorigenesis via activating different intracellular signaling. Also, Imatinib could independently deregulate the other cell survival and proliferation signaling e.g. PI3-K/Akt/mTOR, Wnt/β-catenin and MAPK. Proinflammatory cytokines orchestrated intracellular signaling also involve angiogenic factors to be upregulated during CRC which were also seemed to be independently suppressed by Imatinib. Restoration of physiological apoptosis by increasing the release of intracellular calcium to generate ROS thereby reducing the mitochondrial membrane potential for the release of cytochrome c and activation of caspase-3 was also reported with Imatinib administration. Thus, it may be suggested that Imatinib show synergistic pleiotropy in suppressing the interlinked tumorigenic signaling pathways independently.

Coupling of HSP72 α-Helix Subdomains by the Unexpected Irreversible Targeting of Lysine-56 over Cysteine-17; Coevolution of Covalent Bonding

Covalent inhibition has recently gained a resurgence of interest in several drug discovery areas. The expansion of this approach is based on evidence elucidating the selectivity and potency of covalent inhibitors when bound to particular amino acids of a biological target. The unexpected covalent inhibition of heat shock protein 72 (HSP72) by covalently targeting Lys-56 instead of Cys-17 was an interesting observation. However, the structural basis and conformational changes associated with this preferential coupling to Lys-56 over Cys-17 remain unclear. To resolve this mystery, we employed structural and dynamic analyses to investigate the structural basis and conformational dynamics associated with the unexpected covalent inhibition. Our analyses reveal that the coupling of the irreversible inhibitor to Lys-56 is intrinsically less dynamic than Cys-17. Conformational dynamics analyses further reveal that the coupling of the inhibitor to Lys-56 induced a closed conformation of the nucleotide-binding subdomain (NBD) α-helices, in contrast, an open conformation was observed in the case of Cys-17. The closed conformation maintained the crucial salt-bridge between Glu-268 and Lys-56 residues, which strengthens the interaction affinity of the inhibitor nearly identical to adenosine triphosphate (ADP/Pi) bound to the HSP72-NBD. The outcome of this report provides a substantial shift in the conventional direction for the design of more potent covalent inhibitors.

Identification of Potential Inhibitors of Calcium/Calmodulin-Dependent Protein Kinase IV from Bioactive Phytoconstituents

Calcium/calmodulin-dependent protein kinase IV (CaMKIV) is an upstream regulator of CaMKK-CaMKIV signaling cascade that activates various transcription factors, thereby regulating several cellular activities including, neuronal communication and immune response. Owing to the abnormal expression in cancer and neurodegenerative diseases, the CaMKIV has been considered a potential drug target. In the present study, we checked the binding affinity of plant-derived natural compounds viz., quercetin, ellagic acid (EA), simvastatin, capsaicin, ursolic acid, DL-α-tocopherol acetate, and limonin towards CaMKIV. Molecular docking and fluorescence binding studies showed that EA and quercetin bind to the CaMKIV with a considerable affinity in comparison to other compounds. Enzyme inhibition assay revealed that both EA and quercetin inhibit CaMKIV activity with their IC₅₀ values in the micromolar range. To get atomistic insights into the mode of interactions, inhibition mechanism, and the stability of the CaMKIV-ligand complex, a 100 ns MD simulation analysis was performed. Both EA and quercetin bind to the catalytically important residues of active site pocket of CaMKIV forming enough stabilizing interactions presumably inhibiting enzyme activity. Moreover, no significant structural change in the CaMKIV was observed upon binding of EA and quercetin. In conclusion, this study illustrates the application of phytoconstituents in the development of therapeutic molecules targeting CaMKIV having implications in cancer and neurodegenerative diseases after in vivo validation.

A Mechanistic Study of a Potent and Selective Epidermal Growth Factor Receptor Inhibitor against the L858R/T790M Resistance Mutation

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 (EGFR^DM) relative to the EGFR wild-type form (EGFR^WT). The enhanced therapeutic potency of 20g against EGFR^DM is 263 times greater than that against EGFR^WT, which necessitates a rational explanation for the underlying selective and inhibitory mechanisms. In this work, we investigate the differential binding modes of 20g with EGFR^WT and EGFR^DM 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 EGFR^WT and EGFR^DM. The calculated binding free energies show good correlation with the experimental biological activity. The total binding free energy difference between EGFR^WT-20g and EGFR^DM-20g is -11.47 kcal/mol, implying that 20g binds more strongly to EGFR^DM. This enhanced binding affinity of 20g for EGFR^DM 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 EGFR^DM relative to EGFR^WT.

From genomic variation to protein aberration: Mutational analysis of single nucleotide polymorphism present in ULBP6 gene and implication in immune response

BACKGROUND

Genetic polymorphisms have been identified as one of the underlying factors in disease pathogenesis and drug resistance since they account for protein dysfunctionality, or in some cases, aberrancy. This explains the high degree of inactivity that characterizes the polymorphic variants of ULBP6 binding protein, which in turn disrupts its primary interaction with human Natural Killer Group 2-member D (NKG2D) and accounts for an impediment to immuno-surveillance. The possible identification of deleterious non-synonymous Single Nucleotide Polymorphisms (nsSNPs) present in the ULBP6 gene is essential for the development of novel gene therapies to prevent the translation of dysfunctional protein variants.

METHODS/RESULTS

In this study, for the first time, we employed an SNP-informatics approach (SNPs retrieval, pathogenic/mutational analysis, phenotypic analysis, and structural analysis) and molecular dynamics techniques to identify and characterize undesirable SNPs coupled with their impact on ULBP6 structural activities relative to dysfunctionality. V52F was predictively pathogenic amongst SNPs studied. Conformational and dynamic studies revealed that in comparison to wildtype ULBP6 (ULBP6^wt), pathogenic ULBP6^V52F demonstrated considerable structural inactivity, which could, in turn, impede biological protein-protein interactions. Moreover, ULBP6^V52F showed relatively limited motions in the conformational space as deduced from estimations of structural stability, fluctuations, and principal components.

CONCLUSION

This study provides a workable paradigm for investigating pathological nsSNPS using computational platforms which findings present ULBP6^V52F as a novel and attractive immunotherapeutic target in combatting immune-associated disorders.

Covalent vs. Non-Covalent Inhibition: Tackling Drug Resistance in EGFR - A Thorough Dynamic Perspective

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 (EGFR^DM ) was found to be the most alarming variant. Despite the development of a wide range of inhibitors, none of them could inhibit EGFR^DM effectively. Recently, 11h and 45a, have been found to be potent inhibitors against EGFR^DM 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 (ΔΔG_bind =-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 EGFR^DM by stabilizing its flexible alcohol chain. The 45a binding to EGFR^DM 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.

Exploring the Structural Mechanism of Covalently Bound E3 Ubiquitin Ligase: Catalytic or Allosteric Inhibition?

Covalent inhibition has recently gained a resurgence of interest in several drug discovery areas. The expansion of this approach is based on evidence elucidating the selectivity and potency of covalent inhibitors when bound to particular amino acids of a biological target. The Nedd4-1, an E3 ubiquitin ligase, is characterized by two covalent binding sites, of which catalytic Cys^cat and allosteric Cys^allo are enclosed. This enzyme has demonstrated inhibition at both the above-mentioned binding sites; however, a detailed molecular understanding of the structural mechanism of inhibition upon Cys^cat and Cys^allo binding remains vague. This prompted us to provide the first account of investigating the preferential covalent binding mode and the underlying structural and molecular dynamic implications. Based on the molecular dynamic analyses, it was evident that although both catalytic and allosteric covalent binding led to greater stability of the enzyme, a preferential covalent mechanism of inhibition was seen in the allosteric-targeted system. This was supported by a more favorable binding energy in the allosteric site compared to the catalytic site, in addition to the larger number of residue interactions and stabilizing hydrogen bonds occurring in the allosteric covalent bound complex. The fundamental dynamic analysis presented in this report compliments, as well as adds to previous experimental findings, thus leading to a crucial understanding of the structural mechanism by which Nedd4-1 is inhibited. The findings from this study may assist in the design of more target-specific Nedd4-1 covalent inhibitors exploring the surface-exposed cysteine residues.

The irony of chirality – unveiling the distinct mechanistic binding and activities of 1-(3-(4-amino-5-(7-methoxy-5-methylbenzo[b]thiophen-2-yl)-7H-pyrrolo[2,3-d]pyrimidin-7-yl)pyrrolidin-1-yl)prop-2-en-1-one enantiomers as irreversible covalent FGFR4 inhibitors

Theoretical investigation of the effect of chirality on inhibitors is providing essential insights for drug design.

Targeting HCV polymerase: a structural and dynamic perspective into the mechanism of selective covalent inhibition

Background: Concerns have been raised over the emerging pandemic status of hepatitis C virus (HCV). Current available drugs lack specificity, stability and potency. The HCV NS5B RNA-dependent RNA polymerase (RdRp) is a vital component in viral replication and is often targeted in antiviral therapies. Recent experimental procedures have led to the discovery of a novel covalent RdRp inhibitor, compound 47, which selectively targets cysteine 366 of the HCV RdRp and exhibits promising pharmacokinetic outcomes. Selective covalent inhibition of HCV is, however, a highly neglected subject in the literature, that is reinforced by the lack of efficient structure-based drug design protocols. In this paper, an atomistic insight into a novel selective approach to inhibit HCV RdRp is provided. Methodology/Results: Covalent molecular dynamic analyses revealed the inhibitory mechanism of compound 47 on the RdRp. Inhibitor binding induced distinctive internal movements resulting in the disruption of normal physiological interdomain interactions. Conclusion: Compound 47 stimulates reorganization of key protein elements required for RNA transcription, thus hampering viral replication as well as disrupting the overall conformation of HCV. This study will open new lines of approach for the design of novel selective inhibitors against HCV as well as other viral families.

Covalent simulations of covalent/irreversible enzyme inhibition in drug discovery: a reliable technical protocol

Aim: Irreversible covalent inhibition of biological targets in disease pathogenesis is an emerging field in drug design. Computational techniques have assumed a critical role in understanding covalent enzyme inhibition. However, a gap currently exists with regards to the reliability and reproducibility of currently available protocols available in literature and open scientific forums. Methodology/results: Appropriate ligand and protein target are selected, docked covalently or noncovalently using respective docking tools. Both components are subjected to premolecular dynamic preparations. This was followed by parameterization of the ligand, protein and covalent complex, respectively. The production runs were initiated and the resulting trajectories are saved and analyzed. Conclusion: This protocol is reliable and reproducible, hence would advance the development of irreversible covalent inhibitors toward disease treatment.

An update on the discovery and development of selective heat shock protein inhibitors as anti-cancer therapy

INTRODUCTION

Over the years, not a single HSP inhibitor has progressed into the post-market phase of drug development despite the success recorded in various pre-clinical and clinical studies. The inability of existing drugs to specifically target oncogenic HSPs has majorly accounted for these setbacks. Recent combinatorial strategies that incorporated computer-aided drug design (CADD) techniques are geared towards the development of highly specific HSP inhibitors with increased activities and minimal toxicities. Areas covered: In this review, strategic therapeutic approaches that have recently aided the development of selective HSP inhibitors were highlighted. Also, the significant contributions of CADD techniques over the years were discussed in detail. This article further describes promising computational paradigms and their applications towards the discovery of highly specific inhibitors of oncogenic HSPs. Expert opinion: The recent shift towards highly selective and specific HSP inhibition has shown great promise as evidenced by the development of paralog/isoform-selective HSP drugs. It could be further augmented with computer-aided drug design strategies, which incorporate reliable methods that would greatly enhance the design and optimization of novel inhibitors with improved activities and minimal toxicities.