Conjugated β-Cyclodextrin Enhances the Affinity of Folic Acid towards FRα: Molecular Dynamics Study

Drug targeting is a progressive area of research with folate receptor alpha (FRα) receiving significant attention as a biological marker in cancer drug delivery. The binding affinity of folic acid (FA) to the FRα active site provides a basis for recognition of FRα. In this study, FA was conjugated to beta-cyclodextrin (βCD) and subjected to in silico analysis (molecular docking and molecular dynamics (MD) simulation (100 ns)) to investigate the affinity and stability for the conjugated system compared to unconjugated and apo systems (ligand free). Docking studies revealed that the conjugated FA bound into the active site of FRα with a docking score (free binding energy < −15 kcal/mol), with a similar binding pose to that of unconjugated FA. Subsequent analyses from molecular dynamics (MD) simulations, root mean square deviation (RMSD), root mean square fluctuation (RMSF), and radius of gyration (Rg) demonstrated that FA and FA–βCDs created more dynamically stable systems with FRα than the apo-FRα system. All systems reached equilibrium with stable RMSD values ranging from 1.9–2.4 Å and the average residual fluctuation values of the FRα backbone atoms for all residues (except for terminal residues ARG8, THR9, THR214, and LEU215) were less than 2.1 Å with a consistent Rg value of around 16.8 Å throughout the MD simulation time (0–100 ns). The conjugation with βCD improved the stability and decreased the mobility of all the residues (except residues 149–151) compared to FA–FRα and apo-FRα systems. Further analysis of H-bonds, binding free energy (MM-PBSA), and per residue decomposition energy revealed that besides APS81, residues HIS20, TRP102, HIS135, TRP138, TRP140, and TRP171 were shown to have more favourable energy contributions in the holo systems than in the apo-FRα system, and these residues might have a direct role in increasing the stability of holo systems.     

Effects of Curcumin and Ferulic Acid on the Folding of Amyloid-β Peptide

The polyphenols curcumin (CU) and ferulic acid (FA) are able to inhibit the aggregation of amyloid-β (Aβ) peptide with different strengths. CU is a strong inhibitor while FA is a weaker one. In the present study, we examine the effects of CU and FA on the folding process of an Aβ monomer by 1 µs molecular dynamics (MD) simulations. We found that both inhibitors increase the helical propensity and decrease the non-helical propensity of Aβ peptide. They prevent the formation of a dense bulk core and shorten the average lifetime of intramolecular hydrogen bonds in Aβ. CU makes more and longer-lived hydrogen bonds, hydrophobic, π–π, and cation–π interactions with Aβ peptide than FA does, which is in a good agreement with the observed stronger inhibitory activity of CU on Aβ aggregation.     

Novel and Potential Small Molecule Scaffolds as DYRK1A Inhibitors by Integrated Molecular Docking-Based Virtual Screening and Dynamics Simulation Study

The dual-specificity tyrosine phosphorylation-regulated kinase 1A (DYRK1A) is a novel, promising and emerging biological target for therapeutic intervention in neurodegenerative diseases, especially in Alzheimer’s disease (AD). The molMall database, comprising rare, diverse and unique compounds, was explored for molecular docking-based virtual screening against the DYRK1A protein, in order to find out potential inhibitors. Ligands exhibiting hydrogen bond interactions with key amino acid residues such as Ile165, Lys188 (catalytic), Glu239 (gk+1), Leu241 (gk+3), Ser242, Asn244, and Asp307, of the target protein, were considered potential ligands. Hydrogen bond interactions with Leu241 (gk+3) were considered key determinants for the selection. High scoring structures were also docked by Glide XP docking in the active sites of twelve DYRK1A related protein kinases, viz. DYRK1B, DYRK2, CDK5/p25, CK1, CLK1, CLK3, GSK3β, MAPK2, MAPK10, PIM1, PKA, and PKCα, in order to find selective DYRK1A inhibitors. MM/GBSA binding free energies of selected ligand–protein complexes were also calculated in order to remove false positive hits. Physicochemical and pharmacokinetic properties of the selected six hit ligands were also computed and related with the proposed limits for orally active CNS drugs. The computational toxicity webserver ProTox-II was used to predict the toxicity profile of selected six hits (molmall IDs 9539, 11352, 15938, 19037, 21830 and 21878). The selected six docked ligand–protein systems were exposed to 100 ns molecular dynamics (MD) simulations to validate their mechanism of interactions and stability in the ATP pocket of human DYRK1A kinase. All six ligands were found to be stable in the ATP binding pocket of DYRK1A kinase.     

