Button-Push On-Demand Synthesis for Rapid Optimization of Antiviral Peptidomimetics.

The optimization of hit compounds into drug candidates is a pivotal phase in drug discovery but often hampered by cumbersome manual synthesis of derivatives. While automated organic molecule synthesis has enhanced efficiency, safety, and cost-effectiveness, achieving fully automated multistep synthesis remains a formidable challenge due to issues such as solvent and reagent incompatibilities and the accumulation of side-products. We herein demonstrate an automated solid-phase flow platform for synthesizing α-keto-amides and nitrile peptidomimetics, guided by docking simulations, to identify potent broad-spectrum antiviral leads. A compact parallel synthesizer was built in-house, capable of producing 5 distinct molecules per cycle; 525 reactions could be finished within three months to generate 42 derivatives for a structure-activity relationship (SAR) investigation. Among these, ten derivatives exhibited promising target inhibitory activity (IC50 < 100 nM) including two with antiviral activity (EC50 < 250 nM). The platform, coupled with digital chemical recipe files, offers rapid access to a wide range of peptidomimetics, serving as a valuable reservoir for broad-spectrum antiviral candidates. This automated solid-phase flow synthesis approach expedites the generation of previously difficult complex molecular scaffolds. By integration of SPS-flow synthesis with medicinal chemistry campaign, >10-fold target inhibitory activity was achieved from a small set of derivatives, which indicates the potential to shift the paradigm of drug discovery.

Molecular design of hydroxamic acid-based derivatives as urease inhibitors of Helicobacter pylori.

Helicobacter pylori is the main causative agent of gastric cancer, especially non-cardiac gastric cancers. This bacterium relies on urease producing much ammonia to colonize the host. Herein, the study provides valuable insights into structural patterns driving urease inhibition for high-activity molecules designed via exploring known inhibitors. Firstly, an ensemble model was devised to predict the inhibitory activity of novel compounds in an automated workflow (R2 = 0.761) that combines four machine learning approaches. The dataset was characterized in terms of chemical space, including molecular scaffolds, clustering analysis, distribution for physicochemical properties, and activity cliffs. Through these analyses, the hydroxamic acid group and the benzene ring responsible for distinct activity were highlighted. Activity cliff pairs uncovered substituents of the benzene ring on hydroxamic acid derivatives are key structures for substantial activity enhancement. Moreover, 11 hydroxamic acid derivatives were designed, named mol1-11. Results of molecular dynamic simulations showed that the mol9 exhibited stabilization of the active site flap's closed conformation and are expected to be promising drug candidates for Helicobacter pylori infection and further in vitro, in vivo, and clinical trials to demonstrate in future.

Unraveling the influential mechanism of short-chain branching on the crystallization of trimodal polyethylene by molecular dynamics simulation.

Trimodal polyethylene (PE) has become the focus of research in recent years due to its excellent performance. By means of molecular dynamics simulations, we aim to expound the molecular mechanism of short-chain branching (SCB) in the nucleation process, crystallization process and chain entanglement of trimodal PE. In this study, a series of polyethylene models including different short-chain branching concentrations (SCBCs), short-chain branching lengths (SCBLs), and short-chain branching distributions (SCBDs) were considered. The increase of SCBCs greatly reduces the ability of flipping and movement of PE chains, resulting in more time for nucleation and crystallization and a significant reduction of crystallinity. In contrast, an increase in the SCBL only slightly slows down the diffusion rate of the chain, which leads to a little increase in crystallization time. Most important of all, in the study of SCBD, we find that the distribution of SCBs on a high molecular weight chain, which is the characteristic of trimodal PE, is conducive to the chain entanglement and prevents the occurrence of micro phase separation compared with the case where the SCBs are distributed on a medium molecular weight chain. The mechanism of chain entanglement is proposed to explain the effect of SCBs on tie chain entanglement.

Approaching the Dimerization Mechanism of Small Molecule Inhibitors Targeting PD-L1 with Molecular Simulation.

