Kinetics, in silico docking, molecular dynamics, and MM-GBSA binding studies on prototype indirubins, KT5720, and staurosporine as phosphorylase kinase ATP-binding site inhibitors: the role of water molecules examined

Evidence for the Quercetin Binding Site of Glycogen Phosphorylase as a Target for Liver-Isoform-Selective Inhibitors against Glioblastoma: Investigation of Flavanols Epigallocatechin Gallate and Epigallocatechin

Glycogen phosphorylase (GP) is the rate-determining enzyme in glycogenolysis, and its druggability has been extensively studied over the years for the development of therapeutics against type 2 diabetes (T2D) and, more recently, cancer. However, the conservation of binding sites between the liver and muscle isoforms makes the inhibitor selectivity challenging. Using a combination of kinetic, crystallographic, modeling, and cellular studies, we have probed the binding of dietary flavonoids epigallocatechin gallate (EGCG) and epigallocatechin (EGC) to GP isoforms. The structures of rmGPb-EGCG and rmGPb-EGC complexes were determined by X-ray crystallography, showing binding at the quercetin binding site (QBS) in agreement with kinetic studies that revealed both compounds as noncompetitive inhibitors of GP, with EGCG also causing a significant reduction in cell viability and migration of U87-MG glioblastoma cells. Interestingly, EGCG exhibits different binding modes to GP isoforms, revealing QBS as a promising site for GP targeting, offering new opportunities for the design of liver-selective GP inhibitors.

Phosphorylation-driven epichaperome assembly is a regulator of cellular adaptability and proliferation

The intricate network of protein-chaperone interactions is crucial for maintaining cellular function. Recent discoveries have unveiled the existence of specialized chaperone assemblies, known as epichaperomes, which serve as scaffolding platforms that orchestrate the reconfiguration of protein-protein interaction networks, thereby enhancing cellular adaptability and proliferation. This study explores the structural and regulatory aspects of epichaperomes, with a particular focus on the role of post-translational modifications (PTMs) in their formation and function. A key finding is the identification of specific PTMs on HSP90, particularly at residues Ser226 and Ser255 within an intrinsically disordered region, as critical determinants of epichaperome assembly. Our data demonstrate that phosphorylation of these serine residues enhances HSP90's interactions with other chaperones and co-chaperones, creating a microenvironment conducive to epichaperome formation. Moreover, we establish a direct link between epichaperome function and cellular physiology, particularly in contexts where robust proliferation and adaptive behavior are essential, such as in cancer and pluripotent stem cell maintenance. These findings not only provide mechanistic insights but also hold promise for the development of novel therapeutic strategies targeting chaperone assemblies in diseases characterized by epichaperome dysregulation, thereby bridging the gap between fundamental research and precision medicine.

A computational odyssey: uncovering classical β-lactamase inhibitors in dry fruits

In the antibacterial arsenal, β-lactams have held a prominent position, but increasing resistance due to unauthorized use and genetic factors requires new strategies. Combining β-lactamase inhibitors with broad-spectrum β-lactams proves effective in combating this resistance. ESBL producers demand new inhibitors, leading to the exploration of plant-derived secondary metabolites for potent β-lactam antibiotics or alternative inhibitors. Using virtual screening, molecular docking, ADMET analysis, and molecular dynamic simulation, this study actively analyzed the inhibitory activity of figs, cashews, walnuts, and peanuts against SHV-1, NDM-1, KPC-2, and OXA-48 β-lactamases. Using AutoDock Vina, the docking affinities of various compounds for target enzymes were initially screened, revealing 12 bioactive compounds with higher affinities for the target enzymes compared to Avibactam and Tazobactam. Top-scoring metabolites, including Oleanolic acid, Protocatechuic acid, and Tannin, were subjected to MD simulation studies to further analyze the stability of the docked complexes using WebGro. The simulation coordinates, in terms of RMSD, RMSF, SASA, Rg, and hydrogen bonds formed, showed that these phytocompounds are stable enough to retain in the active sites at various orientations. The PCA and FEL analysis also showed the stability of the dynamic motion of Cα residues of phytochemical-bound enzymes. The pharmacokinetic analysis of the top phytochemicals was performed to analyze their bioavailability and toxicity. This study provides new insights into the therapeutic potential of phytochemicals of selected dry fruits and contributes to future experimental studies to identify βL inhibitors from plants.Communicated by Ramaswamy H. Sarma.

Multidisciplinary docking, kinetics and X-ray crystallography studies of baicalein acting as a glycogen phosphorylase inhibitor and determination of its' potential against glioblastoma in cellular models

Glycogen phosphorylase (GP) is the rate-determining enzyme in the glycogenolysis pathway. Glioblastoma (GBM) is amongst the most aggressive cancers of the central nervous system. The role of GP and glycogen metabolism in the context of cancer cell metabolic reprogramming is recognised, so that GP inhibitors may have potential treatment benefits. Here, baicalein (5,6,7-trihydroxyflavone) is studied as a GP inhibitor, and for its effects on glycogenolysis and GBM at the cellular level. The compound is revealed as a potent GP inhibitor against human brain GPa (Ki = 32.54 μM), human liver GPa (Ki = 8.77 μM) and rabbit muscle GPb (Ki = 5.66 μM) isoforms. It is also an effective inhibitor of glycogenolysis (IC50 = 119.6 μM), measured in HepG2 cells. Most significantly, baicalein demonstrated anti-cancer potential through concentration- and time-dependent decrease in cell viability for three GBM cell-lines (U-251 MG, U-87 MG, T98-G) with IC50 values of ∼20-55 μM (48- and 72-h). Its effectiveness against T98-G suggests potential against GBM with resistance to temozolomide (the first-line therapy) due to a positive O6-methylguanine-DNA methyltransferase (MGMT) status. The solved X-ray structure of rabbit muscle GP-baicalein complex will facilitate structure-based design of GP inhibitors. Further exploration of baicalein and other GP inhibitors with different isoform specificities against GBM is suggested.

3-deazaadenosine: A promising novel p38γ antagonist with potential as a breast cancer therapeutic agent

