Determination of Ligand Binding Modes in Hydrated Viral Ion Channels to Foster Drug Design and Repositioning

Inhibition of xanthine oxidase-catalyzed xanthine and 6-mercaptopurine oxidation by luteolin, naringenin, myricetin, ampelopsin and their conjugated metabolites

Luteolin, naringenin, myricetin, and ampelopsin are abundant flavonoids in nature, and several dietary supplements also contain them at very high doses. After the peroral intake, flavonoids go through extensive presystemic biotransformation; therefore, typically their sulfate/glucuronic acid conjugates reach high concentrations in the circulation. Xanthine oxidase (XO) enzyme is involved in uric acid production, and it also takes part in the elimination of certain drugs (e.g., 6-mercaptopurine). The inhibitory effects of flavonoid aglycones on XO have been widely studied; however, only limited data are available regarding their sulfate and glucuronic acid conjugates. In this study, we examined the impacts of luteolin, naringenin, myricetin, ampelopsin, and their sulfate/glucuronide derivatives on XO-catalyzed xanthine and 6-mercaptopurine oxidations employing in vitro enzyme incubation assays and molecular modeling studies. Our major results/conclusions are the following: (1) Sulfate metabolites were stronger while glucuronic acid derivatives were weaker inhibitors of XO compared to the parent flavonoids. (2) Naringenin, ampelopsin, and their metabolites were weak inhibitors of the enzyme. (3) Luteolin, myricetin, and their sulfates were highly potent inhibitors of XO, and the glucuronides of luteolin showed moderate inhibitory impacts. (4) Conjugated metabolites of luteolin and myricetin can be involved in the inhibitory effects of these flavonoids on XO enzyme.

TWN-FS method: A novel fragment screening method for drug discovery

Fragment-based drug discovery (FBDD) is a well-established and effective method for generating diverse and novel hits in drug design. Kinases are suitable targets for FBDD due to their well-defined structure. Water molecules contribute to structure and function of proteins and also influence the environment within the binding pocket. Water molecules form a variety of hydrogen-bonded cyclic water-ring networks, collectively known as topological water networks (TWNs). Analyzing the TWNs in protein binding sites can provide valuable insights into potential locations and shapes for fragments within the binding site. Here, we introduce TWN-based fragment screening (TWN-FS) method, a novel screening method that suggests fragments through grouped TWN analysis within the protein binding site. We used this method to screen known CDK2, CHK1, IGF1R and ERBB4 inhibitors. Our findings suggest that TWN-FS method has the potential to effectively screen fragments. The TWN-FS method package is available on GitHub at https://github.com/pkj0421/TWN-FS.

The Advances and Limitations of the Determination and Applications of Water Structure in Molecular Engineering

Water is a key actor of various processes of nature and, therefore, molecular engineering has to take the structural and energetic consequences of hydration into account. While the present review focuses on the target-ligand interactions in drug design, with a focus on biomolecules, these methods and applications can be easily adapted to other fields of the molecular engineering of molecular complexes, including solid hydrates. The review starts with the problems and solutions of the determination of water structures. The experimental approaches and theoretical calculations are summarized, including conceptual classifications. The implementations and applications of water models are featured for the calculation of the binding thermodynamics and computational ligand docking. It is concluded that theoretical approaches not only reproduce or complete experimental water structures, but also provide key information on the contribution of individual water molecules and are indispensable tools in molecular engineering.

Interaction of Fumonisin B1, N-Palmitoyl-Fumonisin B1, 5-O-Palmitoyl-Fumonisin B1, and Fumonisin B4 Mycotoxins with Human Serum Albumin and Their Toxic Impacts on Zebrafish Embryos

Fumonisins are frequent food contaminants. The high exposure to fumonisins can cause harmful effects in humans and animals. Fumonisin B1 (FB1) is the most typical member of this group; however, the occurrence of several other derivatives has been reported. Acylated metabolites of FB1 have also been described as possible food contaminants, and the very limited data available suggest their significantly higher toxicity compared to FB1. Furthermore, the physicochemical and toxicokinetic properties (e.g., albumin binding) of acyl-FB1 derivatives may show large differences compared to the parent mycotoxin. Therefore, we tested the interactions of FB1, N-palmitoyl-FB1 (N-pal-FB1), 5-O-palmitoyl-FB1 (5-O-pal-FB1), and fumonisin B4 (FB4) with human serum albumin as well as the toxic effects of these mycotoxins on zebrafish embryos were examined. Based on our results, the most important observations and conclusions are the following: (1) FB1 and FB4 bind to albumin with low affinity, while palmitoyl-FB1 derivatives form highly stable complexes with the protein. (2) N-pal-FB1 and 5-O-pal-FB1 likely occupy more high-affinity binding sites on albumin. (3) Among the mycotoxins tested, N-pal-FB1 showed the most toxic effects on zebrafish, followed by 5-O-pal-FB1, FB4, and FB1. (4) Our study provides the first in vivo toxicity data regarding N-pal-FB1, 5-O-pal-FB1, and FB4.