Inhibitors blocking the PD-1/PD-L1 immune checkpoint demonstrate impressive anti-tumor immunity, and small molecule inhibitors disclosed by the Bristol-Myers Squibb (BMS) company have become a hot topic. In this work, by modifying the carbonyl group of BMS-202 into a hydroxyl group to achieve two enantiomers (MS and MR) with a chiral center, we found that this is an effective way to regulate its hydrophobicity and thus to reduce the negative effect of polar solvation free energy, which enhances the stability of PD-L1 dimer/inhibitor complexes. Moreover, we studied the binding modes of BMS-200 and BMS-202-related small molecule inhibitors by molecular dynamics simulation to explore their inhibitory mechanism targeting PD-L1 dimerization. The results showed that the size exclusion effect of the inhibitors triggered the rearrangement of the residue ATyr56, leading to the formation of an axisymmetric tunnel-shaped pocket, which is an important structural basis for improving the binding affinity of symmetric inhibitors with PD-L1. Furthermore, after inhibitor dissociation, the conformation of ATyr123 and BMet115 rearranged, which blocked the entrance of the binding pocket, while the reverse rearrangements of the same residues occurred when the PD-L1 monomer was complexed with the inhibitors, preparing PD-L1 for dimerization. Overall, this study casts a new light on the inhibitory mechanism of BMS inhibitors targeting PD-L1 dimerization and provides an idea for designing novel small molecule inhibitors for future cancer immunotherapy.

Is the Triggering of PD-L1 Dimerization a Potential Mechanism for Food-Derived Small Molecules in Cancer Immunotherapy? A Study by Molecular Dynamics.

Using small molecules to inhibit the PD-1/PD-L1 pathway is an important approach in cancer immunotherapy. Natural compounds such as capsaicin, zucapsaicin, 6-gingerol and curcumin have been proposed to have anticancer immunologic functions by downregulating the PD-L1 expression. PD-L1 dimerization promoted by small molecules was recently reported to be a potential mechanism to inhibit the PD-1/PD-L1 pathway. To clarify the molecular mechanism of such compounds on PD-L1 dimerization, molecular docking and molecular dynamics simulations were performed. The results evidenced that these compounds could inhibit PD-1/PD-L1 interactions by directly targeting PD-L1 dimerization. Binding free energy calculations showed that capsaicin, zucapsaicin, 6-gingerol and curcumin have strong binding ability with the PD-L1 dimer, where the affinities of them follow the trend of zucapsaicin > capsaicin > 6-gingerol ≈ curcumin. Analysis by residue energy decomposition, contact numbers and nonbonded interactions revealed that these compounds have a tight interaction with the C-sheet, F-sheet and G-sheet fragments of the PD-L1 dimer, which were also involved in the interactions with PD-1. Moreover, non-polar interactions between these compounds and the key residues Ile54, Tyr56, Met115 and Ala121 play a key role in stabilizing the protein−ligand complexes in solution, in which the 4'-hydroxy-3'-methoxyphenyl group and the carbonyl group of zucapsaicin, capsaicin, 6-ginger and curcumin were significant for the complexation of small molecules with the PD-L1 dimer. The conformational variations of these complexes were further analyzed by free energy landscape (FEL) and principal component analysis (PCA) and showed that these small molecules could make the structure of dimers more stable. This work provides a mechanism insight for food-derived small molecules blocking the PD-1/PD-L1 pathway via directly targeting the PD-L1 dimerization and offers theoretical guidance to discover more effective small molecular drugs in cancer immunotherapy.

Effect of Fragment 1 on the Binding of Epigallocatechin Gallate to the PD-L1 Dimer Explored by Molecular Dynamics.

