OPLS4: Improving Force Field Accuracy on Challenging Regimes of Chemical Space

Molecular interactions between a diphenyl scaffold and PED/PEA15: Implications for type II diabetes therapeutics targeting PED/PEA15 - Phospholipase D1 interaction

In a recent study, we have identified BPH03 as a promising scaffold for the development of compounds aimed at modulating the interaction between PED/PEA15 (Phosphoprotein Enriched in Diabetes/Phosphoprotein Enriched in Astrocytes 15) and PLD1 (phospholipase D1), with potential applications in type II diabetes therapy. PED/PEA15 is known to be overexpressed in certain forms of diabetes, where it binds to PLD1, thereby reducing insulin-stimulated glucose transport. The inhibition of this interaction reestablishes basal glucose transport, indicating PED as a potential target of ligands capable to recover glucose tolerance and insulin sensitivity. In this study, we employ computational methods to provide a detailed description of BPH03 interaction with PED, evidencing the presence of a hidden druggable pocket within its PLD1 binding surface. We also elucidate the conformational changes that occur during PED interaction with BPH03. Moreover, we report new NMR data supporting the in-silico findings and indicating that BPH03 disrupts the PED/PLD1 interface displacing PLD1 from its interaction with PED. Our study represents a significant advancement toward the development of potential therapeutics for the treatment of type II diabetes.

Synthesis, carbonic anhydrase inhibition studies and modelling investigations of phthalimide-hydantoin hybrids

A novel series of hydantoins incorporating phthalimides has been synthesised by condensation of activated phthalimides with 1-aminohydantoin and investigated for their inhibitory activity against a panel of human (h) carbonic anhydrase (CA, EC 4.2.1.1): the cytosolic isoforms hCA I, hCA II, and hCA VII, secreted isoform hCA VI, and the transmembrane hCA IX, by a stopped-flow CO2 hydrase assay. Although all newly developed compounds were totally inactive on hCA I and mainly ineffective towards hCA II, they generally exhibited moderate repressing effects on hCA VI, VII, and IX with KIs values in the submicromolar to micromolar ranges. The salts 3a and 3b, followed by derivative 5, displayed the best inhibitory activity of all the evaluated compounds and their binding mode was proposed in silico. These compounds can also be considered interesting starting points for the development of novel pharmacophores for this class of enzyme inhibitors.

Rational engineering of homospermidine synthase for enhanced catalytic efficiency toward spermidine synthesis

Spermidine is a naturally occurring polyamine widely utilized in the prevention and treatment of various diseases. Current spermidine biosynthetic methods have problems such as low efficiency and complex multi-enzyme catalysis. Based on sequence-structure-function relationships, we engineered the widely studied homospermidine synthase from Blastochloris viridis (BvHSS) and obtained mutants that could catalyze the production of spermidine from 1,3-diaminopropane and putrescine. The specific activities of BvHSS and the mutants D361E and E232D + D361E (E232D-D) were 8.72, 46.04 and 48.30 U/mg, respectively. The optimal pH for both mutants was 9.0, and the optimal temperature was 50 °C. Molecular docking and dynamics simulations revealed that mutating aspartic acid at position 361 to glutamic acid narrowed the substrate binding pocket, promoting stable spermidine production. Conversely, mutating glutamic acid at position 232 to aspartic acid enlarged the substrate channel entrance, facilitating substrate entry into the active pocket and enhancing spermidine generation. In whole-cell catalysis lasting 6 h, D361E and E232D-D synthesized 725.3 and 933.5 mg/L of spermidine, respectively. This study offers a practical approach for single-enzyme catalyzed spermidine synthesis and sheds light on the crucial residues influencing homospermidine synthase catalytic activity in spermidine production.

Exploring Asphodelus microcarpus as a source of xanthine oxidase inhibitors: Insights from in silico and in vitro studies

Xanthine oxidase (XO) plays a critical role in purine catabolism, catalyzing the conversion of hypoxanthine to xanthine and xanthine to uric acid, contributing to superoxide anion production. This process is implicated in various human diseases, particularly gout. Traditional XO inhibitors, such as allopurinol and febuxostat, while effective, may present side effects. Our study focuses on Asphodelus microcarpus, a plant renowned for traditional anti-inflammatory uses. Recent investigations into its phenolic-rich flowers, notably abundant in luteolin derivatives, reveal its potential as a natural source of XO inhibitors. In the present research, XO inhibition by an ethanolic flowers extract from A. microcarpus is reported. In silico docking studies have highlighted luteolin derivatives as potential XO inhibitors, and molecular dynamics support that luteolin 7-O-glucoside has the highest binding stability compared to other compounds and controls. In vitro studies confirm that luteolin 7-O-glucoside inhibits XO more effectively than the standard inhibitor allopurinol, with an IC50 value of 4.8 μg/mL compared to 11.5 μg/mL, respectively. These findings underscore the potential therapeutic significance of A. microcarpus in managing conditions related to XO activity. The research contributes valuable insights into the health-promoting properties of A. microcarpus and its potential application in natural medicine, presenting a promising avenue for further exploration in disease management.

In silico modelling of ciprofloxacin specific aptamer for the development of high-performance biosensor

Ciprofloxacin (CFX), a widely used fluoroquinolone antibiotic, is critical in healthcare settings for treating patients. However, improper treatment of wastewater from these facilities can lead to environmental contamination with CFX. This underscores the need for an efficient, straightforward method for early detection. In this study, a DNA aptamer was selected through a hierarchical docking workflow, and the stability and interactions were assessed by Molecular Dynamics (MD) simulation. The aptamer-CFX complex that showed the most promise had a docking score of -8.596 kcal/mol and was further analyzed using MD simulation and MM/PBSA. Based on the overall results, the identified ssDNA sequence length of 60 nt (CAGCGCTAGGGCTTTTAGCGTAATGGGTAGGGTGGTGCGGTGCAGATATCGGAATTGGTG) was immobilized over a gold transducer surface through the self-assembled monolayer (SAM; Au-S-ssDNA) method. The ssDNA-modified surface has demonstrated a high affinity towards CFX, which is confirmed by cyclic voltammogram (CV) and electrochemical impedance spectroscopy measurements (EIS). The DNA-aptamer modified electrode demonstrated a good linear range (10 × 10-9 - 200 × 10-9 M), detection limit (1.0 × 10-9 M), selectivity, reproducibility, and stability. The optimized DNA-aptamer-based CFX sensor was further utilized for the accurate determination of CFX with good recoveries in real samples.

Development of Penicillin-Based Carbonic Anhydrase Inhibitors Targeting Multidrug-Resistant Neisseria gonorrhoeae

The development of antibacterial drugs with new mechanisms of action is crucial in combating the rise of antibiotic-resistant infections. Bacterial carbonic anhydrases (CAs, EC 4.2.1.1) have been validated as promising antibacterial targets against pathogens such as Helicobacter pylori, Neisseria gonorrhoeae, and vancomycin-resistant enterococci. A multitarget strategy is proposed to design penicillin-based CA inhibitor hybrids for tackling resistance by targeting multiple bacterial pathways, thereby resensitizing drug-resistant strains to clinical antibiotics. The sulfonamide derivatives potently inhibited the CAs from N. gonorrhoeae and Escherichia coli with KI values in the range of 7.1-617.2 nM. Computational simulations with the main penicillin-binding protein (PBP) of N. gonorrhoeae indicated that these hybrid derivatives maintained the mechanism of action of the lead β-lactams. A subset of derivatives showed potent PBP-related antigonococcal effects against multidrug-resistant N. gonorrhoeae strains, with several compounds significantly outperforming both the lead β-lactam and CA inhibitor drugs (MIC values in the range 0.25 to 0.5 μg/mL).

Structural basis for urate recognition and apigenin inhibition of human GLUT9

Urate, the physiological form of uric acid and a potent antioxidant in serum, plays a pivotal role in scavenging reactive oxygen species. Yet excessive accumulation of urate, known as hyperuricemia, is the primary risk factor for the development of gout. The high-capacity urate transporter GLUT9 represents a promising target for gout treatment. Here, we present cryo-electron microscopy structures of human GLUT9 in complex with urate or its inhibitor apigenin at overall resolutions of 3.5 Å and 3.3 Å, respectively. In both structures, GLUT9 exhibits an inward open conformation, wherein the substrate binding pocket faces the intracellular side. These structures unveil the molecular basis for GLUT9's substrate preference of urate over glucose, and show that apigenin acts as a competitive inhibitor by occupying the substrate binding site. Our findings provide critical information for the development of specific inhibitors targeting GLUT9 as potential therapeutics for gout and hyperuricemia.

Implementation of Girsanov Reweighting in OpenMM and Deeptime