Design,Synthesis, and Biological Evaluation of Tetra-Substituted Thiophenes as Inhibitors of p38α MAPK

p38α mitogen-activated protein kinase (MAPK) plays a role in several cellular processes and consequently has been a therapeutic target in inflammatory diseases, cancer, and cardiovascular disease. A number of known p38α MAPK inhibitors contain vicinal 4-fluorophenyl/4-pyridyl rings connected to either a 5- or 6-membered heterocycle. In this study, a small library of substituted thiophene-based compounds bearing the vicinal 4-fluorophenyl/4-pyridyl rings was designed using computational docking as a visualisation tool. Compounds were synthesised and evaluated in a fluorescence polarisation binding assay. The synthesised analogues had a higher binding affinity to the active phosphorylated form of p38α MAPK than the inactive nonphosphorylated form of the protein. 4-(2-(4-fluorophenyl)thiophen-3-yl)pyridine had a K_i value of 0.6 μm to active p38α MAPK highlighting that substitution of the core ring to a thiophene retains affinity to the enzyme and can be utilised in p38α MAPK inhibitors. This compound was further elaborated using a substituted phenyl ring in order to probe the second hydrophobic pocket. Many of these analogues exhibited low micromolar affinity to active p38α MAPK. The suppression of neonatal rat fibroblast collagen synthesis was also observed suggesting that further development of these compounds may lead to potential therapeutics having cardioprotective properties.     

Molecular Dynamic Simulations of Bromodomain and Extra-Terminal Protein 4 Bonded to Potent Inhibitors

Bromodomain and extra-terminal domain (BET) subfamily is the most studied subfamily of bromodomain-containing proteins (BCPs) family which can modulate acetylation signal transduction and produce diverse physiological functions. Thus, the BET family can be treated as an alternative strategy for targeting androgen-receptor (AR)-driven cancers. In order to explore the effect of inhibitors binding to BRD4 (the most studied member of BET family), four 150 ns molecular dynamic simulations were performed (free BRD4, Cpd4-BRD4, Cpd9-BRD4 and Cpd19-BRD4). Docking studies showed that Cpd9 and Cpd19 were located at the active pocket, as well as Cpd4. Molecular dynamics (MD) simulations indicated that only Cpd19 binding to BRD4 can induce residue Trp81-Ala89 partly become α-helix during MD simulations. MM-GBSA calculations suggested that Cpd19 had the best binding effect with BRD4 followed by Cpd4 and Cpd9. Computational alanine scanning results indicated that mutations in Phe83 made the greatest effects in Cpd9-BRD4 and Cpd19-BRD4 complexes, showing that Phe83 may play crucial roles in Cpd9 and Cpd19 binding to BRD4. Our results can provide some useful clues for further BCPs family search.     

Computational Simulations Highlight the IL2Rα Binding Potential of Polyphenol Stilbenes from Fenugreek