Blocking the interaction between programmed cell death-1 (PD-1) and programmed cell death-ligand 1 (PD-L1) by directly targeting the PD-L1 dimer has emerged as a hot topic in the field of cancer immunotherapy. Epigallocatechin gallate (EGCG), a natural product, has been demonstrated binding to the PD-L1 dimer in our previous study, but has a weaker binding capacity, moreover, EGCG is located at the end of the binding pocket of the PD-L1 dimer. The inhibitor fragment 1 (FRA) lies at the other end. So, we proposed that the introduction of FRA might be able to improve the binding ability. To illuminate this issue, molecular dynamics (MD) simulation was performed in the present study. Binding free energy calculations show that the binding affinity is significantly increased by 17 kcal/mol upon the introduction of FRA. It may be due to the energy contributions of emerging key residues ATyr56, AMet115, BTyr123, AIle54 and the enhanced contributions of initial key residues ATyr123 and BVal68. Binding mode and non-bonded interaction results indicate that FRA_EGCG (EGCG in combination with FRA) binds to the C-, F- and G-sheet of the PD-L1 dimer. Importantly, the introduction of FRA mainly strengthened the nonpolar interactions. The free energy landscape and secondary structure results further show that FRA_EGCG can interact with the PD-L1 dimer more stably. These data demonstrated here provide the theoretical basis for screening two or more natural products with additive inhibitory effect on this pathway and therefore exerting more effective anticancer immunity.

Effects of Bisphenol A on reproductive toxicity and gut microbiota dysbiosis in male rats.

Bisphenol A (BPA) is an environmental endocrine disruptor. Recent studies have shown an association between decreased spermatogenesis and gut microbiota alteration. However, the potential associations and mechanisms of BPA exposure on spermatogenesis, hormone production, and gut microbiota remain unknown. This study aims to investigate BPA-induced male reproductive toxicity and the potential link with gut microbiota dysbiosis. Male Sprague Dawley rats were exposed to BPA at different doses by oral gavage for thirty consecutive days. The extent of testicular damage was evaluated by basic parameters of body weight and hematoxylin-eosin (H&E) staining. Next, we determined the mRNA levels and protein levels of apoptosis, histone-related factors, and mammalian target of rapamycin (mTOR) pathway in testes. Finally, 16 S rDNA sequencing was used to analyze gut microbiota composition after BPA exposure. BPA exposure damaged testicular histology, significantly decreased sperm count, and increased sperm abnormalities. In addition, BPA exposure caused oxidative stress and cell apoptosis in testes. The levels of histone (H2A, H3) were significantly increased, while ubiquitin histone H2A (ub-H2A) and ubiquitin histone H2B (ub-H2B) were markedly reduced. Furthermore, BPA activated the PI3K and AKT expression, but the protein expressions of mTOR and 4EBP1 in testes were inhibited significantly. Additionally, the relative abundance of class Gammaproteobacteria, and order Betaproteobacteriales was significantly higher when treated with a high dose of BPA compared to the control group, which was negatively correlated with testosterone level. This study highlights the relationship between BPA-induced reproductive toxicity and gut microbiota disorder and provides new insights into the prevention and treatment of BPA-induced reproductive damage.

Exposure to Bisphenol A Caused Hepatoxicity and Intestinal Flora Disorder in Rats.