Classical molecular dynamics (MD) simulations provide invaluable insights into complex molecular systems but face limitations in capturing phenomena occurring on time scales beyond their reach. To bridge this gap, various enhanced sampling techniques have been developed, which are complemented by reweighting techniques to recover the unbiased dynamics. Girsanov reweighting is a reweighting technique that reweights simulation paths, generated by a stochastic MD integrator, without evoking an effective model of the dynamics. Instead, it calculates the relative path probability density at the time resolution of the MD integrator. Efficient implementation of Girsanov reweighting requires that the reweighting factors are calculated on-the-fly during the simulations and thus needs to be implemented within the MD integrator. Here, we present a comprehensive guide for implementing Girsanov reweighting into MD simulations. We demonstrate the implementation in the MD simulation package OpenMM by extending the library openmmtools. Additionally, we implemented a reweighted Markov state model estimator within the time series analysis package Deeptime.

A chemical modification of a peroxisome proliferator-activated receptor pan agonist produced a shift to a new dual alpha/gamma partial agonist endowed with mitochondrial pyruvate carrier inhibition and antidiabetic properties

New analogs of the PPAR pan agonist AL29-26 encompassed ligand (S)-7 showing potent activation of PPARα and -γ subtypes as a partial agonist. In vitro experiments and docking studies in the presence of PPAR antagonists were performed to help interpretation of biological data and investigate the main interactions at the binding sites. Further in vitro experiments showed that (S)-7 induced anti-steatotic effects and enhancement of the glucose uptake. This latter effect could be partially ascribed to a significant inhibition of the mitochondrial pyruvate carrier demonstrating that (S)-7 also acted through insulin-independent mechanisms. In vivo experiments showed that this compound reduced blood glucose and lipid levels in a diabetic mice model displaying no toxicity on bone, kidney, and liver. To our knowledge, this is the first example of dual PPARα/γ partial agonist showing these combined effects representing, therefore, the potential lead of new drugs for treatment of dyslipidemic type 2 diabetes.

Structure-Activity Relationships, Deuteration, and Fluorination of Synthetic Cannabinoid Receptor Agonists Related to AKB48, 5F-AKB-48, and AFUBIATA

Synthetic cannabinoid receptor agonists (SCRAs) are a growing class of new psychoactive substances (NPS) commonly derived from an N-alkylated indole, indazole, or 7-azaindole scaffold. Diversification of this core (at the 3-position) with amide-linked pendant amino acid groups and modular N-alkylation (of the indole/indazole/7-azaindole core) ensures that novel SCRAs continue to enter the illicit drug market rapidly. In response to the large number of SCRAs that have been detected, pharmacological evaluation of this NPS class has become increasingly common. Adamantane-derived SCRAs have consistently appeared throughout the market since 2011, and as such, a systematic set of these derivatives was synthesized and pharmacologically evaluated. Deuterated and fluorinated adamantane derivatives were prepared to evaluate typical hydrogen bioisosteres, as well as evaluation of the newly detected AFUBIATA.

Atomistic description of the OCTN1 recognition mechanism via in silico methods

The Organic Cation Transporter Novel 1 (OCTN1), also known as SLC22A4, is widely expressed in various human tissues, and involved in numerous physiological and pathological processes remains. It facilitates the transport of organic cations, zwitterions, with selectivity for positively charged solutes. Ergothioneine, an antioxidant compound, and acetylcholine (Ach) are among its substrates. Given the lack of experimentally solved structures of this protein, this study aimed at generating a reliable 3D model of OCTN1 to shed light on its substrate-binding preferences and the role of sodium in substrate recognition and transport. A chimeric model was built by grafting the large extracellular loop 1 (EL1) from an AlphaFold-generated model onto a homology model. Molecular dynamics simulations revealed domain-specific mobility, with EL1 exhibiting the highest impact on overall stability. Molecular docking simulations identified cytarabine and verapamil as highest affinity ligands, consistent with their known inhibitory effects on OCTN1. Furthermore, MM/GBSA analysis allowed the categorization of substrates into weak, good, and strong binders, with molecular weight strongly correlating with binding affinity to the recognition site. Key recognition residues, including Tyr211, Glu381, and Arg469, were identified through interaction analysis. Ach demonstrated a low interaction energy, supporting the hypothesis of its one-directional transport towards to outside of the membrane. Regarding the role of sodium, our model suggested the involvement of Glu381 in sodium binding. Molecular dynamics simulations of systems at increasing levels of Na+ concentrations revealed increased sodium occupancy around Glu381, supporting experimental data associating Na+ concentration to molecule transport. In conclusion, this study provides valuable insights into the 3D structure of OCTN1, its substrate-binding preferences, and the role of sodium in the recognition. These findings contribute to the understanding of OCTN1 involvement in various physiological and pathological processes and may have implications for drug development and disease management.

Classical molecular dynamics simulation identifies catechingallate as a promising antiviral polyphenol against MPOX palmitoylated surface protein

Cumulative global prevalence of the emergent monkeypox (MPX) infection in the non-endemic countries has been professed as a global public health predicament. Lack of effective MPX-specific treatments sets the baseline for designing the current study. This research work uncovers the effective use of known antiviral polyphenols against MPX viral infection, and recognises their mode of interaction with the target F13 protein, that plays crucial role in formation of enveloped virions. Herein, we have employed state-of-the-art machine learning based AlphaFold2 to predict the three-dimensional structure of F13 followed by molecular docking and all-atoms molecular dynamics (MD) simulations to investigate the differential mode of F13-polyphenol interactions. Our extensive computational approach identifies six potent polyphenols Rutin, Epicatechingallate, Catechingallate, Quercitrin, Isoquecitrin and Hyperoside exhibiting higher binding affinity towards F13, buried inside a positively charged binding groove. Intermolecular contact analysis of the docked and MD simulated complexes divulges three important residues Asp134, Ser137 and Ser321 that are observed to be involved in ligand binding through hydrogen bonds. Our findings suggest that ligand binding induces minor conformational changes in F13 to affect the conformation of the binding site. Concomitantly, essential dynamics of the six-MD simulated complexes reveals Catechin gallate, a known antiviral agent as a promising polyphenol targeting F13 protein, dominated with a dense network of hydrophobic contacts. However, assessment of biological activities of these polyphenols need to be confirmed through in vitro and in vivo assays, which may pave the way for development of new novel antiviral drugs.

Discovering potential WRN inhibitors from natural product database through computational methods

Microsatellite instability (MSI) is a relatively common feature associated with multiple cancers, and Werner syndrome (WRN) ATP-dependent helicase has been recognized as a novel target for treating MSI cancers, such as colorectal cancer. A small-molecule inhibitor targeting WRN would be a promising strategy for treating colorectal cancer with high MSI expression. In this study, we employed a computer-assisted drug discovery strategy to screen over 30,000 natural product molecules. By using a combination of docking, ligand efficiency, Molecular Mechanics/Generalized Born Surface Area (MM/GBSA), and thermodynamic integration (TI) calculations, we identified MOL008980, MOL010740, MOL011832, T4743, TN1166, and TNP-002173 as potential WRN inhibitors. Subsequent molecular dynamics simulation revealed that these screened natural products possessed better binding dynamic characteristics than ATP substrate and were capable of inhibiting the dynamic process of WRN, making them potential strong ATP competitive inhibitors. In conclusion, our computational approach revealed potential WRN inhibitors from a natural product database, providing a theoretical basis for future research.

An emerging terpolymeric nanoparticle pore former as an internal recrystallization inhibitor of celecoxib in controlled release amorphous solid dispersion beads: Experimental studies and molecular dynamics analysis

Solid oral controlled release formulations feature numerous clinical advantages for drug candidates with adequate solubility and dissolution rate. However, most new chemical entities exhibit poor water solubility, and hence are exempt from such benefits. Although combining drug amorphization with controlled release formulation is promising to elevate drug solubility, like other supersaturating systems, the problem of drug recrystallization has yet to be resolved, particularly within the dosage form. Here, we explored the potential of an emerging, non-leachable terpolymer nanoparticle (TPN) pore former as an internal recrystallization inhibitor within controlled release amorphous solid dispersion (CRASD) beads comprising a poorly soluble drug (celecoxib) reservoir and insoluble polymer (ethylcellulose) membrane. Compared to conventional pore former, polyvinylpyrrolidone (PVP), TPN-containing membranes exhibited superior structural integrity, less crystal formation at the CRASD bead surface, and greater extent of celecoxib release. All-atom molecular dynamics analyses revealed that in the presence of TPN, intra-molecular bonding, crystal formation tendency, diffusion coefficient, and molecular flexibility of celecoxib were reduced, while intermolecular H-bonding was increased as compared to PVP. This work suggests that selection of a pore former that promotes prolonged molecular separation within a nanoporous controlled release membrane structure may serve as an effective strategy to enhance amorphicity preservation inside CRASD.

