A Computational Workflow for the Identification of Novel Fragments Acting as Inhibitors of the Activity of Protein Kinase CK1δ

Fragment-Based Drug Discovery (FBDD) has become, in recent years, a consolidated approach in the drug discovery process, leading to several drug candidates under investigation in clinical trials and some approved drugs. Among these successful applications of the FBDD approach, kinases represent a class of targets where this strategy has demonstrated its real potential with the approved kinase inhibitor Vemurafenib. In the Kinase family, protein kinase CK1 isoform δ (CK1δ) has become a promising target in the treatment of different neurodegenerative diseases such as Alzheimer’s disease, Parkinson’s disease, and amyotrophic lateral sclerosis. In the present work, we set up and applied a computational workflow for the identification of putative fragment binders in large virtual databases. To validate the method, the selected compounds were tested in vitro to assess the CK1δ inhibition.     

Computational Analysis of the Interactions between the S100B Extracellular Chaperone and Its Amyloid β Peptide Client

S100B is an astrocytic extracellular Ca²⁺-binding protein implicated in Alzheimer’s disease, whose role as a holdase-type chaperone delaying Aβ₄₂ aggregation and toxicity was recently uncovered. Here, we employ computational biology approaches to dissect the structural details and dynamics of the interaction between S100B and Aβ₄₂. Driven by previous structural data, we used the Aβ_25–35 segment, which recapitulates key aspects of S100B activity, as a starting guide for the analysis. We used Haddock to establish a preferred binding mode, which was studied with the full length Aβ using long (1 μs) molecular dynamics (MD) simulations to investigate the structural dynamics and obtain representative interaction complexes. From the analysis, Aβ-Lys28 emerged as a key candidate for stabilizing interactions with the S100B binding cleft, in particular involving a triad composed of Met79, Thr82 and Glu86. Binding constant calculations concluded that coulombic interactions, presumably implicating the Lys28(Aβ)/Glu86(S100B) pair, are very relevant for the holdase-type chaperone activity. To confirm this experimentally, we examined the inhibitory effect of S100B over Aβ aggregation at high ionic strength. In agreement with the computational predictions, we observed that electrostatic perturbation of the Aβ-S100B interaction decreases anti-aggregation activity. Altogether, these findings unveil features relevant in the definition of selectivity of the S100B chaperone, with implications in Alzheimer’s disease.     

Molecular Structure of Cefuroxime Axetil Complexes with α-, β-, γ-, and 2-Hydroxypropyl-β-Cyclodextrins: Molecular Simulations and Raman Spectroscopic and Imaging Studies

The formation of cefuroxime axetil+cyclodextrin (CA+CD) complexes increases the aqueous solubility of CA, improves its physico-chemical properties, and facilitates a biomembrane-mediated drug delivery process. In CD-based tablet formulations, it is crucial to investigate the molecular details of complexes in final pharmaceutical preparation. In this study, Raman spectroscopy and mapping were applied for the detection and identification of chemical groups involved in α-, β-, γ-, and 2-hydroxypropyl-β-CD (2-HP- β-CD)+CA complexation process. The experimental studies have been complemented by molecular dynamics-based investigations, providing additional molecular details of CA+CD interactions. It has been demonstrated that CA forms the guest–host type inclusion complexes with all studied CDs; however, the nature of the interactions is slightly different. It seems that both α- and β-CD interact with furanyl and methoxy moieties of CA, γ-CD forms a more diverse pattern of interactions with CA, which are not observed in other CDs, whereas 2HP-β-CD binds CA with the contribution of hydrogen bonding. Apart from supporting this interpretation of the experimental data, molecular dynamics simulations allowed for ordering the CA+CD binding affinities. The obtained results proved that the molecular details of the host–guest complexation can be successfully predicted from the combination of Raman spectroscopy and molecular modeling.     

α4-α5 Helices on Surface of KRAS Can Accommodate Small Compounds That Increase KRAS Signaling While Inducing CRC Cell Death