Bisphenol A (BPA) is a globally utilized industrial chemical and is commonly used as a monomer of polycarbonate plastics and epoxy resins. Recent research reveals that BPA could cause potential adverse biological effects and liver dysfunction. However, the underlying mechanisms of BPA-induced hepatoxicity and gut dysbiosis remain unclear and deserve further study. In this study, male Sprague Dawley rats were exposed to different doses (0, 30, 90, and 270 mg/kg bw) of BPA by gavage for 30 days. The results showed that the high dose of BPA decreased superoxide dismutase (SOD), glutathione (GSH), and increased malondialdehyde (MDA) levels. Moreover, a high dose of BPA caused a significant increase in serum alanine aminotransferase (ALT), aspartate aminotransferase (AST), total cholesterol (TC), and low-density lipoprotein cholesterol (LDL-C), while high-density lipoprotein cholesterol (HDL-C) was significantly decreased in BPA-treated rats. The gene expression of PGC-1α and Nrf1 were decreased in the liver of high doses of BPA-administrated rats, as well as the protein levels of SIRT1, PGC-1α, Nrf2, and TFAM. However, the protein expression of IL-1ß was significantly increased in BPA-treated rats. In addition, BPA weakened the mitochondrial function of hepatocytes and promoted cell apoptosis in the liver by up-regulating the protein levels of Bax, cleaved-Caspase3, and cleaved-PARP1 while down-regulating the Bcl-2 in the liver. More importantly, a high dose of BPA caused a dramatic change in microbiota structure, as characterized at the genus level by increasing the ratio of Firmicutes to Bacteroidetes (F/B), and the relative abundance of Proteobacteria in feces, while decreasing the relative abundance of Prevotella_9 and Ruminococcaceae_UCG-014, which is positively correlated with the content of short-chain fatty acids (SCFAs). In summary, our data indicated that BPA exposure caused hepatoxicity through apoptosis and the SIRT1/PGC-1α pathway. BPA-induced intestinal flora and SCFA changes may be associated with hepatic damage. The results of this study provide a new sight for the understanding of BPA-induced hepatoxicity.

A head-to-head comparison of the inhibitory activities of 15 peptidomimetic SARS-CoV-2 3CLpro inhibitors.

The COVID-19 pandemic caused by SARS-CoV-2 has created an unprecedented global health emergency. As of July 2021, only three antiviral therapies have been approved by the FDA for treating infected patients, highlighting the urgent need for more antiviral drugs. The SARS-CoV-2 3CL protease (3CLpro) is deemed an attractive drug target due to its essential role in viral polyprotein processing and pathogenesis. Indeed, a number of peptidomimetic 3CLpro inhibitors armed with electrophilic warheads have been reported by various research groups that can potentially be developed for treating COVID-19. However, it is currently impossible to compare their relative potencies due to the different assays employed. To solve this, we conducted a head-to-head comparison of fifteen reported peptidomimetic inhibitors in a standard FRET-based SARS-CoV-2 3CLpro inhibition assay to compare and identify potent inhibitors for development. Inhibitor design and the suitability of various warheads are also discussed.

Fragment-based lead discovery of indazole-based compounds as AXL kinase inhibitors.

AXL is a member of the TAM (TYRO3, AXL, MER) subfamily of receptor tyrosine kinases. It is upregulated in a variety of cancers and its overexpression is associated with poor disease prognosis and acquired drug resistance. Utilizing a fragment-based lead discovery approach, a new indazole-based AXL inhibitor was obtained. The indazole fragment hit 11, identified through a high concentration biochemical screen, was expeditiously improved to fragment 24 by screening our in-house expanded library of fragments (ELF) collection. Subsequent fragment optimization guided by docking studies provided potent inhibitor 54 with moderate exposure levels in mice. X-ray crystal structure of analog 50 complexed with the I650M mutated kinase domain of Mer revealed the key binding interactions for the scaffold. The good potency coupled with reasonable kinase selectivity, moderate in vivo exposure levels, and availability of structural information for the series makes it a suitable starting point for further optimization efforts.

Molecular Mechanism of Small-Molecule Inhibitors in Blocking the PD-1/PD-L1 Pathway through PD-L1 Dimerization.