Developing insulin-like peptide 5-based antagonists for the G protein-coupled receptor, RXFP4

Human insulin-like peptide 5 (INSL5) is a gut hormone produced by colonic L-cells, and its biological functions are mediated by Relaxin Family Peptide Receptor 4 (RXFP4). Our preliminary data indicated that RXFP4 agonists are potential drug leads for the treatment of constipation. More recently, we designed and developed a novel RXFP4 antagonist, A13-nR that was shown to block agonist-induced activity in cells and animal models. We showed that A13-nR was able to block agonist-induced increases in colon motility in mice of both genders that express the receptor, RXFP4. Our data also showed that colorectal propulsion induced by intracolonic administration of short-chain fatty acids was antagonized by A13-nR. Therefore, A13-nR is an important research tool and potential drug lead for the treatment of colon motility disorders, such as bacterial diarrhea. However, A13-nR acted as a partial agonist at high concentrations in vitro and demonstrated modest antagonist potency (∼35 nM). Consequently, the primary objective of this study is to pinpoint novel modifications to A13-nR that eliminate partial agonist effects while preserving or augmenting antagonist potency. In this work, we detail the creation of a series of A13-nR-modified analogues, among which analogues 3, 4, and 6 demonstrated significantly improved RXFP4 affinity (∼3 nM) with reduced partial agonist activity, enhanced antagonist potency (∼10 nM) and maximum agonist inhibition (∼80 %) when compared with A13-nR. These compounds have potential as candidates for further preclinical evaluations, marking a significant stride toward innovative therapeutics for colon motility disorders.

Bioactive compounds from Ocimum tenuiflorum and Poria cocos: A novel natural Compound for insomnia treatment based on A computational approach

Insomnia, a widespread public health issue, is associated with substantial distress and daytime functionality impairments and can predispose to depression and cardiovascular disease. Cognitive Behavioral Anti-insomnia therapies including benzodiazepines often face limitations due to patient adherence or potential adverse effects. This study focused on identifying novel bioactive compounds from medicinal plants, aiming to discover and develop new therapeutic agents with low risk-to-benefit ratios using computational drug discovery methods. Through a systematic framework involving compound library preparation, evaluation of drug-likeness and pharmacokinetics, toxicity prediction, molecular docking, and molecular dynamic simulations, two natural compounds such as 2-(4-hydroxy-3-methoxyphenyl)-8-methoxy-6-prop-2-enyl-3,4-dihydro-2H-chromen-3-ol from Ocimum tenuiflorum and 7-(2-hydroxypropan-2-yl)-1,4a-dimethyl-9-oxo-3,4,10,10a-tetrahydro-2H-phenanthrene-1-carboxylic acid from Poria cocos exhibited high binding affinity with orexin receptor type 1 (OX1R) and type 2 (OX2R), surpassing commercial drugs used in insomnia treatment. Additionally, they showed interactions with critical amino acid residues within the receptors that play crucial roles in competitive inhibitor activity, like commercial drugs such as Suvorexant, Lemborexant, and Daridorexant. Further, molecular dynamics simulations of the protein-ligand complexes under conditions that mimic the in vivo environment revealed both compounds' sustained and robust interactions with the OX1R and OX2R, reinforcing their potential as effective therapeutic candidates. Furthermore, upon evaluating both compounds' drug-likeness, pharmacokinetics, and toxicity profiles, it was discerned that they displayed considerable drug-like properties and favorable pharmacokinetics, along with diminished toxicity. The research provides a solid foundation for further exploring and validating these compounds as potential anti-insomnia therapeutics.

Identification of FTY720 and COH29 as novel topoisomerase I catalytic inhibitors by experimental and computational studies

The development of novel topoisomerase I (TOP1) inhibitors is crucial for overcoming the drawbacks and limitations of current TOP1 poisons. Here, we identified two potential TOP1 inhibitors, namely, FTY720 (a sphingosine 1-phosphate antagonist) and COH29 (a ribonucleotide reductase inhibitor), through experimental screening of known active compounds. Biological experiments verified that FTY720 and COH29 were nonintercalative TOP1 catalytic inhibitors that did not induce the formation of DNA-TOP1 covalent complexes. Molecular docking revealed that FTY720 and COH29 interacted favorably with TOP1. Molecular dynamics simulations revealed that FTY720 and COH29 could affect the catalytic domain of TOP1, thus resulting in altered DNA-binding cavity size. The alanine scanning and interaction entropy identified Arg536 as a hotspot residue. In addition, the bioinformatics analysis predicted that FTY720 and COH29 could be effective in treating malignant breast tumors. Biological experiments verified their antitumor activities using MCF-7 breast cancer cells. Their combinatory effects with TOP1 poisons were also investigated. Further, FTY720 and COH29 were found to cause less DNA damage compared with TOP1 poisons. The findings provide reliable lead compounds for the development of novel TOP1 catalytic inhibitors and offer new insights into the potential clinical applications of FTY720 and COH29 in targeting TOP1.

Identifying and Assessing Putative Allosteric Sites and Modulators for CXCR4 Predicted through Network Modeling and Site Identification by Ligand Competitive Saturation

The chemokine receptor CXCR4 is a critical target for the treatment of several cancer types and HIV-1 infections. While orthosteric and allosteric modulators have been developed targeting its extracellular or transmembrane regions, the intramembrane region of CXCR4 may also include allosteric binding sites suitable for the development of allosteric drugs. To investigate this, we apply the Gaussian Network Model (GNM) to the monomeric and dimeric forms of CXCR4 to identify residues essential for its local and global motions located in the hinge regions of the protein. Residue interaction network (RIN) analysis suggests hub residues that participate in allosteric communication throughout the receptor. Mutual residues from the network models reside in regions with a high capacity to alter receptor dynamics upon ligand binding. We then investigate the druggability of these potential allosteric regions using the site identification by ligand competitive saturation (SILCS) approach, revealing two putative allosteric sites on the monomer and three on the homodimer. Two screening campaigns with Glide and SILCS-Monte Carlo docking using FDA-approved drugs suggest 20 putative hit compounds including antifungal drugs, anticancer agents, HIV protease inhibitors, and antimalarial drugs. In vitro assays considering mAB 12G5 and CXCL12 demonstrate both positive and negative allosteric activities of these compounds, supporting our computational approach. However, in vivo functional assays based on the recruitment of β-arrestin to CXCR4 do not show significant agonism and antagonism at a single compound concentration. The present computational pipeline brings a new perspective to computer-aided drug design by combining conformational dynamics based on network analysis and cosolvent analysis based on the SILCS technology to identify putative allosteric binding sites using CXCR4 as a showcase.

Discovery of potent CSK inhibitors through integrated virtual screening and molecular dynamic simulation

Oncogenic overexpression or activation of C-terminal Src kinase (CSK) has been shown to play an important role in triple-negative breast cancer (TNBC) progression, including tumor initiation, growth, metastasis, drug resistance. This revelation has pivoted the focus toward CSK as a potential target for novel treatments. However, until now, there are few inhibitors designed to target the CSK protein. Responding to this, our research has implemented a comprehensive virtual screening protocol. By integrating energy-based screening methods with AI-driven scoring functions, such as Attentive FP, and employing rigorous rescoring methods like Glide docking and molecular mechanics generalized Born surface area (MM/GBSA), we have systematically sought out inhibitors of CSK. This approach led to the discovery of a compound with a potent CSK inhibitory activity, reflected by an IC50 value of 1.6 nM under a homogeneous time-resolved fluorescence (HTRF) bioassay. Subsequently, molecule 2 exhibits strong growth inhibition of MD anderson - metastatic breast (MDA-MB) -231, Hs578T, and SUM159 cells, showing a level of growth inhibition comparable to that observed with dasatinib. Treatment with molecule 2 also induced significant G1 phase accumulation and cell apoptosis. Furthermore, we have explored the explicit binding interactions of the compound with CSK using molecular dynamics simulations, providing valuable insights into its mechanism of action.

Bromo-substituted indirubins for inhibition of protein kinase-mediated signalling involved in inflammatory mediator release in human monocytes