Human p38γ protein kinase, or MAPK12, is a crucial signaling protein that is important in channelizing membrane signals to the nucleus in the MAPK cascade pathway, associated with breast and colorectal cancer, besides other forms of malignancies and atherosclerotic lesions too. P38γ has a significant contribution to the progression of breast carcinoma due to its multifaceted functions. Targeting p38γ for defining potent antagonists against p38γ can turn out to be an attractive and novel means of breast cancer therapeutics. Novel and potent lead molecules were designed utilizing computational drug design methodologies. Using high-throughput virtual screening, 1909 geometrically similar analogs of known inhibitors were generated, primarily using BIRB796, SB202190, ANP, CHEBI: 620708, and CHEBI: 524699. Chemical correctness was ensured using LigPrep for the standalone library, and Prep Wizard for p38γ using Maestro v.11.5. Using the Glide v5.5 flexible docking procedure on a standalone library of p38γ binding sites, we defined 18 potential leads and assessed their ADMET properties. Lead "1", among the proposed four p38γ antagonists with high-scoring and favorable interactions, was considered for 100 ns molecular dynamics simulations. Among the four proposed leads, Lead '1' displayed consistent and stable bonding interactions with p38γ throughout the 100 ns molecular dynamics (MD) simulations. Additionally, it formed water bridges, contributing to its strong association with the protein. Notably, Lead '1' (3-deazaadenosine) exhibited favorable root-mean-square deviation (RMSD) and root-mean-square fluctuation (RMSF) within the acceptable range of pharmacological properties. Thus, 3-deazaadenosine and its mimetic might be promising new directions for developing a novel class of antagonists for breast cancer treatment.

In Silico-Motivated Discovery of Novel Potent Glycogen Synthase-3 Inhibitors: 1-(Alkyl/arylamino)-3H-naphtho[1,2,3-de]quinoline-2,7-dione Identified as a Scaffold for Kinase Inhibitor Development

Glycogen synthase kinase-3 (GSK-3) isoforms α and β have diverse roles within cell biology, and have been linked with multiple diseases that include prominent CNS conditions such as Alzheimer's disease and several psychiatric disorders. In this study, motivated by computation, we aimed to identify novel ATP-binding site inhibitors of GSK-3 with CNS-active potential. A ligand screening (docking) protocol against GSK-3β was first optimized, employing an active/decoy benchmarking set, with the final protocol selected based on statistical performance analysis. The optimized protocol involved pre-filtering of ligands using a three-point 3D-pharmacophore, followed by Glide-SP docking applying hinge region hydrogen bonding constraints. Using this approach, the Biogenic subset of the ZINC15 compound database was screened, focused on compounds with potential for CNS-activity. Twelve compounds (generation I) were selected for experimental validation using in vitro GSK-3β binding assays. Two hit compounds, 1 and 2, with 6-amino-7H-benzo[e]perimidin-7-one and 1-(phenylamino)-3H-naphtho[1,2,3-de]quinoline-2,7-dione type scaffolds were identified with IC₅₀ values of 1.63 µM and 20.55 µM, respectively. Ten analogues of 2 (generation II) were selected for structure activity relationship (SAR) analysis and revealed four low micromolar inhibitors (<10 µM), with 19 (IC₅₀ = 4.1 µM)~five times more potent than initial hit compound 2. Selectivity screening of low micromolar inhibitors 14 and 19 (comparing aryl- and alkyl-substituents) against 10 homologous kinases revealed unique selectivity profiles, with both compounds more potent against the GSK-3α isoform (IC₅₀s~2 µM) and, additionally, inhibitors of PKBβ (IC₅₀s < 25 µM). Compound 14 also inhibited ERK2 and 19, PKCγ, but generally good selectivity for GSK-3 isoforms over the other kinases was observed. The compounds had excellent predicted oral bioavailability and CNS-activity profiles, presenting promising candidates for future testing in cellular models of disease.

Modeling disruption of Apis mellifera (honey bee) odorant-binding protein function with high-affinity binders

Chemical toxins pose a great threat to honey bee health because they affect memory and cognition, diminish immunity, and increase susceptibility to infection, resulting in decreased colony performance, reproduction, and survival. Although the behavioral effects of sub-lethal chemical exposure on honey bees have been intensively studied, how xenobiotics affect olfaction, at the molecular level, still needs to be elucidated. In the present work, in silico tools, such as molecular docking, binding free energy calculations, and molecular dynamics simulations are used to predict if environmental chemicals have stronger binding affinities to honey bee antennal odorant-binding protein 14 (OBP14) than the representative floral odors citralva, eugenol, and the fluorescent probe 1-N-phenylnaphthylamine. Based on structural analysis, 21 chemicals from crop pesticides, household appliances, cosmetics, food, public health-related products, and other sources, many of which are pervasive in the hive environment, have higher binding affinities than the floral odors. These results suggest that chemical exposures are likely to interfere with the honey bee's sense of smell and this disruptive mechanism may be responsible for the lower associative learning and memory based on olfaction found in bees exposed to pesticides. Moreover, bees mainly rely on olfactory cues to perceive their environment and orient themselves as well as to discriminate and identify their food, predators, nestmates, and diseased individuals that need to be removed with hygienic behavior. In summary, sub-lethal exposure to environmental toxins can contribute to colony collapse in several ways from the disruption of proper olfaction functioning.

In silico trials to design potent inhibitors against matrilysin (MMP-7)

The study deals with structure-based rational drug design against the chief zinc-rely endopeptidase called matrilysin (MMP-7) that is involved in inflammatory and metastasis process of several carcinomas. Hyperactivated matrilysin of human was targeted, because of its hydrolytic actions on extracellular matrix (ECM) protein components constitutes fibrillar collagens, gelatins, fibronectins and it also activates zymogen forms of vital matrix metalloproteinases (gelatinase A-MMP-2 and B-MMP-9) responsible for ECM destruction in many cancers. In the present work, e-pharmacophores were generated for the respective five co-crystal structures of human matrilysin by mapping ligand's pharmacophoric features. During the lead-optimization campaign, the five e-pharmacophores-based shape screening against an in-house library of >21 million compounds created a dataset of 5000 structural analogs. The subsequent three different docking strategies, including rigid-receptor docking, quantum-polarized-ligand docking, induced-fit docking and free energy binding calculations resulted four leads as novel and potent MMP-7 binders. These four leads were observed with good pharmacological features and good receiver operating characteristics curve metrics (ROC: 0.93) in post-docking evaluations against five existing co-crystal inhibitors and 1000 decoy molecules with MMP-7. Moreover, stability and dynamics behavior of matrilysin-lead1 complex and matrilysin-cocrystal ligand (TQJ) complex were analyzed in natural physiological milieu of 1000 ns or 1 µs molecular dynamics simulations. Lead1-MMP-7 complex was found with an average Cα root-mean-square deviation (RMSD) of 2.35 Å, average ligand root-mean-square fluctuations (RMSF) of 0.66 Å and the strong metallic interactions with E220, a key residue for proteolytic action thereby hinders ECM proteolysis that in turn can halt metastatic cancerous condition.Communicated by Ramaswamy H. Sarma.