Programmed cell death-1 (PD-1), which is a molecule involved in the inhibitory signal in the immune system and is important due to blocking of the interactions between PD-1 and programmed cell death ligand-1 (PD-L1), has emerged as a promising immunotherapy for treating cancer. In this work, molecular dynamics simulations were performed on complex systems consisting of the PD-L1 dimer with (S)-BMS-200, (R)-BMS-200 and (MOD)-BMS-200 (i.e., S, R and MOD systems) to systematically evaluate the inhibitory mechanism of BMS-200-related small-molecule inhibitors in detail. Among them, (MOD)-BMS-200 was modified from the original (S)-BMS-200 by replacing the hydroxyl group with a carbonyl to remove its chirality. Binding free energy analysis indicates that BMS-200-related inhibitors can promote the dimerization of PD-L1. Meanwhile, no significant differences were observed between the S and MOD systems, though the R system exhibited a slightly higher energy. Residue energy decomposition, nonbonded interaction, and contact number analyses show that the inhibitors mainly bind with the C, F and G regions of the PD-L1 dimer, while nonpolar interactions of key residues Ile54, Tyr56, Met115, Ala121 and Tyr123 on both PD-L1 monomers are the dominant binding-related stability factors. Furthermore, compared with (S)-BMS-200, (R)-BMS-200 is more likely to form hydrogen bonds with charged residues. Finally, free energy landscape and protein-protein interaction analyses show that the key residues of the PD-L1 dimer undergo remarkable conformational changes induced by (S)-BMS-200, which boosts its intimate interactions. This systematic investigation provides a comprehensive molecular insight into the ligand recognition process, which will benefit the design of new small-molecule inhibitors targeting PD-L1 for use in anticancer therapy.

Molecular Mechanism of Food-Derived Polyphenols on PD-L1 Dimerization: A Molecular Dynamics Simulation Study.

In cancer immunotherapy, an emerging approach is to block the interactions of programmed cell death-1 (PD-1) and programmed cell death-ligand 1 (PD-L1) using small-molecule inhibitors. The food-derived polyphenols curcumin (CC), resveratrol (RSV) and epigallocatechin gallate (EGCG) have anticancer immunologic functions, which, recently, have been proposed to act via the downregulation of PD-L1 expression. However, it remains unclear whether they can directly target PD-L1 dimerization and, thus, interrupt the PD-1/PD-L1 pathway. To elucidate the molecular mechanism of such compounds on PD-L1 dimerization, molecular docking and nanosecond molecular dynamics simulations were performed. Binding free energy calculations show that the affinities of CC, RSV and EGCG to the PD-L1 dimer follow a trend of CC > RSV > EGCG. Hence, CC is the most effective inhibitor of the PD-1/PD-L1 pathway. Analysis on contact numbers, nonbonded interactions and residue energy decomposition indicate that such compounds mainly interact with the C-, F- and G-sheet fragments of the PD-L1 dimer, which are involved in interactions with PD-1. More importantly, nonpolar interactions between these compounds and the key residues Ile54, Tyr56, Met115, Ala121 and Tyr123 play a dominant role in binding. Free energy landscape and secondary structure analyses further demonstrate that such compounds can stably interact with the binding domain of the PD-L1 dimer. The results provide evidence that CC, RSV and EGCG can inhibit PD-1/PD-L1 interactions by directly targeting PD-L1 dimerization. This provides a novel approach to discovering food-derived small-molecule inhibitors of the PD-1/PD-L1 pathway with potential applications in cancer immunotherapy.

Strategic Design of Catalytic Lysine-Targeting Reversible Covalent BCR-ABL Inhibitors*.

Targeted covalent inhibitors have re-emerged as validated drugs to overcome acquired resistance in cancer treatment. Herein, by using a carbonyl boronic acid (CBA) warhead, we report the structure-based design of BCR-ABL inhibitors via reversible covalent targeting of the catalytic lysine with improved potency against both wild-type and mutant ABL kinases, especially ABLT315I bearing the gatekeeper residue mutation. We show the evolutionarily conserved lysine can be targeted selectively, and the selectivity depends largely on molecular recognition of the non-covalent pharmacophore in this class of inhibitors, probably due to the moderate reactivity of the warhead. We report the first co-crystal structures of covalent inhibitor-ABL kinase domain complexes, providing insights into the interaction of this warhead with the catalytic lysine. We also employed label-free mass spectrometry to evaluate off-targets of our compounds at proteome-wide level in different mammalian cells.

Mechanistic investigation of zinc-promoted silylation of phenylacetylene and chlorosilane: a combined experimental and computational study.