KRAS is the most frequently mutated oncogene associated with the genesis and progress of pancreatic, lung and colorectal (CRC) tumors. KRAS has always been considered as a therapeutic target in cancer but until now only two compounds that inhibit one specific KRAS mutation have been approved for clinical use. In this work, by molecular dynamics and a docking process, we describe a new compound (P14B) that stably binds to a druggable pocket near the α4-α5 helices of the allosteric domain of KRAS. This region had previously been identified as the binding site for calmodulin (CaM). Using surface plasmon resonance and pulldown analyses, we prove that P14B binds directly to oncogenic KRAS thus competing with CaM. Interestingly, P14B favors oncogenic KRAS interaction with BRAF and phosphorylated C-RAF, and increases downstream Ras signaling in CRC cells expressing oncogenic KRAS. The viability of these cells, but not that of the normal cells, is impaired by P14B treatment. These data support the significance of the α4-α5 helices region of KRAS in the regulation of oncogenic KRAS signaling, and demonstrate that drugs interacting with this site may destine CRC cells to death by increasing oncogenic KRAS downstream signaling.     

Effects of Visfatin on Intracellular Mechanics and Catabolism in Human Primary Chondrocytes through Glycogen Synthase Kinase 3β Inactivation

Osteoarthritis (OA) is still a recalcitrant musculoskeletal disease on account of its complex biochemistry and mechanical stimulations. Apart from stimulation by external mechanical forces, the regulation of intracellular mechanics in chondrocytes has also been linked to OA development. Recently, visfatin has received significant attention because of the clinical finding of the positive correlation between its serum/synovial level and OA progression. However, the precise mechanism involved is still unclear. This study determined the effect of visfatin on intracellular mechanics and catabolism in human primary chondrocytes isolated from patients. The intracellular stiffness of chondrocytes was analyzed by the particle-tracking microrheology method. It was shown that visfatin damages the microtubule and microfilament networks to influence intracellular mechanics to decrease the intracellular elasticity and viscosity via glycogen synthase kinase 3β (GSK3β) inactivation induced by p38 signaling. Further, microtubule network destruction in human primary chondrocytes is predominantly responsible for the catabolic effect of visfatin on the cyclooxygenase 2 upregulation. The present study shows a more comprehensive interpretation of OA development induced by visfatin through biochemical and biophysical perspectives. Finally, the role of GSK3β inactivation, and subsequent regulation of intracellular mechanics, might be considered as theranostic targets for future drug development for OA.     

The Interaction between microRNAs and the Wnt/β-Catenin Signaling Pathway in Osteoarthritis

Osteoarthritis (OA) is a chronic disease affecting the whole joint, which still lacks a disease-modifying treatment. This suggests an incomplete understanding of underlying molecular mechanisms. The Wnt/β-catenin pathway is involved in different pathophysiological processes of OA. Interestingly, both excessive stimulation and suppression of this pathway can contribute to the pathogenesis of OA. microRNAs have been shown to regulate different cellular processes in different diseases, including the metabolic activity of chondrocytes and osteocytes. To bridge these findings, here we attempt to give a conclusive overview of microRNA regulation of the Wnt/β-catenin pathway in bone and cartilage, which may provide insights to advance the development of miRNA-based therapeutics for OA treatment.     

The Effects of One-Point Mutation on the New Delhi Metallo Beta-Lactamase-1 Resistance toward Carbapenem Antibiotics and β-Lactamase Inhibitors: An In Silico Systematic Approach