Widely used in global households, fenugreek is well known for its culinary and medicinal uses. The various reported medicinal properties of fenugreek are by virtue of the different natural phytochemicals present in it. Regarded as a promising target, interleukin 2 receptor subunit alpha (IL2Rα) has been shown to influence immune responses. In the present research, using in silico techniques, we have demonstrated the potential IL2Rα binding properties of three polyphenol stilbenes (desoxyrhaponticin, rhaponticin, rhapontigenin) from fenugreek. As the first step, molecular docking was performed to assess the binding potential of the fenugreek phytochemicals with IL2Rα. All three phytochemicals demonstrated interactions with active site residues. To confirm the reliability of our molecular docking results, 100 ns molecular dynamics simulations studies were undertaken. As discerned by the RMSD and RMSF analyses, IL2Rα in complex with the desoxyrhaponticin, rhaponticin, and rhapontigenin indicated stability. The RMSD analysis of the phytochemicals alone also demonstrated no significant structural changes. Based on the stable molecular interactions and comparatively slightly better MM/PBSA binding free energy, rhaponticin seems promising. Additionally, ADMET analysis performed for the stilbenes indicated that all of them obey the ADMET rules. Our computational study thus supports further in vitro IL2Rα binding studies on these stilbenes, especially rhaponticin.     

Identification of a nanomolar affinity α-synuclein fibril imaging probe by ultra-high throughput in silico screening

Small molecules that bind with high affinity and specificity to fibrils of the α-synuclein (αS) protein have the potential to serve as positron emission tomography (PET) imaging probes to aid in the diagnosis of Parkinson''s disease and related synucleinopathies. To identify such molecules, we employed an ultra-high throughput in silico screening strategy using idealized pseudo-ligands termed exemplars to identify compounds for experimental binding studies. For the top hit from this screen, we used photo-crosslinking to confirm its binding site and studied the structure–activity relationship of its analogs to develop multiple molecules with nanomolar affinity for αS fibrils and moderate specificity for αS over Aβ fibrils. Lastly, we demonstrated the potential of the lead analog as an imaging probe by measuring binding to αS-enriched homogenates from mouse brain tissue using a radiolabeled analog of the identified molecule. This study demonstrates the validity of our powerful new approach to the discovery of PET probes for challenging molecular targets.     

Potential Anti-Alzheimer Agents from Guanidinyl Tryptophan Derivatives with Activities of Membrane Adhesion and Conformational Transition Inhibitions

Guanidinyl tryptophan derivatives TGN1, TGN2, TGN3, and TGN4 were synthesized, and these compounds were shown to possess in vitro inhibitory activity for amyloid aggregation in a previous study. Nevertheless, the influence of the TGN series of compounds on the binding and permeation behaviors of an Aβ monomer to the cell membranes was not elucidated. In this study, we investigated the effect of compounds in the TGN series on the behavior of an Aβ monomer regarding its toxicity toward the bilayer lipid membrane using molecular dynamics (MD) simulation. MD simulations suggest that TGN4 is a potential agent that can interfere with the movement of the Aβ monomer into the membrane. The MM-GBSA result demonstrated that TGN4 exhibits the highest affinity to the Aβ_1–42 monomer but has the lowest affinity to the bilayer. Moreover, TGN4 also contributes to a decrease in the binding affinity between the Aβ_1–42 monomer and the POPC membrane. Regarding the results of the binding mode and conformational analyses, a high number of amino-acid residues were shown to provide the binding interactions between TGN4 and the Aβ_1–42 monomer. TGN4 also reduces the conformational transition of the Aβ_1–42 monomer by means of interacting with the monomer. The present study presents molecular-level insights into how the TGN series of compounds affect the membrane adsorption and the conformational transition of the Aβ_1–42 monomer, which could be valuable for the further development of new anti-Alzheimer agents.     

Identification of a Novel Delta Opioid Receptor Agonist Chemotype with Potential Negative Allosteric Modulator Capabilities