Probing the Interactions of 31 Mycotoxins with Xanthine Oxidase: Alternariol, Alternariol-3-Sulfate, and α-Zearalenol Are Allosteric Inhibitors of the Enzyme

Mycotoxins are frequent toxic contaminants in foods and beverages, causing a significant health threat. Interactions of mycotoxins with biotransformation enzymes (e.g., cytochrome P450 enzymes, sulfotransferases, and uridine 5'-diphospho-glucuronosyltransferases) may be important due to their possible detoxification or toxic activation during enzymatic processes. Furthermore, mycotoxin-induced enzyme inhibition may affect the biotransformation of other molecules. A recent study described the strong inhibitory effects of alternariol and alternariol-9-methylether on the xanthine oxidase (XO) enzyme. Therefore, we aimed to test the impacts of 31 mycotoxins (including the masked/modified derivatives of alternariol and alternariol-9-methylether) on XO-catalyzed uric acid formation. Besides the in vitro enzyme incubation assays, mycotoxin depletion experiments and modeling studies were performed. Among the mycotoxins tested, alternariol, alternariol-3-sulfate, and α-zearalenol showed moderate inhibitory actions on the enzyme, representing more than tenfold weaker impacts compared with the positive control inhibitor allopurinol. In mycotoxin depletion assays, XO did not affect the concentrations of alternariol, alternariol-3-sulfate, and α-zearalenol in the incubates; thus, these compounds are inhibitors but not substrates of the enzyme. Experimental data and modeling studies suggest the reversible, allosteric inhibition of XO by these three mycotoxins. Our results help the better understanding of the toxicokinetic interactions of mycotoxins.

Testing Serum Albumins and Cyclodextrins as Potential Binders of the Mycotoxin Metabolites Alternariol-3-Sulfate, Alternariol-9-Monomethylether and Alternariol-9-Monomethylether-3-Sulfate

Alternaria mycotoxins, including alternariol (AOH), alternariol-9-monomethylether (AME), and their masked/modified derivatives (e.g., sulfates or glycosides), are common food contaminants. Their acute toxicity is relatively low, while chronic exposure can lead to the development of adverse health effects. Masked/modified metabolites can probably release the more toxic parent mycotoxin due to their enzymatic hydrolysis in the intestines. Previously, we demonstrated the complex formation of AOH with serum albumins and cyclodextrins; these interactions were successfully applied for the extraction of AOH from aqueous matrices (including beverages). Therefore, in this study, the interactions of AME, alternariol-3-sulfate (AS), and alternariol-9-monomethylether-3-sulfate (AMS) were investigated with albumins (human, bovine, porcine, and rat) and with cyclodextrins (sulfobutylether-β-cyclodextrin, sugammadex, and cyclodextrin bead polymers). Our major results/conclusions are the following: (1) The stability of mycotoxin-albumin complexes showed only minor species dependent variations. (2) AS and AMS formed highly stable complexes with albumins in a wide pH range, while AME-albumin interactions preferred alkaline conditions. (3) AME formed more stable complexes with the cyclodextrins examined than AS and AMS. (4) Beta-cyclodextrin bead polymer proved to be highly suitable for the extraction of AME, AS, and AMS from aqueous solution. (5) Albumins and cyclodextrins are promising binders of the mycotoxins tested.

Binding Networks Identify Targetable Protein Pockets for Mechanism-Based Drug Design

The human genome codes only a few thousand druggable proteins, mainly receptors and enzymes. While this pool of available drug targets is limited, there is an untapped potential for discovering new drug-binding mechanisms and modes. For example, enzymes with long binding cavities offer numerous prerequisite binding sites that may be visited by an inhibitor during migration from a bulk solution to the destination site. Drug design can use these prerequisite sites as new structural targets. However, identifying these ephemeral sites is challenging. Here, we introduce a new method called NetBinder for the systematic identification and classification of prerequisite binding sites at atomic resolution. NetBinder is based on atomistic simulations of the full inhibitor binding process and provides a networking framework on which to select the most important binding modes and uncover the entire binding mechanism, including previously undiscovered events. NetBinder was validated by a study of the binding mechanism of blebbistatin (a potent inhibitor) to myosin 2 (a promising target for cancer chemotherapy). Myosin 2 is a good test enzyme because, like other potential targets, it has a long internal binding cavity that provides blebbistatin with numerous potential prerequisite binding sites. The mechanism proposed by NetBinder of myosin 2 structural changes during blebbistatin binding shows excellent agreement with experimentally determined binding sites and structural changes. While NetBinder was tested on myosin 2, it may easily be adopted to other proteins with long internal cavities, such as G-protein-coupled receptors or ion channels, the most popular current drug targets. NetBinder provides a new paradigm for drug design by a network-based elucidation of binding mechanisms at an atomic resolution.

Inside and Out of the Pore: Comparing Interactions and Molecular Dynamics of Influenza A M2 Viroporin Complexes in Standard Lipid Bilayers

Ion channels located at viral envelopes (viroporins) have a critical function for the replication of infectious viruses and are important drug targets. Over the last decade, the number and duration of molecular dynamics (MD) simulations of the influenza A M2 ion channel owing to the increased computational efficiency. Here, we aimed to define the system setup and simulation conditions for the correct description of the protein-pore and the protein-lipid interactions for influenza A M2 in comparison with experimental data. We performed numerous MD simulations of the influenza A M2 protein in complex with adamantane blockers in standard lipid bilayers using OPLS2005 and CHARMM36 (C36) force fields. We explored the effect of varying the M2 construct (M2(22-46) and M2(22-62)), the lipid buffer size and type (stiffer DMPC or softer POPC with or without 20% cholesterol), the simulation time, the H37 protonation site (N_δ or Ν_ε), the conformational state of the W41 channel gate, and M2's cholesterol binding sites (BSs). We report that the 200 ns MD with M2(22-62) (having N_ε Η37) in the 20 Å lipid buffer with the C36 force field accurately describe: (a) the M2 pore structure and interactions inside the pore, that is, adamantane channel blocker location, water clathrate structure, and water or chloride anion blockage/passage from the M2 pore in the presence of a channel blocker and (b) interactions between M2 and the membrane environment as reflected by the calculation of the M2 bundle tilt, folding of amphipathic helices, and cholesterol BSs. Strikingly, we also observed that the C36 1 μs MD simulations using M2(22-62) embedded in a 20 Å POPC:cholesterol (5:1) scrambled membrane produced frequent interactions with cholesterol, which when combined with computational kinetic analysis, revealed the experimentally observed BSs of cholesterol and suggested three similarly long-interacting positions in the top leaflet that have previously not been observed experimentally. These findings promise to be useful for other viroporin systems.

Prerequisite Binding Modes Determine the Dynamics of Action of Covalent Agonists of Ion Channel TRPA1

Transient receptor potential ankyrin 1 (TRPA1) is a transmembrane protein channeling the influx of calcium ions. As a polymodal nocisensor, TRPA1 can be activated by thermal, mechanical stimuli and a wide range of chemically damaging molecules including small volatile environmental toxicants and endogenous algogenic lipids. After activation by such compounds, the ion channel opens up, its central pore widens allowing calcium influx into the cytosol inducing signal transduction pathways. Afterwards, the calcium influx desensitizes irritant evoked responses and results in an inactive state of the ion channel. Recent experimental determination of structures of apo and holo forms of TRPA1 opened the way towards the design of new agonists, which can activate the ion channel. The present study is aimed at the elucidation of binding dynamics of agonists using experimental structures of TRPA1-agonist complexes at the atomic level applying molecular docking and dynamics methods accounting for covalent and non-covalent interactions. Following a test of docking methods focused on the final, holo structures, prerequisite binding modes were detected involving the apo forms. It was shown how reversible interactions with prerequisite binding sites contribute to structural changes of TRPA1 leading to covalent bonding of agonists. The proposed dynamics of action allowed a mechanism-based forecast of new, druggable binding sites of potent agonists.