Antibiotic resistance has been becoming more and more critical due to bacteria’s evolving hydrolysis enzymes. The NDM-1 enzyme could hydrolyze not only carbapenems but also most of β-lactam’s antibiotics and inhibitors. In fact, variant strains could impose a high impact on the resistance of bacteria producing NDM-1. Although previous studies showed the effect of some variants toward antibiotics and inhibitors binding, there has been no research systematically evaluating the effects of alternative one-point mutations on the hydrolysis capacity of NDM-1. This study aims to identify which mutants could increase or decrease the effectiveness of antibiotics and β-lactamase inhibitors toward bacteria. Firstly, 35 different variants with a high probability of emergence based on the PAM-1 matrix were constructed and then docked with 5 ligands, namely d-captopril, l-captopril, thiorphan, imipenem, and meropenem. The selected complexes underwent molecular dynamics simulation and free energy binding estimation, with the results showing that the substitutions at residues 122 and 124 most influenced the binding ability of NDM-1 toward inhibitors and antibiotics. The H122R mutant decreases the binding ability between d-captopril and NDM-1 and diminishes the effectiveness of this antibiotic toward Enterobacteriaceae. However, the H122R mutant has a contrary impact on thiorphan, which should be tested in vitro and in vivo in further experiments.     

Molecular Docking and Molecular Dynamics Simulation Studies of Triterpenes from Vernonia patula with the Cannabinoid Type 1 Receptor

A molecular docking approach was employed to evaluate the binding affinity of six triterpenes, namely epifriedelanol, friedelin, α-amyrin, α-amyrin acetate, β-amyrin acetate, and bauerenyl acetate, towards the cannabinoid type 1 receptor (CB1). Molecular docking studies showed that friedelin, α-amyrin, and epifriedelanol had the strongest binding affinity towards CB1. Molecular dynamics simulation studies revealed that friedelin and α-amyrin engaged in stable non-bonding interactions by binding to a pocket close to the active site on the surface of the CB1 target protein. The studied triterpenes showed a good capacity to penetrate the blood–brain barrier. These results help to provide some evidence to justify, at least in part, the previously reported antinociceptive and sedative properties of Vernonia patula.     

Inhibition of Inducible Nitric Oxide Synthase Prevents IL-1β-Induced Mitochondrial Dysfunction in Human Chondrocytes

Interleukin (IL)-1β is an important pro-inflammatory cytokine in the progression of osteoarthritis (OA), which impairs mitochondrial function and induces the production of nitric oxide (NO) in chondrocytes. The aim was to investigate if blockade of NO production prevents IL-1β-induced mitochondrial dysfunction in chondrocytes and whether cAMP and AMP-activated protein kinase (AMPK) affects NO production and mitochondrial function. Isolated human OA chondrocytes were stimulated with IL-1β in combination with/without forskolin, L-NIL, AMPK activator or inhibitor. The release of NO, IL-6, PGE₂, MMP3, and the expression of iNOS were measured by ELISA or Western blot. Parameters of mitochondrial respiration were measured using a seahorse analyzer. IL-1β significantly induced NO release and mitochondrial dysfunction. Inhibition of iNOS by L-NIL prevented IL-1β-induced NO release and mitochondrial dysfunction but not IL-1β-induced release of IL-6, PGE₂, and MMP3. Enhancement of cAMP by forskolin reduced IL-1β-induced NO release and prevented IL-1β-induced mitochondrial impairment. Activation of AMPK increased IL-1β-induced NO production and the negative impact of IL-1β on mitochondrial respiration, whereas inhibition of AMPK had the opposite effects. NO is critically involved in the IL-1β-induced impairment of mitochondrial respiration in human OA chondrocytes. Increased intracellular cAMP or inhibition of AMPK prevented both IL-1β-induced NO release and mitochondrial dysfunction.     

Molecular Modeling and Molecular Dynamics Studies of Hydralazine with Human DNA Methyltransferase 1

A series of DNA methyltransferase 1 (DNMT1) inhibitors were modeled by docking and molecular dynamics studies to rationalize their activity. Our findings will be valuable in guiding research efforts toward the rational design and virtual screening of novel DNMT inhibitors.

     

Synthesis,Biochemistry, and Computational Studies of Brominated Thienyl Chalcones: A New Class of Reversible MAO‐B Inhibitors