Targeting protein kinases that regulate signalling pathways in inflammation is an effective pharmacological approach to alleviate uncontrolled inflammatory diseases. In this context, the natural product indirubin and its 6-bromo-substituted analogue 6-bromoindirubin-3 -glycerol-oxime ether (6BIGOE; 1) were identified as potent inhibitors of glycogen synthase kinase-3β (GSK-3β). These inhibitors suppress the release of pro-inflammatory cytokines and prostaglandins (PG) from human monocytes. However, indirubin derivatives target several protein kinases such as cyclin-dependent kinases (CDKs) which has been a major concern for their application in inflammation therapy. Here, we report on a library of 13 5-bromo-substituted indirubin derivatives that have been designed to improve potency and target selectivity. Side-by-side comparison of reference compound 1 (6BIGOE) with 5-bromo derivatives revealed its isomer 2 (5BIGOE), as the most potent derivative able to supress pro-inflammatory cytokine and PG release in lipopolysaccharide-stimulated human monocytes. Analysis of protein kinase inhibition in intact monocytes, supported by our in silico findings, proposed higher selectivity of 1 for GSK-3β inhibition with lesser potency against CDKs 8 and 9. In contrast, 2 supressed the activity of these CDKs with higher effectiveness than GSK-3β, representing additional targets of indirubins within the inflammatory response. Encapsulation of 1 and 2 into polymer-based nanoparticles (NP) improved their pharmacological potential. In conclusion, the 5- and 6-brominated indirubins 1 and 2 as dual GSK-3β and CDK8/9 inhibitors represent a novel concept for intervention with inflammatory disorders.

Discovery of a cell-active chikungunya virus nsP2 protease inhibitor using a covalent fragment-based screening approach

Chikungunya virus (CHIKV) is a mosquito-borne alphavirus that has been responsible for numerous large-scale outbreaks in the last twenty years. Currently, there are no FDA-approved therapeutics for any alphavirus infection. CHIKV non-structural protein 2 (nsP2), which contains a cysteine protease domain, is essential for viral replication, making it an attractive target for a drug discovery campaign. Here, we optimized a CHIKV nsP2 protease (nsP2pro) biochemical assay for the screening of a 6,120-compound cysteine-directed covalent fragment library. Using a 50% inhibition threshold, we identified 153 hits (2.5% hit rate). In dose-response follow up, RA-0002034, a covalent fragment that contains a vinyl sulfone warhead, inhibited CHIKV nsP2pro with an IC 50 of 58 ± 17 nM, and further analysis with time-dependent inhibition studies yielded a k inact /K I of 6.4 x 10 3 M -1 s -1 . LC-MS/MS analysis determined that RA-0002034 covalently modified the catalytic cysteine in a site-specific manner. Additionally, RA-0002034 showed no significant off-target reactivity against a panel of cysteine proteases. In addition to the potent biochemical inhibition of CHIKV nsP2pro activity and exceptional selectivity, RA-0002034 was tested in cellular models of alphavirus infection and effectively inhibited viral replication of both CHIKV and related alphaviruses. This study highlights the discovery and characterization of the chemical probe RA-0002034 as a promising hit compound from covalent fragment-based screening for development toward a CHIKV or pan-alphavirus therapeutic.

Significance Statement

Chikungunya virus is one of the most prominent and widespread alphaviruses and has caused explosive outbreaks of arthritic disease. Currently, there are no FDA-approved drugs to treat disease caused by chikungunya virus or any other alphavirus-caused infection. Here, we report the discovery of a covalent small molecule inhibitor of chikungunya virus nsP2 protease activity and viral replication of four diverse alphaviruses. This finding highlights the utility of covalent fragment screening for inhibitor discovery and represents a starting point towards the development of alphavirus therapeutics targeting nsP2 protease.

A CARM1 Inhibitor Potently Suppresses Breast Cancer Both In Vitro and In Vivo

CARM1, belonging to the protein arginine methyltransferase (PRMT) family, is intricately associated with the progression of cancer and is viewed as a promising target for both cancer diagnosis and therapy. However, the number of specific and potent CARM1 inhibitors is limited. We herein discovered a CARM1 inhibitor, iCARM1, that showed better specificity and activity toward CARM1 compared to the known CARM1 inhibitors, EZM2302 and TP-064. Similar to CARM1 knockdown, iCARM1 suppressed the expression of oncogenic estrogen/ERα-target genes, whereas activated type I interferon (IFN) and IFN-induced genes (ISGs) in breast cancer cells. Consequently, iCARM1 potently suppressed breast cancer cell growth both in vitro and in vivo. The combination of iCARM1 with either endocrine therapy drugs or etoposide demonstrated synergistic effects in inhibiting the growth of breast tumors. In summary, targeting CARM1 by iCARM1 effectively suppresses breast tumor growth, offering a promising therapeutic approach for managing breast cancers in clinical settings.

Identification of 1,3,8-Triazaspiro[4.5]Decane-2,4-Dione Derivatives as a Novel δ Opioid Receptor-Selective Agonist Chemotype

δ opioid receptors (DORs) hold potential as a target for neurologic and psychiatric disorders, yet no DOR agonist has proven efficacious in critical phase II clinical trials. The exact reasons for the failure to produce quality drug candidates for the DOR are unclear. However, it is known that certain DOR agonists can induce seizures and exhibit tachyphylaxis. Several studies have suggested that those adverse effects are more prevalent in delta agonists that share the (+)-4-[(αR)-α-((2S,5R)-4-allyl-2,5-dimethyl-1-piperazinyl)-3-methoxybenzyl]-N,N-diethylbenzamide (SNC80)/4-[(αR*)-α-((2S*,5R*)-4-allyl-2,5-dimethyl-1-piperazinyl)-3-hydroxybenzyl]-N,N-diethylbenzamide chemotype. There is a need to find novel lead candidates for drug development that have improved pharmacological properties to differentiate them from the current failed delta agonists. Our objective in this study was to identify novel DOR agonists. We used a β-arrestin assay to screen a small G-protein coupled receptors (GPCR)-focused chemical library. We identified a novel chemotype of DOR agonists that appears to bind to the orthosteric site based of docking and molecular dynamic simulation. The most potent agonist hit compound is selective for the DOR over a panel of 167 other GPCRs, is slightly biased toward G-protein signaling and has anti-allodynic efficacy in a complete Freund's adjuvant model of inflammatory pain in C57BL/6 male and female mice. The newly discovered chemotype contrasts with molecules like SNC80 that are highly efficacious β-arrestin recruiters and may suggest this novel class of DOR agonists could be expanded on to develop a clinical candidate drug. SIGNIFICANCE STATEMENT: δ opioid receptors are a clinical target for various neurological disorders, including migraine and chronic pain. Many of the clinically tested delta opioid agonists share a single chemotype, which carries risks during drug development. Through a small-scale high-throughput screening assay, this study identified a novel δ opioid receptor agonist chemotype, which may serve as alternative for the current analgesic clinical candidates.

Identification and Characterization of a Bacterial Periplasmic Solute Binding Protein That Binds l-Amino Acid Amides

Periplasmic solute-binding proteins (SBPs) are key ligand recognition components of bacterial ATP-binding cassette (ABC) transporters that allow bacteria to import nutrients and metabolic precursors from the environment. Periplasmic SBPs comprise a large and diverse family of proteins, of which only a small number have been empirically characterized. In this work, we identify a set of 610 unique uncharacterized proteins within the SBP_bac_5 family that are found in conserved operons comprising genes encoding (i) ABC transport systems and (ii) putative amidases from the FmdA_AmdA family. From these uncharacterized SBP_bac_5 proteins, we characterize a representative periplasmic SBP from Mesorhizobium sp. A09 (MeAmi_SBP) and show that MeAmi_SBP binds l-amino acid amides but not the corresponding l-amino acids. An X-ray crystal structure of MeAmi_SBP bound to l-serinamide highlights the residues that impart distinct specificity for l-amino acid amides and reveals a structural Ca2+ binding site within one of the lobes of the protein. We show that the residues involved in ligand and Ca2+ binding are conserved among the 610 SBPs from experimentally uncharacterized FmdA_AmdA amidase-associated ABC transporter systems, suggesting these homologous systems are also likely to be involved in the sensing, uptake, and metabolism of l-amino acid amides across many Gram-negative nitrogen-fixing soil bacteria. We propose that MeAmi_SBP is involved in the uptake of such solutes to supplement pathways such as the citric acid cycle and the glutamine synthetase-glutamate synthase pathway. This work expands our currently limited understanding of microbial interactions with l-amino acid amides and bacterial nitrogen utilization.

Structure-based design of multitargeting ChEs-MAO B inhibitors based on phenyl ring bioisosteres: AChE/BChE selectivity switch and drug-like characterization