The δ-opioid receptor (δOR) holds great potential as a therapeutic target. Yet, clinical drug development, which has focused on δOR agonists that mimic the potent and selective tool compound SNC80 have largely failed. It has increasingly become apparent that the SNC80 scaffold carries with it potent and efficacious β-arrestin recruitment. Here, we screened a relatively small (5120 molecules) physical drug library to identify δOR agonists that underrecruit β-arrestin, as it has been suggested that compounds that efficaciously recruit β-arrestin are proconvulsant. The screen identified a hit compound and further characterization using cellular binding and signaling assays revealed that this molecule (R995045, compound 1) exhibited ten-fold selectivity over µ- and κ-opioid receptors. Compound 1 represents a novel chemotype at the δOR. A subsequent characterization of fourteen analogs of compound 1, however did not identify a more potent δOR agonist. Computational modeling and in vitro characterization of compound 1 in the presence of the endogenous agonist leu-enkephalin suggest compound 1 may also bind allosterically and negatively modulate the potency of Leu-enkephalin to inhibit cAMP, acting as a ‘NAM-agonist’ in this assay. The potential physiological utility of such a class of compounds will need to be assessed in future in vivo assays.     

Molecular Dynamic Simulation Analysis on the Inclusion Complexation of Plumbagin with β-Cyclodextrin Derivatives in Aqueous Solution

Stable encapsulation of medically active compounds can lead to longer storage life and facilitate the slow-release mechanism. In this work, the dynamic and molecular interactions between plumbagin molecule with β-cyclodextrin (BCD) and its two derivatives, which are dimethyl-β-cyclodextrin (MBCD), and 2-O-monohydroxypropyl-β-cyclodextrin (HPBCD) were investigated. Molecular dynamics simulations (MD) with GLYCAM-06 and AMBER force fields were used to simulate the inclusion complex systems under storage temperature (4 °C) in an aqueous solution. The simulation results suggested that HPBCD is the best encapsulation agent to produce stable host–guest binding with plumbagin. Moreover, the observation of the plumbagin dynamic inside the binding cavity revealed that it tends to orient the methyl group toward the wider rim of HPBCD. Therefore, HPBCD is a decent candidate for the preservation of plumbagin with a promising longer storage life and presents the opportunity to facilitate the slow-release mechanism.     

Inhibitory Effect of Fisetin on α-Glucosidase Activity: Kinetic and Molecular Docking Studies

The inhibition of α-glucosidase is a clinical strategy for the treatment of type 2 diabetes mellitus (T2DM), and many natural plant ingredients have been reported to be effective in alleviating hyperglycemia by inhibiting α-glucosidase. In this study, the α-glucosidase inhibitory activity of fisetin extracted from Cotinus coggygria Scop. was evaluated in vitro. The results showed that fisetin exhibited strong inhibitory activity with an IC₅₀ value of 4.099 × 10⁻⁴ mM. Enzyme kinetic analysis revealed that fisetin is a non-competitive inhibitor of α-glucosidase, with an inhibition constant value of 0.01065 ± 0.003255 mM. Moreover, fluorescence spectrometric measurements indicated the presence of only one binding site between fisetin and α-glucosidase, with a binding constant (lgKa) of 5.896 L·mol⁻¹. Further molecular docking studies were performed to evaluate the interaction of fisetin with several residues close to the inactive site of α-glucosidase. These studies showed that the structure of the complex was maintained by Pi-Sigma and Pi-Pi stacked interactions. These findings illustrate that fisetin extracted from Cotinus coggygria Scop. is a promising therapeutic agent for the treatment of T2DM.     

Theoretical studies on the structural rearrangement of ligand binding pocket in human vitamin D receptor