A series of (2E)‐1‐(5‐bromothiophen‐2‐yl)‐3‐(para‐substituted phenyl)prop‐2‐en‐1‐ones ( TB1 – TB11 ) was synthesized and tested for inhibitory activity toward human monoamine oxidase (hMAO). All compounds were found to be competitive, selective, and reversible toward hMAO‐B except (2E)‐1‐(5‐bromothiophen‐2‐yl)‐3‐(4‐nitrophenyl)prop‐2‐en‐1‐one ( TB7 ) and (2E)‐1‐(5‐bromothiophen‐2‐yl)‐3‐(4‐chlorophenyl)prop‐2‐en‐1‐one ( TB8 ), which were selective inhibitors of hMAO‐A. The most potent compound, (2E)‐1‐(5‐bromothiophen‐2‐yl)‐3‐[4‐(dimethylamino)phenyl]prop‐2‐en‐1‐one ( TB5 ), showed the best inhibitory activity and higher selectivity toward hMAO‐B, with K_i and SI values of 0.11±0.01 μm and 13.18, respectively. PAMPA assays for all compounds were carried out in order to evaluate the capacity of the compounds to cross the blood–brain barrier. Moreover, the most potent MAO‐B inhibitor, TB5 , was found to be nontoxic at 5 and 25 μm , with 95.75 and 84.59 % viability among cells, respectively. Molecular docking simulations were carried out to understand the crucial interactions responsible for selectivity and potency.     

Co-Incubation with PPARβ/δ Agonists and Antagonists Modeled Using Computational Chemistry: Effect on LPS Induced Inflammatory Markers in Pulmonary Artery

Peroxisome proliferator activated receptor beta/delta (PPARβ/δ) is a nuclear receptor ubiquitously expressed in cells, whose signaling controls inflammation. There are large discrepancies in understanding the complex role of PPARβ/δ in disease, having both anti- and pro-effects on inflammation. After ligand activation, PPARβ/δ regulates genes by two different mechanisms; induction and transrepression, the effects of which are difficult to differentiate directly. We studied the PPARβ/δ-regulation of lipopolysaccharide (LPS) induced inflammation (indicated by release of nitrite and IL-6) of rat pulmonary artery, using different combinations of agonists (GW0742 or L−165402) and antagonists (GSK3787 or GSK0660). LPS induced release of NO and IL-6 is not significantly reduced by incubation with PPARβ/δ ligands (either agonist or antagonist), however, co-incubation with an agonist and antagonist significantly reduces LPS-induced nitrite production and Nos2 mRNA expression. In contrast, incubation with LPS and PPARβ/δ agonists leads to a significant increase in Pdk−4 and Angptl−4 mRNA expression, which is significantly decreased in the presence of PPARβ/δ antagonists. Docking using computational chemistry methods indicates that PPARβ/δ agonists form polar bonds with His287, His413 and Tyr437, while antagonists are more promiscuous about which amino acids they bind to, although they are very prone to bind Thr252 and Asn307. Dual binding in the PPARβ/δ binding pocket indicates the ligands retain similar binding energies, which suggests that co-incubation with both agonist and antagonist does not prevent the specific binding of each other to the large PPARβ/δ binding pocket. To our knowledge, this is the first time that the possibility of binding two ligands simultaneously into the PPARβ/δ binding pocket has been explored. Agonist binding followed by antagonist simultaneously switches the PPARβ/δ mode of action from induction to transrepression, which is linked with an increase in Nos2 mRNA expression and nitrite production.     

Serotonin Promotes Serum Albumin Interaction with the Monomeric Amyloid β Peptide

Prevention of amyloid β peptide (Aβ) deposition via facilitation of Aβ binding to its natural depot, human serum albumin (HSA), is a promising approach to preclude Alzheimer’s disease (AD) onset and progression. Previously, we demonstrated the ability of natural HSA ligands, fatty acids, to improve the affinity of this protein to monomeric Aβ by a factor of 3 (BBRC, 510(2), 248–253). Using plasmon resonance spectroscopy, we show here that another HSA ligand related to AD pathogenesis, serotonin (SRO), increases the affinity of the Aβ monomer to HSA by a factor of 7/17 for Aβ₄₀/Aβ₄₂, respectively. Meanwhile, the structurally homologous SRO precursor, tryptophan (TRP), does not affect HSA’s affinity to monomeric Aβ, despite slowdown of the association and dissociation processes. Crosslinking with glutaraldehyde and dynamic light scattering experiments reveal that, compared with the TRP-induced effects, SRO binding causes more marked changes in the quaternary structure of HSA. Furthermore, molecular docking reveals distinct structural differences between SRO/TRP complexes with HSA. The disintegration of the serotonergic system during AD pathogenesis may contribute to Aβ release from HSA in the central nervous system due to impairment of the SRO-mediated Aβ trapping by HSA.     