Hydrogen bonding penalty used for virtual screening to discover potent inhibitors for Papain-Like cysteine proteases of SARS-CoV-2

The Papain-Like proteases (PLpro) of SARS-CoV-2 play a crucial role in viral replication and the formation of nonstructural proteins. To find available inhibitors, the 3D structure of PLpro of SARS2 was obtained by homologous modelling, and we used this structure as a target to search for inhibitors through molecular docking and MM/GBSA binding free energy rescoring. A novel hydrogen bonding penalty was applied to the screening process, which meanwhile took desolvation into account. Finally, 61 compounds were acquired and 4 of them with IC50 at micromolar level tested in vitro enzyme activity assay, which includes clinical drugs tegaserod. Considering the importance of crystal water molecules, the 4 compounds were re-docked and considered bound waters in the active site as a part of PLpro. The binding modes of these 4 compounds were further explored with metadynamics simulations. The hits will provide a starting point for future key interactions identified and lead optimization targetting PLpro.

The interaction analysis between human serum albumin with chlorpyrifos and its derivatives through sub-atomic docking and molecular dynamics simulation techniques

Chlorpyrifos (CPF) is an extensively used organophosphate pesticide for crop protection. However, there are concerns about it contaminating the environment and human health, with estimated three lakh deaths annually. The molecular modeling protocol was assisted in redesigning thirteen well-known CPF linkers and inserting them at five selectable CPF (R1-R5) positions of CPF to get 258 CPF derivatives. CPF and its derivatives were optimized using LigPrep and docked to a grid centralized on Trp214 using extra precision glide docking. The Binding free energy of complexes was calculated using molecular mechanics/generalized born surface area (MM-GBSA). CPF and CPFD-225 have glide scores of - 3.08 and - 6.152 kcal/mol, respectively, with human serum albumin and ΔG bind for CPF (- 33.041817 kcal/mol) (- 52.825 kcal/mol) for CPF-D225. The top ten CPF derivatives showed at least ninefold better binding free energy than the CPF proposed for polyclonal antibody production. Subsequently, molecular docking studies revealed that CPF and its derivatives could bind to human serum albumin (HSA). Furthermore, using the Desmond package, a 100-ns molecular dynamics (MD) simulation was performed on the potential binding site. The final systems of CPF-HSA and CPF-222D complexes consist of 76,014 and 76,026 atoms, respectively. The physical stability of both the systems (CPF-HSA and CPF-222D) was analyzed by considering the overall potential energy, RMSF, RMSD, Hydrophobic interactions, and water-mediated patterns, which showed total energy of - 141,610 kcal/mol and - 140,150 kcal/mol, respectively.

Supplementary Information

The online version contains supplementary material available at 10.1007/s13205-022-03344-7.

Homology Modelling, Docking-based Virtual Screening, ADME Properties, and Molecular Dynamics Simulation for Identification of Probable Type II Inhibitors of AXL Kinase

Background:

AXL kinase is an important member of the TAM family for kinases which is involved in most cancers. Considering its role in different cancers due to its pro-tumorigenic effects and its involvement in the resistance, it has gained importance recently. Majority of research carried out is on Type I inhibitors and limited studies have been carried out for Type II inhibitors. Taking this into consideration, we have attempted to build Homology models to identify the Type II inhibitors for the AXL kinase.

Methods:

Homology Models for DFG-out C-helix-in/out state were developed using SWISS Model, PRIMO, and Prime. These models were validated by different methods and further evaluated for stability by molecular dynamics simulation using Desmond software. Selected models PED1-EB and PEDI1-EB were used for the docking-based virtual screening of four compound libraries using Glide software. The hits identified were subjected to interaction analysis and shortlisted compounds were subjected to Prime MM-GBSA studies for energy calculation. These compounds were also docked in the DFG-in state to check for binding and elimination of any compounds that may not be Type II inhibitors. The Prime energies were calculated for these complexes as well and some compounds were eliminated. ADMET studies were carried out using Qikprop. Some selected compounds were subjected to molecular dynamics simulation using Desmond for evaluating the stability of the complexes.

Results:

Out of 78 models inclusive of both DFG-out C-helix-in and DFG-out C-helix-out, 5 models were identified after different types of evaluation as well as validation studies. 1 model representing each type (PED1-EB and PEDI1-EB) was selected for the screening studies. The screening studies resulted in the identification of 29 compounds from the screen on PED1-EB and 10 compounds from the screen on PEDI1-EB. Hydrogen bonding interactions with Pro621, Met623, and Asp690 were observed for these compounds primarily. In some compounds, hydrogen bonding with Leu542, Glu544, Lys567, and Asn677 as well as pi-pi stacking interactions with either Phe622 or Phe691 were also seen. 4 compounds identified from PED1-EB screen were subjected to molecular dynamics simulation and their interactions were found to be consistent during the simulation. 2 compounds identified from PEDI1-EB screen were also subjected to the simulation studies, however, their interactions with Asp690 were not observed for a significant time and in both cases differed from the docked pose.

Conclusion:

Multiple models of DFG-out conformations of AXL kinase were built, validated and used for virtual screening. Different compounds were identified in the virtual screening, which may possibly act as Type II inhibitors for AXL kinase. Some more experimental studies can be done to validate these findings in future. This study will play a guiding role in the further development of the newer Type II inhibitors of the AXL kinase for the probable treatment of cancer.

A Multi-target Drug Designing for BTK, MMP9, Proteasome and TAK1 for the Clinical Treatment of Mantle Cell Lymphoma

BACKGROUND