A structure-based drug design approach was focused on incorporating phenyl ring heterocyclic bioisosteres into coumarin derivative 1, previously reported as potent dual AChE-MAO B inhibitor, with the aim of improving drug-like features. Structure-activity relationships highlighted that bioisosteric rings were tolerated by hMAO B enzymatic cleft more than hAChE. Interestingly, linker homologation at the basic nitrogen enabled selectivity to switch from hAChE to hBChE. In the present work, we identified thiophene-based isosteres 7 and 15 as dual AChE-MAO B (IC50 = 261 and 15 nM, respectively) and BChE-MAO B (IC50 = 375 and 20 nM, respectively) inhibitors, respectively. Both 7 and 15 were moderately water-soluble and membrane-permeant agents by passive diffusion (PAMPA-HDM). Moreover, they were able to counteract oxidative damage induced by both H2O2 and 6-OHDA in SH-SY5Y cells and predicted to penetrate into CNS in a cell-based model mimicking blood-brain barrier. Molecular dynamics (MD) simulations shed light on key differences in AChE and BChE recognition processes promoted by the basic chain homologation from 7 to 15.

Chemical manipulation of an activation/inhibition switch in the nuclear receptor PXR

Nuclear receptors are ligand-activated transcription factors that can often be useful drug targets. Unfortunately, ligand promiscuity leads to two-thirds of receptors remaining clinically untargeted. PXR is a nuclear receptor that can be activated by diverse compounds to elevate metabolism, negatively impacting drug efficacy and safety. This presents a barrier to drug development because compounds designed to target other proteins must avoid PXR activation while retaining potency for the desired target. This problem could be avoided by using PXR antagonists, but these compounds are rare, and their molecular mechanisms remain unknown. Here, we report structurally related PXR-selective agonists and antagonists and their corresponding co-crystal structures to describe mechanisms of antagonism and selectivity. Structural and computational approaches show that antagonists induce PXR conformational changes incompatible with transcriptional coactivator recruitment. These results guide the design of compounds with predictable agonist/antagonist activities and bolster efforts to generate antagonists to prevent PXR activation interfering with other drugs.

Comprehensive analysis of two hotspot codons in the TUBB4B gene and associated phenotypes

Our purpose was to elucidate the genotype and ophthalmological and audiological phenotype in TUBB4B-associated inherited retinal dystrophy (IRD) and sensorineural hearing loss (SNHL), and to model the effects of all possible amino acid substitutions at the hotspot codons Arg390 and Arg391. Six patients from five families with heterozygous missense variants in TUBB4B were included in this observational study. Ophthalmological testing included best-corrected visual acuity, fundus examination, optical coherence tomography, fundus autofluorescence imaging, and full-field electroretinography (ERG). Audiological examination included pure-tone and speech audiometry in adult patients and auditory brainstem response testing in a child. Genetic testing was performed by disease gene panel analysis based on genome sequencing. The molecular consequences of the substitutions of residues 390 and 391 on TUBB4B and its interaction with α-tubulin were predicted in silico on its three-dimensional structure obtained by homology modelling. Two independent patients had amino acid exchanges at position 391 (p.(Arg391His) or p.(Arg391Cys)) of the TUBB4B protein. Both had a distinct IRD phenotype with peripheral round yellowish lesions with pigmented spots and mild or moderate SNHL, respectively. Yet the phenotype was milder with a sectorial pattern of bone spicules in one patient, likely due to a genetically confirmed mosaicism for p.(Arg391His). Three patients were heterozygous for an amino acid exchange at position 390 (p.(Arg390Gln) or p.(Arg390Trp)) and presented with another distinct retinal phenotype with well demarcated pericentral retinitis pigmentosa. All showed SNHL ranging from mild to severe. One additional patient showed a variant distinct from codon 390 or 391 (p.(Tyr310His)), and presented with congenital profound hearing loss and reduced responses in ERG. Variants at codon positions 390 and 391 were predicted to decrease the structural stability of TUBB4B and its complex with α-tubulin, as well as the complex affinity. In conclusion, the twofold larger reduction in heterodimer affinity exhibited by Arg391 substitutions suggested an association with the more severe retinal phenotype, compared to the substitution at Arg390.

Identification of Druggable Binding Sites and Small Molecules as Modulators of TMC1

Our ability to hear and maintain balance relies on the proper functioning of inner ear sensory hair cells, which translate mechanical stimuli into electrical signals via mechano-electrical transducer (MET) channels, composed of TMC1/2 proteins. However, the therapeutic use of ototoxic drugs, such as aminoglycosides and cisplatin, which can enter hair cells through MET channels, often leads to profound auditory and vestibular dysfunction. Despite extensive research on otoprotective compounds targeting MET channels, our understanding of how small molecule modulators interact with these channels remains limited, hampering the discovery of novel compounds. Here, we propose a structure-based screening approach, integrating 3D-pharmacophore modeling, molecular simulations, and experimental validation. Our pipeline successfully identified several novel compounds and FDA-approved drugs that reduced dye uptake in cultured cochlear explants, indicating MET modulation activity. Molecular docking and free-energy estimations for binding allowed us to identify three potential drug binding sites within the channel pore, phospholipids, and key amino acids involved in modulator interactions. We also identified shared ligand-binding features between TMC and structurally related TMEM16 protein families, providing novel insights into their distinct inhibition, while potentially guiding the rational design of MET-channel-specific modulators. Our pipeline offers a broad application to discover small molecule modulators for a wide spectrum of mechanosensitive ion channels.

Probing benzenesulfonamide-thiazolidinone hybrids as multitarget directed ligands for efficient control of type 2 diabetes mellitus through targeting the enzymes: α-glucosidase and carbonic anhydrase II

Diabetes mellitus is a chronic metabolic disorder characterized by improper expression/function of a number of key enzymes that can be regarded as targets for anti-diabetic drug design. Herein, we report the design, synthesis, and biological assessment of two series of thiazolidinone-based sulfonamides 4a-l and 5a-c as multitarget directed ligands (MTDLs) with potential anti-diabetic activity through targeting the enzymes: α-glucosidase and human carbonic anhydrase (hCA) II. The synthesized sulfonamides were evaluated for their inhibitory activity against α-glucosidase where most of the compounds showed good to potent activities. Compounds 4d and 4e showed potent inhibitory activities (IC50 = 0.440 and 0.3456 μM), comparable with that of the positive control (acarbose; IC50 = 0.420 μM). All the synthesized derivatives were also tested for their inhibitory activities against hCA I, II, IX, and XII. They exhibited different levels of inhibition against these isoforms. Compound 4d outstood as the most potent one against hCA II with Ki equals to 7.0 nM, more potent than the reference standard (acetazolamide; Ki = 12.0 nM). In silico studies for the most active compounds within the active sites of α-glucosidase and hCA II revealed good binding modes that can explain their biological activities. MM-GBSA refinements and molecular dynamic simulations were performed on the top-ranking docking pose of the most potent compound 4d to confirm the formation of stable complex with both targets. Compound 4d was screened for its in vivo antihyperglycemic efficacy by using the oral glucose tolerance test. Compound 4d decreased blood glucose level to 217 mg/dl, better than the standard acarbose (234 mg/dl). Hence, this revealed its synergistic mode of action on post prandial hyperglycemia and hepatic gluconeogenesis. Thus, these benzenesulfonamide thiazolidinone hybrids could be considered as promising multi-target candidates for the treatment of type II diabetes mellitus.

Structural investigation of Trypanosoma cruzi Akt-like kinase as drug target against Chagas disease

According to the World Health Organization, Chagas disease (CD) is the most prevalent poverty-promoting neglected tropical disease. Alarmingly, climate change is accelerating the geographical spreading of CD causative parasite, Trypanosoma cruzi, which additionally increases infection rates. Still, CD treatment remains challenging due to a lack of safe and efficient drugs. In this work, we analyze the viability of T. cruzi Akt-like kinase (TcAkt) as drug target against CD including primary structural and functional information about a parasitic Akt protein. Nuclear Magnetic Resonance derived information in combination with Molecular Dynamics simulations offer detailed insights into structural properties of the pleckstrin homology (PH) domain of TcAkt and its binding to phosphatidylinositol phosphate ligands (PIP). Experimental data combined with Alpha Fold proposes a model for the mechanism of action of TcAkt involving a PIP-induced disruption of the intramolecular interface between the kinase and the PH domain resulting in an open conformation enabling TcAkt kinase activity. Further docking experiments reveal that TcAkt is recognized by human inhibitors PIT-1 and capivasertib, and TcAkt inhibition by UBMC-4 and UBMC-6 is achieved via binding to TcAkt kinase domain. Our in-depth structural analysis of TcAkt reveals potential sites for drug development against CD, located at activity essential regions.

Evolutionary sequence and structural basis for the distinct conformational landscapes of Tyr and Ser/Thr kinases