Promotion and Inhibition of Amyloid-β Peptide Aggregation: Molecular Dynamics Studies

Aggregates of amyloid-β (Aβ) peptides are known to be related to Alzheimer’s disease. Their aggregation is enhanced at hydrophilic–hydrophobic interfaces, such as a cell membrane surface and air-water interface, and is inhibited by polyphenols, such as myricetin and rosmarinic acid. We review molecular dynamics (MD) simulation approaches of a full-length Aβ peptide, Aβ40, and Aβ(16–22) fragments in these environments. Since these peptides have both hydrophilic and hydrophobic amino acid residues, they tend to exist at the interfaces. The high concentration of the peptides accelerates the aggregation there. In addition, Aβ40 forms a β-hairpin structure, and this structure accelerates the aggregation. We also describe the inhibition mechanism of the Aβ(16–22) aggregation by polyphenols. The aggregation of Aβ(16–22) fragments is caused mainly by the electrostatic attraction between charged amino acid residues known as Lys16 and Glu22. Since polyphenols form hydrogen bonds between their hydroxy and carboxyl groups and these charged amino acid residues, they inhibit the aggregation.     

Biochemical Characterization of a Carboxylesterase from the Archaeon Pyrobaculum sp. 1860 and a Rational Explanation of Its Substrate Specificity and Thermostability

In this work, genome mining was used to identify esterase/lipase genes in the archaeon Pyrobaculum sp. 1860. A gene was cloned and functionally expressed in Escherichia coli as His-tagged protein. The recombinant enzyme (rP186_1588) was verified by western blotting and peptide mass fingerprinting. Biochemical characterization revealed that rP186_1588 exhibited optimum activity at pH 9.0 and 80 °C towards p-nitrophenyl acetate (K_m: 0.35 mM, k_cat: 11.65 s⁻¹). Interestingly, the purified rP186_1588 exhibited high thermostability retaining 70% relative activity after incubation at 90 °C for 6 h. Circular dichroism results indicated that rP186_1588 showed slight structure alteration from 60 to 90 °C. Structural modeling showed P186_1588 possessed a typical α/β hydrolase’s fold with the catalytic triad consisting of Ser₉₇, Asp₁₄₇ and His₁₇₂, and was further confirmed by site-directed mutagenesis. Comparative molecular simulations at different temperatures (300, 353, 373 and 473 K) revealed that its thermostability was associated with its conformational rigidity. The binding free energy analysis by MM-PBSA method revealed that the van der Waals interaction played a major role in p-NP ester binding for P186_1588. Our data provide insights into the molecular structures of this archaeal esterase, and may help to its further protein engineering for industrial applications.     

Disclosing the Interaction of Gold Nanoparticles with Aβ(1–40) Monomers through Replica Exchange Molecular Dynamics Simulations

Amyloid-β aggregation is one of the principal causes of amyloidogenic diseases that lead to the loss of neuronal cells and to cognitive impairments. The use of gold nanoparticles treating amyloidogenic diseases is a promising approach, because the chemistry of the gold surface can be tuned in order to have a specific binding, obtaining effective tools to control the aggregation. In this paper, we show, by means of Replica Exchange Solute Tempering Molecular Simulations, how electrostatic interactions drive the absorption of Amyloid-β monomers onto citrates-capped gold nanoparticles. Importantly, upon binding, amyloid monomers show a reduced propensity in forming β-sheets secondary structures that are characteristics of mature amyloid fibrils.