Mantle cell lymphoma (MCL) is a type of non-Hodgkin lymphoma characterized by the mutation and overexpression of the cyclin D1 protein by the reciprocal chromosomal translocation t(11;14)(q13:q32).

AIM

The present study aims to identify potential inhibition of MMP9, Proteasome, BTK, and TAK1 and determine the most suitable and effective protein target for the MCL.

METHODOLOGY

Nine known inhibitors for MMP9, 24 for proteasome, 15 for BTK and 14 for TAK1 were screened. SB-3CT (PubChem ID: 9883002), oprozomib (PubChem ID: 25067547), zanubrutinib (PubChem ID: 135565884) and TAK1 inhibitor (PubChem ID: 66760355) were recognized as drugs with high binding capacity with their respective protein receptors. 41, 72, 102 and 3 virtual screened compounds were obtained after the similarity search with compound (PubChem ID:102173753), PubChem compound SCHEMBL15569297 (PubChem ID:72374403), PubChem compound SCHEMBL17075298 (PubChem ID:136970120) and compound CID: 71814473 with best virtual screened compounds.

RESULT

MMP9 inhibitors show commendable affinity and good interaction profile of compound holding PubChem ID:102173753 over the most effective established inhibitor SB-3CT. The pharmacophore study of the best virtual screened compound reveals its high efficacy based on various interactions. The virtual screened compound's better affinity with the target MMP9 protein was deduced using toxicity and integration profile studies.

CONCLUSION

Based on the ADMET profile, the compound (PubChem ID: 102173753) could be a potent drug for MCL treatment. Similar to the established SB-3CT, the compound was non-toxic with LD50 values for both the compounds lying in the same range.

Oximes: Novel Therapeutics with Anticancer and Anti-Inflammatory Potential

Oximes have been studied for decades because of their significant roles as acetylcholinesterase reactivators. Over the last twenty years, a large number of oximes have been reported with useful pharmaceutical properties, including compounds with antibacterial, anticancer, anti-arthritis, and anti-stroke activities. Many oximes are kinase inhibitors and have been shown to inhibit over 40 different kinases, including AMP-activated protein kinase (AMPK), phosphatidylinositol 3-kinase (PI3K), cyclin-dependent kinase (CDK), serine/threonine kinases glycogen synthase kinase 3 α/β (GSK-3α/β), Aurora A, B-Raf, Chk1, death-associated protein-kinase-related 2 (DRAK2), phosphorylase kinase (PhK), serum and glucocorticoid-regulated kinase (SGK), Janus tyrosine kinase (JAK), and multiple receptor and non-receptor tyrosine kinases. Some oximes are inhibitors of lipoxygenase 5, human neutrophil elastase, and proteinase 3. The oxime group contains two H-bond acceptors (nitrogen and oxygen atoms) and one H-bond donor (OH group), versus only one H-bond acceptor present in carbonyl groups. This feature, together with the high polarity of oxime groups, may lead to a significantly different mode of interaction with receptor binding sites compared to corresponding carbonyl compounds, despite small changes in the total size and shape of the compound. In addition, oximes can generate nitric oxide. This review is focused on oximes as kinase inhibitors with anticancer and anti-inflammatory activities. Oximes with non-kinase targets or mechanisms of anti-inflammatory activity are also discussed.

Matairesinol, an active constituent of HC9 polyherbal formulation, exhibits HDAC8 inhibitory and anticancer activity

Histone deacetylase 8 (HDAC8) has emerged as a promising drug target for cancer therapeutics development. HDAC8 has been reported to regulate cancer cell proliferation, invasion and promote metastasis through modulation of cell cycle associated proteins. Of late, phytocompounds have been demonstrated to exhibit anticancer and anti-HDAC8 activity. Here, we have shown the HDAC8 inhibitory potential of an active phytocompound from HC9 (herbal composition-9), a polyherbal anticancer formulation based on the traditional Ayurvedic drug, Stanya Shodhan Kashaya. HC9 was recently reported to exhibit anticancer activity against breast cancer cells through induction of cell cycle arrest, decrease in migration and invasion as well as regulation of inflammation and chromatin modulators. In silico studies such as molecular docking, molecular dynamics (MD) simulation and binding free energy analyses showed greater binding energy values and interaction stability of MA with HDAC8 compared to other phytocompounds of HC9. Interestingly, in vitro validation confirmed the anti-HDAC8 activity of MA. Further, in vitro studies showed that MA significantly decreased the viability of breast and prostate cancer cell lines, thereby confirming its anticancer potential.

Structure-Based Design of Potent Selective Nanomolar Type-II Inhibitors of Glycogen Synthase Kinase-3β

For the first time, the in silico design, screening, and in vitro validation of potent GSK-3β type-II inhibitors are presented. In the absence of crystallographic evidence for a DFG-out GSK-3β activation loop conformation, computational models were designed using an adapted DOLPHIN approach and a method consisting of Prime loop refinement, induced-fit docking, and molecular dynamics. Virtual screening of the Biogenics subset from the ZINC database led to an initial selection of 20 Phase I compounds revealing two low micromolar inhibitors in an isolated enzyme assay. Twenty more analogues (Phase II compounds) related to the hit [pyrimidin-2-yl]amino-furo[3,2-b]furyl-urea scaffold were selected for structure-activity relationship analysis. The Phase II studies led to five highly potent nanomolar inhibitors, with compound 23 (IC₅₀ =0.087 μM) > 100 times more potent than the best Phase I inhibitor, and selectivity for GSK-3β inhibition compared to homologous kinases was observed. Ex vivo experiments (SH-SY5Y cell lines) for tau hyperphosphorylation revealed promising neuroprotective effects at low micromolar concentrations. The type-II inhibitor design has been unraveled as a potential route toward more clinically effective GSK-3β inhibitors.

Affinity Crystallography Reveals Binding of Pomegranate Juice Anthocyanins at the Inhibitor Site of Glycogen Phosphorylase: The Contribution of a Sugar Moiety to Potency and Its Implications to the Binding Mode