The structural rearrangement of the ligand binding domain (LBD) of human Vitamin D receptor (hVDR) complexed with 1α, 25‐dihydroxyvitamin D₃ (natural ligand) and its analogues (denoted as b and c ) was studied by molecular dynamics (MD) simulations. MD simulations revealed that these ligands could induce different structural changes of LBD, in which 1α, 25‐dihydroxyvitamin D₃ only led to a minute change, suggesting that LBD adopted its canonical active conformation upon binding the natural ligand, while b and c could provoke a clear structural rearrangement of the LBD. In complex of hVDR‐LBD/ b , it is found that helix 6 (H6) and subsequent loop 6‐7 shift outward and the last turn of H11 shifts away from H12, which generate a new cavity at the bottom of binding pocket to accommodate the extra butyl group on the side chain of ligand b . As for hVDR‐LBD/ c , the steric exclusion of the second side chain of ligand c makes the N‐terminal of H7 move outsides and C‐terminal of H11 close to H12, expanding the bottom of the pocket. These calculation results agree well the experimental observations. © 2010 Wiley Periodicals, Inc. Int J Quantum Chem, 2011     

<Emphasis Type="Italic">In-silico</Emphasis> Screening using Flexible Ligand Binding Pockets: A Molecular Dynamics-based Approach

Summary In-silico screening of flexible ligands against flexible ligand binding pockets (LBP) is an emerging approach in structure-based drug discovery. Here, we describe a molecular dynamics (MD) based docking approach to investigate the influence on the high-throughput in-silico screening of small molecules against flexible ligand binding pockets. In our approach, an ensemble of 51 energetically favorable structures of the LBP of human estrogen receptor α (hERα) were collected from 3 ns MD simulations. In-silico screening of 3500 endocrine disrupting compounds against these flexible ligand binding pockets resulted in thousands of ER–ligand complexes of which 582 compounds were unique. Detailed analysis of MD generated structures showed that only 17 of the LBP residues significantly contribute to the overall binding pocket flexibility. Using the flexible LBP conformations generated, we have identified 32 compounds that bind better to the flexible ligand-binding pockets compared to the crystal structure. These compounds, though chemically divergent, are structurally similar to the natural hormone. Our MD-based approach in conjunction with grid–based distributed computing could be applied routinely for in-silico screening of large databases against any given target.     

Virtual Screening and Molecular Dynamics Simulation Study of Influenza Polymerase PB2 Inhibitors

Influenza A virus is the main cause of worldwide epidemics and annual influenza outbreaks in humans. In this study, a virtual screen was performed to identify compounds that interact with the PB2 cap-binding domain (CBD) of influenza A polymerase. A virtual screening workflow based on Glide docking was used to screen an internal database containing 8417 molecules, and then the output compounds were selected based on solubility, absorbance, and structural fingerprints. Of the 16 compounds selected for biological evaluation, six compounds were identified that rescued cells from H1N1 virus-mediated death at non-cytotoxic concentrations, with EC50 values ranging from 2.5–55.43 μM, and that could bind to the PB2 CBD of H1N1, with Kd values ranging from 0.081–1.53 μM. Molecular dynamics (MD) simulations of the docking complexes of our active compounds revealed that each compound had its own binding characteristics that differed from those of VX-787. Our active compounds have novel structures and unique binding modes with PB2 proteins, and are suitable to serve as lead compounds for the development of PB2 inhibitors. An analysis of the MD simulation also helped us to identify the dominant amino acid residues that play a key role in binding the ligand to PB2, suggesting that we should focus on increasing and enhancing the interaction between inhibitors and these major amino acids during lead compound optimization to obtain more active PB2 inhibitors.     

Discovery of Amide-Functionalized Benzimidazolium Salts as Potent α-Glucosidase Inhibitors