Protein kinases are molecular machines with rich sequence variation that distinguishes the two main evolutionary branches - tyrosine kinases (TKs) from serine/threonine kinases (STKs). Using a sequence co-variation Potts statistical energy model we previously concluded that TK catalytic domains are more likely than STKs to adopt an inactive conformation with the activation loop in an autoinhibitory "folded" conformation, due to intrinsic sequence effects. Here we investigated the structural basis for this phenomenon by integrating the sequence-based model with structure-based molecular dynamics (MD) to determine the effects of mutations on the free energy difference between active and inactive conformations, using a novel thermodynamic cycle involving many (n=108) protein-mutation free energy perturbation (FEP) simulations in the active and inactive conformations. The sequence and structure-based results are consistent and support the hypothesis that the inactive conformation "DFG-out Activation Loop Folded", is a functional regulatory state that has been stabilized in TKs relative to STKs over the course of their evolution via the accumulation of residue substitutions in the activation loop and catalytic loop that facilitate distinct substrate binding modes in trans and additional modes of regulation in cis for TKs.

Amyloid Accelerator Polyphosphate Implicated as the Mystery Density in α-Synuclein Fibrils

Aberrant aggregation of α-Synuclein is the pathological hallmark of a set of neurodegenerative diseases termed synucleinopathies. Recent advances in cryo-electron microscopy have led to the structural determination of the first synucleinopathy-derived α-Synuclein fibrils, which contain a non-proteinaceous, "mystery density" at the core of the protofilaments, hypothesized to be highly negatively charged. Guided by previous studies that demonstrated that polyphosphate (polyP), a universally conserved polyanion, significantly accelerates α-Synuclein fibril formation, we conducted blind docking and molecular dynamics simulation experiments to model the polyP binding site in α-Synuclein fibrils. Here we demonstrate that our models uniformly place polyP into the lysine-rich pocket, which coordinates the mystery density in patient-derived fibrils. Subsequent in vitro studies and experiments in cells revealed that substitution of the two critical lysine residues K43 and K45 leads to a loss of all previously reported effects of polyP binding on α-Synuclein, including stimulation of fibril formation, change in filament conformation and stability as well as alleviation of cytotoxicity. In summary, our study demonstrates that polyP fits the unknown electron density present in in vivo α-Synuclein fibrils and suggests that polyP exerts its functions by neutralizing charge repulsion between neighboring lysine residues.

Structure-based virtual screening of novel USP5 inhibitors targeting the zinc finger ubiquitin-binding domain

The equilibrium of cellular protein levels is pivotal for maintaining normal physiological functions. USP5 belongs to the deubiquitination enzyme (DUBs) family, controlling protein degradation and preserving cellular protein homeostasis. Aberrant expression of USP5 is implicated in a variety of diseases, including cancer, neurodegenerative diseases, and inflammatory diseases. In this paper, a multi-level virtual screening (VS) approach was employed to target the zinc finger ubiquitin-binding domain (ZnF-UBD) of USP5, leading to the identification of a highly promising candidate compound 0456-0049. Molecular dynamics (MD) simulations were then employed to assess the stability of complex binding and predict hotspot residues in interactions. The results indicated that the candidate stably binds to the ZnF-UBD of USP5 through crucial interactions with residues ARG221, TRP209, GLY220, ASN207, TYR261, TYR259, and MET266. Binding free energy calculations, along with umbrella sampling (US) simulations, underscored a superior binding affinity of the candidate relative to known inhibitors. Moreover, US simulations revealed conformational changes of USP5 during ligand dissociation. These insights provide a valuable foundation for the development of novel inhibitors targeting USP5.

Bioactivity and molecular docking studies of selected plant compounds

BACKGROUND

The shortcomings of synthetic pesticides are becoming increasingly evident. In the pursuit of looking for alternative, environmentally friendly pesticides, the potential use of two naturally occurring sesquiterpenes, α-bisabolol and bisabolol oxide A, as bactericides or acaricides, was investigated against three gram-positive and three gram-negative bacteria, and the eggs and larvae of the cattle tick Hyalomma scupense, respectively. The primary targeted site of synthetic insecticides is acetylcholinesterase (AChE), so the inhibitory activity of both sesquiterpenes was evaluated. Then, molecular docking and molecular mechanics generalized Born surface area (MM-GBSA) studies were performed to gain insight into their binding mode in the AChE active site.

RESULTS

Bisabolol oxide A showed good activity against Staphylococcus aureus, Escherichia coli, and Salmonella enteritidis with inhibition zones ranging from 16 to 34 mm and minimal inhibitory concentration (CIM) of 3.125  mg mL-1. Ovicidal and larvicidal assays revealed promising activity of α-bisabolol and bisabolol oxide A against H. scupense, with bisabolol oxide A being more effective against larval ticks with lethal concentration (LC50) value of 0.78%. Both compounds also displayed potent anti-AChE activity with inhibition concentration (IC50) values of 37.09 and 28.14 μg mL-1, respectively. Furthermore, α-bisabolol and bisabolol oxide A exhibited good and comparable docking scores (-7.289 and -7.139 Kcal mol-1, respectively) and were found to accommodate in the active-site gorge of AChE via hydrogen bonding and hydrophobic interactions.

CONCLUSION

Bisabolol oxide A and, to a lesser degree, α-bisabolol are active against bacteria and ectoparasites and may represent an economical and sustainable alternative to toxic synthetic pesticides to control pathogens. © 2024 Society of Chemical Industry.

Equisetum arvense standardized dried extract hinders age-related osteosarcopenia

Age-associated osteosarcopenia is an unresolved syndrome characterized by the concomitant loss of bone (osteopenia) and skeletal muscle (sarcopenia) tissues increasing falls, immobility, morbidity, and mortality. Unbalanced resorption of bone in the remodeling process and excessive protein breakdown, especially fast type II myosin heavy chain (MyHC-II) isoform and myofiber metabolic shift, are the leading causes of bone and muscle deterioration in the elderly, respectively. Equisetum arvense (EQ) is a plant traditionally recommended for many pathological conditions due to its anti-inflammatory properties. Thus, considering that a chronic low-grade inflammatory state predisposes to both osteoporosis and sarcopenia, we tested a standardized hydroalcoholic extract of EQ in in vitro models of muscle atrophy [C2C12 myotubes treated with proinflammatory cytokines (TNFα/IFNγ), excess glucocorticoids (dexamethasone), or the osteokine, receptor activator of nuclear factor kappa-B ligand (RANKL)] and osteoclastogenesis (RAW 264.7 cells treated with RANKL). We found that EQ counteracted myotube atrophy, blunting the activity of several pathways depending on the applied stimulus, and reduced osteoclast formation and activity. By in silico target fishing, IKKB-dependent nuclear factor kappa-B (NF-κB) inhibition emerges as a potential common mechanism underlying EQ's anti-atrophic effects. Consumption of EQ (500 mg/kg/day) by pre-geriatric C57BL/6 mice for 3 months translated into: i) maintenance of muscle mass and performance; ii) restrained myofiber oxidative shift; iii) slowed down age-related modifications in osteoporotic bone, significantly preserving trabecular connectivity density; iv) reduced muscle- and spleen-related inflammation. EQ can preserve muscle functionality and bone remodeling during aging, potentially valuable as a natural treatment for osteosarcopenia.

In vitro and in silico investigation of FDA-approved drugs to be repurposed against Alzheimer's disease

Alzheimer's disease (AD), one of the main causes of dementia, is a neurodegenerative disorder. Cholinesterase inhibitors are used in the treatment of AD, but prolonged use of these drugs can lead to serious side effects. Drug repurposing is an approach that aims to reveal the effectiveness of drugs in different diseases beyond their clinical uses. In this work, we investigated in vitro and in silico inhibitory effects of 11 different drugs on cholinesterases. The results showed that trimebutine, theophylline, and levamisole had the highest acetylcholinesterase inhibitory actions among the tested drugs, and these drugs inhibited by 68.70 ± 0.46, 53.25 ± 3.40, and 44.03 ± 1.20%, respectively at 1000 µM. In addition, these drugs are bound to acetylcholinesterase via competitive manner. Molecular modeling predicted good fitness in acetylcholinesterase active site for these drugs and possible central nervous system action for trimebutine. All of these results demonstrated that trimebutine was determined to be the drug with the highest potential for use in AD.

Vasorelaxant Activity of (2S)-Sakuranetin and Other Flavonoids Isolated from the Green Propolis of the Caatinga Mimosa tenuiflora