Anthocyanins (ACNs) are dietary phytochemicals with an acknowledged therapeutic significance. Pomegranate juice (PJ) is a rich source of ACNs with potential applications in nutraceutical development. Glycogen phosphorylase (GP) catalyzes the first step of glycogenolysis and is a molecular target for the development of antihyperglycemics. The inhibitory potential of the ACN fraction of PJ is assessed through a combination of in vitro assays, ex vivo investigation in hepatic cells, and X-ray crystallography studies. The ACN extract potently inhibits muscle and liver isoforms of GP. Affinity crystallography reveals the structural basis of inhibition through the binding of pelargonidin-3-O-glucoside at the GP inhibitor site. The glucopyranose moiety is revealed as a major determinant of potency as it promotes a structural binding mode different from that observed for other flavonoids. This inhibitory effect of the ACN scaffold and its binding mode at the GP inhibitor binding site may have significant implications for future structure-based drug design endeavors.

The crystal structure of Arabidopsis BON1 provides insights into the copine protein family

The Arabidopsis thaliana BON1 gene product is a member of the evolutionary conserved eukaryotic calcium-dependent membrane-binding protein family. The copine protein is composed of two C2 domains (C2A and C2B) followed by a vWA domain. The BON1 protein is localized on the plasma membrane, and is known to suppress the expression of immune receptor genes and to positively regulate stomatal closure. The first structure of this protein family has been determined to 2.5-Å resolution and shows the structural features of the three conserved domains C2A, C2B and vWA. The structure reveals the third Ca²⁺ -binding region in C2A domain is longer than classical C2 domains and a novel Ca²⁺ binding site in the vWA domain. The structure of BON1 bound to Mn²⁺ is also presented. The binding of the C2 domains to phospholipid (PSF) has been modeled and provides an insight into the lipid-binding mechanism of the copine proteins. Furthermore, the selectivity of the separate C2A and C2B domains and intact BON1 to bind to different phospholipids has been investigated, and we demonstrated that BON1 could mediate aggregation of liposomes in response to Ca²⁺ . These studies have formed the basis of further investigations into the important role that the copine proteins play in vivo.

Estrogen receptor modulation of some polyphenols extracted from Daucus carota as a probable mechanism for antifertility effect: An in silico study

The seeds of Daucus carota, traditionally used by women in many countries to prevent conception, were proved to have negative impact on reproductive hormone levels as well as on the estrous cycles in albino mice affecting the fertility status. This study is an attempt to investigate the possible role of polyphenols present in the seeds in hampering the reproductive processes. The Molecular Docking, Molecular Dynamics (MD) simulation and binding free energy calculation studies reveal that six polyphenols present in the seeds can bind with the active sites of human Estrogen Receptor (ER) and may interfere in the estrogen signaling in human. These polyphenols were found to bind to a conservative pocket of ER[Formula: see text], which is comprised of residues 343–388, 421–428 and 525–540. Docking studies indicated the presence of strong hydrogen bonding, pi–pi interactions and numerous hydrophobic interactions that stabilize the ER[Formula: see text]-polyphenol complexes.

The docked complexes were further subjected to MM/GBSA analysis to calculate binding free energies. Molecular dynamic simulation studies carried out for a period of 20[Formula: see text]ns revealed low RMS deviation values suggesting high accuracy of the docking poses and stability of the complexes. Out of the six polyphenols, catechin and epicatechin have shown highest binding affinity towards the ER[Formula: see text] receptor. These findings will help in identifying ER modulators of plant origin targeting ER alpha and predicting their effects on the reproductive hormone homeostasis. Moreover, this study may form preliminary basis for further identification of potential herbal antifertility agents.

Identification of Phosphoinositide-3 Kinases Delta and Gamma Dual Inhibitors Based on the p110δ/γ Crystal Structure

Background:

Phosphoinositide-3 kinases (PI3Ks) are key signaling molecules that affect a diverse array of biological processes in cells, including proliferation, differentiation, survival, and metabolism. The abnormal activity of PI3K signals is closely related to the occurrence of many diseases, which has become a very promising drug target, especially for the treatment of cancer. PI3Kδ/γ inhibitors can reduce toxicity concerns for chronic indications such as asthma and rheumatoid arthritis compared with pan PI3Ks inhibitors.

Methods:

With the aim of finding more effective PI3Kδ/γ dual inhibitors, virtual screening, ADMET prediction Molecular Dynamics (MD) simulations and MM-GBSA were executed based on the known p110δ/γ crystal structure. Compound ZINC28564067 with high docking score and low toxicity was obtained.

Results:

By MD simulations and MM-GBSA, we could observe that ZINC28564067 had more favorable conformation binding to the PI3Kδ/γ than the original ligands.

Conclusion:

The results provided a rapid approach for the discovery of novel PI3Kδ/γ dual inhibitors which might be a potential anti-tumor lead compound.

Understanding the binding affinity of noscapines with protease of SARS-CoV-2 for COVID-19 using MD simulations at different temperatures

The current outbreak of a novel coronavirus, named as SARS-CoV-2 causing COVID-19 occurred in 2019, is in dire need of finding potential therapeutic agents. Recently, ongoing viral epidemic due to coronavirus (SARS-CoV-2) primarily affected mainland China that now threatened to spread to populations in most countries of the world. In spite of this, there is currently no antiviral drug/ vaccine available against coronavirus infection, COVID-19. In the present study, computer-aided drug design-based screening to find out promising inhibitors against the coronavirus (SARS-CoV-2) leads to infection, COVID-19. The lead therapeutic molecule was investigated through docking and molecular dynamics simulations. In this, binding affinity of noscapines(23B)-protease of SARS-CoV-2 complex was evaluated through MD simulations at different temperatures. Our research group has established that noscapine is a chemotherapeutic agent for the treatment of drug resistant cancers; however, noscapine was also being used as anti-malarial, anti-stroke and cough-suppressant. This study suggests for the first time that noscapine exerts its antiviral effects by inhibiting viral protein synthesis.

Identification of C-β-d-Glucopyranosyl Azole-Type Inhibitors of Glycogen Phosphorylase That Reduce Glycogenolysis in Hepatocytes: In Silico Design, Synthesis, in Vitro Kinetics, and ex Vivo Studies