The zinc-promoted silylation method is of great importance to synthesize high-performance silicon-containing arylacetylene (PSA) resins in the industry. However, it is difficult to eliminate the accompanied by-product of terminal alkenes due to the lack of mechanistic understanding of the silylation. The initiation of zinc-promoted silylation is facilitated by the interaction between zinc and phenylacetylene. Our DFT calculations indicated that the intermolecular hydrogen transfer of phenylacetylene follows an ionic pathway, which generates a phenylacetylene anion and the corresponding alkene moieties on the zinc surface. The styrene by-product is observed in this stage, with its alkene moieties desorbing as radicals into the solvent under the high reaction temperature. Three possible intermediates of surface phenylacetylene anions were proposed including PhC[triple bond, length as m-dash]C-Zn, PhC[triple bond, length as m-dash]CZnCl, and (PhC[triple bond, length as m-dash]C)2Zn. These carbanion-zinc intermediates undergo an SN2 reaction with Me3SiCl to afford the alkynylsilane on the zinc surface, which is calculated to be the rate-determining step for the zinc-promoted silylation reaction.

Molecular dynamics simulation of shish-kebab crystallization of polyethylene: Unraveling the effects of molecular weight distribution.

By means of molecular dynamics simulations, extensional flow was performed on five polyethylene models with different molecular weight distributions (MWDs) precisely designed in view of Grubbs, metallocene, Ziegler-Natta, and chromium-based catalysts, while ignoring the sequence distributions of short branches to shed light on the molecular mechanism of MWD on shish-kebab formation. The formation of shish-kebab crystallites can be divided into three stages: the emergence of precursors, evolution from precursors to shish nuclei, and the formation of lamellar crystallites. The results demonstrated that the precursors initiated from trans-rich segments with local order and minor crystallinity grew into large shish nuclei and eventually evolved into lamellae. There were more inconsecutively trans-state bonds occurring in long chains rather than in short chains, which promoted an easier transformation from precursors to shish nuclei. Therefore, broader MWDs make positive contributions to the formation of shish nuclei, increase the crystallization speed, and the generation of a more regular, compact, and thicker lamella with less tie molecule fractions, while the final crystallinity is independent of MWD.

A unified model of Grignard reagent formation.

Grignard reagents are among the most fundamental reagents in organic synthesis, yet studies have hitherto failed to fully explain the selectivity and kinetics of Grignard reagent formation (GRF). The present study provides new insights into the intermediates and pathways of GRF using density functional theory (DFT) calculations. Potential energy surfaces of RX dissociation along different directions reveal the origin of configuration retention of alkenyl and aromatic halides. Radical intermediates participate solely in the dissociation stage, and depend on the geometry of the reactant halide. Dissociation of organic halides yields stabilized surface anions, and the rest of the reaction is ionic in nature. MgX+/RMg+ were proposed as the key intermediates of Mg leaving from the surface in the self-activation of GRF, which explains the accelerated kinetics upon addition of RMgX or MgX2. The intermediacy of the cations was supported by a simple electrochemical experiment. To the best of our knowledge, this is the first unified ionic model (I-model) developed for resolving the controversial issues of GRF.

Escherichia coli Topoisomerase IV E Subunit and an Inhibitor Binding Mode Revealed by NMR Spectroscopy.

Bacterial topoisomerases are attractive antibacterial drug targets because of their importance in bacterial growth and low homology with other human topoisomerases. Structure-based drug design has been a proven approach of efficiently developing new antibiotics against these targets. Past studies have focused on developing lead compounds against the ATP binding pockets of both DNA gyrase and topoisomerase IV. A detailed understanding of the interactions between ligand and target in a solution state will provide valuable information for further developing drugs against topoisomerase IV targets. Here we describe a detailed characterization of a known potent inhibitor containing a 9H-pyrimido[4,5-b]indole scaffold against the N-terminal domain of the topoisomerase IV E subunit from Escherichia coli (eParE). Using a series of biophysical and biochemical experiments, it has been demonstrated that this inhibitor forms a tight complex with eParE. NMR studies revealed the exact protein residues responsible for inhibitor binding. Through comparative studies of two inhibitors of markedly varied potencies, it is hypothesized that gaining molecular interactions with residues in the α4 and residues close to the loop of ß1-α2 and residues in the loop of ß3-ß4 might improve the inhibitor potency.