α-Glucosidase inhibitors (AGIs) are used as medicines for the treatment of diabetes mellitus. The α-Glucosidase enzyme is present in the small intestine and is responsible for the breakdown of carbohydrates into sugars. The process results in an increase in blood sugar levels. AGIs slow down the digestion of carbohydrates that is helpful in controlling the sugar levels in the blood after meals. Among heterocyclic compounds, benzimidazole moiety is recognized as a potent bioactive scaffold for its wide range of biologically active derivatives. The aim of this study is to explore the α-glucosidase inhibition ability of benzimidazolium salts. In this study, two novel series of benzimidazolium salts, i.e., 1-benzyl-3-{2-(substituted) amino-2-oxoethyl}-1H-benzo[d]imidazol-3-ium bromide 9a–m and 1-benzyl-3-{2-substituted) amino-2-oxoethyl}-2-methyl-1H-benzo[d] imidazol-3-ium bromide 10a–m were screened for their in vitro α-glucosidase inhibitory potential. These compounds were synthesized through a multistep procedure and were characterized by ¹H-NMR, ¹³C-NMR, and EI-MS techniques. Compound 10d was identified as the potent α-glucosidase inhibitor among the series with an IC₅₀ value of 14 ± 0.013 μM, which is 4-fold higher than the standard drug, acarbose. In addition, compounds 10a, 10e, 10h, 10g, 10k, 10l, and 10m also exhibited pronounced potential for α-glucosidase inhibition with IC₅₀ value ranging from 15 ± 0.037 to 32.27 ± 0.050 µM when compared with the reference drug acarbose (IC₅₀ = 58.8 ± 0.12 μM). A molecular docking study was performed to rationalize the binding interactions of potent inhibitors with the active site of the α-glucosidase enzyme.     

Holospiniferoside: A New Antitumor Cerebroside from The Red Sea Cucumber Holothuria spinifera: In Vitro and In Silico Studies

Chemical investigation of the methanolic extract of the Red Sea cucumber Holothuria spinifera led to the isolation of a new cerebroside, holospiniferoside (1), together with thymidine (2), methyl-α-d-glucopyranoside (3), a new triacylglycerol (4), and cholesterol (5). Their chemical structures were established by NMR and mass spectrometric analysis, including gas chromatography–mass spectrometry (GC–MS) and high-resolution mass spectrometry (HRMS). All the isolated compounds are reported in this species for the first time. Moreover, compound 1 exhibited promising in vitro antiproliferative effect on the human breast cancer cell line (MCF-7) with IC₅₀ of 20.6 µM compared to the IC₅₀ of 15.3 µM for the drug cisplatin. To predict the possible mechanism underlying the cytotoxicity of compound 1, a docking study was performed to elucidate its binding interactions with the active site of the protein Mdm2–p53. Compound 1 displayed an apoptotic activity via strong interaction with the active site of the target protein. This study highlights the importance of marine natural products in the design of new anticancer agents.     

Synthesis,Pharmacological Evaluation,and Computational Studies of Cyclic Opioid Peptidomimetics Containing β3-Lysine

Our formerly described pentapeptide opioid analog Tyr-c[D-Lys-Phe-Phe-Asp]NH₂ (designated RP-170), showing high affinity for the mu (MOR) and kappa (KOR) opioid receptors, was much more stable than endomorphine-2 (EM-2) in the rat brain homogenate and displayed remarkable antinociceptive activity after central (intracerebroventricular) and peripheral (intravenous ) administration. In this report, we describe the further modification of this analog, which includes the incorporation of a β³-amino acid, (R)- and (S)-β³-Lys, instead of D-Lys in position 2. The influence of such replacement on the biological properties of the obtained analogs, Tyr-c[(R)-β³-Lys-Phe-Phe-Asp]NH₂ (RP-171) and Tyr-c[(S)-β³-Lys-Phe-Phe-Asp]NH₂, (RP-172), was investigated in vitro. Receptor radiolabeled displacement and functional calcium mobilization assays were performed to measure binding affinity and receptor activation of the new analogs. The obtained data revealed that only one of the diastereoisomeric peptides, RP-171, was able to selectively bind and activate MOR. Molecular modeling (docking and molecular dynamics (MD) simulations) suggests that both compounds should be accommodated in the MOR binding site. However, in the case of the inactive isomer RP-172, fewer hydrogen bonds, as well as instability of the canonical ionic interaction to Asp¹⁴⁷, could explain its very low MOR affinity.