Some in vitro and in vivo evidence is consistent with the cardiovascular beneficial activity of propolis. As the single actors responsible for this effect have never been identified, an in-depth investigation of flavonoids isolated from the green propolis of the Caatinga Mimosa tenuiflora was performed and their mechanism of action was described. A comprehensive electrophysiology, functional, and molecular docking approach was applied. Most flavanones and flavones were effective CaV1.2 channel blockers with a potency order of (2S)-sakuranetin > eriodictyol-7,3'-methyl ether > quercetin 3-methyl ether > 5,4'-dihydroxy-6,7-dimethoxyflavanone > santin > axillarin > penduletin > kumatakenin, ermanin and viscosine being weak or modest stimulators. Except for eriodictyol 5-O-methyl ether, all the flavonoids were also effective spasmolytic agents of vascular rings, kumatakenin and viscosine also showing an endothelium-dependent activity. (2S)-Sakuranetin also stimulated KCa1.1 channels both in single myocytes and vascular rings. In silico analysis provided interesting insights into the mode of action of (2S)-sakuranetin within both CaV1.2 and KCa1.1 channels. The green propolis of the Caatinga Mimosa tenuiflora is a valuable source of multi-target vasoactive flavonoids: this evidence reinforces its nutraceutical value in the cardiovascular disease prevention arena.

Extracting binding energies and binding modes from biomolecular simulations of fragment binding to endothiapepsin

Fragment-based drug discovery (FBDD) aims to discover a set of small binding fragments that may be subsequently linked together. Therefore, in-depth knowledge of the individual fragments' structural and energetic binding properties is essential. In addition to experimental techniques, the direct simulation of fragment binding by molecular dynamics (MD) simulations became popular to characterize fragment binding. However, former studies showed that long simulation times and high computational demands per fragment are needed, which limits applicability in FBDD. Here, we performed short, unbiased MD simulations of direct fragment binding to endothiapepsin, a well-characterized model system of pepsin-like aspartic proteases. To evaluate the strengths and limitations of short MD simulations for the structural and energetic characterization of fragment binding, we predicted the fragments' absolute free energies and binding poses based on the direct simulations of fragment binding and compared the predictions to experimental data. The predicted absolute free energies are in fair agreement with the experiment. Combining the MD data with binding mode predictions from molecular docking approaches helped to correctly identify the most promising fragments for further chemical optimization. Importantly, all computations and predictions were done within 5 days, suggesting that MD simulations may become a viable tool in FBDD projects.

Synthesis and Functional Evaluation of Synthetic Cannabinoid Receptor Agonists Related to ADB-HEXINACA

ADB-HEXINACA has been recently reported as a synthetic cannabinoid receptor agonist (SCRA), one of the largest classes of new psychoactive substances (NPSs). This compound marks the entry of the n-hexyl tail group into the SCRA landscape, which has continued in the market with recent, newly detected SCRAs. As such, a proactive characterization campaign was undertaken, including the synthesis, characterization, and pharmacological evaluation of ADB-HEXINACA and a library of 41 closely related analogues. Two in vitro functional assays were employed to assess activity at CB1 and CB2 cannabinoid receptors, measuring Gβγ-coupled agonism through a fluorescence-based membrane potential assay (MPA) and β-arrestin 2 (βarr2) recruitment via a live cell-based nanoluciferase complementation reporter assay. ADB-HEXINACA was a potent and efficacious CB1 agonist (CB1 MPA pEC50 = 7.87 ± 0.12 M; Emax = 124 ± 5%; βarr2 pEC50 = 8.27 ± 0.14 M; Emax = 793 ± 42.5), as were most compounds assessed. Isolation of the heterocyclic core and alkyl tails allowed for the comprehensive characterization of structure-activity relationships in this compound class, which were rationalized in silico via induced fit docking experiments. Overall, most compounds assessed are possibly emerging NPSs.

Identification of Riluzole derivatives as novel calmodulin inhibitors with neuroprotective activity by a joint synthesis, biosensor, and computational guided strategy

The development of new molecules for the treatment of calmodulin related cardiovascular or neurodegenerative diseases is an interesting goal. In this work, we introduce a novel strategy with four main steps: (1) chemical synthesis of target molecules, (2) Förster Resonance Energy Transfer (FRET) biosensor development and in vitro biological assay of new derivatives, (3) Cheminformatics models development and in vivo activity prediction, and (4) Docking studies. This strategy is illustrated with a case study. Firstly, a series of 4-substituted Riluzole derivatives 1-3 were synthetized through a strategy that involves the construction of the 4-bromoriluzole framework and its further functionalization via palladium catalysis or organolithium chemistry. Next, a FRET biosensor for monitoring Ca2+-dependent CaM-ligands interactions has been developed and used for the in vitro assay of Riluzole derivatives. In particular, the best inhibition (80%) was observed for 4-methoxyphenylriluzole 2b. Besides, we trained and validated a new Networks Invariant, Information Fusion, Perturbation Theory, and Machine Learning (NIFPTML) model for predicting probability profiles of in vivo biological activity parameters in different regions of the brain. Next, we used this model to predict the in vivo activity of the compounds experimentally studied in vitro. Last, docking study conducted on Riluzole and its derivatives has provided valuable insights into their binding conformations with the target protein, involving calmodulin and the SK4 channel. This new combined strategy may be useful to reduce assay costs (animals, materials, time, and human resources) in the drug discovery process of calmodulin inhibitors.

Machine learning and classical MD simulation to identify inhibitors against the P37 envelope protein of monkeypox virus

Monkeypox virus (MPXV) outbreak is a serious public health concern that requires international attention. P37 of MPXV plays a pivotal role in DNA replication and acts as one of the promising targets for antiviral drug design. In this study, we intent to screen potential analogs of existing FDA approved drugs of MPXV against P37 using state-of-the-art machine learning and computational biophysical techniques. AlphaFold2 guided all-atoms molecular dynamics simulations optimized P37 structure is used for molecular docking and binding free energy calculations. Similar to members of Phospholipase-D family , the predicted P37 structure also adopts a β-α-β-α-β sandwich fold, harbouring strongly conserved HxKxxxxD motif. The binding pocket comprises of Tyr48, Lys86, His115, Lys117, Ser130, Asn132, Trp280, Asn240, His325, Lys327 and Tyr346 forming strong hydrogen bonds and dense hydrophobic contacts with the screened analogs and is surrounded by positively charged patches. Loops connecting the two domains and C-terminal region exhibit high degree of flexibility. In some structural ensembles, the partial disorderness in the C-terminal region is presumed to be due to its low confidence score, acquired during structure prediction. Transition from loop to β-strands (244-254 aa) in P37-Cidofovir and its analog complexes advocates the need for further investigations. MD simulations support the accuracy of the molecular docking results, indicating the potential of analogs as potent binders of P37. Taken together, our results provide preferable understanding of molecular recognition and dynamics of ligand-bound states of P37, offering opportunities for development of new antivirals against MPXV. However, the need of in vitro and in vivo assays for confirmation of these results still persists.Communicated by Ramaswamy H. Sarma.

Development of Aptamer-DNAzyme based metal-nucleic acid frameworks for gastric cancer therapy

The metal-nucleic acid nanocomposites, first termed metal-nucleic acid frameworks (MNFs) in this work, show extraordinary potential as functional nanomaterials. However, thus far, realized MNFs face limitations including harsh synthesis conditions, instability, and non-targeting. Herein, we discover that longer oligonucleotides can enhance the synthesis efficiency and stability of MNFs by increasing oligonucleotide folding and entanglement probabilities during the reaction. Besides, longer oligonucleotides provide upgraded metal ions binding conditions, facilitating MNFs to load macromolecular protein drugs at room temperature. Furthermore, longer oligonucleotides facilitate functional expansion of nucleotide sequences, enabling disease-targeted MNFs. As a proof-of-concept, we build an interferon regulatory factor-1(IRF-1) loaded Ca2+/(aptamer-deoxyribozyme) MNF to target regulate glucose transporter (GLUT-1) expression in human epidermal growth factor receptor-2 (HER-2) positive gastric cancer cells. This MNF nanodevice disrupts GSH/ROS homeostasis, suppresses DNA repair, and augments ROS-mediated DNA damage therapy, with tumor inhibition rate up to 90%. Our work signifies a significant advancement towards an era of universal MNF application.

Construction of High Accuracy Machine Learning Interatomic Potential for Surface/Interface of Nanomaterials-A Review

The inherent discontinuity and unique dimensional attributes of nanomaterial surfaces and interfaces bestow them with various exceptional properties. These properties, however, also introduce difficulties for both experimental and computational studies. The advent of machine learning interatomic potential (MLIP) addresses some of the limitations associated with empirical force fields, presenting a valuable avenue for accurate simulations of these surfaces/interfaces of nanomaterials. Central to this approach is the idea of capturing the relationship between system configuration and potential energy, leveraging the proficiency of machine learning (ML) to precisely approximate high-dimensional functions. This review offers an in-depth examination of MLIP principles and their execution and elaborates on their applications in the realm of nanomaterial surface and interface systems. The prevailing challenges faced by this potent methodology are also discussed.