Biophysical Studies of Bacterial Topoisomerases Substantiate Their Binding Modes to an Inhibitor.

Bacterial DNA topoisomerases are essential for bacterial growth and are attractive, important targets for developing antibacterial drugs. Consequently, different potent inhibitors that target bacterial topoisomerases have been developed. However, the development of potent broad-spectrum inhibitors against both Gram-positive (G(+)) and Gram-negative (G(-)) bacteria has proven challenging. In this study, we carried out biophysical studies to better understand the molecular interactions between a potent bis-pyridylurea inhibitor and the active domains of the E-subunits of topoisomerase IV (ParE) from a G(+) strain (Streptococcus pneumoniae (sParE)) and a G(-) strain (Pseudomonas aeruginosa (pParE)). NMR results demonstrated that the inhibitor forms a tight complex with ParEs and the resulting complexes adopt structural conformations similar to those observed for free ParEs in solution. Further chemical-shift perturbation experiments and NOE analyses indicated that there are four regions in ParE that are important for inhibitor binding, namely, α2, the loop between ß2 and α3, and the ß2 and ß6 strands. Surface plasmon resonance showed that this inhibitor binds to sParE with a higher KD than pParE. Point mutations in α2 of ParE, such as A52S (sParE), affected its binding affinity with the inhibitor. Taken together, these results provide a better understanding of the development of broad-spectrum antibacterial agents.

Probing the binding mechanism of Mnk inhibitors by docking and molecular dynamics simulations.

Mitogen-activated protein kinases-interacting kinase 1 and 2 (Mnk1/2) activate the oncogene eukaryotic initiation factor 4E (eIF4E) by phosphorylation. High level of phosphorylated eIF4E is associated with various types of cancers. Inhibition of Mnk prevents eIF4E phosphorylation, making them potential therapeutic targets for cancer. Recently, we have designed and synthesized a series of novel imidazopyridine and imidazopyrazine derivatives that inhibit Mnk1/2 kinases with a potency in the nanomolar to micromolar range. In the current work we model the inhibition of Mnk kinase activity by these inhibitors using various computational approaches. Combining homology modeling, docking, molecular dynamics simulations, and free energy calculations, we find that all compounds bind similarly to the active sites of both kinases with their imidazopyridine and imidazopyrazine cores anchored to the hinge regions of the kinases through hydrogen bonds. In addition, hydrogen bond interactions between the inhibitors and the catalytic Lys78 (Mnk1), Lys113 (Mnk2) and Ser131 (Mnk1), Ser166 (Mnk2) appear to be important for the potency and stability of the bound conformations of the inhibitors. The computed binding free energies (ΔGPred) of these inhibitors are in accord with experimental bioactivity data (pIC50) with correlation coefficients (r(2)) of 0.70 and 0.68 for Mnk1 and Mnk2 respectively. van der Waals energies and entropic effects appear to dominate the binding free energy (ΔGPred) for each Mnk-inhibitor complex studied. The models suggest that the activities of these small molecule inhibitors arise from interactions with multiple residues in the active sites, particularly with the hydrophobic residues.

Characterization of the interaction between Escherichia coli topoisomerase IV E subunit and an ATP competitive inhibitor.

Bacterial topoisomerase IV (ParE) is essential for DNA replication and serves as an attractive target for antibacterial drug development. The X-ray structure of the N-terminal 24 kDa ParE, responsible for ATP binding has been solved. Due to the accessibility of structural information of ParE, many potent ParE inhibitors have been discovered. In this study, a pyridylurea lead molecule against ParE of Escherichia coli (eParE) was characterized with a series of biochemical and biophysical techniques. More importantly, solution NMR analysis of compound binding to eParE provides better understanding of the molecular interactions between the inhibitor and eParE.