Several C-β-d-glucopyranosyl azoles have recently been uncovered as among the most potent glycogen phosphorylase (GP) catalytic site inhibitors discovered to date. Toward further exploring their translational potential, ex vivo experiments have been performed for their effectiveness in reduction of glycogenolysis in hepatocytes. New compounds for these experiments were predicted in silico where, for the first time, effective ranking of GP catalytic site inhibitor potencies using the molecular mechanics-generalized Born surface area (MM-GBSA) method has been demonstrated. For a congeneric training set of 27 ligands, excellent statistics in terms of Pearson (R_P) and Spearman (R_S) correlations (both 0.98), predictive index (PI = 0.99), and area under the receiver operating characteristic curve (AU-ROC = 0.99) for predicted versus experimental binding affinities were obtained, with ligand tautomeric/ionization states additionally considered using density functional theory (DFT). Seven 2-aryl-4(5)-(β-d-glucopyranosyl)-imidazoles and 2-aryl-4-(β-d-glucopyranosyl)-thiazoles were subsequently synthesized, and kinetics experiments against rabbit muscle GPb revealed new potent inhibitors with best K_i values in the low micromolar range (5c = 1.97 μM; 13b = 4.58 μM). Ten C-β-d-glucopyranosyl azoles were then tested ex vivo in mouse primary hepatocytes. Four of these (5a-c and 9d) demonstrated significant reduction of glucagon stimulated glycogenolysis (IC₅₀ = 30-60 μM). Structural and predicted physicochemical properties associated with their effectiveness were analyzed with permeability related parameters identified as crucial factors. The most effective ligand series 5 contained an imidazole ring, and the calculated pK_a (Epik: 6.2; Jaguar 5.5) for protonated imidazole suggests that cellular permeation through the neutral state is favored, while within the cell, there is predicted more favorable binding to GP in the protonated form.

High Consistency of Structure-Based Design and X-Ray Crystallography: Design, Synthesis, Kinetic Evaluation and Crystallographic Binding Mode Determination of Biphenyl-N-acyl-β-d-Glucopyranosylamines as Glycogen Phosphorylase Inhibitors

Structure-based design and synthesis of two biphenyl-N-acyl-β-d-glucopyranosylamine derivatives as well as their assessment as inhibitors of human liver glycogen phosphorylase (hlGPa, a pharmaceutical target for type 2 diabetes) is presented. X-ray crystallography revealed the importance of structural water molecules and that the inhibitory efficacy correlates with the degree of disturbance caused by the inhibitor binding to a loop crucial for the catalytic mechanism. The in silico-derived models of the binding mode generated during the design process corresponded very well with the crystallographic data.

Preliminary pharmacokinetic study of the anticancer 6BIO in mice using an UHPLC-MS/MS approach

Indirubins represent a group of natural and synthetic products with bio-activities against numerous human cancer cell lines acting by inhibiting protein kinases. The natural sources of indirubins are plants of Isatis sp., Indigofera sp., and Polygonum sp., recombinant bacteria, mammalian urine and some marine mollusks. Specifically, the halogenated derivative 6-bromo indirubin-3'-oxime (6BIO) possesses increased selectivity against GSK-3. However, to our knowledge, no analytical method to determine 6BIO in biological fluids has been developed till now. Therefore, a rapid, sensitive and high throughput UHPLC-MS/MS methods were developed and validated to evaluate the concentrations of 6BIO in mice plasma. Plasma samples were pre-treated by protein precipation using cold mixture of methanol: acetonitrile (9:1, v/v) and separations were carried out on a Hypersil Gold C18 column (50 × 2.1 mm i.d.; 1.9 μm p.s.) using 0.1% acetic acid and methanol as mobile phase at a flow rate of 500 mL/min in a gradient mode. For quantitation, a hybrid LTQ-Orbitrap MS equipped with an electro-spray ionization source was used applying a selected reaction monitoring (SRM) option. The monitored transitions were m/z 354.0 → 324.0 for 6BIO and 297.1 → 282.1 for afromorsin (used as the internal standard) in the negative mode. Following the EMA, ICH and FDA guidelines for validation of analytical procedures, the assay method was fully validated in terms of selectivity, linearity, recovery, matrix effect, accuracy, precision, stability, and robustness. The validated methods were successfully applied to the pharmacokinetic studies of 6BIO following an oral administration to mice at the dose of 50 mg/kg. The results indicated that 6BIO possesses a T_max of 30 min, a half-life of 1 h, and low plasma bioavailability.

Targeting Nucleotide Binding Domain of Multidrug Resistance-associated Protein-1 (MRP1) for the Reversal of Multi Drug Resistance in Cancer

Multidrug resistance (MDR) is the major cause, by which cancer cells expel the drugs out, developing a challenge against the current chemotherapeutic drugs regime. This mechanism is attributed to the over expression of ABC transporters like MRP1 on the surface of cells. Since nucleotide binding domains (NBD) of ABC transporters are the site of ATP binding and hydrolysis, thereby in this study we have targeted NBD1 of MRP1using molecular docking and molecular dynamic simulations (MDS). The compounds present in the FDA approved library were docked against NBD1 of the human multidrug resistance associated protein 1 (PDB ID: 2CBZ). For the docking studies, Standard Precision and Extra Precision methods were employed. After the EP docking studies, ligands showed an extremely low docking score that was indicative of very high binding affinity of the ligands to the NBD. Apart from the low docking score, another short listing criterion in simulation studies was the interaction of incoming ligand with the desired conserved residues of NDB involved in ATP binding and hydrolysis. Based on these measures, potassium citrate (DB09125) and technetium Tc-99m medronate (DB09138) were chosen and subjected to 100 ns simulation studies. From the MDS study we concluded that between these two compounds, potassium citrate is a better candidate for targeting MRP1.

Evidence for Novel Action at the Cell-Binding Site of Human Angiogenin Revealed by Heteronuclear NMR Spectroscopy, in silico and in vivo Studies

A member of the ribonuclease A superfamily, human angiogenin (hAng) is a potent angiogenic factor. Heteronuclear NMR spectroscopy combined with induced-fit docking revealed a dual binding mode for the most antiangiogenic compound of a series of ribofuranosyl pyrimidine nucleosides that strongly inhibit hAng's angiogenic activity in vivo. While modeling suggests the potential for simultaneous binding of the inhibitors at the active and cell-binding sites, NMR studies indicate greater affinity for the cell-binding site than for the active site. Additionally, molecular dynamics simulations at 100 ns confirmed the stability of binding at the cell-binding site with the predicted protein-ligand interactions, in excellent agreement with the NMR data. This is the first time that a nucleoside inhibitor is reported to completely inhibit the angiogenic activity of hAng in vivo by exerting dual inhibitory activity on hAng, blocking both the entrance of hAng into the cell and its ribonucleolytic activity.