Sulfonyl thioureas with a benzo[d]thiazole ring as dual acetylcholinesterase/butyrylcholinesterase and human monoamine oxidase A and B inhibitors: An in vitro and in silico study

A series of sulfonyl thioureas 6a-q containing a benzo[d]thiazole ring with an ester functional group was synthesized from corresponding substituted 2-aminobenzo[d]thiazoles 3a-q and p-toluenesulfonyl isothiocyanate. They had remarkable inhibitory activity against acetylcholinesterase (AChE), butyrylcholinesterase (BChE), monoamine oxidase (MAO)-A, and MAO-B. Among thioureas, several compounds had notable activity in the order of 6k > 6 h > 6c (AChE), 6j > 6g > 6k (BChE), 6k > 6g > 6f (MAO-A), and 6i > 6k > 6h (MAO-B). Compound 6k was an inhibitor of interest due to its potent or good activity against all studied enzymes, with IC50 values of 0.027 ± 0.008 μM (AChE), 0.043 ± 0.004 μM (BChE), 0.353 ± 0.01 μM (MAO-A), and 0.716 ± 0.02 μM (MAO-B). This inhibitory capacity was comparable to that of the reference drugs for each enzyme. Kinetic studies of two compounds with potential activity, 6k (against AChE) and 6j (against BChE), had shown that both 6k and 6j followed competitive-type enzyme inhibition, with Ki constants of 24.49 and 12.16 nM, respectively. Induced fit docking studies for enzymes 4EY7, 7BO4, 2BXR, and 2BYB showed active interactions between sulfonyl thioureas of benzo[d]thiazoles and the residues in the active pocket with ligands 6k, 6i, and 6j, respectively. The stability of the ligand-protein complexes while each ligand entered the active site of each enzyme (4EY7, 7BO4, 2BXR, or 2BYB) was confirmed by molecular dynamics simulations.

Modular preparation of biphenyl triazoles via click chemistry as non-competitive hyaluronidase inhibitors

Hyaluronidase is a promising target in drug discovery, given its overexpression in a range of physiological and pathological processes, including tumor migration, skin aging, sagging, and wrinkling, as well as inflammation and bacterial infections. In this study, to identify novel hyaluronidase inhibitors, we applied click chemistry for the modular synthesis of 370 triazoles in 96-well plates, starting with biphenyl azide. Utilizing an optimized turbidimetric screening assay in microplates, we identified Fmoc-containing triazoles 5 and 6, as well as quinoline-containing triazoles 15 and 16, as highly effective hyaluronidase inhibitors. Subsequent research indicated that these triazoles potentially interact with a novel binding site of hyaluronidase. Notably, these inhibitors displayed minimal cytotoxicity and showed promising anti-inflammatory effects in LPS-stimulated macrophages. Remarkably, compound 6 significantly reduced NO release by 74 % at a concentration of 20 μM.

Binding properties of selective inhibitors of P323L mutated RdRp of SARS-CoV-2: a combined molecular screening, docking and dynamics simulation study

Since 2019 the SARS-CoV-2 and its variants caused COVID-19, such incidents brought the world in pandemic situation. This happened due to furious mutations in SARS-CoV-2, in which some variants had high transmissibility and infective, this led the virus emerged as virulent and worsened the COVID-19 situation. Among the variants, P323L is one of the important mutants of RdRp in SARS-CoV-2. To inhibit the erroneous function of this mutated RdRp, we have screened 943 molecules against the P323L mutated RdRp with the criteria that the molecules with 90% similar to the structure of remdesivir (control drug) resulted nine molecules. Further, these molecules were evaluated by induced fit docking (IFD) identified two molecules (M2 & M4) which are forming strong intermolecular interactions with the key residues of mutated RdRp and has high binding affinity. Docking score of the M2 and M4 molecules with mutated RdRp are -9.24 and -11.87 kcal/mol, respectively. Further, to understand the intermolecular interactions, conformational stability, the molecular dynamics simulation and binding free energy calculations were performed. The binding free energy values of M2 and M4 molecules with the P323L mutated RdRp complexes are -81.60 and -83.07 kcal/mol, respectively. The results of this in silico study confirm that M4 is a potential molecule; hence, it may be considered as the potential inhibitor of P323L mutated RdRp to treat COVID-19 after clinical investigation.Communicated by Ramaswamy H. Sarma.

Do electrostatic interactions make a difference in physics-based AutoDock4 scoring function?

Physics-based scoring function AutoDock4 is one of the most successfully applied tools in the area of structure-based drug design. However, current scoring functions are still far from being perfect. In a recent work highlighting the strengths and deficiencies of current scoring functions, we discovered that the residual error of ΔGbind predictions made by AutoDock4 is highly correlated to the presence of formally charged fragments in a ligand. In this work, we study how the use of the high-quality atomic charges, applied for contemporary force fields calculation, affects the quality of the experimental ΔGbind prediction by means of AutoDock4. We initially expected that the previously found discrepancy could be attributed to the Gasteiger charges used within AutoDock4. We show that AutoDock4 is, surprisingly, not sensitive to the charges used, and the use of QC-derived atomic charges does not lead to any statistical improvements. We also briefly discuss the role of the explicit empirical hydrogen bond term along with the electrostatic term.

Flexible scaffold-based cheminformatics approach for polypharmacological drug design

Effective treatments for complex central nervous system (CNS) disorders require drugs with polypharmacology and multifunctionality, yet designing such drugs remains a challenge. Here, we present a flexible scaffold-based cheminformatics approach (FSCA) for the rational design of polypharmacological drugs. FSCA involves fitting a flexible scaffold to different receptors using different binding poses, as exemplified by IHCH-7179, which adopted a "bending-down" binding pose at 5-HT2AR to act as an antagonist and a "stretching-up" binding pose at 5-HT1AR to function as an agonist. IHCH-7179 demonstrated promising results in alleviating cognitive deficits and psychoactive symptoms in mice by blocking 5-HT2AR for psychoactive symptoms and activating 5-HT1AR to alleviate cognitive deficits. By analyzing aminergic receptor structures, we identified two featured motifs, the "agonist filter" and "conformation shaper," which determine ligand binding pose and predict activity at aminergic receptors. With these motifs, FSCA can be applied to the design of polypharmacological ligands at other receptors.

Robust prediction of relative binding energies for protein-protein complex mutations using free energy perturbation calculations

Computational free energy-based methods have the potential to significantly improve throughput and decrease costs of protein design efforts. Such methods must reach a high level of reliability, accuracy, and automation to be effectively deployed in practical industrial settings in a way that impacts protein design projects. Here, we present a benchmark study for the calculation of relative changes in protein-protein binding affinity for single point mutations across a variety of systems from the literature, using free energy perturbation (FEP+) calculations. We describe a method for robust treatment of alternate protonation states for titratable amino acids, which yields improved correlation with and reduced error compared to experimental binding free energies. Following careful analysis of the largest outlier cases in our dataset, we assess limitations of the default FEP+ protocols and introduce an automated script which identifies probable outlier cases that may require additional scrutiny and calculates an empirical correction for a subset of charge-related outliers. Through a series of three additional case study systems, we discuss how protein FEP+ can be applied to real-world protein design projects, and suggest areas of further study.

Computational Design of Novel Cyclic Peptides Endowed with Autophagy-Inhibiting Activity on Cancer Cell Lines

(1) Autophagy plays a significant role in development and cell proliferation. This process is mainly accomplished by the LC3 protein, which, after maturation, builds the nascent autophagosomes. The inhibition of LC3 maturation results in the interference of autophagy activation. (2) In this study, starting from the structure of a known LC3B binder (LIR2-RavZ peptide), we identified new LC3B ligands by applying an in silico drug design strategy. The most promising peptides were synthesized, biophysically assayed, and biologically evaluated to ascertain their potential antiproliferative activity on five humans cell lines. (3) A cyclic peptide (named Pep6), endowed with high conformational stability (due to the presence of a disulfide bridge), displayed a Kd value on LC3B in the nanomolar range. Assays accomplished on PC3, MCF-7, and A549 cancer cell lines proved that Pep6 exhibited cytotoxic effects comparable to those of the peptide LIR2-RavZ, a reference LC3B ligand. Furthermore, it was ineffective on both normal prostatic epithelium PNT2 and autophagy-defective prostate cancer DU145 cells. (4) Pep6 can be considered a new autophagy inhibitor that can be employed as a pharmacological tool or even as a template for the rational design of new small molecules endowed with autophagy inhibitory activity.