Design, sythesis and evaluation of a series of 3- or 4-alkoxy substituted phenoxy derivatives as PPARs agonists

Peroxisome proliferators-activated receptors (PPARα, γ and δ) are potentially effective targets for Type 2 diabetes mellitus therapy. The severe effects of known glitazones and the successfully approved agents (saroglitazar and lobeglitazone) motivated us to study novelly potent PPARs drugs with improved safety profile. In this work, we received 15 carboxylic acids based on the combination principle to integrate the polar head of bezafibrate with the hydrophobic tail of pioglitazone. Another 12 tetrazoles based on the bioisosterism principle were obtained accordingly. Furthermore, in vitro PPARs transactivation assays on these 3- or 4-alkoxy substituted phenoxy derivatives afforded six compounds. Interactions and binding stability from the docking analysis and 20 ns molecular dynamic simulations confirmed the representative compounds to be suitable and plausible for PPARs pockets. The above-mentioned results demonstrated that the compounds may be used as reference for further optimization for enhanced PPARs activities and wide safety range.

Interactions of omeprazole-based analogues with cytochrome P450 2C19: a computational study

Cytochrome P450 2C19 (CYP2C19) is one of 57 drug metabolizing enzymes in humans and is responsible for the metabolism of ∼7-10% of drugs in clinical use. Recently omeprazole-based analogues were reported to be the potent inhibitors of CYP2C19 and have the potential to be used as the tool compounds for studying the substrate selectivity of CYP2C19. However, the binding modes of these compounds with CYP2C19 remain to be elucidated. In this study, a combination of molecular docking, molecular dynamics (MD), and MM/GBSA calculations was employed to systematically investigate the interactions between these compounds and CYP2C19. The binding modes of these analogues were analyzed in detail. The results indicated that the inclusion of explicit active site water molecules could improve binding energy prediction when the water molecules formed a hydrogen bonding network between the ligand and protein. We also found that the effect of active site water molecules on binding free energy prediction was dependent on the ligand binding modes. Our results unravel the interactions of these omeprazole-based analogues with CYP2C19 and might be helpful for the future design of potent CYP2C19 inhibitors with improved metabolic properties.

Molecular Dynamics Simulations for Biological Systems

Molecular dynamics simulation is an important tool to capture the dynamicity of biological molecule and the atomistic insights. These insights are helpful to explore biological functions. Molecular dynamics simulation from femto seconds to milli seconds scale give a large ensemble of conformations that can reveal many biological mysteries. The main focus of the chapter is to throw light on theories, requirement of molecular dynamics for biological studies and application of molecular dynamics simulations. Molecular dynamics simulations are widely used to study protein-protein interaction, protein-ligand docking, effects of mutation on interactions, protein folding and flexibility of the biological molecules. This chapter also deals with various methods/algorithms of protein tertiary structure prediction, their strengths and weaknesses.

Identification of dual ligands targeting angiotensin II type 1 receptor and peroxisome proliferator-activated receptor-γ by core hopping of telmisartan

It has been reported previously that some angiotensin II receptor blockers not only antagonize angiotensin II type 1 receptor (AT₁R), but also exert stimulation in peroxisome proliferator-activated receptor γ (PPARγ) partial activation, among which telmisartan displays the best. Telmisartan has been tested as a bifunctional ligand with antihypertensive and hypoglycemic activity. Aiming at more potent leads with selective AT₁R antagonism and PPARγ partial agonism, the three parts of telmisartan including the distal benzimidazole ring, the biphenyl moiety, and the carboxylic acid group experienced modification by core hopping method in our study. The central benzimidazole ring, however, remained intact considering its great affinity toward AT₁R and PPARγ. We utilized computational techniques for the sake of details on the binding interactions and conformational stability. Standard precision docking analysis and absorption, distribution, metabolism, excretion, and toxicity prediction received 10 molecules with higher Glide scores, similar interactions, and improved pharmacokinetic profiles compared to telmisartan. Comp#91 with highest scores for AT₁R (-11.92 kcal/mol) and PPARγ (-13.88 kcal/mol) exhibited excellent binding modes and pharmacokinetic parameters. Molecular dynamics trajectories on best docking pose of comp#91 confirmed the docking results and verified the conformational stability with both receptors throughout the course of 20-ns simulations. Thus, comp#91 could be identified as a promising lead in the development of dual AT₁R antagonist and PPARγ partial agonist against hypertension and type 2 diabetes.

E-pharmacophore-based virtual screening to identify GSK-3β inhibitors

Glycogen synthase kinase-3β (GSK-3β) is a serine/threonine kinase which has attracted significant attention during recent years in drug design studies. The deregulation of GSK-3β increased the loss of hippocampal neurons by triggering apoptosis-mediating production of neurofibrillary tangles and alleviates memory deficits in Alzheimer's disease (AD). Given its role in the formation of neurofibrillary tangles leading to AD, it has been a major therapeutic target for intervention in AD, hence was targeted in the present study. Twenty crystal structures were refined to generate pharmacophore models based on energy involvement in binding co-crystal ligands. Four common e-pharmacophore models were optimized from the 20 pharmacophore models. Shape-based screening of four e-pharmacophore models against nine established small molecule databases using Phase v3.9 had resulted in 1800 compounds having similar pharmacophore features. Rigid receptor docking (RRD) was performed for 1800 compounds and 20 co-crystal ligands with GSK-3β to generate dock complexes. Interactions of the best scoring lead obtained through RRD were further studied with quantum polarized ligand docking (QPLD), induced fit docking (IFD) and molecular mechanics/generalized Born surface area. Comparing the obtained leads to 20 co-crystal ligands resulted in 18 leads among them, lead1 had the lowest docking score, lower binding free energy and better binding orientation toward GSK-3β. The 50 ns MD simulations run confirmed the stable nature of GSK-3β-lead1 docking complex. The results from RRD, QPLD, IFD and MD simulations confirmed that lead1 might be used as a potent antagonist for GSK-3β.

Structural signatures of DRD4 mutants revealed using molecular dynamics simulations: Implications for drug targeting

Human Dopamine Receptor D4 (DRD4) orchestrates several neurological functions and represents a target for many psychological disorders. Here, we examined two rare variants in DRD4; V194G and R237L, which elicit functional alterations leading to disruption of ligand binding and G protein coupling, respectively. Using atomistic molecular dynamics (MD) simulations, we provide in-depth analysis to reveal structural signatures of wild and mutant complexes with their bound agonist and antagonist ligands. We constructed intra-protein network graphs to discriminate the global conformational changes induced by mutations. The simulations also allowed us to elucidate the local side-chain dynamical variations in ligand-bound mutant receptors. The data suggest that the mutation in transmembrane V (V194G) drastically disrupts the organization of ligand binding site and causes disorder in the native helical arrangement. Interestingly, the R237L mutation leads to significant rewiring of side-chain contacts in the intracellular loop 3 (site of mutation) and also affects the distant transmembrane topology. Additionally, these mutations lead to compact ICL3 region compared to the wild type, indicating that the receptor would be inaccessible for G protein coupling. Our findings thus reveal unreported structural determinants of the mutated DRD4 receptor and provide a robust framework for design of effective novel drugs.

Molecular dynamics simulation studies of GSK-3β ATP competitive inhibitors: understanding the factors contributing to selectivity

Glycogen synthase kinase-3 is a constitutively acting, multifunctional serine threonine kinase, the role of which has been implicated in several physiological pathways and has emerged as a promising target for the treatment of type-II diabetes and Alzheimer's disease. In order to provide a detailed understanding of the origin of selectivity determinants of ATP competitive inhibitors, molecular dynamics simulations in combination with MM-PBSA binding energy calculations were performed using crystal structures of GSK-3β and CDK-2 in complex with 12 ATP competitive inhibitors. Analysis of energy contributions indicate that electrostatic interaction energy dictates the selectivity of ATP competitive inhibitors against CDK-2. Key interactions as well as residues that potentially make a major contribution to the binding free energy were identified at the ATP binding site. This analysis stresses the need for the inhibitors to interact with Lys85, Thr138, and Arg141 in the binding site of GSK-3β to show selectivity. The residue-wise energy decomposition analysis further suggested the additional role of Gln185 in determining the selectivity of maleimides. The results obtained in this study can be utilized to design new selective GSK-3 ATP competitive inhibitors.

A systematic methodology for large scale compound screening: A case study on the discovery of novel S1PL inhibitors

Decrease in sphingosine 1-phosphate (S1P) concentration induces migration of pathogenic T cells to the blood stream, disrupts the CNS and it is implicated in multiple sclerosis (MS), a progressive inflammatory disorder of the central nervous system (CNS), and Alzheimer's disease (AD). A promising treatment alternative for MS and AD is inhibition of the activity of the microsomal enzyme sphingosine 1-phosphate lyase (S1PL), which degrades intracellular S1P. This report describes an integrated systematic approach comprising virtual screening, molecular docking, substructure search and molecular dynamics simulation to discover novel S1PL inhibitors. Virtual screening of the ZINC database via ligand-based and structure-based pharmacophore models yielded 10000 hits. After molecular docking, common substructures of the top ranking hits were identified. The ligand binding poses were optimized by induced fit docking. MD simulations were performed on the complex structures to determine the stability of the S1PL-ligand complex and to calculate the binding free energy. Selectivity of the selected molecules was examined by docking them to hERG and cytochrome P450 receptors. As a final outcome, 15 compounds from different chemotypes were proposed as potential S1PL inhibitors. These molecules may guide future medicinal chemistry efforts in the discovery of new compounds against the destructive action of pathogenic T cells.

Structure guided design of potent and selective ponatinib-based hybrid inhibitors for RIPK1

RIPK1 and RIPK3, two closely related RIPK family members, have emerged as important regulators of pathologic cell death and inflammation. In the current work, we report that the Bcr-Abl inhibitor and anti-leukemia agent ponatinib is also a first-in-class dual inhibitor of RIPK1 and RIPK3. Ponatinib potently inhibited multiple paradigms of RIPK1- and RIPK3-dependent cell death and inflammatory tumor necrosis factor alpha (TNF-α) gene transcription. We further describe design strategies that utilize the ponatinib scaffold to develop two classes of inhibitors (CS and PN series), each with greatly improved selectivity for RIPK1. In particular, we detail the development of PN10, a highly potent and selective "hybrid" RIPK1 inhibitor, capturing the best properties of two different allosteric RIPK1 inhibitors, ponatinib and necrostatin-1. Finally, we show that RIPK1 inhibitors from both classes are powerful blockers of TNF-induced injury in vivo. Altogether, these findings outline promising candidate molecules and design approaches for targeting RIPK1- and RIPK3-driven inflammatory pathologies.

Newly discovered chroman-2,4-diones neutralize the in vivo DNA damage induced by alkylation through the inhibition of Topoisomerase IIα: A story behind the molecular modeling approach

Eight chroman-2,4-diones, namely 2a-h, previously investigated as anticoagulants, of which 2a and 2f as the most active, were evaluated as in vivo genotoxic agents in Wistar rat livers and kidneys using the comet assay. Compounds 2a, 2b, and 2f without genotoxic activity were applied prior to ethyl methanesulfonate (EMS) and diminished EMS-induced DNA damage according to the total score and percentage of reduction. EMS produce harmful O(6)-ethylguanine lesion which is incorporated in aberrant genotoxic GT and TG pairing after ATP-dependent DNA strand breaks have been catalyzed by rat Topoisomerase IIα (rTopIIα, EC 5.99.1.3). Therefore, the mechanism of 2a, 2b, and 2f antigenotoxic activity was investigated on the enzyme level using molecular docking and molecular dynamics simulations insamuch as it had been determined that compounds do not intercalate DNA but instead inhibit the ATPase activity. Calculations predicted that compounds inhibit ATP hydrolysis before the DNA-EMS cleavage is being catalyzed by rTopIIα, prevent EMS mutagenic and carcinogenic effects, and beside anticoagulant activity can even be applied in the cancer treatment to control the rate of anticancer alkylation drugs.